Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computational inelasticity”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Computing inelastic neutron scattering spectra from molecular dynamics trajectories

Abstract Inelastic neutron scattering (INS) provides a weighted density of phonon modes. Currently, INS spectra can only be interpreted for perfectly crystalline materials because of high computational cost for electronic simulations. INS has the potential to provide detailed morphological information if sufficiently large volumes and appropriate structural variety are simulated. Here, we propose a method that allows direct comparison between INS data with molecular dynamics simulations, a simulation method that is frequently used to simulate semicrystalline/amorphous materials. We illustrate the technique by analyzing spectra of a well-studied conjugated polymer, poly(3-hexylthiophene-2,5-diyl) (P3HT) and conclude that our technique provides improved volume and structural variety, but that the classical force field requires improvement before the morphology can be accurately interpreted.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Capturing Coupled Structural and Electronic Motions During Excited-State Intramolecular Proton Transfer via Computational Multiedge Resonant Inelastic X-ray Scattering

Proton-transfer processes build the foundation of many chemical processes. In Excited State Intramolecular Proton Transfer (ESIPT) processes the proton transfer process is impulsively started through light. Here, in this study, we explore the changes in coupled atomic and electronic motions during and following ESIPT through computational time-resolved Resonant Inelastic X-ray Scattering (trRIXS). Excited-state Ab Initio Molecular Dynamics (AIMD) simulations combined with Time-Dependent Density Functional Theory (TDDFT) calculations were performed for 10-hydroxybenzo[h]quinoline to obtain trRIXS signatures. RIXS at both the nitrogen and oxygen K-edges were computed to resolve the dynamics electronic structure from both the proton donor and acceptor perspective. The results show how RIXS can reveal the local electronic structure, the coupling between different electronic states and how electronic structure and coupling change during the proton transfer process. Additionally, we observe a strong correlation between spectral changes and structural changes during ESIPT.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Using photoelectron spectroscopy to measure resonant inelastic X-ray scattering: a computational investigation

Resonant inelastic X-ray scattering (RIXS) has become an important scientific tool. Nonetheless, conventional high-resolution (few hundred meV or less) RIXS measurements, especially in the soft X-ray range, require low-throughput grating spectrometers, which limits measurement accuracy. Here, the performance of a different method for measuring RIXS, i.e. photoelectron spectrometry for analysis of X-rays (PAX), is computationally investigated. This method transforms the X-ray measurement problem of RIXS to an electron measurement problem, enabling use of high-throughput, compact electron spectrometers. X-rays to be measured are incident on a converter material and the energy distribution of the resultant photoelectrons, the PAX spectrum, is measured with an electron spectrometer. A deconvolution algorithm for analysis of such PAX data is proposed. It is shown that the deconvolution algorithm works well on data recorded with ~0.5 eV resolution. Additional simulations show the potential of PAX for estimation of RIXS features with smaller widths. For simulations using the 3d levels of Ag as a converter material, and with 10 5 simulated detected electrons, it is estimated that features with a few hundred meV width can be accurately estimated in a model RIXS spectrum. For simulations using a sharp Fermi edge to encode RIXS spectra, it is estimated that one can accurately distinguish 100 meV FWHM peaks separated by 45 meV with 10 5 simulated detected electrons that were photoemitted from within 0.4 eV of the Fermi level.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Phonon second harmonic generation in NaBr studied by inelastic neutron scattering and computer simulation

The phenomenon of second harmonic generation (SHG) was found for phonons in anharmonic NaBr by inelastic neutron scattering. The temperature dependence of this phonon SHG was measured from 300 K to 650 K. At 300 K the second harmonic (SH) is seen as a high-energy branch around 33 meV, nearly independent of $\overrightarrow{Q}$. The temperature effective potential (TDEP) method and classical molecular dynamics (MD) simulation with machine learning interatomic potential were able to reproduce the SH, and showed that SHG occurs with the flat transverse optical (TO) phonon branch. A classical model of a nonlinear medium explains the intensity and lifetime of the SH, compared to those of the TO modes. Also successful was a quantum model based on the Heisenberg-Langevin equation for interacting phonons coupled to a thermal bath, which also predicts a spectral distribution of the SH. In conclusion, the measured temperature dependence of the intensity of the second harmonic showed that it follows the Planck distribution of a one-phonon quasiparticle, and not two TO phonons.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Parallel projection—An improved return mapping algorithm for finite element modeling of shape memory alloys

Here, we present a novel finite element analysis of inelastic structures containing Shape Memory Alloys (SMAs). Phenomenological constitutive models for SMAs lead to material nonlinearities, that require substantial computational effort to resolve. Finite element analysis methods, which rely on Gauss quadrature integration schemes, must solve two sets of coupled differential equations: one at the global level and the other at the local, i.e. Gauss point level. In contrast to the conventional return mapping algorithm, which solves these two sets of coupled differential equations separately using a nested Newton procedure, we propose a scheme to solve the local and global differential equations simultaneously. In the process we also derive closed-form expressions used to update the internal/constitutive state variables, and unify the popular closest-point and cutting plane methods with our formulas. Numerical testing indicates that our method allows for larger thermomechanical loading steps and provides increased computational efficiency, over the standard return mapping algorithm.

42 ENGINEERING↗
DOE OSTI2023

Improvements to Contributions from Neutron Inelastic Scattering for Next-Event Estimators in MCNP ® Software

An improvement to the handling of contributions from neutron inelastic scattering to next-event estimators has been implemented in the MCNP6 ® software for release with version 6.3.1. The kinematic equations that govern the outgoing energy of inelastic neutron scattering contributions to next-event estimators have two roots. Historically, the implementation in the MCNP coding has only used the upper root to the quadratic equation. This includes all versions predating version 6.3.1 including all versions of MCNP5 and MCNPX software. However, a review of the neutron next-event estimator physics has shown that this does not reproduce the track-length estimator results at low energies. Several examples are presented that test single neutron inelastic scattering reaction types, Level Scattering (Law 3), Tabulated Energy Angle (Law 61), and Kalbach-Mann distribution (Law 44). The test problems compare the track-length estimator (f4 tally) in the MCNP software, with the existing implementation of the neutron next-event estimator (f5 tally), and the modified changes to the neutron next-event estimator implementation. As the MCNP implementation ignores the lower root, the unmodified neutron next-event estimators will generally underestimate the lower energy contribution. However, a second issue with the Kalbach-Mann distribution (Law 44) implementation allows contributions to backward scattering in the center-of-mass frame that is not kinematically possible, thus overestimating backward scattering contributions. A third issue with the way the MCNP implementation handles floating point comparison for scattering directly ahead or directly backward in the center-of-mass frame generally leads to underestimation (except for backwards scattering for Law 44).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Study of Anharmonicity in Zirconium Hydrides Using Inelastic Neutron Scattering and Ab-Initio Computer Modeling

The anharmonic phonon behavior in zirconium hydrides and deuterides, including ϵ-ZrH 2 , γ-ZrH, and γ-ZrD, has been investigated from aspects of inelastic neutron scattering (INS) and lattice dynamics calculations within the framework of density functional theory (DFT). The harmonic model failed to reproduce the spectral features observed in the experimental data, indicating the existence of anharmonicity in those materials and the necessity of further explanations. Here, we present a detailed study on the anharmonicity in zirconium hydrides/deuterides by exploring the 2D potential energy surface of hydrogen/deuterium atoms and solving the corresponding 2D single-particle Schrödinger equation to obtain the eigenfrequencies, which are then convoluted with the instrument resolution. The convoluted INS spectra qualitatively describe the anharmonic peaks in the experimental INS spectra and demonstrate that the anharmonicity originates from the deviations of hydrogen potentials from quadratic behavior in certain directions; the effects are apparent for the higher-order excited vibrational states, but small for the ground and first excited states.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Soft integrals and soft anomalous dimensions at $\mathrm{N^3LO}$ and beyond

We calculate soft phase-space and loop master integrals for the computation of color-singlet cross sections through N 3 LO in perturbative QCD. Our results are functions of homogeneous transcendental weight and include the first nine terms in the expansion in the dimensional regulator ϵ. We discuss the application of our results to the computation of deeply-inelastic scattering and e + e - annihilation processes. We use these results to compute the perturbative coefficient functions for the Drell-Yan and gluon-fusion Higgs boson production cross sections to higher orders in ϵ through N 3 LO in QCD in the limit where only soft partons are produced on top of the colorless final state. Furthermore, we extract the anomalous dimension of the inclusive threshold soft function and of the N-Jettiness beam and jet functions to N 4 LO in perturbative QCD.

higher-order perturbative calculations↗
DOE OSTI2025

Polarization-dependent resonant inelastic x-ray scattering of β-Ga2O3: An experimental and computational study

The bulk electronic structure of β-Ga2O3 single crystals was investigated using oxygen K-edge x-ray emission and absorption spectroscopy as well as resonant inelastic soft x-ray scattering. Spectra were obtained for different orientations of the crystal planes with respect to the polarization vector of the incident x-ray beam. The spectra are analyzed with calculations based on density functional theory and using the Bethe–Salpeter equation in the ocean code to take the core–hole interaction into account. These calculations correctly capture all the main features in the experimental spectra, demonstrating the potential of the approach to predict the electronic properties of similar compounds. We find a pronounced anisotropy as a function of the excitation polarization vector and significant differences in the spectral contributions from the inequivalent oxygen atoms of β-Ga2O3.

Pyatenko, Elizaveta↗
DOE OSTI2021

Comparison of MCNP thermal scattering laws to inelastic neutron scattering data

Neutron transport computations employ thermal scattering laws (TSLs) to account for the influence of material structure and collective excitations at low energy (≲ 10 eV). Recent work has been undertaken to improve TSL accuracy and expand the number of available materials in the Evaluated Nuclear Data File (ENDF). However, comparison of TSL to measured data is often sparse. In this work, we present a method to evaluate TSL accuracy in the Monte Carlo N-Particle (MCNP) code by simulating inelastic neutron scattering (INS) measurement and comparing to data. The method is briefly described with example application to the recently released ENDF-VIII TSL library for polyethylene as well as a cryogenic PE TSL we are separately developing.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Large change of interlayer vibrational coupling with stacking in Mo 1-x W x Te 2

Stacking variations in quasi-two-dimensional materials can have an important influence on material properties, such as changing the topology of the band structure. Unfortunately, the weakness of van der Waals (vdW) interactions makes it difficult to compute the stacking dependence of properties, and even in a material as simple as graphite the stacking energetics remain unclear. Mo 1 - x W x Te 2 is a material in which three differently stacked phases are conveniently accessible by temperature changes: 1 T ' , T d * , and the reported Weyl semimetal phase T d . Additionally, the transitions proceed via layer sliding, and the corresponding interlayer shear mode (ISM) is relevant not just for the stacking energetics but also for understanding the relationship between Weyl physics and structural changes. However, the interlayer interactions of Mo 1 - x W x Te 2 are not well understood, with wide variation in computed properties. We report inelastic neutron scattering of the ISM in a Mo 0.91 W 0.09 Te 2 crystal. The ISM energies are generally consistent with the linear chain model, as expected given the weak interlayer interaction, though there are some discrepancies from predicted intensities. However, the interlayer force constants K x in the T d * and 1 T ' phases are substantially weaker than that of T d at 75(3) and 83(3)%, respectively. Considering that the relative positioning of atoms in neighboring layers is approximately the same regardless of overall stacking, our results suggest that longer-range influences, such as stacking-induced electronic band-structure changes, may be responsible for the substantial change in the interlayer vibrational coupling and thus the C 55 elastic constant. These findings should elucidate the stacking energetics of Mo 1 - x W x Te 2 and other vdW layered materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Back-to-Back Inclusive Dijets in Deep Inelastic Scattering at Small x : Complete NLO Results and Predictions

We compute the back-to-back dijet cross section in deep inelastic scattering at small x to next-to-leading order (NLO) in the color glass condensate effective field theory. Our result can be factorized into a convolution of the Weizsäcker-Williams gluon transverse-momentum-dependent distribution function (WW gluon TMD) with a universal soft factor and an NLO coefficient function. The soft factor includes both double and single logarithms in the ratio of the relative transverse momentum P ⊥ of the dijet pair to the dijet momentum imbalance q ⊥ ; its renormalization group (RG) evolution is resummed into the Sudakov factor. Likewise, the WW TMD obeys a nonlinear RG equation in x that is kinematically constrained to satisfy both the lifetime and rapidity ordering of the projectile. Exact analytical expressions are obtained for the NLO coefficient function of transversely and longitudinally polarized photons. Our results allow for the first quantitative separation of the dynamics of Sudakov suppression from that of gluon saturation. They can be extended to other final states and provide a framework for precision tests of novel QCD many-body dynamics at the Electron-Ion Collider. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Scalable quantum simulations of scattering in scalar field theory on 120 qubits

Simulations of collisions of fundamental particles on a quantum computer are expected to have an exponential advantage over classical methods and promise to enhance searches for new physics. Furthermore, scattering in scalar field theory has been shown to be bounded-error quantum polynomial time (BQP) complete, making it a representative problem for which quantum computation is efficient. As a step toward large-scale quantum simulations of collision processes, scattering of wave packets in one-dimensional scalar field theory is simulated using 120 qubits of IBM’s Heron superconducting quantum computer ibm_fez. Variational circuits compressing vacuum preparation, wave packet initialization, and time evolution are determined using classical resources. By leveraging physical properties of states in the theory, such as symmetries and locality, the variational quantum algorithm constructs scalable circuits that can be used to simulate arbitrarily large system sizes. A new strategy is introduced to mitigate errors in quantum simulations, which enables the extraction of meaningful results from circuits with up to 4924 two-qubit gates and two-qubit gate depths of 103. The effect of interactions is clearly seen, and is found to be in agreement with classical matrix product state simulations. Finally, the developments that will be necessary to simulate high-energy inelastic collisions on a quantum computer are discussed.

quantum circuits↗
DOE OSTI2021

Dijet impact factor in DIS at next-to-leading order in the Color Glass Condensate

In this work, we compute the next-to-leading order impact factor for inclusive dijet production in deeply inelastic electron-nucleus scattering at small x Bj . Our computation, performed in the framework of the Color Glass Condensate effective field theory, includes all real and virtual contributions in the gluon shock wave background of all-twist lightlike Wilson line correlators. We demonstrate explicitly that the rapidity evolution of these correlators, to leading logarithmic accuracy, is described by the JIMWLK Hamiltonian. When combined with the next-to-leading order JIMWLK Hamiltonian, our results for the impact factor improve the accuracy of the inclusive dijet cross-section to \( \mathcal{O} \) ( \( {\alpha}_s^2 \) ln( x f /x Bj )), where x f is a rapidity factorization scale. These results are an essential ingredient in assessing the discovery potential of inclusive dijets to uncover the physics of gluon saturation at the Electron-Ion Collider.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Explainable AI classification for parton density theory

Quantitatively connecting properties of parton distribution functions (PDFs, or parton densities) to the theoretical assumptions made within the QCD analyses which produce them has been a longstanding problem in HEP phenomenology. To confront this challenge, we introduce an ML-based explainability framework, XAI4PDF, to classify PDFs by parton flavor or underlying theoretical model using ResNet-like neural networks (NNs). By leveraging the differentiable nature of ResNet models, this approach deploys guided backpropagation to dissect relevant features of fitted PDFs, identifying x-dependent signatures of PDFs important to the ML model classifications. By applying our framework, we are able to sort PDFs according to the analysis which produced them while constructing quantitative, human-readable maps locating the x regions most affected by the internal theory assumptions going into each analysis. This technique expands the toolkit available to PDF analysis and adjacent particle phenomenology while pointing to promising generalizations.

Artificial Intelligence↗
DOE OSTI2021

First-principles characterization of the magnetic properties of Cu 2 (OH) 3 Br

Low-dimensional spin-1/2 transition metal oxides and oxyhalides continue to be at the forefront of research investigating nonclassical phases such as quantum spin liquids. In this work, we examine the magnetic properties of the oxyhalide Cu 2 (OH) 3 Br in the botallackite structure using first-principles density functional theory, linear spin-wave theory, and exact diagonalization calculations. This quasi-two-dimensional system consists of Cu 2+ S = 1/2 moments arranged on a distorted triangular lattice. Our exact diagonalization calculations, which rely on a first-principles-based magnetic model, generate spectral functions consistent with inelastic neutron scattering data. By performing computational experiments to disentangle the chemical and steric effects of the halide ions, we find that the dominant effect of the halogen ions is steric in the Cu 2 (OH) 3 X series of compounds.

36 MATERIALS SCIENCE↗