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,621 records · Page 90

The Impact of Nonequilibrium Ionization on SDO/AIA and Hinode/EIS Observations of Impulsively Heated Plasmas

Most plasma diagnostics assume the emitting material is in a state of ionization equilibrium. For example, the AIA temperature response functions have been derived on this basis. The assumption is reasonable whenever the plasma is evolving slowly or is very dense, but these are not the conditions that apply during impulsive heating events. It is now widely believed that many coronal loops are bundles of unresolved strands that are heated quasi-randomly by nanoflares. Full blown flares are thought to have similar sub-structure. We have studied the importance of nonequilibrium effects in these circumstances by modeling nanoflare-heated loops and simulating their observation by AIA and the EIS spectrometer on Hinode. We find that the intensities of hot emission lines can be highly suppressed and that the net emission from the loop tends to be dominated by strands that have entered a slow cooling phase, well after the impulsive energy release has ended. The hottest strands are relatively invisible, both because they are tenuous and because they cool rapidly by thermal conduction. Thus, AIA channels that are normally thought of as being sensitive to hot plasma, such 131 and 94, are in fact frequently not able to detect the hot plasma that is present. The magnitude of the effect is case dependent. Great care must be exercised when using the standard temperature response functions in situations where nonequilibrium ionization is likely to be important.

Klimchuk, James A.↗

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↗

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↗

Near Hartree-Fock quality GTO basis sets for the second-row atoms

Energy optimized, near Hartree-Fock quality Gaussian basis sets ranging in size from (17s12p) to (20s15p) are presented for the ground states of the second-row atoms for Na(2P), Na(+), Na(-), Mg(3P), P(-), S(-), and Cl(-). In addition, optimized supplementary functions are given for the ground state basis sets to describe the negative ions, and the excited Na(2P) and Mg(3P) atomic states. The ratios of successive orbital exponents describing the inner part of the 1s and 2p orbitals are found to be nearly independent of both nuclear charge and basis set size. This provides a method of obtaining good starting estimates for other basis set optimizations.

Partridge, Harry↗

Near Hartree-Fock quality GTO basis sets for the second-row atoms

Energy optimized, near Hartree-Fock quality Gaussian basis sets ranging in size from (17s12p) to (20s15p) are presented for the ground states of the second-row atoms for Na(2P), Na(+), Na(-), Mg(3P), P(-), S(-), and Cl(-). In addition, optimized supplementary functions are given for the ground state basis sets to describe the negative ions, and the excited Na(2P) and Mg(3P) atomic states. The ratios of successive orbital exponents describing the inner part of the 1s and 2p orbitals are found to be nearly independent of both nuclear charge and basis set size. This provides a method of obtaining good starting estimates for other basis set optimizations.

Partridge, Harry↗

TiCl, TiH and TiH+ Bond Energies, a Test of a Correlation Consistent Ti Basis Set

Correlation consistent basis sets are developed for Ti atom. The polarization functions are optimized for the average of the 3F and 5F states. One series of correlation consistent basis sets is for 3d and 4s correlation, while the second series includes 3s and 3p correlation as well as 3d and 4s correlation. These basis sets are tested using the Ti 3F-5F separation and the dissociation energies of TiCl X4Phi, TiH X4Phi, and TiH(+) X3Phi. The CCSD(T) complete basis set limit values are determined by extrapolation. The Douglas Kroll approach is used to compute the scalar relativistic effect. Spin-orbit effects are taken from experiment and/or computed at the CASSCF level. The Ti 3F-5F separation is in excellent agreement with experiment, while the TiCl, TiH, and TiH(+) bond energies are in good agreement with experiment. Extrapolation with the valence basis set is consistent with other atoms, while including 3s and 3p correlation appears to make extrapolation.

Bauschlicher, Charles W., Jr.↗

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↗

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↗

Finite-field Calculations of Molecular Polarizabilities Using Field-Induced Polarization Functions: Second- and Third-order Perturbation Correlation Corrections to the Coupled Hartree-Fock Polarizability of H2O

Ordinary Rayleigh-Schroudinger perturbation theory with Moller-Plesset (RSMP) partitioning is used to calculate second- and third-order correlation corrections to the CHF polarizability and dipole moment of the water molecule by a finite-field procedure. Pade approximants are found to be useful in accelerating the convergence of the property perturbation expansions. Field-induced polarization functions suitable for polarizability calculations are determined. The average polarizability calculated, neglecting vibrational averaging, with Dunning's (9s5p/4s-4s2p/2s) contracted GTO basis set augmented by field-induced lslp2d/lp polarization functions is within 3 per cent of the experimental result. Correlation corrections to the dipole moment and polarizability of the water molecule calculated by the finite-field RSMP and single + double excitation CI(SDCI) methods for the same basis set are found to be in close agreement. The RSMP approach has the advantages of being size-consistent and of being capable of greater efficiency than the SCDI method. Comparative calculations carried out using Epstein-Nesbet partitioning show that through third order RSEN correlation perturbation expansions for the dipole moment and polarizability are less rapidly convergent than RSMP expansions. However, reasonable accord with RSMP results can be achieved by using Pade approximants to accelerate the convergence of RSEN energy perturbation expansions. The convergence of RSMP property correlation expansions based on the zeroth-order uncoupled-Hartree-Fock (UCHF) and coupled-Hartree-Fock (CHF) approximations are compared through third order. Whereas the CHF + RSMP expansions are for practical purposes fully converged, the UCHF + RSMP expansions are not adequately converged.

Langhoff, S. R.↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

3,5,4′-Trimethoxystilbene, a natural methoxylated analog of resveratrol, inhibits breast cancer cell invasiveness by downregulation of PI3K/Akt and Wnt/β-catenin signaling cascades and reversal of epithelial–mesenchymal transition

The molecular basis of epithelial–mesenchymal transition (EMT) functions as a potential therapeutic target for breast cancer because EMT may endow breast tumor-initiating cells with stem-like characteristics and enable the dissemination of breast cancer cells. We have recently verified the antitumor activity of 3,5,4′-trimethoxystilbene (MR-3), a naturally methoxylated derivative of resveratrol, in colorectal cancer xenografts via an induction of apoptosis. The effect of MR-3 on EMT and the invasiveness of human MCF-7 breast adenocarcinoma cell line were also explored. We found that MR-3 significantly increased epithelial marker E-cadherin expression and triggered a cobblestone-like morphology of MCF-7 cells, while reciprocally decreasing the expression of mesenchymal markers, such as snail, slug, and vimentin. In parallel with EMT reversal, MR-3 downregulated the invasion and migration of MCF-7 cells. Exploring the action mechanism of MR-3 on the suppression of EMT and invasion indicates that MR-3 markedly reduced the expression and nuclear translocation of β-catenin, accompanied with the downregulation of β-catenin target genes and the increment of membrane-bound β-catenin. These results suggest the involvement of Wnt/β-catenin signaling in the MR-3-induced EMT reversion of MCF-7 cells. Notably, MR-3 restored glycogen synthase kinase-3β activity by inhibiting the phosphorylation of Akt, the event required for β-catenin destruction via a proteasome-mediated system. Overall, these findings indicate that the anti-invasive activity of MR-3 on MCF-7 cells may result from the suppression of EMT via down-regulating phosphatidylinositol 3-kinase (PI3K)/AKT signaling, and consequently, β-catenin nuclear translocation. These occurrences ultimately lead to the blockage of EMT and the invasion of breast cancer cells. - Highlights: • MR-3 blocked MCF-7 cell invasion by inducing a reversal of EMT. • Wnt/β-catenin signaling is involved in MR-3-induced EMT reversion of MCF-7 cells. • Knockdown of β-catenin was sufficient to restore epithelial marker E-cadherin levels. • MR-3 recovered the function of GSK-3β that inhibits β-catenin nuclear translocation.

60 APPLIED LIFE SCIENCES↗