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,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Multi-petascale highly efficient parallel supercomputer

A Multi-Petascale Highly Efficient Parallel Supercomputer of 100 petaOPS-scale computing, at decreased cost, power and footprint, and that allows for a maximum packaging density of processing nodes from an interconnect point of view. The Supercomputer exploits technological advances in VLSI that enables a computing model where many processors can be integrated into a single Application Specific Integrated Circuit (ASIC). Each ASIC computing node comprises a system-on-chip ASIC utilizing four or more processors integrated into one die, with each having full access to all system resources and enabling adaptive partitioning of the processors to functions such as compute or messaging I/O on an application by application basis, and preferably, enable adaptive partitioning of functions in accordance with various algorithmic phases within an application, or if I/O or other processors are underutilized, then can participate in computation or communication nodes are interconnected by a five dimensional torus network with DMA that optimally maximize the throughput of packet communications between nodes and minimize latency.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2014

Conformational effect of dicyclo-hexano-18-crown-6 on isotopic fractionation of zinc: DFT approach

Generalized gradient approximated BP86 density functional employing triple zeta valence plus polarization (TZVP) basis set has been used to compute the reduced partition function ratio and isotopic separation factor for zinc isotopes. The isotopic separation factor was found to be in good agreement with the experimental results. The isotopic separation factor was found to depend on the conformation of the crown ether ligand. The trans-trans conformation shows the highest fractionation compared to cis-cis conformer. The present theoretical results can thus be used to plan the isotope separation experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Deformation properties of lead isotopes

The deformation properties of a long lead isotopic chain up to the neutron drip line are analyzed on the basis of the energy density functional (EDF) in the FaNDF{sup 0} Fayans form. The question of whether the ground state of neutron-deficient lead isotopes can have a stable deformation is studied in detail. The prediction of this deformation is contained in the results obtained on the basis of the HFB-17 and HFB-27 Skyrme EDF versions and reported on Internet. The present analysis reveals that this is at odds with experimental data on charge radii and magnetic moments of odd lead isotopes. The Fayans EDF version predicts a spherical ground state for all light lead isotopes, but some of them (for example, {sup 180}Pb and {sup 184}Pb) prove to be very soft—that is, close to the point of a phase transition to a deformed state. Also, the results obtained in our present study are compared with the predictions of some other Skyrme EDF versions, including SKM*, SLy4, SLy6, and UNE1. By and large, their predictions are closer to the results arising upon the application of the Fayans functional. For example, the SLy4 functional predicts, in just the same way as the FaNDF{sup 0} functional, a spherical shape for all nuclei of this region. The remaining three Skyrme EDF versions lead to a deformation of some light lead isotopes, but their number is substantially smaller than that in the case of the HFB-17 and HFB-27 functionals. Moreover, the respective deformation energy is substantially lower, which gives grounds to hope for the restoration of a spherical shape upon going beyond the mean-field approximation, which we use here. Also, the deformation properties of neutron-rich lead isotopes are studied up to the neutron drip line. Here, the results obtained with the FaNDF{sup 0} functional are compared with the predictions of the HFB-17, HFB-27, SKM*, and SLy4 Skyrme EDF versions. All of the EDF versions considered here predict the existence of a region where neutron-rich lead isotopes undergo deformations, but the size of this region is substantially different for the different functionals being considered. Once again, it is maximal for the HFB-17 and HFB-27 functionals, is substantially narrower for the FaNDF{sup 0} functional, and is still narrower for the SKM* and SLy4 functionals. The two-neutron drip line proved to be A{sub drip}{sup 2n} = 266 for all of the EDF versions considered here, with the exception of SKM*, for which it is shifted to A{sub drip}{sup 2n}(SKM*) = 272.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Huzinaga projection embedding for efficient and accurate energies of systems with localized spin-densities

We demonstrate the accuracy and efficiency of the restricted open-shell and unrestricted formulation of the absolutely localized Huzinaga projection operator embedding method. Restricted open-shell and unrestricted Huzinaga projection embedding in the full system basis is formally exact to restricted open-shell and unrestricted Kohn–Sham density functional theory, respectively. By utilizing the absolutely localized basis, we significantly improve the efficiency of the method while maintaining high accuracy. Furthermore, the absolutely localized basis allows for high accuracy open-shell wave function methods to be embedded into a closed-shell density functional theory environment. The open-shell embedding method is shown to calculate electronic energies of a variety of systems to within 1 kcal/mol accuracy of the full system wave function result. For certain highly localized reactions, such as spin transition energies on transition metals, we find that very few atoms are necessary to include in the wave function region in order to achieve the desired accuracy. This extension further broadens the applicability of our absolutely localized Huzinaga level-shift projection operator method to include open-shell species. Here, we apply our method to several representative examples, such as spin splitting energies, catalysis on transition metals, and radical reactions.

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

Understanding the conductance switching of permethyloligosilanes: A theoretical approach

On the basis of ab initio density functional theory coupled with non-equilibrium Green’s function technique, we have presented a molecular level understanding on the stereoelectronic switching of conducting properties of oligosilane molecules. Su et al. [Nat. Chem. 7, 215–220 (2015)] demonstrated that these types of oligosilane molecules exhibit three stereoconformers which show two distinct conducting profiles. On the basis of break-junction technique, the authors show that manipulating a specific dihedral angle and thereby controlling the length of the molecular contact, it is possible to switch the conducting states minutely. However, their discussions scarce the proper energy level alignment upon which the molecule-lead tunneling amplitude depends. On the basis of electronic structure and non-equilibrium electron transport calculations, we interpret such switching behavior and thus quantify the switching parameter demonstrating how the metal-molecule contact geometry along with the electronic energy level alignment is responsible for such kind of junction process. We also provide the variation of switching parameter and the type of majority carrier with the conjugation length of the oligosilanes.

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

Squeezed states and Hermite polynomials in a complex variable

Following the lines of the recent paper of J.-P. Gazeau and F. H. Szafraniec [J. Phys. A: Math. Theor. 44, 495201 (2011)], we construct here three types of coherent states, related to the Hermite polynomials in a complex variable which are orthogonal with respect to a non-rotationally invariant measure. We investigate relations between these coherent states and obtain the relationship between them and the squeezed states of quantum optics. We also obtain a second realization of the canonical coherent states in the Bargmann space of analytic functions, in terms of a squeezed basis. All this is done in the flavor of the classical approach of V. Bargmann [Commun. Pure Appl. Math. 14, 187 (1961)].

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

Quantifying local exciton, charge resonance, and multiexciton character in correlated wave functions of multichromophoric systems

A new method for quantifying the contributions of local excitation, charge resonance, and multiexciton configurations in correlated wave functions of multichromophoric systems is presented. The approach relies on fragment-localized orbitals and employs spin correlators. Its utility is illustrated by calculations on model clusters of hydrogen, ethylene, and tetracene molecules using adiabatic restricted-active-space configuration interaction wave functions. In addition to the wave function analysis, this approach provides a basis for a simple state-specific energy correction accounting for insufficient description of electron correlation. The decomposition scheme also allows one to compute energies of the diabatic states of the local excitonic, charge-resonance, and multi-excitonic character. The new method provides insight into electronic structure of multichromophoric systems and delivers valuable reference data for validating excitonic models.

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

Complexity and Structure of Circuits for Parity Functions

The paper is devoted to circuits implementing parity functions. A review of results establishing exact values of the complexity of parity functions is given. The structure of optimal circuits implementing parity functions is described for some bases. For one infinite basis, an upper bound for the complexity of parity functions is given.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2011

Quasiparticle-phonon interaction in the theory of finite Fermi systems

Within the Green's function method and on the basis of the method developed by V.A. Khodel for analyzing anharmonic effects, effects of quasiparticle-phonon interaction in the second order in the amplitude of phonon production are studied in two problems as a natural development of A.B. Migdal's theory of finite Fermi systems. Transitions between excited states and static moments of magic and nonmagic nuclei in excited states, each of which is described in the random-phase approximation, are considered. The results for this problem are found to differ considerably from those in the quasiparticle random-phase approximation. The inclusion of all second-order anharmonic effects in the extended theory of finite Fermi systems that extends the standard theory of finite Fermi systems to the case of taking into account quasiparticle-phonon interaction in order to describe excited states, but which does not take into account all such effects, is also considered. They are taken into account at a level that makes it possible to calculate static moments of odd nuclei-more precisely, the respective equation for the vertex function, which, in the theory of finite Fermi systems, is a basic ingredient that describes the interaction of a nucleus with an external field, is derived. Some numerical results obtained within the recently implemented self-consistent version of the extended theory of finite Fermi systems are also presented for 15 stable and unstable tin isotopes. These results give sufficient grounds to conclude that phenomenological systematics are inapplicable to giant dipole resonances in neutron-rich isotopes. The cross sections for radiative neutron capture that are calculated by usingmicroscopic strength functions for the neutron-rich isotopes 132Sn and 150Sn differ strongly from the cross sections calculated on the basis of a phenomenological description of giant dipole resonances. These results are of paramount importance for astrophysics and for the theory of nuclear data for reactors since they highlight the inapplicability of phenomenological systematics to giant dipole resonances and the need for self-consistent calculations for neutron-rich nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2014

A divide-conquer-recombine algorithmic paradigm for large spatiotemporal quantum molecular dynamics simulations

We introduce an extension of the divide-and-conquer (DC) algorithmic paradigm called divide-conquer-recombine (DCR) to perform large quantum molecular dynamics (QMD) simulations on massively parallel supercomputers, in which interatomic forces are computed quantum mechanically in the framework of density functional theory (DFT). In DCR, the DC phase constructs globally informed, overlapping local-domain solutions, which in the recombine phase are synthesized into a global solution encompassing large spatiotemporal scales. For the DC phase, we design a lean divide-and-conquer (LDC) DFT algorithm, which significantly reduces the prefactor of the O(N) computational cost for N electrons by applying a density-adaptive boundary condition at the peripheries of the DC domains. Our globally scalable and locally efficient solver is based on a hybrid real-reciprocal space approach that combines: (1) a highly scalable real-space multigrid to represent the global charge density; and (2) a numerically efficient plane-wave basis for local electronic wave functions and charge density within each domain. Hybrid space-band decomposition is used to implement the LDC-DFT algorithm on parallel computers. A benchmark test on an IBM Blue Gene/Q computer exhibits an isogranular parallel efficiency of 0.984 on 786 432 cores for a 50.3 × 10{sup 6}-atom SiC system. As a test of production runs, LDC-DFT-based QMD simulation involving 16 661 atoms is performed on the Blue Gene/Q to study on-demand production of hydrogen gas from water using LiAl alloy particles. As an example of the recombine phase, LDC-DFT electronic structures are used as a basis set to describe global photoexcitation dynamics with nonadiabatic QMD (NAQMD) and kinetic Monte Carlo (KMC) methods. The NAQMD simulations are based on the linear response time-dependent density functional theory to describe electronic excited states and a surface-hopping approach to describe transitions between the excited states. A series of techniques are employed for efficiently calculating the long-range exact exchange correction and excited-state forces. The NAQMD trajectories are analyzed to extract the rates of various excitonic processes, which are then used in KMC simulation to study the dynamics of the global exciton flow network. This has allowed the study of large-scale photoexcitation dynamics in 6400-atom amorphous molecular solid, reaching the experimental time scales.

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

Influence of the Basset force on the resonant behavior of an oscillator with fluctuating frequency

The influence of hydrodynamic interactions, such as Stokes and Basset forces, on the dynamics of a harmonically trapped Brownian tracer is considered. A generalized Langevin equation is used to describe the tracer’s response to an external periodic force and to dichotomous fluctuations of the stiffness of the trapping potential. Relying on the Shapiro-Loginov formula, exact expressions for the complex susceptibility and for the response function are presented. On the basis of these exact formulas, it is demonstrated that interplay of a multiplicative colored noise and the Basset force induced memory effects can generate a variety of cooperation effects, such as multiresonance versus the driving frequency, as well as stochastic resonance versus noise parameters. In particular, in certain parameter regions the response function exhibits a resonance-like enhancement at intermediate values of the intensity of the Basset force. Conditions for the appearance of these effects are also discussed.

42 ENGINEERING↗
DOE OSTI2015

Model for transport and reaction of defects and carriers within displacement cascades in gallium arsenide

A model is presented for recombination of charge carriers at evolving displacement damage in gallium arsenide, which includes clustering of the defects in atomic displacement cascades produced by neutron or ion irradiation. The carrier recombination model is based on an atomistic description of capture and emission of carriers by the defects with time evolution resulting from the migration and reaction of the defects. The physics and equations on which the model is based are presented, along with the details of the numerical methods used for their solution. The model uses a continuum description of diffusion, field-drift and reaction of carriers, and defects within a representative spherically symmetric cluster of defects. The initial radial defect profiles within the cluster were determined through pair-correlation-function analysis of the spatial distribution of defects obtained from the binary-collision code MARLOWE, using recoil energies for fission neutrons. Properties of the defects are discussed and values for their parameters are given, many of which were obtained from density functional theory. The model provides a basis for predicting the transient response of III-V heterojunction bipolar transistors to displacement damage from energetic particle irradiation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

First principles treatment of structural, optical, and thermoelectric properties of Li{sub 7}MnN{sub 4} as electrode for a Li secondary battery

The electronic structure, electronic charge density and linear optical properties of the metallic Li{sub 7}MnN{sub 4} compound, having cubic symmetry, are calculated using the full potential linearized augmented plane wave (FP-LAPW) method. The calculated band structure and density of states using the local density, generalized gradient and Engel–Vosko approximations, depict the metallic nature of the cubic Li{sub 7}MnN{sub 4} compound. The bands crossing the Fermi level in the calculated band structure are mainly from the Mn-d states with small support of N-p states. In addition, the Mn-d states at the Fermi level enhance the density of states, which is very useful for the electronic transport properties. The valence electronic charge density depicts strong covalent bond between Mn and two N atoms and polar covalent bond between Mn and Li atoms. The frequency dependent linear optical properties like real and imaginary part of the dielectric function, optical conductivity, reflectivity and energy loss function are calculated on the basis of the computed band structure. Both intra-band and inter-band transitions contribute to the calculated optical parameters. Using the BoltzTraP code, the thermoelectric properties like electrical and thermal conductivity, Seebeck coefficient, power coefficient and heat capacity of the Li{sub 7}MnN{sub 4} are also calculated as a function of temperature and studied.

36 MATERIALS SCIENCE↗
DOE OSTI2022

KSSOLV 2.0: An efficient MATLAB toolbox for solving the Kohn-Sham equations with plane-wave basis set

KSSOLV (Kohn-Sham Solver) is a MATLAB toolbox for performing Kohn-Sham density functional theory (DFT) calculations with a plane-wave basis set. KSSOLV 2.0 preserves the design features of the original KSSOLV software to allow users and developers to easily set up a problem and perform ground-state calculations as well as to prototype and test new algorithms. Furthermore, it includes new functionalities such as new iterative diagonalization algorithms, k-point sampling for electron band structures, geometry optimization and advanced algorithms for performing DFT calculations with local, semi-local, and hybrid exchange-correlation functionals. It can be used to study the electronic structures of both molecules and solids. So, we describe these new capabilities in this work through a few use cases. We also demonstrate the numerical accuracy and computational efficiency of KSSOLV on a variety of examples.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2013

Long-range interactions and the sign of natural amplitudes in two-electron systems

In singlet two-electron systems, the natural occupation numbers of the one-particle reduced density matrix are given as squares of the natural amplitudes which are defined as the expansion coefficients of the two-electron wave function in a natural orbital basis. In this work, we relate the sign of the natural amplitudes to the nature of the two-body interaction. We show that long-range Coulomb-type interactions are responsible for the appearance of positive amplitudes and give both analytical and numerical examples that illustrate how the long-distance structure of the wave function affects these amplitudes. We further demonstrate that the amplitudes show an avoided crossing behavior as function of a parameter in the Hamiltonian and use this feature to show that these amplitudes never become zero, except for special interactions in which infinitely many of them can become zero simultaneously when changing the interaction strength. This mechanism of avoided crossings provides an alternative argument for the non-vanishing of the natural occupation numbers in Coulomb systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

All-electron ab initio Bethe-Salpeter equation approach to neutral excitations in molecules with numeric atom-centered orbitals

The Bethe-Salpeter equation (BSE) based on GW quasiparticle levels is a successful approach for calculating the optical gaps and spectra of solids and also for predicting the neutral excitations of small molecules. We here present an all-electron implementation of the GW+BSE formalism for molecules, using numeric atom-centered orbital (NAO) basis sets. We present benchmarks for low-lying excitation energies for a set of small organic molecules, denoted in the literature as “Thiel’s set.” Literature reference data based on Gaussian-type orbitals are reproduced to about one millielectron-volt precision for the molecular benchmark set, when using the same GW quasiparticle energies and basis sets as the input to the BSE calculations. For valence correlation consistent NAO basis sets, as well as for standard NAO basis sets for ground state density-functional theory with extended augmentation functions, we demonstrate excellent convergence of the predicted low-lying excitations to the complete basis set limit. A simple and affordable augmented NAO basis set denoted “tier2+aug2” is recommended as a particularly efficient formulation for production calculations. We finally demonstrate that the same convergence properties also apply to linear-response time-dependent density functional theory within the NAO formalism.

Liu, Chi↗
DOE OSTI2019

Understanding Active Sites in Pyrolyzed Fe–N–C Catalysts for Fuel Cell Cathodes by Bridging Density Functional Theory Calculations and 57 Fe Mössbauer Spectroscopy

Pyrolyzed Fe–N–C materials are promising platinum-group-metal-free catalysts for proton-exchange membrane fuel cell cathodes. However, the detailed structure, oxidation, and spin states of their active sites are still undetermined. As such, 57 Fe Mössbauer spectroscopy has identified FeN x moieties as the most active sites, with their fingerprint being a doublet in low-temperature Mössbauer spectra. However, the interpretation of the doublets for such materials has lacked theoretical basis. Here, we applied density functional theory to calculate the quadrupole splitting energy of doublets (Δ E QS ) for a range of FeN x structures in different oxidation and spin states. The calculated and experimental values are then compared for a reference Fe–N–C catalyst, whereas further information on the Fe oxidation and spin states was obtained from electron paramagnetic resonance, superconducting quantum interference device, and 57 Fe Mössbauer spectroscopy under external magnetic field. The combined theoretical and experimental results identify the main presence of FeN x moieties in Fe(II) low-spin and Fe(III) high-spin states, whereas a minor fraction of sites could exist in the Fe(II) S = 1 state. From the analysis of the 57 Fe Mössbauer spectrum under the external magnetic field and the comparison of calculated and measured Δ E QS values, we assign the experimental doublet D1 with a mean Δ E QS value of around 0.9 mm·s –1 to Fe(III)N 4 C 12 moieties in high-spin state and the experimental doublet D2 with a mean Δ E QS value of around 2.3 mm·s –1 to Fe(II)N 4 C 10 moieties in low and medium spin. These conclusions indicate that D1 corresponds to surface-exposed sites, whereas D2 may correspond either to bulk sites that are inaccessible to O 2 or to surface sites that bind O 2 weaker than D1.

fuel cell↗