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,747 records · Page 97
All sourcesNASA NTRSDOE OSTI
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↗
DOE OSTI2022

DMC-ICE13 : Ambient and high pressure polymorphs of ice from diffusion Monte Carlo and density functional theory

Ice is one of the most important and interesting molecular crystals, exhibiting a rich and evolving phase diagram. Recent discoveries mean that there are now 20 distinct polymorphs; a structural diversity that arises from a delicate interplay of hydrogen bonding and van der Waals dispersion forces. This wealth of structures provides a stern test of electronic structure theories, with Density Functional Theory (DFT) often not able to accurately characterize the relative energies of the various ice polymorphs. Thanks to recent advances that enable the accurate and efficient treatment of molecular crystals with Diffusion Monte Carlo (DMC), we present here the DMC-ICE13 dataset; a dataset of lattice energies of 13 ice polymorphs. This dataset encompasses the full structural complexity found in the ambient and high-pressure molecular ice polymorphs, and when experimental reference energies are available, our DMC results deliver sub-chemical accuracy. Using this dataset, we then perform an extensive benchmark of a broad range of DFT functionals. Of the functionals considered, revPBE-D3 and RSCAN reproduce reference absolute lattice energies with the smallest error, while optB86b-vdW and SCAN+rVV10 have the best performance on the relative lattice energies. Our results suggest that a single functional achieving reliable performance for all phases is still missing, and that care is needed in the selection of the most appropriate functional for the desired application. The insights obtained here may also be relevant to liquid water and other hydrogen-bonded and dispersion-bonded molecular crystals.

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

Learning physics-based reduced-order models from data using nonlinear manifolds

Here we present a novel method for learning reduced-order models of dynamical systems using nonlinear manifolds. First, we learn the manifold by identifying nonlinear structure in the data through a general representation learning problem. The proposed approach is driven by embeddings of low-order polynomial form. A projection onto the nonlinear manifold reveals the algebraic structure of the reduced-space system that governs the problem of interest. The matrix operators of the reduced-order model are then inferred from the data using operator inference. Numerical experiments on a number of nonlinear problems demonstrate the generalizability of the methodology and the increase in accuracy that can be obtained over reduced-order modeling methods that employ a linear subspace approximation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2011

Transmission characteristics of circular metallic waveguides for terahertz radiation

Transmission characteristics of oversized circular metallic waveguides excited by linearly polarised Gaussian laser beams in the terahertz range (4 - 28 THz) are studied theoretically and experimentally. Calculating the transmission characteristics, we have determined the conditions of applicability of the method of the eigenoscillations in the approximations of a real metal by an ideal metal or dielectric, depending on the transmitted radiation frequency. The existence of the transition region is established in the behaviour of the electrodynamic properties of metallic waveguides in the frequency range of 7.5 - 15 THz.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Quenching of highly rotationally excited HCl in collisions with He

We report rotational quenching cross sections and rate coefficients of HCl due to collisions with He. The close-coupling method and the coupled-states approximation are applied in quantum-mechanical scattering calculations of state-to-state cross sections for HCl with initial rotational levels up to j = 20 for kinetic energies from 10{sup –5} to 15,000 cm{sup –1}. State-to-state rate coefficients for temperatures between 0.1 and 3000 K are also presented. Comparison of the present rate coefficients with previous results reported in the literature for lowly excited rotational levels shows reasonable agreement. Small differences are attributed to the differences in the interaction potential energy surfaces. The uncertainty in the computed cross sections and rate coefficients is estimated by varying the potential well depth. Applications of current results to astrophysical systems are also briefly discussed.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2012

BAYESIAN ANGULAR POWER SPECTRUM ANALYSIS OF INTERFEROMETRIC DATA

We present a Bayesian angular power spectrum and signal map inference engine which can be adapted to interferometric observations of anisotropies in the cosmic microwave background (CMB), 21 cm emission line mapping of galactic brightness fluctuations, or 21 cm absorption line mapping of neutral hydrogen in the dark ages. The method uses Gibbs sampling to generate a sampled representation of the angular power spectrum posterior and the posterior of signal maps given a set of measured visibilities in the uv-plane. We use a mock interferometric CMB observation to demonstrate the validity of this method in the flat-sky approximation when adapted to take into account arbitrary coverage of the uv-plane, mode-mode correlations due to observations on a finite patch, and heteroschedastic visibility errors. The computational requirements scale as O(n{sub p} log n{sub p}) where n{sub p} measures the ratio of the size of the detector array to the inter-detector spacing, meaning that Gibbs sampling is a promising technique for meeting the data analysis requirements of future cosmology missions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

S 8 ring content in liquid sulfur determined from x-ray diffraction

The S 8 ring content in liquid sulfur has been determined from the 4.47 Å peak in the x-ray pair distribution function over the temperature range −55 °C –330 °C. This peak represents the characteristic 3rd and 4th nearest atom-atom distances across an S 8 -ring. Three procedures were used to calculate the S 8 -ring fraction; scaling to crystalline β-S, a linear combination fit and a direct integration method. All methods gave similar values, indicating that approximately one third of the S 8 -rings break up upon melting into a distribution of short chains with an average length of 5 ± 2 atoms. Our results contradict established models of low temperature sulfur that indicate below the λ-transition the liquid comprises almost entirely of S8-rings. In addition, we find the S 8 -ring fraction increases to 76 ± 10% at −30 °C in the supercooled liquid, and decreases to 23% ± 3% at 330 °C in the high temperature melt. First principle molecular dynamics calculations confirm the onset of dynamic ring distortion via broad bond angle distributions, associated with the interstitial region between the first and second peaks in the pair distribution function. Our results are consistent with the presence of a propensity of short chains in the low temperature liquid, that are required for the bond switching mechanism responsible for the polymerization process to occur at the λ-transition temperature, T λ .

lambda transition↗
DOE OSTI2019

Modeling the Kondo effect of a magnetic atom adsorbed on graphene

The Kondo effect due to a magnetic atom adsorbed on graphene is investigated using two theoretical approaches, including the non-equilibrium Green's function method within the slave-boson mean-field approximation and a newly developed tight-binding (TB) model with an effectively infinite Hubbard U term. Both methods reveal the presence of a Kondo peak in the local density of state (LDOS) of graphene near the Fermi level. A sharp Kondo peak is only observed in the odd-neighboring C sites of the C atom directly below the magnetic atom placed in top position. The peak intensity is found to decay quickly with respect to the distance between the C site and the magnetic atom. A theoretical model of scanning tunneling microscopy (STM) of graphene is presented as a means to identify the manifestation of the Kondo effect via direct topographic STM measurements. STM simulations show that in the proximity of the magnetic atom, only the C atoms of one sublattice are visible at low bias potential and thus a triangular lattice can be seen, in agreement with recent experiments. Further away from the magnetic atom, the Kondo effect dies down and eventually both sublattices are visible, leading to the recovery of the hexagonal lattice. Furthermore, this work not only provides a new TB method to study the Kondo effect in graphene, but also presents both scanning tunneling spectroscopy and microscopy fingerprints of the Kondo effect to facilitate its experimental verification.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Inelastic rate coefficients for collisions of C 4 H - with H 2

Carbon-chain anions were recently detected in the interstellar medium. These very reactive species are used as tracers of the physical and chemical conditions in a variety of astrophysical environments. However, the local thermodynamic equilibrium conditions are generally not fulfilled in these environments. Therefore, collisional as well as radiative rates are needed to accurately model the observed emission lines. We determine in this work the state-to-state rate coefficients of C 4 H- in collision with both ortho- and para-H 2 . A new ab initio 4D potential energy surface was computed using explicitly correlated coupled-cluster procedures. This surface was then employed to determine rotational excitation and de-excitation cross-sections and rate coefficients for the first 21 rotational levels (up to rotational level j 1 = 20) using the close-coupling method, while the coupled-state approximation was used to extend the calculations up to j 1 = 30. State-to-state rate coefficients were obtained for the temperature range 2–100K. The differences between the ortho- and para-H 2 rate coefficients are found to be small.

79 ASTRONOMY AND ASTROPHYSICS↗