Search NASA⌕ Search

SEARCH · Search NASA

Results for “theoretical calculations”

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 19 records

CompF2: Theoretical Calculations and Simulation Topical Group Report

This report summarizes the work of the Computational Frontier topical group on theoretical calculations and simulation for Snowmass 2021. We discuss the challenges, potential solutions, and needs facing six diverse but related topical areas that span the subject of theoretical calculations and simulation in high energy physics (HEP): cosmic calculations, particle accelerator modeling, detector simulation, event generators, perturbative calculations, and lattice QCD (quantum chromodynamics). The challenges arise from the next generations of HEP experiments, which will include more complex instruments, provide larger data volumes, and perform more precise measurements. Calculations and simulations will need to keep up with these increased requirements. The other aspect of the challenge is the evolution of computing landscape away from general-purpose computing on CPUs and toward special-purpose accelerators and coprocessors such as GPUs and FPGAs. These newer devices can provide substantial improvements for certain categories of algorithms, at the expense of more specialized programming and memory and data access patterns.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Unveiling the electronic structure and chemical bonding of the deprotonated cisplatin anion [(NH 3 )(NH 2 )PtCl 2 ] − via low-temperature photoelectron spectroscopy and theoretical calculations

The dehydrogenated cisplatin anion, [(NH 3 )(NH 2 )PtCl 2 ] − , was investigated via low-temperature photoelectron spectroscopy and theoretical calculations. Seven and four spectral peaks are respectively resolved at 193 and 266 nm, yielding rich electronic structure information for both the anion and neutral. From the threshold and maximum of the lowest electron binding energy band, the experimental adiabatic (ADE) and vertical detachment energies (VDE) are determined to be 3.3 ± 0.1 and 3.525 ± 0.025 eV, respectively. Theoretical calculations indicate the dominant isomer adopting a cis-geometry, in which the platinum center is coplanar with two chlorine and two nitrogen ligands. The calculated VDE of 3.54 eV based on this structure agrees well with the experimental value. Charge analyses reveal that the excess electron in the anion is primarily localized on the Pt and Cl atoms. A suite of theoretical analysis tools was employed to elucidate the bonding characteristics and interaction strength between Pt and its ligands.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Theoretical Calculation Study of Nanomaterials

In recent years, theoretical simulations using computational physics and chemistry models have been used widely to predict the physical, chemical, and structural properties of nanomaterials. To reflect such rapid progress, two Special Issues on the theoretical study of nanomaterials were published. The 36 articles collected cover broad re-search fields and applications. In this Editorial, we provide an overview of these pub-lications and summarize the main findings. Based on these research outcomes, we lay out some future research directions.

First-principles calculations↗

A theoretical calculation of stacking fault energy of Ni alloys: The effects of temperature and composition

Combining cluster expansion (CE) method with one dimensional axial Ising model, this work investigated the effects of alloying elements and configurational variations due to temperature on the stacking fault energy (SFE) of FCC Ni binary alloys. Ensembles of large numbers of atomistic structures, each with more than ~400 atoms, were generated to consider sufficient long-range chemical disorder and temperature effects due to configurational entropy. Here, a Monte Carlo Metropolis algorithm was used to generate these structures, whose energies were then evaluated based on the effective cluster interactions obtained from CE. As a baseline, this work had shown the SFE of pure Ni and Al to be 127 mJ/m 2 and 137 mJ/m 2 , respectively, which agreed with the experimental values of 125 mJ/m 2 and 150 mJ/m 2 reported in the literature. Additions of Al, Ti, Cr and Co to pure Ni were found to decrease the SFE to different extents. Although temperature does not strongly influence the SFE of the FCC Ni-Al and Ni-Cr binary alloys, it can lead to significant changes to the SFE of the FCC Ni-Ti and Ni-Co alloys. While effects of temperature and composition on SFE observed in this work were calculated from binary Ni alloys, the general trends are nonetheless expected to be valid in the γ phase of multicomponent Ni superalloys.

36 MATERIALS SCIENCE↗

Threshold Ionization Spectroscopy and Theoretical Calculations of LnO (Ln = La and Ce)

Mass-analyzed threshold ionization (MATI) spectroscopy was used to measure the vibronic spectra of LnO (Ln = La and Ce). Single-reference coupled cluster and relativistic multireference configuration calculations were carried out to compare with the measured vibronic energies. The spectrum of LaO displays a single vibronic band system, while that of CeO shows multiple ones. The ionization energies of LaO and CeO are measured as 5.2446 (6) and 5.3332(6) eV, respectively, which are a hundredfold improvement over the literature values. The vibrational energies of the neutral molecule and corresponding ion reveal the charge effect on the metal-oxygen bond of both species. The single band system in the spectrum of LaO arises from the transition of the ground state of the neutral molecule with the La(6s 1 )O(2p 6 ) valence configuration to the ground state of the singly charged ion with the La(6s o )O(2p 6 ) configuration. Here, the multiple band systems in the spectrum of CeO are attributed to the spin-orbit coupling for the Ce(4f 1 6s 1 )O(2p 6 ) configuration of the neutral molecule and an excited state for the Ce(4f 1 )O(2p 6 ) configuration of the ion.

74 ATOMIC AND MOLECULAR PHYSICS↗

Informing air–carbon ablation modeling with theoretical calculations of atomic oxygen and nitrogen interacting with carbon surfaces

To understand the gas–surface chemistry above the thermal protection system of a hypersonic vehicle, it is necessary to map out the kinetics of key elementary reaction steps. In this work, extensive periodic density functional theory (DFT) calculations are performed to elucidate the interaction of atomic oxygen and nitrogen with both the basal plane and edge sites of highly oriented pyrolytic graphite (HOPG). Reaction energies and barriers are determined for adsorption, desorption, diffusion, recombination, and several reactions. These DFT results are compared with the most recent finite-rate model for air–carbon ablation. Our DFT results corroborated some of the parameters used in the model but suggest that further refinement may be necessary for others. The calculations reported here will help to establish a predictive kinetic model for the complex reaction network present under hypersonic flight conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Experimental determination of the magnetic anisotropy in five-coordinated Co( II ) field-induced single molecule magnets

Magnetic anisotropy of the central metal atom is a crucial property of single molecule magnets (SMMs). Small structural changes can alter the magnetic properties, and accurate experimental methods to investigate magnetic anisotropy are therefore critical. Here, we investigate two five-coordinated Co( II ) SMMs, [CoCl 2 Cltpy] (1) and [CoBr 2 Cltpy] (2) (Cltpy = 4′-chloro-2,2′:6′,2′′-terpyridine), through multiple techniques. Ab initio theoretical calculations performed on the two compounds show that both possess axial magnetic anisotropy with the magnetic easy axis pointing towards one of the terminal halogen atoms. Theoretical calculations on SMMs are typically done on isolated molecular species, and to validate this approximation the magnetic anisotropy was further studied through experimental techniques. EPR measurements confirm an axial anisotropy of 1, and magnetic measurements provide experimental Zero-Field Splitting (ZFS) parameters, showing that the values from theoretical calculations are slightly overestimated. The X-ray electron density determined from 20 K single-crystal synchrotron X-ray diffraction data provides estimated d-orbital populations also suggesting axial magnetic anisotropy in both systems, and furthermore suggesting a more pronounced axiality in 1 compared to 2. This is in good agreement with the results obtained from both magnetic measurements and theoretical calculations. The magnetic anisotropy of 1 is quantified experimentally through polarized powder neutron diffraction via the site susceptibility method, confirming an axial magnetic anisotropy of the compound. A slight deviation in the easy axis direction is observed between experimental and theoretical results. This, together with the overestimation of the ZFS parameters from theoretical calculations, shows that experimental investigation of the magnetic anisotropy of SMMs is of high relevance. Magnetic anisotropy of the central metal atom is a crucial property of single molecule magnets (SMMs).

Slavensky, Hannah H. [Aarhus Univ. (Denmark)] (ORC↗

Collapse of Nuclear Collectivity along the 𝑁 = 𝑍 Line

The lifetime of the 𝐽 𝜋 = 2 + state in the self-conjugate 88 Ru nucleus has been determined in an experiment performed using rare-isotope beams provided by the new Facility for Rare Isotope Beams. This is the heaviest 𝑁 = 𝑍 nucleus for which such a measurement has been achieved. 88 Ru was populated by both one-neutron knockout and charge-exchange reactions, and the lifetime of 14.3$^{+2.5}_{−3.4}$ ps was determined using the triple-foil plunger technique. The extracted electromagnetic transition strength shows that the quadrupole collectivity has dropped significantly compared with the highly deformed 𝑁 = 𝑍 region around 𝐴 ∼ 80. These results are compared with state-of-the-art large-scale shell-model and discrete nonorthogonal shell-model calculations. The theoretical calculations indicate a moderate triaxial deformation and suggest that low-lying states in this nucleus are no longer dominated by strong many-particle many-hole excitations, unlike the lighter, highly deformed 𝑁 = 𝑍 nuclei nearby.

59 ≤ A ≤ 89↗

Theoretical and calculable dependent variables and their covariance in nuclear data libraries [Slides]

This presentation begins by defining observables, including theoretical observable, calculable observable, and measured observable. It also provides definitions and examples of experimental effects. Additionally, uncertainty in evaluated libraries and a metric to probe evaluated uncertainty is presented. In conclusion, being related to theoretical quantities such as resonance parameters, the uncertainty in ENDF libraries (e.g., ENDF/B-VIII.0) can be overestimated because it is usually evaluated under the guidance of experimental uncertainty: lim Δσexp.corr.→0 Δσ = Δσ theoretical model . The presentation states that a clear distinction of the uncertainty between nuclear theoretical models and experimental corrections should be revisited and this is important since (e.g., transport) simulations need theoretical quantities convoluted with specific operational parameters and material configurations. It also states that strong coupling between uncertainty quantification methodologies and optimization procedures exists, and it is necessary to develop methodologies to obtain uncertainty on theoretical models from physical and mathematical constraints.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Measurements of Gamow-Teller transitions from 59 Co via the 59 Co ⁢(𝑡, 3 He +𝛾) charge-exchange reaction and its application to the stellar electron-capture rates

Electron-capture reactions on iron-group nuclei play a crucial role in the late stages of massive star evolution. Since stellar evolution simulations depend on accurate electron-capture rates—which are highly sensitive to the detailed Gamow-Teller (GT) strength distributions—reliable theoretical models are essential. However, experimental data on GT strength distributions are scarce. High-resolution measurements are therefore vital for benchmarking and improving these theoretical calculations. To provide high-resolution data on Gamow-Teller strength distributions of iron-group nuclei and to compare these results with theoretical calculations within this mass region. Differential cross sections for the 59 Co ⁢(𝑡, 3 He)⁢ 59 Fe charge-exchange reaction at 115 MeV/u were measured using the S800 spectrometer. Furthermore, to resolve individual levels that are not distinguishable in the S800 particle singles data, coincident 𝛾 rays from the 59 Fe residual nucleus were detected by using the Gamma-Ray Energy Tracking In-beam Nuclear Array 𝛾-ray tracking array. Here, the Gamow-Teller transition strength distribution from the ground state of 59 Co to 59 Fe was extracted up to an excitation energy of 10 MeV. Additionally, transition strengths for several low-lying states were determined from coincident 𝛾-ray measurements. Electron-capture rates calculated using the present data indicate that these low-lying states contribute significantly to the overall rates in relevant stellar environments. The experimental results show reasonable agreement with theoretical predictions based on both shell-model and projected shell-model calculations. High-resolution data on Gamow-Teller strength distributions—particularly for individual low-lying states—are essential for accurately determining electron-capture rates in iron-group nuclei. Coincident 𝛾-ray measurements provide a powerful tool for obtaining such detailed information. While the present work demonstrates that shell-model calculations successfully reproduce the experimental results, such comparisons are scarce and more experimental data are desirable.

59 ≤ A ≤ 89↗

Measurement of ω meson production in pp collisions at $\sqrt{s}$ = 13 TeV

The p T -differential cross section of ω meson production in pp collisions at $\sqrt{s}$ =13 TeV at midrapidity (|y| < 0.5) was measured with the ALICE detector at the LHC, covering an unprecedented transverse-momentum range of 1.6 < p T < 50 GeV/c. The meson is reconstructed via the ω → π + π − π 0 decay channel. The results are compared with various theoretical calculations: PYTHIA8.2 with the Monash 2013 tune overestimates the data by up to 50%, whereas good agreement is observed with Next-to-Leading Order (NLO) calculations incorporating ω fragmentation using a broken SU(3) model. The ω/π 0 ratio is presented and compared with theoretical calculations and the available measurements at lower collision energies. The presented data triples the p T ranges of previously available measurements. A constant ratio of C ω/π0 = 0.578 ± 0.006 (stat.) ± 0.013 (syst.) is found above a transverse momentum of 4 GeV/c, which is in agreement with previous findings at lower collision energies within the systematic and statistical uncertainties.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Understanding Biases in Sample Preparation Techniques for Coupled Scanning Electron Microscopy and MAMA PuO 2 Morphological Analysis

In this project, the scanning electron microscopy (SEM) sampling method used during the statistical design study (SDS) was investigated to determine if any sampling biases were present in the analyzed data. Using standard particle size distribution powders from the National Institute of Standards and Technology (NIST 1984 standard reference material) with the origin wet dispersion method, it was determined that a bias to smaller particles was present. This was supported by theoretical calculations using Stokes’ law to determine the settling rate of spherical particles of roughly the same size and mass as those found in the SDS. Based on the theoretical calculations, it was determined that the settling rate for each of the 76 powder sets in the SDS could be unique based on specific particle shape and mass distributions, making a universal correction factor/formula not applicable. Therefore, priority shifted to developing an improved wet dispersion method that significantly reduced the particle settling rate for all particle size and shapes. This was achieved by replacing the original solvent (isopropyl alcohol) with a heavy liquid (lithium heteropolytungstates), which dramatically slowed the settling rate and allowed for the capture of a suitable homogeneous aliquot. SEM imaging and Morphological Analysis for Material Attribution (MAMA) software analysis were conducted on the NIST standard, and the SEM/MAMA data were compared to data captured by a dynamic image analysis particle size analyzer. The resulting data confirmed that the new wet dispersion method does indeed deliver an improved representative aliquot to the SEM stub. For instance, in the NIST certificate, the average particle size is ~17.1 µm ± 2.2 µm with a normal distribution. The initial wet dispersion method resulted in a drastically reduced average particle size of 6.1 µm in addition to a non-representative heavy bi-modal distribution whereas the improved LST wet dispersion method resulting in an average particle size that was much closer to the NIST certificate (12.7 µm) with a similar normal distribution. Although the improved method was still short of the NIST certificate average, atomic force microscopy analysis determined that the resulting ~20-25% reduction in size was due to particles sinking into the carbon sticky tape used for SEM imaging. It is believed that that this bias can be calibrated in a much more predicable manner than the original settling rate bias. In addition, the matching normal distribution curves between the NIST certificate and the heavy liquid method indicate a much-improved representative aliquot has been sampled and imaged. A surrogate CeO 2 powder was used to reflect PuO 2 more accurately and to aid in implementing radiological controls and shielding. The resulting data sets from the SEM/MAMA method and the particle size analyzer give almost identical average particle sizes and particle distribution statistics. Future work will re-analyze several select runs from the SDS to determine if morphological signatures can be found with the improved sampling method.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Electronic structure of a nodal line semimetal candidate TbSbTe

The 𝐿⁢𝑛⁢SbTe (𝐿⁢𝑛 = Lanthanides) family, like isostructural ZrSiS-type compounds, has emerged as a fertile playground for exploring the interaction of electronic correlations and magnetic ordering with the nodal line band topology. Here, we report on a detailed electronic band structure investigation of TbSbTe, corroborated by electrical transport, thermodynamic, and magnetic studies. Temperature-dependent magnetic susceptibility and thermodynamic transport studies indicate the onset of antiferromagnetic ordering below 𝑇 𝑁 ∼ 5.8K. The electronic band structure study, carried out with high-resolution angle-resolved photoemission spectroscopy measurements aided with density functional theory–based first-principles calculations reveal presence of a nonsymmorphic symmetry-protected Dirac crossing in the Γ-X high-symmetry (HS) direction, which is part of a nodal line along the X-R HS direction. Another Dirac crossing occurs along the Γ-X direction at a relatively higher binding energy, which occurs from $\tilde{𝒞}$ 2⁢𝜈⁢ 𝒫 symmetry which is gapped in the theoretical calculations with the effect of spin-orbit coupling considered. Parallel to this direction, our theoretical calculations and experimental results exhibit strongly momentum-dependent surface bands. In this study, we open an avenue to further uncover the intricate interplay among symmetry-protected topological band structure, spin-orbit coupling, and magnetism in this material and the 𝐿⁢𝑛⁢SbTe family, in general.

Angle-resolved photoemission spectroscopy↗

Spectroscopic Investigation into P(NDI2OD-T2) Charge Localization

P(NDI2OD-T2), commonly referred to as N2200, stands out as a promising electron-transporting (n-type) polymer for low-cost, flexible (photo)electrochemical applications due to its reversible two-electron reduction and high electron mobility. UV–vis spectroelectrochemistry in the tetrabutyl ammonium hexafluorophosphate/acetonitrile electrolyte shows two sets of chemically reversible redox signals in the cyclic voltammetry corresponding to the reduction of the neutral polymer film to polaronic and bipolaronic species. These electrochemical signatures suggest a distinct electronic reorganization upon reduction (polaron/bipolaron formation), highlighting the need for molecular-level insights into how charges are accommodated within the polymer backbone. While it has been previously hypothesized that charge predominantly localizes on the naphthalene diimide (NDI) unit during reductive charging, specific changes in atomic environments that confirm this localization have not been characterized in n-type polymers. Herein, we use near-edge X-ray absorption fine structure (NEXAFS) spectroscopy to probe electronic transitions in an electrochemically charged polymer to deduce charge localization. The O K-edge (1s) spectra exhibit two distinct π* peaks; the intensity of the lower-energy π* a peak that corresponds to an excitation to a largely localized carbonyl state decreases with reductive potentials relative to the higher-energy π* b peak. We corroborate this with Raman spectroscopy at different potentials, which shows a decrease in intensity on the C–C/C═C and C═O stretching bands of NDI as well as a red shift of the carbonyl band due to the formation of a polaron on the NDI. Additionally, new Raman active NDI signals associated with elongated C═O and C═C bonds are observed at lower energy during the formation of charged states. Together with theoretical calculations, these findings show that the injected charge spatially localizes on the NDI units and is dominantly distributed on the carbonyl groups. In conclusion, the combination of NEXAFS, optical and vibrational spectroscopies, and theoretical calculations is generalizable to other pi-conjugated polymers and can identify charge localization for the further development of organic semiconductors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

New Measurements of the Beam-Normal Single Spin Asymmetry in Elastic Electron Scattering Over a Range of Spin-0 Nuclei

We report precision determinations of the beam-normal single spin asymmetries (A n ) in the elastic scattering of 0.95 and 2.18 GeV electrons off 12 C, 40 Ca, 48 Ca, and 208 Pb at very forward angles where the most detailed theoretical calculations have been performed. The first measurements of A n for 40 Ca and 48 Ca are found to be similar to that of 12 C, consistent with expectations and thus demonstrating the validity of theoretical calculations for nuclei with Z ≤ 20. We also report A n for 208 Pb at two new momentum transfers (Q 2 ) extending the previous measurement. Our new data confirm the surprising result previously reported, with all three data points showing significant disagreement with the results from the Z ≤ 20 nuclei. These data confirm our basic understanding of the underlying dynamics that govern A n for nuclei containing ≲ 50 nucleons, but point to the need for further investigation to understand the unusual A n behavior discovered for scattering off 208 Pb.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

The role of collisional energy transfer in the thermal and prompt dissociation of 1-methyl allyl

Collisional energy transfer plays a pivotal role in a-priori calculations of the pressure-dependent rate constants obtained from a master equation. However, accurate determinations of collisional energy transfer parameters are rare and so most kinetic studies rely on best-fits of these parameters to experimental measurements (if available) or estimates guided from literature studies. In this work, we have quantified the effect of the uncertainty in energy transfer parameters on the thermal and prompt dissociation kinetics of a resonance stabilized radical, 1-methyl allyl (1MA), of relevance to the combustion of 1- and 2-butene isomers. Simulations using literature kinetics models were performed to assess the impact of these uncertainties on flame propagation and speciation data in laminar flames of 1- and 2-butene. Analyses of the uncertainty propagated by the energy transfer parameters for 1MA dissociation to the flame simulations, in particular the laminar flame speed, indicate an intricate coupling between the kinetics of 1MA dissociation (chain propagation) and its reaction with H-atoms (chain-termination). Ab-initio based theoretical calculations were also performed to obtain pressure-dependent kinetics for the reaction of 1MA with H-atoms. Lastly, theoretically calculated energy transfer parameters were used to best characterize the kinetics and branching between the chain propagating 1MA dissociation and the chain-terminating reaction 1MA + H.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗