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,423 records · Page 79

Fundamental Vibrational Analyses of the HCN Monomer, Dimer and Associated Isotopologues

In this work we provide high level ab initio treatments of the structures, vibrational frequencies, and electronic energies of the HCN monomer and dimer systems along with several isotopologues. Additionally, the plethora of information related to this system within the literature is summarized and serves as a basis for comparison with the results of this paper. The geometry of the dimer and monomer are reported at the all electroncoupled-cluster singles, doubles, and perturbative triples level of theory [AE-CCSD(T)] with the correlation consistent quadruple-zeta quality basis sets with extra core functions (cc-pCVQZ) from Dunning. The theoretical geometries and electronic structures are further analyzed through the use of the Natural Bond Orbital (NBO) method and Natural Resonance Theory (NRT). At the AE-CCSD(T)/cc-pCVQZ level of theory, the full cubic with semi-diagonal quartic force field for nine dimer and four monomer isotopologues (the parent isotopologue along with 15N, 13C, and D derivatives) were obtained to treat the anharmonicity of the vibrations via second order vibrational perturbation theory (VPT2). Lastly, the enthalpy change associated with the formation of the dimer from two monomer units was determined using the focal point analysis. Computations including coupled-cluster through perturbative quadruples as well as basis sets up to six-zeta quality, including core functions (cc-pCVXZ, X=D,T,Q,5,6) were used to extrapolate to the AE-CCSDT(Q)/CBS energy associated with this hydrogen-bond forming process. After appending anharmonic zero-point vibrational, relativistic, and diagonal Born-Oppenheimer corrections, we report a value of -3.93 kcal mol -1 for the enthalpy of formation. To our knowledge, each set of results (geometries, vibrational frequencies, and energetics) reported in this study represents the highest-level and most reliable theoretical predictions reported for this system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ab Initio Thermodynamics of Iridium Surface Oxidation and Oxygen Evolution Reaction

Iridium-based materials are considered as state-of-the-art electrocatalysts for oxygen evolution reaction (OER), however, their stability and catalytic activity greatly depend on surface-state changes induced by electrochemical cycling. To better understand the behavior of the low-index Ir surfaces in an electrochemical environment, we perform a systematic thermodynamic analysis by means of the density functional theory (DFT) calculations. On the basis of computed surface energies of the Ir (111), (110) and (100) facets as a function of applied electrode potential and coverage of adsorbed water species we determine stability maps and predict equilibrium shapes of Ir nanoparticles. Our calculations also show that metastable oxide precursors formed at the initial stages of Ir surface oxidation are responsible for enhanced catalytic activity toward OER as compared to metal surfaces covered by oxygen adsorbates and thick-oxide films. Such enhancement occurs not only due to the modified thermodynamic stability of OER intermediates, but also because thin-oxide layers may display the more energetically favorable I2M (interaction of two M–O units) rather than WNA (water nucleophilic attack) OER mechanism.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Properties of a potential energy matrix in oscillator basis

Matrix elements of potential energy are examined in detail. We consider a model problem — a particle in a central potential. The most popular forms of central potential are taken up, namely, square-well potential, Gaussian, Yukawa and exponential potentials. We study eigenvalues and eigenfunctions of the potential energy matrix constructed with oscillator functions. It is demonstrated that eigenvalues coincide with the potential energy in coordinate space at some specific discrete points. We establish approximate values for these points. It is also shown that the eigenfunctions of the potential energy matrix are the expansion coefficients of the spherical Bessel functions in a harmonic oscillator basis. We also demonstrate a close relation between the separable approximation and L{sup 2} basis (J-matrix) method for the quantum theory of scattering.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A NEW TOOL FOR IMAGE ANALYSIS BASED ON CHEBYSHEV RATIONAL FUNCTIONS: CHEF FUNCTIONS

We introduce a new approach to the modeling of the light distribution of galaxies, an orthonormal polar basis formed by a combination of Chebyshev rational functions and Fourier polynomials that we call CHEF functions, or CHEFs. We have developed an orthonormalization process to apply this basis to pixelized images, and implemented the method as a Python pipeline. The new basis displays remarkable flexibility, being able to accurately fit all kinds of galaxy shapes, including irregulars, spirals, ellipticals, highly compact, and highly elongated galaxies. It does this while using fewer components than similar methods, as shapelets, and without producing artifacts, due to the efficiency of the rational Chebyshev polynomials to fit quickly decaying functions like galaxy profiles. The method is linear and very stable, and therefore is capable of processing large numbers of galaxies in a fast and automated way. Due to the high quality of the fits in the central parts of the galaxies, and the efficiency of the CHEF basis modeling galaxy profiles up to very large distances, the method provides highly accurate estimates of total galaxy fluxes and ellipticities. Future papers will explore in more detail the application of the method to perform multiband photometry, morphological classification, and weak shear measurements.

79 ASTRONOMY AND ASTROPHYSICS↗

Theoretical study of the low-lying bound states of O2

It is demonstrated that a complete-active-space self-consistent-field (CASSCF) (2p)/MRCI + Q (multireference configuration interaction with a Davidson correction) description in a (13s8p6d 4f2g)/((5s4p3d 2f1g) atomic natural orbits (ANO) basis set supplemented with diffuse functions provides a quantitative description of the six lowest states of O2. The calculated potentials are within 0.05 eV (1.2 kilocal/mol) of accurate experimental results. The importance of substantially expanding the primitive basis set has been investigated, and it is demonstrated that such expansions yield insignificant improvement in the spectroscopic constants. Potential energy curves have also been reported for the weakly bound states of O2. The 5Pi(g) state is estimated to have a D(e) of 0.16 +/- 0.03 eV. The upper bound of D(e) is found to be sufficiently large that the importance of this state as a precursor for the formation of O2 (b 1Sigma(t)(+)) and O(1S) should be reconsidered.

Partridge, Harry↗

First-Principles Modeling of Non-Covalent Interactions in Supramolecular Systems: The Role of Many-Body Effects

Supramolecular host–guest systems play an important role for a wide range of applications in chemistry and biology. The prediction of the stability of host–guest complexes represents a great challenge to first-principles calculations due to an interplay of a wide variety of covalent and noncovalent interactions in these systems. In particular, van der Waals (vdW) dispersion interactions frequently play a prominent role in determining the structure, stability, and function of supramolecular systems. On the basis of the widely used benchmark case of the buckyball catcher complex (C60@C60H28), we assess the feasibility of computing the binding energy of supramolecular host–guest complexes from first principles. Large-scale diffusion Monte Carlo (DMC) calculations are carried out to accurately determine the binding energy for the C60@C60H28 complex (26 ± 2 kcal/mol). On the basis of the DMC reference, we assess the accuracy of widely used and efficient density-functional theory (DFT) methods with dispersion interactions. The inclusion of vdW dispersion interactions in DFT leads to a large stabilization of the C60@C60H28 complex. However, DFT methods including pairwise vdW interactions overestimate the stability of this complex by 9–17 kcal/mol compared to the DMC reference and the extrapolated experimental data. A significant part of this overestimation (9 kcal/mol) stems from the lack of dynamical dielectric screening effects in the description of the molecular polarizability in pairwise dispersion energy approaches. The remaining overstabilization arises from the isotropic treatment of atomic polarizability tensors and the lack of many-body dispersion interactions. A further assessment of a different supramolecular system – glycine anhydride interacting with an amide macrocycle – demonstrates that both the dynamical screening and the many-body dispersion energy are complex contributions that are very sensitive to the underlying molecular geometry and type of bonding. We discuss the required improvements in theoretical methods for achieving “chemical accuracy” in the first-principles modeling of supramolecular systems.

Chemistry↗

Exogenic and endogenic albedo and color patterns on Europa

New global and high-resolution multispectral mosaics of Europa have been produced from the Voyager imaging data. Photometric normalizations are based on multiple-image techniques that explicitly account for intrinsic albedo variations through pixel-by-pixel solutions. The exogenic color and albedo pattern on Europa is described by a second-order function of the cosine of the angular distance from the apex of orbital motion. On the basis of this second-order function and of color trends that are different on the leading and trailing hemispheres, the exogenic pattern is interpreted as being due to equilibrium between two dominant processes: (1) impact gardening and (2) magnetospheric interactions, including sulfur-ion implantation and sputtering redistribution. Removal of the model exogenic pattern in the mosaics reveals the endogenic variations, consisting of only two major units: darker (redder) and bright materials. Therefore Europa's visual spectral reflectivity is simple, having one continuous exogenic pattern and two discrete endogenic units.

Mcewen, A. S.↗

Tunable-fidelity wave functions for the ab initio description of scattering and reactions

Here, the no-core shell model (NCSM) is an ab initio method that solves the nuclear many-body problem by expanding the many-particle wave function into a (typically) harmonic oscillator basis and minimizing the energy to obtain the expansion coefficients. Extensions of the NCSM, such as its coupling with microscopic-cluster basis states, further allow for an ab initio treatment of light-ion nuclear reactions of interest for both astrophysics and nuclear technology applications. A downside of the method is the exponential scaling of the basis size with increasing number of nucleons and excitation quanta, which limits its applicability to mass A ≲ 16 nuclei, except for variants where the basis is further down-selected via some truncation scheme. We consider a basis selection method for the NCSM that was first introduced in the context of the large-scale shell model and captures the essential degrees of freedom of the nuclear wave function leading to a favorable complexity scaling for calculations and enabling ab initio reaction calculations in sd-shell nuclei. The particle configurations within the NCSM basis are ordered based on their contribution to the first moment of the Hamiltonian matrix that results from the projection onto the many-body basis. The truncation scheme then consists in retaining only the lowest-first-moment configurations, which typically contain only few many-body basis states (Slater determinants). As the energy threshold above which configurations are disregarded is increased, the size of the basis becomes an almost-continuous variable, allowing for tunable fidelity in the obtained wave functions. The resulting wave functions can then be used directly in ab initio reaction calculations. We present calculations for 7 Li and n + 12 C scattering using nucleon-nucleon interactions derived from chiral effective field theory and softened using the similarity renormalization group method. The obtained energy levels invariably demonstrate exponential convergence with the size of the basis, and we find improved convergence in scattering calculations. To demonstrate the possibilities enabled by the approach, we also present a first calculation for the scattering of neutrons from 24 Mg. The method presented in this work appears promising for future studies of nuclei with mass A > 16, opening multiple future research directions impacting both nuclear astrophysics and nuclear technology applications.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Min-Max Spaces and Complexity Reduction in Min-Max Expansions

Idempotent methods have been found to be extremely helpful in the numerical solution of certain classes of nonlinear control problems. In those methods, one uses the fact that the value function lies in the space of semiconvex functions (in the case of maximizing controllers), and approximates this value using a truncated max-plus basis expansion. In some classes, the value function is actually convex, and then one specifically approximates with suprema (i.e., max-plus sums) of affine functions. Note that the space of convex functions is a max-plus linear space, or moduloid. In extending those concepts to game problems, one finds a different function space, and different algebra, to be appropriate. Here we consider functions which may be represented using infima (i.e., min-max sums) of max-plus affine functions. It is natural to refer to the class of functions so represented as the min-max linear space (or moduloid) of max-plus hypo-convex functions. We examine this space, the associated notion of duality and min-max basis expansions. In using these methods for solution of control problems, and now games, a critical step is complexity-reduction. In particular, one needs to find reduced-complexity expansions which approximate the function as well as possible. We obtain a solution to this complexity-reduction problem in the case of min-max expansions.

97 MATHEMATICS AND COMPUTING↗

Hybrid preconditioning for iterative diagonalization of ill-conditioned generalized eigenvalue problems in electronic structure calculations

The iterative diagonalization of a sequence of large ill-conditioned generalized eigenvalue problems is a computational bottleneck in quantum mechanical methods employing a nonorthogonal basis for ab initio electronic structure calculations. We propose a hybrid preconditioning scheme to effectively combine global and locally accelerated preconditioners for rapid iterative diagonalization of such eigenvalue problems. In partition-of-unity finite-element (PUFE) pseudopotential density-functional calculations, employing a nonorthogonal basis, we show that the hybrid preconditioned block steepest descent method is a cost-effective eigensolver, outperforming current state-of-the-art global preconditioning schemes, and comparably efficient for the ill-conditioned generalized eigenvalue problems produced by PUFE as the locally optimal block preconditioned conjugate-gradient method for the well-conditioned standard eigenvalue problems produced by planewave methods.

74 ATOMIC AND MOLECULAR PHYSICS↗

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↗

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↗

Preliminary Altitude Performance Characteristics of the J57-P-1 Turbojet Engine with Fixed-area Exhaust Nozzle

An investigation to determine the altitude performance of the J57-P-1 turbojet engine and components was conducted at the NACA Lewis altitude wind tunnel. Data were obtained over a corrected inboard rotor speed range from 56 to 106 percent of rated speed, with intercompressor bleeds both open and closed, at altitudes from 15,000 to 50,000 feet and at a flight Mach number of 0.81. The corresponding range of Reynolds number indices was from 0.858 to 0.213. All data presented were obtained with a fixed-area exhaust nozzle sized according to the manufacturer's specification. Over-all engine performance parameters are presented as functions of inboard rotor speed corrected on the basis of engine inlet temperature. Component parameters are presented as functions of their respective corrected rotor speeds. A tabulation of all performance data is included in addition to the graphical presentation. Corrected net thrust is unusually sensitive to changes in corrected inboard rotor speed in the high speed region. A change of 1 percent in speed, at sated speed, produced a change of 6 percent in corrected net thrust . At rated engine speed, increasing the altitude from 15,000 to 50,000 feet at a constant flight Mach number of 0.81 increased the specific fuel consumption 13 percent but did not affect corrected net thrust.

Bloomer, Harry E.↗

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↗

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↗

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↗

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↗