Engineering Papers⌕ Search

SEARCH · Engineering Papers

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,729 records · Page 96
All sourcesNASA NTRSDOE OSTI
DOE OSTI2011

{Lambda}K*{Lambda}(1116) photoproduction and nucleon resonances

In this presentation, we report our recent studies on the {Lambda}K*{Lambda}(1116) photoproduction off the proton target, using the tree-level Born approximation, via the effective Lagrangian approach. In addition, we include the nine (three- or four-star confirmed) nucleon resonances below the threshold {radical}(s{sub th}){approx_equal}2008 MeV, to interpret the discrepancy between the experiment and previous theoretical studies, in the vicinity of the threshold region. From the numerical studies, we observe that the S{sub 11}(1535) and S11(1650) play an important role for the cross-section enhancement near the {radical}(s{sub th}). It also turns out that, in order to reproduce the data, we have the vector coupling constants gK*S{sub 11}(1535){Lambda} = (7.0{approx}9.0) and gK*S{sub 11}(1650){Lambda} (5.0{approx}6.0).

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2011

Spin-orbit relativistic long-range corrected time-dependent density functional theory for investigating spin-forbidden transitions in photochemical reactions

A long-range corrected (LC) time-dependent density functional theory (TDDFT) incorporating relativistic effects with spin-orbit couplings is presented. The relativistic effects are based on the two-component zeroth-order regular approximation Hamiltonian. Before calculating the electronic excitations, we calculated the ionization potentials (IPs) of alkaline metal, alkaline-earth metal, group 12 transition metal, and rare gas atoms as the minus orbital (spinor) energies on the basis of Koopmans' theorem. We found that both long-range exchange and spin-orbit coupling effects are required to obtain Koopmans' IPs, i.e., the orbital (spinor) energies, quantitatively in DFT calculations even for first-row transition metals and systems containing large short-range exchange effects. We then calculated the valence excitations of group 12 transition metal atoms and the Rydberg excitations of rare gas atoms using spin-orbit relativistic LC-TDDFT. We found that the long-range exchange and spin-orbit coupling effects significantly contribute to the electronic spectra of even light atoms if the atoms have low-lying excitations between orbital spinors of quite different electron distributions.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2012

Origin-independent calculation of quadrupole intensities in X-ray spectroscopy

For electronic excitations in the ultraviolet and visible range of the electromagnetic spectrum, the intensities are usually calculated within the dipole approximation, which assumes that the oscillating electric field is constant over the length scale of the transition. For the short wavelengths used in hard X-ray spectroscopy, the dipole approximation may not be adequate. In particular, for metal K-edge X-ray absorption spectroscopy (XAS), it becomes necessary to include higher-order contributions. In quantum-chemical approaches to X-ray spectroscopy, these so-called quadrupole intensities have so far been calculated by including contributions depending on the square of the electric-quadrupole and magnetic-dipole transition moments. However, the resulting quadrupole intensities depend on the choice of the origin of the coordinate system. Here, we show that for obtaining an origin-independent theory, one has to include all contributions that are of the same order in the wave vector consistently. This leads to two additional contributions depending on products of the electric-dipole and electric-octupole and of the electric-dipole and magnetic-quadrupole transition moments, respectively. We have implemented such an origin-independent calculation of quadrupole intensities in XAS within time-dependent density-functional theory, and demonstrate its usefulness for the calculation of metal and ligand K-edge XAS spectra of transition metal complexes.

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

On the alternatives for bath correlators and spectral densities from mixed quantum-classical simulations

We investigate on the procedure of extracting a 'spectral density' from mixed QM/MM calculations and employing it in open quantum systems models. In particular, we study the connection between the energy gap correlation function extracted from ground state QM/MM and the bath spectral density used as input in open quantum system approaches. We introduce a simple model which can give intuition on when the ground state QM/MM propagation will give the correct energy gap. We also discuss the role of higher order correlators of the energy-gap fluctuations which can provide useful information on the bath. Further, various semiclassical corrections to the spectral density, are applied and investigated. Finally, we apply our considerations to the photosynthetic Fenna-Matthews-Olson complex. For this system, our results suggest the use of the Harmonic prefactor for the spectral density rather than the Standard one, which was employed in the simulations of the system carried out to date.

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

Irradiation-induced composition patterns in binary solid solutions

A theoretical/computational model for the irradiation-driven compositional instabilities in binary solid solutions has been developed. The model is suitable for investigating the behavior of structural alloys and metallic nuclear fuels in a reactor environment as well as the response of alloy thin films to ion beam irradiation. The model is based on a set of reaction-diffusion equations for the dynamics of vacancies, interstitials, and lattice atoms under irradiation. The dynamics of these species includes the stochastic generation of defects by collision cascades as well as the defect reactions and diffusion. The atomic fluxes in this model are derived based on the transitions of lattice defects. The set of reaction-diffusion equations are stiff, hence a stiffly stable method, also known as the Gear method, has been used to numerically approximate the equations. For the Cu-Au alloy in the solid solution regime, the model results demonstrate the formation of compositional patterns under high-temperature particle irradiation, with Fourier space properties (Fourier spectrum, average wavelength, and wavevector) depending on the cascade damage characteristics, average composition, and irradiation temperature.

36 MATERIALS SCIENCE↗
DOE OSTI2013

Heteroepitaxial growth and surface structure of L1 0 -MnGa(111) ultra-thin films on GaN(0001)

L1 0 -structured MnGa(111) ultra-thin films were heteroepitaxially grown on GaN(0001) under lightly Mn-rich conditions using molecular beam epitaxy. Room-temperature scanning tunneling microscopy (STM) investigations reveal smooth terraces and angular step edges, with the surface structure consisting primarily of a 2x2 reconstruction along with small patches of 1x2. Theoretical calculations were undertaken using density functional theory, and the simulated STM images were calculated using the Tersoff-Hamann approximation, demonstrating that a stoichiometric 1x2 and a Mn-rich 2x2 surface structure give the best agreement with the observed experimental images.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

First principles molecular dynamics without self-consistent field optimization

We present a first principles molecular dynamics approach that is based on time-reversible extended Lagrangian Born-Oppenheimer molecular dynamics [A. M. N. Niklasson, Phys. Rev. Lett. 100, 123004 (2008)] in the limit of vanishing self-consistent field optimization. The optimization-free dynamics keeps the computational cost to a minimum and typically provides molecular trajectories that closely follow the exact Born-Oppenheimer potential energy surface. Only one single diagonalization and Hamiltonian (or Fockian) construction are required in each integration time step. The proposed dynamics is derived for a general free-energy potential surface valid at finite electronic temperatures within hybrid density functional theory. Even in the event of irregular functional behavior that may cause a dynamical instability, the optimization-free limit represents a natural starting guess for force calculations that may require a more elaborate iterative electronic ground state optimization. Our optimization-free dynamics thus represents a flexible theoretical framework for a broad and general class of ab initio molecular dynamics simulations.

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

Energy density functionals from the strong-coupling limit applied to the anions of the He isoelectronic series

Anions and radicals are important for many applications including environmental chemistry, semiconductors, and charge transfer, but are poorly described by the available approximate energy density functionals. Here we test an approximate exchange-correlation functional based on the exact strong-coupling limit of the Hohenberg-Kohn functional on the prototypical case of the He isoelectronic series with varying nuclear charge Z < 2, which includes weakly bound negative ions and a quantum phase transition at a critical value of Z, representing a big challenge for density functional theory. We use accurate wavefunction calculations to validate our results, comparing energies and Kohn-Sham potentials, thus also providing useful reference data close to and at the quantum phase transition. We show that our functional is able to bind H{sup −} and to capture in general the physics of loosely bound anions, with a tendency to strongly overbind that can be proven mathematically. We also include corrections based on the uniform electron gas which improve the results.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
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↗