Search NASA⌕ Search

SEARCH · Search NASA

Results for “Basis function”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,153 records · Page 64
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Indirect Validation of Probe Speed Data on Arterial Corridors

This study aimed to estimate the accuracy of probe speed data on arterial corridors on the basis of roadway geometric attributes and functional classification. It was assumed that functional class (medium and low) along with other road characteristics (such as weighted average of the annual average daily traffic, average signal density, average access point density, and average speed) were available as correlation factors to estimate the accuracy of probe traffic data. This study tested these factors as predictors of the fidelity of probe traffic data by using the results of an extensive validation exercise. This study showed strong correlations between these geometric attributes and the accuracy of probe data when they were assessed by using average absolute speed error. Linear models were regressed to existing data to estimate appropriate models for medium- and low-type arterial corridors. The proposed models for medium- and low-type arterials were validated further on the basis of the results of a slowdown analysis. These models can be used to predict the accuracy of probe data indirectly in medium and low types of arterial corridors.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2013

Functional Bounding Content Envelope for Gamma and Neutron Emitting Isotopes in the Model 9977 Packaging

The contents approved for shipment in a Type B radioactive material transportation package have historically been descriptions of discrete items, or groupings of well-defined similar items, in the package safety basis documentation. The need for a comprehensive functional content envelope of both gamma and neutron emitting nuclides compliant with regulatory limits has become necessary as the DOE complex requires shipments of unique mixtures of radionuclides and impurities. Recent publications have presented a calculational model and a corresponding content envelope intended to be compliant with federal regulatory external radiation limits as well as design decay heat limits based on the Model 9977 Packaging. The methodology used to develop this content envelope consisted of determining the external radiation dose rates based on one gram of a given isotope combined, in the case of actinides, with various levels of light element impurities and determining the allowable mass to the regulatory limits based on these dose rates. The method of ratioing from these calculations to the masses that meet the regulatory limits fails in the case of some combinations of actinides and light element impurities because of the effect of subcritical multiplication, resulting in non-conservative and non-compliant dose rates in some cases. This is particularly true for many actinides combined with beryllium, boron, fluorine, lithium, and sodium. In this study, the source was modeled as a sphere with the appropriate dimensions based on the actinide or gamma source density and placed at the bottom of the containment vessel for maximum conservatism. The results presented in this report rectify the above-noted deficiencies by adjusting the previously determined mass limits so that the resulting dose rates are compliant with regulatory limits. Some amount of iteration was required in some instances to adjust the masses for dose compliance. In some instances the masses were increased since they were too conservative. In addition, the method of ratioing the masses when applied to several gamma emitting isotopes also resulted in either non-conservative or overly conservative data. The mass imits presented here are for bare sources and are complaint with the design decay heat limit of 19W. The masses are based on package surface dose rates (the limiting dose rate) that fall in the range 185-195 mrem/h, thus giving an additional margin of between 7.5% and 2.5% to the regulatory limit of 200 mrem/h. This report presents a revised set of mass limits (i.e., content envelope) for several neutron emitting actinides with varying levels of light element impurities compliant with both external radiation and design decay heat limits. In addition, revised limits for gamma sources are also presented. The neutron emitter mass limits are in Tables 1 through 9 and the gamma emitter mass limits are in Table 10. The revised content envelope covering a wide range of materials present in the DOE complex, was developed for the Model 9977 Packaging, providing conservative limits for other commonly-used Type B Packagings such as the Models 9975 and 9978. It should be noted that for the Model 9977 package, shipment is also possible with three shielded containers – lead or tungsten for gamma sources and polyethylene for neutron emitters. The design decay heat limits for these are as follows: 19W for tungsten, 6W for lead and 3W for polyethylene. The limits presented in this report would need to be checked against the lower decay heat limits for the lead and polyethylene and adjusted as required for compliance. Criticality safety is not factored into these limits and should be evaluated on a case by case basis.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2012

Description of deuteron-nucleus scattering in terms of a microscopic optical potential based on Skyrme forces

A model of a deuteron-nucleus (dA) optical potential is developed. The real and imaginary parts of this potential were determined from the calculation of the mass operator in the respective Green's function by perturbation theory on the basis of effective nuclear-density-dependent nucleon-nucleon (NN) Skyrme forces and with allowance for the rearrangement potential. Radial and energy dependences of various components of the dA potential were studied in the calculations with various versions of Skyrme forces. The presentmodel was applied to analyzing differential cross sections for elastic deuteron scattering on nuclei over a broad range of mass numbers at various energies and total cross sections for dA reactions. A reasonable description of respective experimental data was obtained.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2016

Raynal–Revai coefficients for a general kinematic rotation

In a three-body system, transitions between different sets of normalized Jacobi coordinates are described as general kinematic transformations that include an orthogonal or a pseudoorthogonal rotation. For such rotations, the Raynal–Revai coefficients execute a unitary transformation between three-body hyperspherical functions. Recurrence relations that make it possible to calculate the Raynal–Revai coefficients for arbitrary angular momenta are derived on the basis of linearized representations of products of hyperspherical functions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Improvements in the orbitalwise scaling down of Perdew–Zunger self-interaction correction in many-electron regions

The Perdew-Zunger (PZ) method offers a way to remove the self-interaction (SI) error from density functional approximations on an orbital by orbital basis. The PZ method provides significant improvements for the properties such as barrier heights or dissociation energies but results in over-correcting the properties well described by SI-uncorrected semi-local functional. One cure to rectify the over-correcting tendency is to scale down the magnitude of SI-correction of each orbital in the many electron region. We have implemented the orbitalwise scaled down SI-correction (OSIC) scheme of Vydrov et al. [J. Chem. Phys. 124, 094108 (2006)] using the Fermi-L owdin SI-correction method. After validating the OSIC implementation with previously reported OSICLSDA results, we examine its performance with the most successful non-empirical SCAN meta-GGA functional. Using different forms of scaling factors to identify one-electron regions, we assess the performance of OSIC-SCAN for a wide range of properties: total energies, ionization potentials and electron affinities for atoms, atomization energies, dissociation and reaction energies, and reaction barrier heights of molecules. Our results show that OSIC-SCAN provides superior results than the previously reported OSIC-LSDA, -PBE, and -TPSS results. Furthermore, we propose selective scaling of OSIC (SOSIC) to remove its major shortcoming that destroys the -1/r asymptotic behavior of the potentials. The SOSIC method gives the highest occupied orbital eigenvalues practically identical to those in PZSIC and unlike OSIC provides bound atomic anions even with larger powers of scaling factors. Furthermore, SOSIC compared to PZSIC or OSIC provides more balanced description of total energies and barrier heights.

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

Blind test of density-functional-based methods on intermolecular interaction energies

In the past decade, a number of approaches have been developed to fix the failure of (semi)local density-functional theory (DFT) in describing intermolecular interactions. In this work, the performance of several such approaches with respect to highly accurate benchmarks is compared on a set of separation-dependent interaction energies for ten dimers. Since the benchmarks were unknown before the DFT-based results were collected, this comparison constitutes a blind test of these methods.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Molecular dynamics simulation studies of 1,3-dimethyl imidazolium nitrate ionic liquid with water

The fundamental understanding of intermolecular interactions of ionic liquids (ILs) with water is essential in predicting IL–water thermodynamic properties. In this study, intermolecular or noncovalent interactions were studied for 1,3-dimethyl imidazolium [DMIM]+ cation and nitrate [NO 3 ] – anion with water, employing quantum mechanics and molecular dynamics simulations. Molecular dynamics simulations were performed using a revised multipolar polarizable force field. The effect of water on ionic liquids was evaluated in terms of thermodynamic and dynamic properties. Thermodynamic properties included liquid densities ρ, excess molar volumes ΔV E , and liquid structures g(r)⁠. Dynamic properties included self-diffusion coefficients D of mixture constituents as a function of water concentration. The density of ionic liquid–water mixtures monotonically decrease with increasing concentration of water. A negative excess volume was obtained for low and high water concentrations, demonstrating strong intermolecular interactions of water with ionic liquid components. Liquid structures of ionic liquid–water mixtures revealed a tendency for anions to interact with cations at shorter intermolecular distances when the water concentration is increased. Diffusion rates were found to increase for all mixture components with increase in the fraction of water. A significant change in the diffusion rate was found at ~0.3 weight fraction of water. However, the water self-diffusion coefficient was dominant at all concentrations. Furthermore, the ratio of water/anion and anion/cation self-diffusion coefficients was found to decrease linearly with increasing concentration of water molecules.

36 MATERIALS SCIENCE↗
DOE OSTI2020

QMCPACK: Advances in the development, efficiency, and application of auxiliary field and real-space variational and diffusion quantum Monte Carlo

We review recent advances in the capabilities of the open source ab initio Quantum Monte Carlo (QMC) package QMCPACK and the workflow tool Nexus used for greater efficiency and reproducibility. The auxiliary field QMC (AFQMC) implementation has been greatly expanded to include k-point symmetries, tensor-hypercontraction, and accelerated graphical processing unit (GPU) support. These scaling and memory reductions greatly increase the number of orbitals that can practically be included in AFQMC calculations, increasing the accuracy. Advances in real space methods include techniques for accurate computation of bandgaps and for systematically improving the nodal surface of ground state wavefunctions. Results of these calculations can be used to validate application of more approximate electronic structure methods, including GW and density functional based techniques. To provide an improved foundation for these calculations, we utilize a new set of correlation-consistent effective core potentials (pseudopotentials) that are more accurate than previous sets; these can also be applied in quantum-chemical and other many-body applications, not only QMC. These advances increase the efficiency, accuracy, and range of properties that can be studied in both molecules and materials with QMC and QMCPACK.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Full dimensional Franck-Condon factors for the acetylene A{sup ~} {sup 1}A{sub u}—X{sup ~1}Σ{sup +}{sub g} transition. I. Method for calculating polyatomic linear—bent vibrational intensity factors and evaluation of calculated intensities for the gerade vibrational modes in acetylene

Franck-Condon vibrational overlap integrals for the A{sup ~1}A{sub u}—X{sup ~1}Σ{sup +}{sub g} transition in acetylene have been calculated in full dimension in the harmonic normal mode basis. The calculation uses the method of generating functions first developed for polyatomic Franck-Condon factors by Sharp and Rosenstock [J. Chem. Phys. 41(11), 3453–3463 (1964)], and previously applied to acetylene by Watson [J. Mol. Spectrosc. 207(2), 276–284 (2001)] in a reduced-dimension calculation. Because the transition involves a large change in the equilibrium geometry of the electronic states, two different types of corrections to the coordinate transformation are considered to first order: corrections for axis-switching between the Cartesian molecular frames and corrections for the curvilinear nature of the normal modes at large amplitude. The angular factor in the wave function for the out-of-plane component of the trans bending mode, ν{sub 4}{sup ″}, is treated as a rotation, which results in an Eckart constraint on the polar coordinates of the bending modes. To simplify the calculation, the other degenerate bending mode, ν{sub 5}{sup ″}, is integrated in the Cartesian basis and later transformed to the constrained polar coordinate basis, restoring the conventional v and l quantum numbers. An updated A{sup ~}-state harmonic force field obtained recently in the R. W. Field research group is evaluated. The results for transitions involving the gerade vibrational modes are in qualitative agreement with experiment. Calculated results for transitions involving ungerade modes are presented in Paper II of this series [G. B. Park, J. H. Baraban, and R. W. Field, “Full dimensional Franck–Condon factors for the acetylene A{sup ~1}A{sub u}—X{sup ~1}Σ{sup +}{sub g} transition. II. Vibrational overlap factors for levels involving excitation in ungerade modes,” J. Chem. Phys. 141, 134305 (2014)].

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

An electronic criterion for assessing intrinsic brittleness of metallic glasses

Bulk metallic glasses (BMGs) are characterized by a number of remarkable physical and mechanical properties. Unfortunately, these same materials are often intrinsically brittle, which limits their utility. Consequently, considerable effort has been expended searching for correlations between the phenomenologically complex mechanical properties of metallic glasses and more basic properties, such correlations might provide insight into the structure and bonding controlling the deformation properties of BMGs. While conducting such a search, we uncovered a weak correlation between a BMG’s work function and its susceptibility to brittle behavior. We argue that the basis for this correlation is a consequence of a component of the work function – the surface dipole – and a fundamental bond property related to the shape of the charge density at a bond critical point. Together these observations suggest that simple first principle calculations might be useful in the search for tougher BMGs.

36 MATERIALS SCIENCE↗
DOE OSTI2020

On the light-front wave functions of quarkonia

The light-front wave functions of hadrons allow us to calculate a wide range of physical observables; however, the wave functions themselves cannot be measured. We discuss recent results for quarkonia obtained in basis light-front quantization using an effective Hamiltonian with a confining model in both the transverse and longitudinal directions and with explicit one-gluon exchange. In particular, we focus on the numerical convergence of the basis expansion, as well as the asymptotic behavior of the light-front wave functions. Finally, we also illustrate that, for mesons with unequal quark masses, the maxima of the light-front wave functions depend in a non-trivial way on the valence quark-mass difference.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

Norm-conserving pseudopotentials and basis sets to explore Lanthanide chemistry in complex environments

A complete set of pseudopotentials and accompanying basis sets for all lanthanide elements are presented based on the relativistic, norm-conserving, separable, dual-space Gaussian-type pseudopotential protocol of Goedecker, Teter and Hutter (GTH) within the generalized gradient approximation (GGA) and the exchange-correlation functional of Perdew, Burke and Ernzerhof (PBE). The corresponding basis sets have been molecularly optimized (MOLOPT) using a contracted form with a single set of Gaussian exponents for s, p and d states. The f states are uncontracted explicitly with Gaussian exponents. Moreover, the Hubbard U values for each lanthanide element, to be used in DFT+U calculations, are also tabulated, allowing for the proper treatment of the strong on-site Coulomb interactions of localized 4f electrons. The accuracy and reliability of our GTH pseudopotentials and companion basis sets optimized for lanthanides is illustrated by a series of test calculations on lanthanide-centered molecules, and solid-state systems, with the most common oxidation states. We anticipate that these pseudopotentials and basis sets will enable larger-scale density functional theory calculations and ab initio molecular dynamics simulations of lanthanide molecules in either gas or condensed phases, as well as of solid state lanthanide-containing materials, allowing to further explore the chemical and physical properties of lanthanide systems.

Lu, Jun-Bo↗
DOE OSTI2022

Magnetic exchange interactions in binuclear and tetranuclear iron(III) complexes described by spin‐flip DFT and Heisenberg effective Hamiltonians

Abstract Low‐energy spectra of single‐molecule magnets (SMMs) are often described by Heisenberg Hamiltonians. Within this formalism, exchange interactions between magnetic centers determine the ground‐state multiplicity and energy separation between the ground and excited states. In this contribution, we extract exchange coupling constants ( J ) for a set of iron (III) binuclear and tetranuclear complexes from all‐electron calculations using non‐collinear spin‐flip time‐dependent density functional theory (NC‐SF‐TDDFT). For 12 binuclear complexes with J ‐values ranging from −6 to −132 cm −1 , our benchmark calculations using the short‐range hybrid ω PBEh functional and 6‐31G(d,p) basis set agree well with the experimentally derived values (mean absolute error of 4.7 cm −1 ). For the tetranuclear SMMs, the computed J constants are within 6 cm −1 from the experimentally derived values. We explore the range of applicability of the Heisenberg model by analyzing bonding patterns in these Fe(III) complexes using natural orbitals (NO), their occupations, and the number of effectively unpaired electrons. The results illustrate the efficiency of the spin‐flip protocol for computing the exchange couplings and the utility of the NO analysis in assessing the validity of effective spin Hamiltonians.

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

Leading-color two-loop amplitudes for four partons and a $\textit{W}$ boson in QCD

We present the leading-color two-loop QCD corrections for the scattering of four partons and a W boson, including its leptonic decay. The amplitudes are assembled from the planar two-loop helicity amplitudes for four partons and a vector boson decaying to a lepton pair, which are also used to determine the planar two-loop amplitudes for four partons and a $Z/γ^\ast$ boson with a leptonic decay. The analytic expressions are obtained by setting up a dedicated Ansatz and constraining the free parameters from numerical samples obtained within the framework of numerical unitarity. The large linear systems that must be solved to determine the analytic expressions are constructed to be in Vandermonde form. Such systems can be very efficiently solved, bypassing the bottleneck of Gaussian elimination. Our results are expressed in a basis of one-mass pentagon functions, which opens the possibility of their efficient numerical evaluation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

Finite deformation cohesive polygonal finite elements for modeling pervasive fracture

In this work, we introduce a framework for modeling dynamic fracture problems using cohesive polygonal finite elements. Random polygonal meshes provide a robust, efficient method for generating an unbiased network of fracture surfaces. Further, these meshes have more facets per element than standard triangle or quadrilateral meshes, providing more possible facets per element to insert cohesive surfaces. This property of polygonal meshes is advantageous for the modeling of pervasive fracture. We use both Wachspress and maximum entropy shape functions to form a finite element basis over the polygons. Fracture surfaces are captured through dynamically inserted cohesive zone elements at facets between the polygons in the mesh. Contact is enforced through a penalty method that is applied to both closed cohesive surfaces and general interpenetration of two polygonal elements. Finally, several numerical examples are presented that illustrate the capabilities of the method and demonstrate convergence of solutions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Decisive role of magnetism in the interaction of chromium and nickel solute atoms with 1/2$\langle$111$\rangle$-screw dislocation core in body-centered cubic iron

Resistance to swelling under irradiation and a low rate of corrosion in high temperature environments make Fe-Cr and Fe-Cr-Ni alloys promising structural materials for energy technologies. In this paper we report the results obtained using a combination of density functional theory (DFT) techniques: plane wave basis set solutions for pseudo-potentials and multiple scattering solutions for all electron potentials. We have found a very strong role of magnetism in the stability of screw dislocation cores in pure Fe and their interaction with Cr and Ni magnetic impurities. In particular, the screw dislocation quadrupole in Fe is stabilized only in the presence of ferromagnetism. In addition, Ni atoms, who's magnetic moment is oriented along the magnetization direction of the Fe matrix, prefer to occupy in core positions whereas Cr atoms, which couple anti-ferromagnetically with the Fe matrix, prefer out of the dislocation core positions. In effect, Ni impurities are attracted to, while Cr impurities are repelled by the dislocation core. Moreover, we demonstrate that this contrasting behavior can be explained only by the nature of magnetic coupling of the impurities to the Fe matrix. In addition, Cr interaction with the dislocation core mirrors that of Ni if the Cr magnetic moment is constrained to be along the direction of Fe matrix magnetization. In addition, we have shown that the magnetic contribution can affect the impurity-impurity interaction at distances up to a few Burgers vectors. In particular, the distance between Cr atoms in Fe matrix should be at least 3–4 lattice parameters in order to eliminate finite size effects.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A case study of high-temperature polyetherimide film capacitor fabrication

Capacitors, as the basic building block of power electronics and electrical systems, are the major constraint of the increasingly integrated power systems that require new capable polymer dielectric films operating at higher temperatures >125 °C. Despite tremendous research efforts on the lab scale, there remains a considerable barrier and knowledge basis for converting polymeric films to functioning capacitors desired in actual applications. Here, the authors tackled the polyetherimide (PEI) film scale up issues and developed various engineering processes for film de-wrinkling, optimal metallization, static elimination, and capacitor fabrication improvement. The authors discovered various fabrication challenges such as the electrostatic charge, metallization scheme, winding tension, and end electrode disconnects. To fabricate a high yield of PEI capacitor bobbins, it is necessary to utilize static eliminators (radioactive ionizers), moderately thick aluminum metallization (15–30 Ω/sq), and winding tension (60–100 g). It is also effective to evaluate capacitance, dielectric loss, equivalent series resistance, and thermal cycling stability of capacitors. This work sheds some light on converting PEI films and other polymeric films to capacitors on a scale up fabrication effort. The practical learning from film handling to capacitor fabrication in this work provided the necessary knowledge for manufacturing high-temperature polar film capacitors.

36 MATERIALS SCIENCE↗