Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Collisional excitation of C2H− by H2: New interaction potential and scattering calculations

Interstellar anions play an important role in astrochemistry as being tracers of the physical and chemical conditions in cold molecular clouds and circumstellar gas. The local thermodynamic equilibrium is generally not fulfilled in media where anions are detected and radiative and collisional data are required to model the observed lines. The C2H− anion has not yet been detected in the interstellar medium; however, collisional data could be used for non-LTE models that would help in identifying the most intense lines. For this purpose, we have computed the first 4D potential energy surface (PES) of the C2H−–H2 complex using an explicitly correlated coupled-cluster approach. The PES is characterized by a single deep minimum with a well-depth of 924.96 cm−1. From this interaction potential, we derived excitation cross sections and rate coefficients of C2H− induced by collisions with para- and ortho-H2. The results obtained for collisions with para-H2 are compared to previous calculations performed using a 2D-PES obtained from an average over H2 rotations.

Chemistry↗
DOE OSTI2023

Accelerating and stabilizing the convergence of vibrational self-consistent field calculations via the direct inversion of the iterative subspace (vDIIS) algorithm

The vibrational self-consistent field (VSCF) method yields anharmonic states and spectra for molecular vibrations, and it serves as the starting point for more sophisticated correlated-vibration methods. Convergence of the iterative, non-linear optimization in VSCF calculations can be erratic or altogether unsuccessful, particularly for chemical systems involving low-frequency motions. In this work, a vibrational formulation of the Direct Inversion of the Iterative Subspace method of Pulay is presented and investigated. This formulation accounts for distinct attributes of the vibrational and electronic cases, including the expansion of each single-mode vibrational wavefunction in its own basis set. The resulting Direct Inversion of the Iterative Subspace method is shown to substantially accelerate VSCF convergence in all convergent cases as well as rectify many cases where Roothaan-based methods fail. Performance across systems ranging from small, rigid molecules to weakly bound molecular clusters is investigated in this analysis.

Chemistry↗
DOE OSTI2023

Ab initio calculations of low-energy quasiparticle lifetimes in bilayer graphene

Motivated by recent experimental results we calculate from first-principles the lifetime of low-energy quasiparticles in bilayer graphene (BLG). Here, we take into account the scattering rate arising from electron-electron interactions within the GW approximation for the electron self-energy and consider several p-type doping levels ranging from 0 to ρ ≈ 2.4 × 10 12 holes/cm 2 . In the undoped case we find that the average inverse lifetime scales linearly with energy away from the charge neutrality point, with values in good agreement with experiments. The decay rate is approximately three times larger than in monolayer graphene, a consequence of the enhanced screening in BLG. In the doped case, the dependence of the inverse lifetime on quasiparticle energy acquires a non-linear component due to the opening of an additional decay channel mediated by acoustic plasmons.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Dynamic and relativistic effects on Pt–Pt indirect spin–spin coupling in aqueous solution studied by ab initio molecular dynamics and two- vs four-component density functional NMR calculations

Treating 195Pt nuclear magnetic resonance parameters in solution remains a considerable challenge from a quantum chemistry point of view, requiring a high level of theory that simultaneously takes into account the relativistic effects, the dynamic treatment of the solvent–solute system, and the dynamic electron correlation. A combination of Car-Parrinello molecular dynamics (CPMD) and relativistic calculations based on two-component zeroth order regular approximation spin–orbit Kohn–Sham (2c-ZKS) and four-component Dirac–Kohn–Sham (4c-DKS) Hamiltonians is performed to address the solvent effect (water) on the conformational changes and JPtPt1 coupling. A series of bridged PtIII dinuclear complexes [L1–Pt2(NH3)4(Am)2–L2]n+ (Am = α–pyrrolidonate and pivalamidate; L = H2O, Cl−, and Br−) are studied. The computed Pt–Pt coupling is strongly dependent on the conformational dynamics of the complexes, which, in turn, is correlated with the trans influence among axial ligands and with the angle N–C–O from the bridging ligands. The J-coupling is decomposed in terms of dynamic contributions. The decomposition reveals that the vibrational and explicit solvation contributions reduce JPtPt1 of diaquo complexes (L1 = L2 = H2O) in comparison to the static gas-phase magnitude, whereas the implicit solvation and bulk contributions correspond to an increase in JPtPt1 in dihalo (L1 = L2 = X−) and aquahalo (L1 = H2O; L2 = X−) complexes. Relativistic treatment combined with CPMD shows that the 2c-ZKS Hamiltonian performs as well as 4c-DKS for the JPtPt1 coupling.

Chemistry↗
DOE OSTI2024

Reconciling calculations and measurements of inverse bremsstrahlung absorption

It was recently shown that the use of Coulomb logarithms appropriate for bremsstrahlung radiation (rather than transport processes) along with corrections for the Langdon effect and ion screening reproduced measurements of collisional absorption in well-characterized underdense plasmas [D. Turnbull et al., Phys. Rev. Lett. 130, 145103 (2023)]. However, it was recognized at the time that the use of the standard absorption-reduction factor from Langdon's seminal paper was inconsistent with the use of Coulomb logarithms that are thermally averaged over a Maxwellian electron-velocity distribution function. A more accurate approach would be to average over the expected super-Gaussian distribution function while accounting for the Gaunt factor's velocity dependence, which somewhat mitigates the Langdon effect; however, at that time, this theory matched the data less well. This conflict is now eliminated with the additional insight that the ionization state of our mid-Z ion species (when present) was lower than had been assumed, as evidenced by the Thomson-scattering data and time-dependent Cretin simulations. We are now able to show that an improved treatment of the Langdon effect provides the best match to data. Otherwise, the prior conclusions remain unchanged. We also show an example of the substantial expected impact to the absorption rate in calculations of indirect-drive hohlraums.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

PyLRO: A Python calculator for analyzing long-range structural order

We present PyLRO, an open-source Python calculator designed to detect, quantify, and display long-range order in periodic structures. The program’s design methodology, workflow, and approach to order quantification are described and demonstrated using a simple toy model. Additionally, we apply PyLRO to a series of metastable AlPO 4 structural intermediates from a prior high-pressure study, demonstrating how to compute and visualize structural order in all directions on a Miller sphere. We further highlight the program’s capabilities through a high-throughput analysis of structural patterns in the pressure-induced amorphization of AlPO 4 , revealing atomistic insights into specific energy regions of massive amorphous structures. These results suggest that PyLRO can be a valuable tool for investigating crystal–amorphous transition in materials research.

36 MATERIALS SCIENCE↗
DOE OSTI2025

GPU acceleration of hybrid functional calculations in the SPARC electronic structure code

We present a Graphics Processing Unit (GPU)-accelerated version of the real-space SPARC electronic structure code for performing hybrid functional calculations in generalized Kohn–Sham density functional theory. In particular, we develop a batch variant of the recently formulated Kronecker product-based linear solver for the simultaneous solution of multiple linear systems. We then develop a modular, math kernel based implementation for hybrid functionals on NVIDIA architectures, where computationally intensive operations are offloaded to the GPUs, while the remaining workload is handled by the central processing units (CPUs). Considering bulk and slab examples, we demonstrate that GPUs enable up to 8× speedup in node-hours and 80× in core-hours compared to CPU-only execution, reducing the time to solution on V100 GPUs to around 300 s for a metallic system with over 6000 electrons, and significantly reducing the computational resources required for a given wall time.

Kohn-Sham density functional theory↗
DOE OSTI2026

pyRMG: A framework for high-throughput, large-cell DFT calculations on supercomputers

Exascale computing delivers the raw power to simulate ever larger and more chemically realistic systems, but realizing this potential requires codes that can efficiently use thousands of processors. Our real-space multigrid (RMG) density functional theory (DFT) code’s grid-decomposition approach scales nearly linearly with the number of graphics processing units (GPUs), even for simulations exceeding thousands of atoms. This scalability makes RMG a compelling tool for high-throughput DFT studies of materials that would otherwise be bottlenecked in other codes (for example, by global fast Fourier transforms in plane-wave DFT). However, the limited workflow infrastructure for RMG has thus far constrained its adoption to a small user community. In this work, we present pyRMG, a Python package designed to streamline the setup and execution of RMG DFT calculations. Built on the pymatgen and ASE (Atomic Simulation Environment) computational materials science Python packages, pyRMG automates input generation and convergence checking, and it integrates with modern job schedulers (e.g., Flux) on leadership-class platforms such as Frontier and Perlmutter. Here, we demonstrate pyRMG for a high-throughput study of strain effects in 2D 2L-Bi 2 Se 3 /2L-NbSe 2 heterostructures, which offers chemical insights into this system and shows that RMG-based workflows can converge with limited user intervention.

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

Calculating shock Hugoniot and isentropes using multiphase equation of state tables and application to shock and release of diamond ablators in inertial confinement fusion implosions

Advances in shock and ramp compression techniques now allow experimental access to unprecedented extreme conditions of pressure and temperature, providing a means to test theoretical models. Here, we describe a simple methodology to compute multi-phase shock Hugoniot and isentropes using multiphase equation of state tables. We treat explicitly the phase coexistence along the phase boundary to reveal the evolution of the sample as it undergoes the phase transformation in adiabatic conditions. We illustrate the method by calculating the predicted shock and shock-and-release behavior of diamond at conditions relevant for the initial stage of inertial confinement fusion implosions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

SHarmonic: A fast and accurate implementation of spherical harmonics for electronic-structure calculations

The authors present SHarmonic, a new implementation of the spherical harmonics targeted for electronic-structure calculations. Their approach is to use explicit formulas for the harmonics written in terms of normalized Cartesian coordinates. This approach results in a code that is as precise as other implementations while being at least one order of magnitude more computationally efficient. The library can run on graphics processing units as well, achieving an additional order of magnitude in execution speed. This new implementation is simple to use and is provided under an open-source license; it can be readily used by other codes to avoid the error-prone and cumbersome implementation of the spherical harmonics.

Mathematics and Computing↗
DOE OSTI2026

Reduced fast-ion transport calculations of m = n = 1 fishbone-like instabilities in MAST-U

Fast-ion transport associated with an m = n = 1 fishbone-like burst in MAST-U discharge 47128 is investigated using a reduced guiding-center-based transport model (ORBIT-Kick) constrained by multi-diagnostic measurements. The two-dimensional beam-emission spectroscopy system provides measurements of the core poloidal mode structure and fluctuation amplitude, while EFIT++ reconstructions constrained by the motional Stark effect diagnostic indicate a flat q-profile with q 0 > 1⁠, indicating the absence of a resonant q = 1 surface and supporting a pressure-driven infernal-mode interpretation. Analytic m = n = 1 displacement profiles consistent with the measured core mode structure and equilibrium constraints are used as the mode structure inputs to ORBIT-Kick. The calculations show that the dominant resonances occur between the mode and co-passing fast ions, producing redistribution localized near the magnetic axis. Synthetic neutron camera signals from TRANSP-Kick recover up to 90% of the experimentally observed neutron deficit at the time of peak mode amplitude, indicating that the measured m = n = 1 mode is a dominant contributor to core fast-ion transport. However, the synthetic neutron signals recover rapidly, whereas the measured neutron emission continues to decrease after the peak amplitude. In conclusion, the remaining discrepancy may arise from contributions not included in the present single-harmonic model, including higher-m and higher-n harmonics, multi-harmonic interactions, and additional transport mechanisms, motivating future diagnostic development and modeling efforts to resolve and incorporate these additional contributions.

Wong, Henry H. [University of California, Los Ange↗
DOE OSTI2021

Neutronics Calculation Advances at Los Alamos: Manhattan Project to Monte Carlo

The history and advances of neutronics calculations at Los Alamos during the Manhattan Project through the present are reviewed. Substantial improvements to neutron diffusion methods and the invention of both the Monte Carlo neutron transport methods in 1947 and deterministic discrete ordinates Sn in 1953 were all made at Los Alamos just after the Manhattan Project. We briefly summarize early simpler and more approximate neutronics methods and then describe the need to better predict neutronics behavior through consideration of theoretical equations, models and algorithms, experimental measurements, and available computing capabilities and their limitations. This paper briefly covers key advances in deterministic methods during the Manhattan Project. These capabilities, coupled with increasing postwar defense needs and the invention of electronic computing with the Electronic Numeric Integrator and Computer, known as ENIAC, and the Mathematical Analyzer Numerical Integrator and Automatic Computer Model, known as MANIAC, led to the creation of Monte Carlo and deterministic discrete ordinates neutronics transport methods. We note the important role that the scientific comradery between the Los Alamos scientists played in the process. This paper briefly covers the early methods, algorithms, computers, and electronic and women pioneers that enabled Monte Carlo to spread to all areas of science. We focus heavily on these early developments and the subsequent creation of the MCNP® code, advances in its associated nuclear data, and its applications to problems of national defense at Los Alamos.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Multigroup Constant Calculation with Static α-Eigenvalue Monte Carlo for Time-Dependent Neutron Transport Simulations

Herein, multigroup constants for deterministic methods that preserve the time-dependent physics of the neutron transport equations are derived. Alternative multigroup constant weighting spectra are discussed: (1) the fundamental k-eigenfunction, (2) the fundamental α-eigenfunction, and (3) a composite of several α-modes. To generate the fundamental α-eigenfunction for calculating the multigroup constants, a static fundamental α-eigenvalue method is implemented into the open source Monte Carlo code OpenMC. Several static and kinetic problems are devised to verify the implementations and to investigate the relative performance of the alternative multigroup constant weighting spectra. Results emphasize that as a multigroup constant weighting spectrum, the fundamental α-eigenfunction offers physical characteristics that make it advantageous (in producing accurate solutions) over the typically used fundamental k-eigenfunction.

42 ENGINEERING↗
DOE OSTI2023

Deep Learning for Multigroup Cross-Section Representation in Two-Step Core Calculations

Here we investigate using deep learning, a type of machine-learning algorithm employing multiple layers of artificial neurons, for the mathematical representation of multigroup cross sections for use in the Griffin reactor multiphysics code for two-step deterministic neutronics calculations. A three-dimensional fuel element typical of a high-temperature gas reactor as well as a two-dimensional sodium-cooled fast reactor lattice are modeled using the Serpent Monte Carlo code, and multigroup macroscopic cross sections are generated for various state parameters to produce a training data set and a separate validation data set. A fully connected, feedforward neural network is trained using the open-source PyTorch machine-learning framework, and its accuracy is compared against the standard piecewise linear interpolation model. Additionally, we provide in this work a generic technique for propagating the cross-section model errors up to the k eff using sensitivity coefficients with the first-order uncertainty propagation rule. Quantifying the eigenvalue error due to the cross-section regression errors is especially practical for appropriately selecting the mathematical representation of the cross sections. We demonstrate that the artificial neural network model produces lower errors and therefore enables better accuracy relative to the piecewise linear model when the cross sections exhibit nonlinear dependencies; especially when a coarse grid is employed, where the errors can be halved by the artificial neural network. However, for linearly dependent multigroup cross sections as found for the sodium-cooled fast reactor case, a simpler linear regression outperforms deeper networks.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

SCALE Shielding Calculations for Advanced Reactor Accident Scenarios

The study results presented in this paper demonstrate the capabilities of the SCALE computer code for non–light water reactor (non-LWR) radiation source term and shielding calculations in support of the U.S. Nuclear Regulatory Commission non-LWR fuel cycle demonstration project. Representative non-LWR types, including the sodium-cooled fast reactor (SFR), the molten salt reactor (MSR), and the heat pipe microreactor (HPMR), were analyzed to evaluate dose rates associated with postulated accident scenarios (SFR and MSR), reactor operation (MSR and HPMR), and fuel self-protecting characteristics (HPMR). New features were implemented in SCALE depletion codes to better simulate MSR operation.

MSR↗
DOE OSTI2025

Development and validation of a humidified CAPS-PMSSA with an improved methodology to calculate the truncation correction factor

Absorbing aerosols play an important role in the Earth’s radiation budget, and current uncertainties associated with their radiative impacts remain substantial. An element of this uncertainty derives from insufficient understanding of the extent to which water, as a coating on absorbing particles, enhances the absorption. Here, we describe and characterize the performance of a cavity-attenuated phase shift single scatter albedo spectrometer (CAPS-PMSSA) system that has been modified to enable concurrent aerosol extinction and scattering measurements, from which aerosol absorption is derived, at relative humidities up to ~90%. The modified instrument is referred to as the humidified CAPS (HCAPS). Additionally, we develop and validate an update to a method to calculate the truncation correction factor, used to correct the scattering measurements, that can account for the effects of particle growth and compositional changes resulting from water uptake. The method can be applied to experiments that use either size-selected or polydisperse aerosol samples of either known or unknown composition. Overall, we show that measurements made with the HCAPS, when processed using the updated truncation correction method, enable accurate measurement of the influence of water uptake on light absorption.

Dal Porto, Rachael [University of California, Davi↗
DOE OSTI2020

Benchmark calculations of electron impact electronic excitation of the hydrogen molecule

We present benchmark integrated and differential cross-sections for electron collisions with H 2 using two different theoretical approaches, namely, the R-matrix and molecular convergent close-coupling. This is similar to comparative studies conducted on electron–atom collisions for H, He and Mg. Electron impact excitation to the ${\mathrm{b}}^{3}{{\Sigma}}_{\mathrm{u}}^{+}$, ${\mathrm{a}}^{3}{{\Sigma}}_{\mathrm{g}}^{+}$, ${\mathrm{B}}^{1}{{\Sigma}}_{\mathrm{u}}^{+}$, c3Πu, ${\mathrm{E}\mathrm{F}}^{1}{{\Sigma}}_{\mathrm{g}}^{+}$, C1Πu, ${\mathrm{e}}^{3}{{\Sigma}}_{\mathrm{u}}^{+}$, ${\mathrm{h}}^{3}{{\Sigma}}_{\mathrm{g}}^{+}$, ${\mathrm{B}}^{\prime 1}{{\Sigma}}_{\mathrm{u}}^{+}$ and d3Πu excited electronic states are considered. Calculations are presented in both the fixed nuclei and adiabatic nuclei approximations, where the latter is shown only for the ${\mathrm{b}}^{3}{{\Sigma}}_{\mathrm{u}}^{+}$ state. Good agreement is found for all transitions presented. Where available, we compare with existing experimental and recommended data.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

First-principles calculation of Hubbard U for Terbium metal under high pressure

Abstract Using density functional theory (DFT) and linear response approaches, we compute the on-site Hubbard interaction U of elemental Terbium (Tb) metal in the pressure range ∼ 0–65 GPa. The resulting first-principles U values with experimental crystal structures enable us to examine the magnetic properties of Tb using a DFT+U method. The lowest-energy magnetic states in our calculations for different high-pressure Tb phases—including hcp, α -Sm, and dhcp—are found to be compatible with the corresponding magnetic ordering vectors reported in experiments. The result shows that the inclusion of Hubbard U substantially improves the accuracy and efficiency in modeling correlated rare-earth materials. Our study also provides the necessary U information for other quantum many-body techniques to study Tb under extreme pressure conditions.

36 MATERIALS SCIENCE↗