Search NASA⌕ Search

SEARCH · Search NASA

Results for “Approximation methods”

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,603 records · Page 89
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Increasing the applicability of density functional theory. IV. Consequences of ionization-potential improved exchange-correlation potentials

This paper's objective is to create a “consistent” mean-field based Kohn-Sham (KS) density functional theory (DFT) meaning the functional should not only provide good total energy properties, but also the corresponding KS eigenvalues should be accurate approximations to the vertical ionization potentials (VIPs) of the molecule, as the latter condition attests to the viability of the exchange-correlation potential (V{sub XC}). None of the prominently used DFT approaches show these properties: the optimized effective potential V{sub XC} based ab initio dft does. A local, range-separated hybrid potential cam-QTP-00 is introduced as the basis for a “consistent” KS DFT approach. The computed VIPs as the negative of KS eigenvalue have a mean absolute error of 0.8 eV for an extensive set of molecule's electron ionizations, including the core. Barrier heights, equilibrium geometries, and magnetic properties obtained from the potential are in good agreement with experiment. A similar accuracy with less computational efforts can be achieved by using a non-variational global hybrid variant of the QTP-00 approach.

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

Graphics processing units accelerated semiclassical initial value representation molecular dynamics

This paper presents a Graphics Processing Units (GPUs) implementation of the Semiclassical Initial Value Representation (SC-IVR) propagator for vibrational molecular spectroscopy calculations. The time-averaging formulation of the SC-IVR for power spectrum calculations is employed. Details about the GPU implementation of the semiclassical code are provided. Four molecules with an increasing number of atoms are considered and the GPU-calculated vibrational frequencies perfectly match the benchmark values. The computational time scaling of two GPUs (NVIDIA Tesla C2075 and Kepler K20), respectively, versus two CPUs (Intel Core i5 and Intel Xeon E5-2687W) and the critical issues related to the GPU implementation are discussed. The resulting reduction in computational time and power consumption is significant and semiclassical GPU calculations are shown to be environment friendly.

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

Structural, electronic, mechanical, and dynamical properties of graphene oxides: A first principles study

We report the results of a theoretical study on the structural, electronic, mechanical, and vibrational properties of some graphene oxide models (GDO, a-GMO, z-GMO, ep-GMO and mix-GMO) at ambient pressure. The calculations are based on the ab-initio plane-wave pseudo potential density functional theory, within the generalized gradient approximations for the exchange and correlation functional. The calculated values of lattice parameters, bulk modulus, and its first order pressure derivative are in good agreement with other reports. A linear response approach to the density functional theory is used to derive the phonon frequencies. We discuss the contribution of the phonons in the dynamical stability of graphene oxides and detailed analysis of zone centre phonon modes in all the above mentioned models. Our study demonstrates a wide range of energy gap available in the considered models of graphene oxide and hence the possibility of their use in nanodevices.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2014

Structural and orientation effects on electronic energy transfer between silicon quantum dots with dopants and with silver adsorbates

Starting from the atomic structure of silicon quantum dots (QDs), and utilizing ab initio electronic structure calculations within the Förster resonance energy transfer (FRET) treatment, a model has been developed to characterize electronic excitation energy transfer between QDs. Electronic energy transfer rates, K{sub EET}, between selected identical pairs of crystalline silicon quantum dots systems, either bare, doped with Al or P, or adsorbed with Ag and Ag{sub 3}, have been calculated and analyzed to extend previous work on light absorption by QDs. The effects of their size and relative orientation on energy transfer rates for each system have also been considered. Using time-dependent density functional theory and the hybrid functional HSE06, the FRET treatment was employed to model electronic energy transfer rates within the dipole-dipole interaction approximation. Calculations with adsorbed Ag show that: (a) addition of Ag increases rates up to 100 times, (b) addition of Ag{sub 3} increases rates up to 1000 times, (c) collinear alignment of permanent dipoles increases transfer rates by an order of magnitude compared to parallel orientation, and (d) smaller QD-size increases transfer due to greater electronic orbitals overlap. Calculations with dopants show that: (a) p-type and n-type dopants enhance energy transfer up to two orders of magnitude, (b) surface-doping with P and center-doping with Al show the greatest rates, and (c) K{sub EET} is largest for collinear permanent dipoles when the dopant is on the outer surface and for parallel permanent dipoles when the dopant is inside the QD.

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

Practical auxiliary basis implementation of Rung 3.5 functionals

Approximate exchange-correlation functionals for Kohn-Sham density functional theory often benefit from incorporating exact exchange. Exact exchange is constructed from the noninteracting reference system's nonlocal one-particle density matrix γ(r{sup -vector},r{sup -vector}′). Rung 3.5 functionals attempt to balance the strengths and limitations of exact exchange using a new ingredient, a projection of γ(r{sup -vector},r{sup -vector} ′) onto a semilocal model density matrix γ{sub SL}(ρ(r{sup -vector}),∇ρ(r{sup -vector}),r{sup -vector}−r{sup -vector} ′). γ{sub SL} depends on the electron density ρ(r{sup -vector}) at reference point r{sup -vector}, and is closely related to semilocal model exchange holes. We present a practical implementation of Rung 3.5 functionals, expanding the r{sup -vector}−r{sup -vector} ′ dependence of γ{sub SL} in an auxiliary basis set. Energies and energy derivatives are obtained from 3D numerical integration as in standard semilocal functionals. We also present numerical tests of a range of properties, including molecular thermochemistry and kinetics, geometries and vibrational frequencies, and bandgaps and excitation energies. Rung 3.5 functionals typically provide accuracy intermediate between semilocal and hybrid approximations. Nonlocal potential contributions from γ{sub SL} yield interesting successes and failures for band structures and excitation energies. The results enable and motivate continued exploration of Rung 3.5 functional forms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

High angular momentum states of lithium atom, studied with symmetry-projected explicitly correlated Gaussian lobe functions

Method of construction of wave functions approximating eigenfunctions of the L{sup ^2} operator is proposed for high angular momentum states of few-electron atoms. Basis functions are explicitly correlated Gaussian lobes, projected onto irreducible representations of finite point groups. Variational calculations have been carried out for the lowest states of lithium atom, with quantum number L in the range from 1 to 8. Nonrelativistic energies accurate to several dozens of nanohartree have been obtained. For 2{sup 2}P, 3{sup 2}D, and 4{sup 2}F states they agree well with the reference results. Transition frequencies have been computed and compared with available experimental data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Transport and optical properties of warm dense aluminum in the two-temperature regime: Ab initio calculation and semiempirical approximation

This work is devoted to the investigation of transport and optical properties of liquid aluminum in the two-temperature case. At first optical properties, static electrical, and thermal conductivities were obtained in the ab initio calculation which is based on the quantum molecular dynamics, density functional theory, and the Kubo-Greenwood formula. Then the semiempirical approximation was constructed based on the results of our simulation. This approximation yields the dependences σ{sub 1{sub D{sub C}}}∝1/T{sub i}{sup 0.25} and K∝T{sub e}/T{sub i}{sup 0.25} for the static electrical conductivity and thermal conductivity, respectively, for liquid aluminum at ρ = 2.70 g/cm{sup 3}, 3 kK ≤ T{sub i} ≤ T{sub e} ≤ 20 kK. Our results are well described by the Drude model with the effective relaxation time τ∝T{sub i}{sup −0.25}. We have considered a number of other models for the static electrical and thermal conductivities of aluminum, they are all reduced in the low-temperature limit to the Drude model with different expressions for the relaxation time τ. Our results are not consistent with the models in which τ∝T{sub i}{sup −1} and support the models which use the expressions with the slower decrease of the relaxation time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Classical nuclear motion coupled to electronic non-adiabatic transitions

Based on the exact factorization of the electron-nuclear wave function, we have recently proposed a mixed quantum-classical scheme [A. Abedi, F. Agostini, and E. K. U. Gross, Europhys. Lett. 106, 33001 (2014)] to deal with non-adiabatic processes. Here we present a comprehensive description of the formalism, including the full derivation of the equations of motion. Numerical results are presented for a model system for non-adiabatic charge transfer in order to test the performance of the method and to validate the underlying approximations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Band-gap corrected density functional theory calculations for InAs/GaSb type II superlattices

We performed pseudopotential based density functional theory (DFT) calculations for GaSb/InAs type II superlattices (T2SLs), with bandgap errors from the local density approximation mitigated by applying an empirical method to correct the bulk bandgaps. Specifically, this work (1) compared the calculated bandgaps with experimental data and non-self-consistent atomistic methods; (2) calculated the T2SL band structures with varying structural parameters; (3) investigated the interfacial effects associated with the no-common-atom heterostructure; and (4) studied the strain effect due to lattice mismatch between the two components. This work demonstrates the feasibility of applying the DFT method to more exotic heterostructures and defect problems related to this material system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Metal cluster structures and properties from Born-Oppenheimer molecular dynamics

Density functional theory (DFT) Born-Oppenheimer molecular dynamics (BOMD) simulations of metal clusters are presented. The calculations have been performed with the deMon2k [1] code employing all-electron basis sets and local and non-local functionals. The capability to perform reasonable long (∼ 100 ps) first-principle BOMD simulations in order to explore potential energy landscape of metallic clusters will be presented [2,3]. The evolution of the cluster structures and properties, such as polarizability and heat capacity, with temperature is discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Interference enhanced thermoelectricity in quinoid type structures

Quantum interference (QI) effects in molecular junctions may be used to obtain large thermoelectric responses. We study the electrical conductance G and the thermoelectric response of a series of molecules featuring a quinoid core using density functional theory, as well as a semi-empirical interacting model Hamiltonian describing the π-system of the molecule which we treat in the GW approximation. Molecules with a quinoid type structure are shown to have two distinct destructive QI features close to the frontier orbital energies. These manifest themselves as two dips in the transmission, that remain separated, even when either electron donating or withdrawing side groups are added. We find that the position of the dips in the transmission and the frontier molecular levels can be chemically controlled by varying the electron donating or withdrawing character of the side groups as well as the conjugation length inside the molecule. This feature results in a very high thermoelectric power factor S{sup 2}G and figure of merit ZT, where S is the Seebeck coefficient, making quinoid type molecules potential candidates for efficient thermoelectric devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Structural, electronic and optical properties of La{sub x}Sc{sub 1-x}Sb alloys

We present calculations of the structural, electronic and optic properties of LaxSc1-xSb ternary alloys for 0≤x≤1, by using the first principle full potential linear muffin-tin orbital (FPLMTO) method based on the local density approximation (LDA). the lattice constant, bulk modulus, electronic band structures, density of state and optical properties such as dielectric functions, refractive index and extinction coefficient are calculated and discussed for (x=0.25, 0.5 and 0.75). Our results agree well with the available data in the literature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

An ab initio study of multiple phonon scattering resonances in silicon germanium alloys

We have computed here phonon scattering rates and density of states in silicon germanium alloys using Green's function calculations and density functional theory. This method contrasts with the virtual crystal approximation (VCA) used in conjunction with Fermi's golden rule, which cannot capture resonance states occurring through the interaction of substitutional impurities with the host lattice. These resonances are demonstrated by density of states and scattering rate calculations in the dilute limit and show broadening as the concentration increases. Although these deviations become significant from the VCA at high frequencies, the relaxation times obtained for these phonon modes are small in both the full scattering theory and the VCA, resulting in their negligible contribution to thermal transport.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Enhanced carrier mobility and direct tunneling probability of biaxially strained Ge{sub 1−x}Sn{sub x} alloys for field-effect transistors applications

The carrier transport and tunneling capabilities of biaxially strained Ge{sub 1−x}Sn{sub x} alloys with (001), (110), and (111) orientations were comprehensively investigated and compared. The electron band structures of biaxially strained Ge{sub 1−x}Sn{sub x} alloys were calculated by the nonlocal empirical pseudopotential method and the modified virtual crystal approximation was adopted in the calculation. The electron and hole effective masses at the band edges were extracted using a parabolic line fit. It is shown that the applied biaxial strain and the high Sn composition are both helpful for the reduction of carrier effective masses, which leads to the enhanced carrier mobility and the boosted direct band-to-band-tunneling probability. Furthermore, the strain induced valance band splitting reduces the hole interband scattering, and the splitting also results in the significantly enhanced direct tunneling rate along the out-of-plane direction compared with that along the in-plane direction. The biaxially strained (111) Ge{sub 1−x}Sn{sub x} alloys exhibit the smallest band gaps compared with (001) and (110) orientations, leading to the highest in-plane and out-of-plane direct tunneling probabilities. The small effective masses on (110) and (111) planes in some strained conditions also contribute to the enhanced carrier mobility and tunneling probability. Therefore, the biaxially strained (110) and (111) Ge{sub 1−x}Sn{sub x} alloys have the potential to outperform the corresponding (001) Ge{sub 1−x}Sn{sub x} devices. It is important to optimize the applied biaxial strain, the Sn composition, and the substrate orientation for the design of high performance Ge{sub 1−x}Sn{sub x} field-effect transistors.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Discrete photon implementation for plasma simulations

The self-produced light emission from pulsed plasma discharges and its impact on plasma development are challenging to characterize through simulation and modeling, chiefly due to the large number of radiating species and limited computer memory. Often, photo-processes, such as photo-ionization or photo-emission of electrons, are implemented through over-simplifying approximations or neglected altogether. Here, a method applicable to plasma simulations is implemented in a Particle-in-Cell /Monte Carlo Collision model, which is capable of discretely tracking photons and their corresponding wavelengths. Combined with the appropriate cross sections or quantum yields, a wavelength dependent model for photo-ionization or photo-emission may be implemented. Additionally, by resolving the wavelengths of each photon, an emission spectrum for a region of interest may be generated. Simulations for a pure nitrogen environment reveal that the calculated emission profile of the second positive system agrees well with the experimental spectrum of a pulsed, nanosecond discharge in the same spectral region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Ab-initio study of thermal expansion in pure graphene

Graphene is a zero band gap semiconductor with exceptionally high thermal conductivity. The electronic properties having been studied, therole of phonon in contributing to thermal expansion, thermal conductivity and other thermodynamic properties, is required to be investigated. This paper focuses more on thermal expansion. Some others results like phonon dispersion, Grüneisenparameters and bulk modulus,which are essential to estimation of thermal expansion, are also presented. The dynamical matrix was calculated using VASP code using both DFT and DFPT and the phonon frequencies were calculated using phonopy code under harmonic approximation. The linear thermal expansion coefficient of graphene is found to be strongly dependent on temperature but remains negative upto 470 K and positive thereafter, with a room temperature value of −1.44×10{sup −6}. The negative expansion coefficient is very interesting and is found to be in conformity with experimental as well as with recent theoretical estimates. There is only qualitative agreement of our results with experimental data and motivates further investigation, primarily on the high negative values of Grüneisen parameters.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2016

On the existence of the optimal order for wavefunction extrapolation in Born-Oppenheimer molecular dynamics

Wavefunction extrapolation greatly reduces the number of self-consistent field (SCF) iterations and thus the overall computational cost of Born-Oppenheimer molecular dynamics (BOMD) that is based on the Kohn–Sham density functional theory. Going against the intuition that the higher order of extrapolation possesses a better accuracy, we demonstrate, from both theoretical and numerical perspectives, that the extrapolation accuracy firstly increases and then decreases with respect to the order, and an optimal extrapolation order in terms of minimal number of SCF iterations always exists. We also prove that the optimal order tends to be larger when using larger MD time steps or more strict SCF convergence criteria. By example BOMD simulations of a solid copper system, we show that the optimal extrapolation order covers a broad range when varying the MD time step or the SCF convergence criterion. Therefore, we suggest the necessity for BOMD simulation packages to open the user interface and to provide more choices on the extrapolation order. Another factor that may influence the extrapolation accuracy is the alignment scheme that eliminates the discontinuity in the wavefunctions with respect to the atomic or cell variables. We prove the equivalence between the two existing schemes, thus the implementation of either of them does not lead to essential difference in the extrapolation accuracy.

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

Note: A pure-sampling quantum Monte Carlo algorithm with independent Metropolis

Recently, Ospadov and Rothstein published a pure-sampling quantum Monte Carlo algorithm (PSQMC) that features an auxiliary Path Z that connects the midpoints of the current and proposed Paths X and Y, respectively. When sufficiently long, Path Z provides statistical independence of Paths X and Y. Under those conditions, the Metropolis decision used in PSQMC is done without any approximation, i.e., not requiring microscopic reversibility and without having to introduce any G(x → x′; τ) factors into its decision function. This is a unique feature that contrasts with all competing reptation algorithms in the literature. An example illustrates that dependence of Paths X and Y has adverse consequences for pure sampling.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗