Search NASA⌕ Search

SEARCH · Search NASA

Results for “Higher order 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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Measurements of Lund subjet multiplicities in 13 TeV proton-proton collisions with the ATLAS detector

This Letter presents a differential cross-section measurement of Lund subjet multiplicities, suitable for testing current and future parton shower Monte Carlo algorithms. This measurement is made in dijet events in 140 fb -1 of $\sqrt{s}$ =13 TeV proton–proton collision data collected with the ATLAS detector at CERN's Large Hadron Collider. The data are unfolded to account for acceptance and detector-related effects, and are then compared with several Monte Carlo models and to recent resummed analytical calculations. The experimental precision achieved in the measurement allows tests of higher-order effects in QCD predictions. Most predictions fail to accurately describe the measured data, particularly at large values of jet transverse momentum accessible at the Large Hadron Collider, indicating the measurement's utility as an input to future parton shower developments and other studies probing fundamental properties of QCD and the production of hadronic final states up to the TeV-scale.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Prospects for strong coupling measurement at hadron colliders using soft-drop jet mass

We compute the soft-drop jet-mass distribution from pp collisions to NNLL accuracy while including nonperturbative corrections through a field-theory based formalism. Using these calculations, we assess the theoretical uncertainties on an αs precision measurement due to higher order perturbative effects, nonperturbative corrections, and PDF uncertainty. We identify which soft-drop parameters are well-suited for measuring αs, and find that higher-logarithmic resummation has a qualitatively important effect on the shape of the jet-mass distribution. We find that quark jets and gluon jets have similar sensitivity to αs, and emphasize that experimentally distinguishing quark and gluon jets is not required for an αs measurement. We conclude that measuring αs to the 10% level is feasible now, and with improvements in theory a 5% level measurement is possible. Getting down to the 1% level to be competitive with other state-of-the-art measurements will be challenging.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

The influence of thermochemistry on the reactivity of propane, the pentane isomers and $n$-heptane in the low temperature regime

In this study, the influence of thermochemistry on the reactivity of fuels at low temperatures (600–1000 K) is studied here. Specifically, the effect of different sets of thermochemistry on chemical model predictions is explored, where various sets are calculated at different levels of theory in addition to recently updated group additivity values. Experimentally measured ignition delay times for propane, the pentane isomers and n-heptane are simulated using NUIGMech1.2 and replacing the thermochemistry of the low-temperature species with the calculated values. For propane, three different thermochemistry sets were calculated, namely CCSD(T)-F12/TZ-F12//B2PLYP-D3/TZ//B2PLYP-D3/TZ (QM1), CCSD(T)-F12/TZ-F12//B2PLYP-D3/TZ//ωB97X-D/TZ (QM2) and B2PLYP-D3/TZ/ωB97X-D/6-31G*//ωB97X-D/6-31G* (QM3). The QM2 results provide parameters to optimize new group additivity (NGA) values which are used to calculate the fourth set of thermochemistry. The model predictions using these four sets are compared to those using NUIGMech1.2 for propane. As the QM1 and QM2 calculations are expensive, the thermochemistry calculated from the QM3 and NGA calculations are used in the pentane isomer and n-heptane models. For all of the models, it is found that the thermochemistry of the species involved in the low-temperature reaction sequence ($RH$, $\dot{R}$, $RO_2H$, $R\dot{O}_2$, $\dot{Q}OOH$ and $\dot{O}_2QOOH$ species) significantly affect fuel reactivity. The NGA values were developed based on all of these species except $\dot{O}_2QOOH$ radicals. The thermochemistry of $\dot{O}_2QOOH$ species cannot be accurately calculated with the NGA representations due to the importance of non-next-nearest neighbor interactions of –OOH substitution. Further development of the NGA method to capture such interactions is in progress. Overall, the model developed using the NGA thermochemistry shows better agreement with experimental data than the model using thermochemistry from affordable and prominent QM methods, such as QM3. Based on the results presented for propane, the pentane isomers and n-heptane, the thermochemistry calculated using the NGA method can be used to model the oxidation of higher order hydrocarbons at low temperatures.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2021

Reducing negative weights in Monte Carlo event generation with Sherpa

An increase in theoretical precision of Monte Carlo event generators is typically accompanied by an increased need for computational resources. One major obstacle are negative weighted events, which appear in Monte Carlo simulations with higher perturbative accuracy. While they can be handled somewhat easily in fixed-order calculations, they are a major concern for particle level event simulations. In this article, the origin of negative weights in the S-MC@NLO method is reviewed and mechanisms to reduce the negative weight fraction in simulations with the Sherpa event generator are presented, with a focus on V+jets and tt+jets simulations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Measurement of the Sixth-Order Cumulant of Net-Proton Multiplicity Distributions in Au + Au Collisions at s NN = 27 , 54.4, and 200 GeV at RHIC

According to first principle Lattice QCD calculations, the transition from quark-gluon plasma to hadronic matter is a smooth crossover in the region µB ≤ Tc. As a result, higher-order cumulants and their ratios are predicted to be negative, C 6 /C 2 < 0, for example. In this paper, we report the first measurement of the midrapidity net-proton C 6 /C 2 from 27, 54.4 and 200 GeV Au+Au collisions at RHIC. The dependence on collision centrality and kinematic acceptance in (p T, y ) are analyzed. While for 27 and 54.4 GeV collisions the C 6 /C 2 values are close to zero within uncertainties, it is observed that for 200 GeV collisions, the C 6 /C 2 ratio becomes progressively negative from peripheral to central collisions. Transport model calculations without critical dynamics predict values around zero. These observations seem to favor a smooth crossover in the high energy nuclear collisions at RHIC.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Understanding Lithium Local Environments in LiMn 0.5 Ni 0.5 O 2 Cathodes: A DFT-Supported 6 Li Solid-State NMR Study

LiMn 0.5 Ni 0.5 O 2 has been plagued by the structural failure driven by the irreversible phase change upon subsequent cycling and the structural degradation caused by the seemingly unavoidable Li/Ni exchange. This begs for a better understanding of the correlation between lithium local structures, synthesis conditions, and overall electrochemical properties of materials. In this study, we use density functional theory (DFT) to assist experimental 6 Li solid-state magic angle spinning (MAS) nuclear magnetic resonance (NMR) spectroscopy to understand the nature and change in lithium local structures with extended annealing times (168 vs 15 h) synthesized at 900 °C. To correlate the NMR peak changes after different synthesis conditions, we calculated the NMR spectra of representative low-energy model configurations of LiMn 0.5 Ni 0.5 O 2 composition–zigzag, row, and flower structures–as well as the influence of Ni/Mn ordering and Li/Ni exchange on each spectrum using a supercell that is large enough to allow us to consider the effect of different relative configurations of adjacent transition metal (TM) layers. The analysis based on a combination of the calculated NMR and the experimental spectra suggests that most configurations contributing to the experimental spectra are in fact composed of a blend of configurations higher in energy than the classical well-ordered ground-state structures. Further, the extended annealing of LiMn 0.5 Ni 0.5 O 2 slightly enhances Ni/Mn and Li/Mn orderings while reducing the Li/Ni mixing ratio. This study shows how DFT calculations are crucial in providing a better understanding of the experimental characterization data and therefore local environments and domain structures in oxide materials, such as Li-ion battery cathodes.

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

Very massive tracers and higher derivative biases

Most of the upcoming cosmological information will come from analyzing the clustering of the Large Scale Structures (LSS) of the universe through LSS or CMB observations. It is therefore essential to be able to understand their behavior with exquisite precision. The Effective Field Theory of Large Scale Structures (EFTofLSS) provides a consistent framework to make predictions for LSS observables in the mildly non-linear regime. Herein we focus on biased tracers. We argue that in calculations at a given order in the dark matter perturbations, highly biased tracers will underperform because of their larger higher derivative biases. A natural prediction of the EFTofLSS is therefore that by simply adding higher derivative biases, highly massive tracers should perform comparably well. We implement this prediction for the halo-halo and the halo-matter power spectra at one loop, and the halo-halo-halo, halo-halo-matter, and halo-matter-matter bispectra at tree-level, and compare with simulations. We find good agreement with the prediction: at $z$ = 0, for all tracers, we are able to match the power spectra up to $k ≃$ 0.28$h$ Mpc -1 , as well as a small set of about 10 2 bispectra triangles up to $k ≃$ 0.17$h$ Mpc -1 . We also discuss the limitations of our study and some avenues to pursue to further establish these findings.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Hydraulic Response to Thermal Stimulation Efforts at Raft River Based on Stepped Rate Injection Testing

The injection well stimulation project at the Raft River geothermal field tests the effect of long-term cold water injection and high pressure injection on well injectivity, improvements to which could reduce operating costs. The primary data for analysis and interpretation of the injection test are step-rate flow tests run before each new phase of the injection. These tests were analyzed using a combination of standard pump-test analytical solution methods and methods developed expressly for the observed conditions. The stepped rate injection tests, combined with long-term flow and pressure response data suggest that the well is located within a fractured formation of low transmissivity but high storativity. These calculated parameters appeared to increase with pressure during the first injection test and the higher values were reproduced during the second stepped rate test. Calculated transmissivity and storativity are on the order of 4E-5 m cm 2 and 1E-4 m Pa -1 , respectively. The apparent pressure dependence of fitted hydraulic parameters may reflect near-well fracture compliance that increased the effective radius of the wellbore during the first test. While the type curve fit analysis also suggests that the reservoir behaves as a uniformly fractured reservoir with a radial flow regime, the hydraulic parameters indicate that condition may exist only a very limited distance (<10 m) from the well. Longer-term pressure response suggests that flow in the system effectively reaches steady state in a period of less than a day, which may reflect pressure stabilization resulting from pressure-dependent permeability or a region of much higher permeability located with a few meters of the well. The transmissivity estimates obtained from this analysis, converted to approximate fracture density and aperture, provide useful constraints on the distance to which the thermal front may migrate from the well during the cold water injection phase of the stimulation project. We estimate that the cooling front will migrate less than a tenth of a kilometer over an approximately one-year injection period. Here, the effects of that cooling, however, may be substantial, because increases in permeability have maximum effect nearest the well.

cold water injection↗
DOE OSTI2020

First-principles theory of the Dirac semimetal Cd 3 As 2 under Zeeman magnetic field

Time-reversal broken Weyl semimetals have attracted much attention recently, but certain aspects of their behavior, including the evolution of their Fermi surface topology and anomalous Hall conductivity with Fermi-level position, have remained underexplored. A promising route to obtain such materials may be to start with a nonmagnetic Dirac semimetal and break time-reversal symmetry via magnetic doping or magnetic proximity. Here we explore this scenario in the case of the Dirac semimetal Cd 3 As 2 , based on first-principles density-functional calculations and subsequent low-energy modeling of Cd 3 As 2 in the presence of a Zeeman field applied along the symmetry axis. We clarify how each four-fold degenerate Dirac node splits into four Weyl nodes, two with chirality ±1 and two higher-order nodes with chirality ±2. Using a minimal k · p model Hamiltonian whose parameters are fit to the first-principles calculations, we detail the evolution of the Fermi surfaces and their Chern numbers as the Fermi energy is scanned across the region of the Weyl nodes at fixed Zeeman field. We also compute the intrinsic anomalous Hall conductivity as a function of Fermi-level position, finding a characteristic inverted-dome structure. Cd 3 As 2 is especially well suited to such a study because of its high mobility, but the qualitative behavior revealed here should be applicable to other Dirac semimetals as well.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Harmonic Infrared and Raman Spectra in Molecular Environments Using the Polarizable Embedding Model

We present a fully analytic approach to calculate infrared (IR) and Raman spectra of molecules embedded in complex molecular environments modeled using the fragment-based polarizable embedding (PE) model. We provide the theory for the calculation of analytic second-order geometric derivatives of molecular energies and first-order geometric derivatives of electric dipole moments and dipole–dipole polarizabilities within the PE model. The derivatives are implemented using a general open-ended response theory framework, thus allowing for an extension to higher-order derivatives. The embedding-potential parameters used to describe the environment in the PE model are derived through first-principles calculations, thus allowing a wide variety of systems to be modeled, including solvents, proteins, and other large and complex molecular environments. Here, we present proof-of-principle calculations of IR and Raman spectra of acetone in different solvents. This work is an important step toward calculating accurate vibrational spectra of molecules embedded in realistic environments.

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

High temperature stability and transport characteristics of hydrogen in alumina via multiscale computation

Here, the impact of hydrogen charge states on the stability and transport characteristics of hydrogen interstitials in alumina polymorphs is evaluated by multiscale computational methods including density functional theory (DFT), ab initio molecular dynamics (AIMD) and machine learned force fields. Thermodynamic calculations show that the protonic H i +1 interstitial is the most stable defect species for most values of the electronic bandgap in both and amorphous alumina (Al 2 O 3 ). Further, active learned Gaussian approximation potentials (GAP) were developed using AIMD data to study temperature dependent long time proton diffusion in alumina. Diffusivity calculations from GAP-MD simulations are found to be comparable with of the AIMD data, while being ~340 times faster and scalable to larger systems. Comparisons with diffusivity values for other interstitial charge states (H i 0 and H i -1 ) and published experimental literature indicate that H i +1 diffusion is the likely mechanism of hydrogen transport. A good agreement is obtained between H i +1 diffusivity calculated in α-Al 2 O 3 from DFT: 5.05 10 -3 exp(-0.81 eV/k B /T) cm 2 /s and reported experiment: 9.7X10 -4 exp(-0.83 eV/k B /T) cm 2 /s. Computationally and experimentally calculated energy barriers (0.81 and 0.83 eV respectively) only differ by 2.5%. Similarly, the pre-exponential diffusion coefficients only differ by 0.5 orders of magnitude. Moreover, the diffusivity of H i +1 in amorphous Al 2 O 3 in the 1000–2000 K range is calculated to be 2.53X10 -2 exp(-0.89 eV/k B /T), just one order of magnitude higher than the corresponding value in α-Al 2 O 3 . This suggests that local structural disorder does not significantly affect the energy landscape and diffusion behavior of H i +1 in Al 2 O 3 . Overall, these results show promise for the application of alumina polymorphs as hydrogen permeation barriers.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Staged Z-pinch modeling of high and low atomic number liners compressing deuterium targets using parameters of the Z pulsed power facility

The staged Z-pinch is a potentially transformative magneto-inertial fusion energy concept where a high atomic number liner implodes on a deuterium or deuterium-tritium target using multi-MA pulsed current. Over the past several years, this concept was studied experimentally on 1 MA facilities with argon or krypton gas puffs compressing a magnetized deuterium target. Consistent thermonuclear neutron yield of 1010 per shot was measured with krypton liners. In this paper, we investigate the fusion performance of deuterium targets of varying density undergoing compression with low (beryllium) and high atomic number liners (silver, tantalum) using parameters of the Z pulsed power facility. Silver and tantalum liners create strong shocks that preheat the target plasma above 100 eV and pile up liner material at the liner-target interface. The increased mass density at the interface creates strong ram pressure just before the pinch stagnation time. The target plasma is heated to 2–4 keV, in contrast to the < 0.5 kV temperature calculated for the beryllium case. The high atomic number liners produce neutron yield orders of magnitude higher than the yield from the low atomic number liner.

Ruskov, E. (ORCID:0000000336724580)↗
DOE OSTI2021

Minimal two-body quantum absorption refrigerator

Here, we study the phenomenon of absorption refrigeration, where refrigeration is achieved by heating instead of work, in two different setups: a minimal set up based on coupled qubits and two nonlinearly coupled resonators. Considering ZZ interaction between the two qubits, we outline the basic ingredients required to achieve cooling. Using local as well as global master equations, we observe that inclusion of an XX type term in the qubit-qubit coupling is detrimental to cooling. We compare the cooling effect obtained in the qubit case with that of nonlinearly coupled resonators (multilevel system) where the ZZ interaction translates to a Kerr-type nonlinearity. For small to intermediate strengths of nonlinearity, we observe that multilevel quantum systems, for example qutrits, give better cooling effect compared to the qubits. Using Keldysh nonequilibrium Green's function formalism, we go beyond first order sequential tunneling processes and study the effect of higher order processes on refrigeration. We find a reduced cooling effect compared to the master equation calculations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

QCD splitting functions beyond kinematical limits

In this talk, I will present a new decomposition of QCD splitting functions, carried out systematically up to second order in the strong coupling. The core idea is to separate the splitting functions into two components: scalar dipole radiator functions and pure splitting remainders. While these pieces include terms that are formally subleading in the soft or collinear limits, they turn out to be universal and well understood. Unlike conventional approaches, our construction does not rely on any soft or collinear approximations. This framework not only offers a clearer understanding of the structure of higher-order QCD corrections but also paves the way for more robust subtraction schemes in precision calculations.

Knobbe, Max [Fermilab]↗
DOE OSTI2022

High-order accurate solutions of the point kinetics equations with the spectral deferred correction method

Solving initial value problems with higher-order methods can improve the accuracy of the simulation results or the efficiency of the calculation. In this paper, we apply the spectral deferred correction (SDC) method to solve the initial value problem of the point kinetics equations (PKE). SDC is a stable, robust, and efficient high-order time-integration method capable of an arbitrary order of accuracy. For our implementation we show that it is A-stable for orders up to 8 and the order of accuracy is verified for PKE problems with a range of different reactivities. A 5.-order SDC method was then implemented to solve the exact PKE (EPKE) in the Transient Multilevel (TML) method of MPACT. The error from solutions of the EPKE is shown to be negligible. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

A Block-Structured Adaptive Mesh Framework to Solve Radiation Transfer Equation in Irregular Embedded Geometries

Radiation transport arises in various scientific, industrial, and medical fields, and understanding its effect in applications is needed to make accurate predictions, safety assessments and performance optimizations. Solving the Radiation Transport Equation (RTE) is challenging due to its integro-differential nature, which involves both differential and integral terms. The differential term describes the change in radiation intensity due to absorption and emission, while the integral term accounts for scattering. The accurate modeling of radiation is further complicated in many applications due to the complex, irregular geometries. Various methods exist for solving the RTE, including the zonal, Monte Carlo, spherical harmonics, discrete ordinates, and finite volume methods. Traditional mesh-based approaches, which rely on structured or unstructured meshes, struggle with irregular geometries due to: a) the difficulty of conforming structured grids to irregular domains, b) challenges in enforcing boundary conditions correctly, and c) the additional computational cost of unstructured mesh methods. This work presents a second-order accurate method for solving the RTE in irregular geometries. The radiation intensity is discretized using the finite-volume method in both spatial and angular directions on regular Cartesian grid blocks. Leveraging the block-structured adaptive mesh refinement (AMR) framework provided by AMReX, our method refines the grid locally to reduce spatial discretization error, ensuring a converged numerical solution while minimizing computational costs elsewhere. A two-stage deferred correction approach is employed: First, a first-order discretization on grid blocks is solved using an algebraic multigrid method in HYPRE. Second, a correction term is applied explicitly to achieve second-order accuracy. The correction term is calculated by approximating the radiation flux on cell faces using a Total Variation Diminishing (TVD) scheme. This approach ensures quick convergence of the multigrid method while preserving higher-order accuracy of the numerical solution. Irregular geometries are resolved as embedded boundaries (EB), resulting in both cut cells and regular cells. In cut cells, we modify the fluxes using face fractions and incorporate additional contributions from EB boundary conditions. To ensure higher-order convergence near the EB interface, the correction term is modified by interpolating the radiation intensity to fictitious ghost points. The implementation takes advantage of modern supercomputers by leveraging AMReX’sMPI/X parallelization strategy where X can be MPI or a GPU accelerator including CUDA, HIP and DPC++. We validate our solver using classical test cases, both with and without EB, demonstrating accuracy and efficiency. Additionally, we analyze the impact of adaptive mesh refinement on solution accuracy and computational cost, highlighting the advantages of our approach for high-resolution radiation transport simulations.

computational fluid dynamics (CFD)↗
DOE OSTI2024

Anion photoelectron spectroscopy and chemical bonding of ThS 2 – and ThSO –

Anion photoelectron spectra of ThSO – and ThS 2 – were recorded using the third (355 nm) harmonic of an Nd-YAG laser; these provided the measured vertical detachment energies of each anion. The experiments are supported by extensive coupled cluster calculations on ThSO, ThSO – , ThS 2 , and ThS 2 – , as well as the oxygen congeners ThO 2 and ThO 2 – . The ab initio calculations, which included complete basis set extrapolations, spin–orbit effects using four-component coupled cluster, and higher-order correlation contributions through CCSDT(Q), yielded an adiabatic electron affinity for ThO 2 – that was within 0.02 eV of the previously determined experimental value. The singly occupied molecular orbital (SOMO) in all three anions corresponds primarily to the 7s orbital on Th. Successive substitution of S for each O in ThO 2 leads to larger electron affinities and smaller bond angles in the neutral molecules, but larger angles in the anions. As demonstrated by Franck–Condon simulations of the spectra using the CCSD(T) spectroscopic constants, substitution of O by S significantly complicates the resulting detachment spectra due to the lower vibrational frequencies in the sulfur species. Altogether the calculated vertical detachment energies are in very good agreement with the experiment. In addition to the adiabatic electron affinities of each species, atomization energies and heats of formation have also been determined via the FPD approach with expected uncertainties of 1–2 kcal/mol.

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

Modeling phase transformations in Mn-rich disordered rocksalt cathodes with machine-learning interatomic potentials

Mn-rich disordered rocksalt (DRX) cathode materials exhibit a phase transformation from a disordered to a partially disordered spinel-like structure (δ-phase) during electrochemical cycling. Here, in this computational study, we use charge-informed molecular dynamics with a fine-tuned CHGNet foundation potential to investigate the phase transformation in LixMn 0.8 Ti 0.1 O 1.9 F 0.1 . Our results indicate that transition metal migration occurs and reorders to form the spinel-like ordering in an FCC anion framework. The transformed structure contains a higher concentration of nontransition metal (0-TM) face-sharing channels, which are known to improve Li transport kinetics. Analysis of the Mn valence distribution suggests that the appearance of tetrahedral Mn 2+ is a consequence of spinel-like ordering, rather than the trigger for cation migration as previously suggested. Calculated equilibrium intercalation voltage profiles demonstrate that the δ-phase, unlike the ordered spinel, exhibits solid-solution signatures at low voltage. A higher Li capacity is obtained than in the DRX phase. This study provides atomic insights into solid-state phase transformation and its relation to experimental electrochemistry, highlighting the potential of machine-learning interatomic potentials for understanding complex oxide materials.

Zhong, Peichen [University of California, Berkeley↗