SEARCH · Search NASA

Results for “EXTRAPOLATION”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 919 records · Page 51

All sources NASA NTRS DOE OSTI

DOE OSTI2014

Multilevel characteristics and memory mechanisms for nonvolatile memory devices based on CuInS{sub 2} quantum dot-polymethylmethacrylate nanocomposites

Nonvolatile memory devices based on CuInS{sub 2} (CIS) quantum dots (QDs) embedded in a polymethylmethacrylate (PMMA) layer were fabricated using spin-coating method. The memory window widths of the capacitance-voltage (C-V) curves for the Al/CIS QDs embedded in PMMA layer/p-Si devices were 0.3, 0.6, and 1.0 V for sweep voltages of ±3, ±5, and ±7 V, respectively. Capacitance-cycle data demonstrated that the charge-trapping capability of the devices with an ON/OFF ratio value of 2.81 × 10{sup −10} was maintained for 8 × 10{sup 3} cycles without significant degradation and that the extrapolation of the ON/OFF ratio value to 1 × 10{sup 6} cycles converged to 2.40 × 10{sup −10}, indicative of the good stability of the devices. The memory mechanisms for the devices are described on the basis of the C-V curves and the energy-band diagrams.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

DOE OSTI2014

Modeling of molecular nitrogen collisions and dissociation processes for direct simulation Monte Carlo

The Direct Simulation Monte Carlo (DSMC) method typically used for simulating hypersonic Earth re-entry flows requires accurate total collision cross sections and reaction probabilities. However, total cross sections are often determined from extrapolations of relatively low-temperature viscosity data, so their reliability is unknown for the high temperatures observed in hypersonic flows. Existing DSMC reaction models accurately reproduce experimental equilibrium reaction rates, but the applicability of these rates to the strong thermal nonequilibrium observed in hypersonic shocks is unknown. For hypersonic flows, these modeling issues are particularly relevant for nitrogen, the dominant species of air. To rectify this deficiency, the Molecular Dynamics/Quasi-Classical Trajectories (MD/QCT) method is used to accurately compute collision and reaction cross sections for the N{sub 2}({sup 1}Σ{sub g}{sup +})-N{sub 2}({sup 1}Σ{sub g}{sup +}) collision pair for conditions expected in hypersonic shocks using a new potential energy surface developed using a ReaxFF fit to recent advanced ab initio calculations. The MD/QCT-computed reaction probabilities were found to exhibit better physical behavior and predict less dissociation than the baseline total collision energy reaction model for strong nonequilibrium conditions expected in a shock. The MD/QCT reaction model compared well with computed equilibrium reaction rates and shock-tube data. In addition, the MD/QCT-computed total cross sections were found to agree well with established variable hard sphere total cross sections.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2014

In search of the X{sub 2}BO and X{sub 2}BS (X = H, F) free radicals: Ab initio studies of their spectroscopic signatures

The F{sub 2}BO free radical is a known, although little studied, species but similar X{sub 2}BY (X = H, D, F; Y = O, S) molecules are largely unknown. High level ab initio methods have been used to predict the molecular structures, vibrational frequencies (in cm{sup −1}), and relative energies of the ground and first two excited electronic states of these free radicals, as an aid to their eventual spectroscopic identification. The chosen theoretical methods and basis sets were tested on F{sub 2}BO and found to give good agreement with the known experimental quantities. In particular, complete basis set extrapolations of coupled-cluster single and doubles with perturbative triple excitations/aug-cc-pVXZ (X = 3, 4, 5) energies gave excellent electronic term values, due to small changes in geometry between states and the lack of significant multireference character in the wavefunctions. The radicals are found to have planar C{sub 2v} geometries in the X{sup ~2}B{sub 2} ground state, the low-lying A{sup ~2}B{sub 1} first excited state, and the higher B{sup ~2}A{sub 1} state. Some of these radicals have very small ground state dipole moments hindering microwave measurements. Infrared studies in matrices or in the gas phase may be possible although the fundamentals of H{sub 2}BO and H{sub 2}BS are quite weak. The most promising method of identifying these species in the gas phase appears to be absorption or laser-induced fluorescence spectroscopy through the allowed B{sup ~}-X{sup ~} transitions which occur in the visible-near UV region of the electromagnetic spectrum. The ab initio results have been used to calculate the Franck-Condon profiles of the absorption and emission spectra, and the rotational structure of the B{sup ~}-X{sup ~}0{sub 0}{sup 0} bands has been simulated. The calculated single vibronic level emission spectra provide a unique, readily recognizable fingerprint of each particular radical, facilitating the experimental identification of new X{sub 2}BY species in the gas phase.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2014

Unbiased reduced density matrices and electronic properties from full configuration interaction quantum Monte Carlo

Properties that are necessarily formulated within pure (symmetric) expectation values are difficult to calculate for projector quantum Monte Carlo approaches, but are critical in order to compute many of the important observable properties of electronic systems. Here, we investigate an approach for the sampling of unbiased reduced density matrices within the full configuration interaction quantum Monte Carlo dynamic, which requires only small computational overheads. This is achieved via an independent replica population of walkers in the dynamic, sampled alongside the original population. The resulting reduced density matrices are free from systematic error (beyond those present via constraints on the dynamic itself) and can be used to compute a variety of expectation values and properties, with rapid convergence to an exact limit. A quasi-variational energy estimate derived from these density matrices is proposed as an accurate alternative to the projected estimator for multiconfigurational wavefunctions, while its variational property could potentially lend itself to accurate extrapolation approaches in larger systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2015

A general non-Abelian density matrix renormalization group algorithm with application to the C{sub 2} dimer

We extend our previous work [S. Sharma and G. K.-L. Chan, J. Chem. Phys. 136, 124121 (2012)], which described a spin-adapted (SU(2) symmetry) density matrix renormalization group algorithm, to additionally utilize general non-Abelian point group symmetries. A key strength of the present formulation is that the requisite tensor operators are not hard-coded for each symmetry group, but are instead generated on the fly using the appropriate Clebsch-Gordan coefficients. This allows our single implementation to easily enable (or disable) any non-Abelian point group symmetry (including SU(2) spin symmetry). We use our implementation to compute the ground state potential energy curve of the C{sub 2} dimer in the cc-pVQZ basis set (with a frozen-core), corresponding to a Hilbert space dimension of 10{sup 12} many-body states. While our calculated energy lies within the 0.3 mE{sub h} error bound of previous initiator full configuration interaction quantum Monte Carlo and correlation energy extrapolation by intrinsic scaling calculations, our estimated residual error is only 0.01 mE{sub h}, much more accurate than these previous estimates. Due to the additional efficiency afforded by the algorithm, the excitation energies (T{sub e}) of eight lowest lying excited states: a{sup 3}Π{sub u}, b{sup 3}Σ{sub g}{sup −}, A{sup 1}Π{sub u}, c{sup 3}Σ{sub u}{sup +}, B{sup 1}Δ{sub g}, B{sup ′1}Σ{sub g}{sup +}, d{sup 3}Π{sub g}, and C{sup 1}Π{sub g} are calculated, which agree with experimentally derived values to better than 0.06 eV. In addition, we also compute the potential energy curves of twelve states: the three lowest levels for each of the irreducible representations {sup 1}Σ{sub g}{sup +}, {sup 1}Σ{sub u}{sup +}, {sup 1}Σ{sub g}{sup −}, and {sup 1}Σ{sub u}{sup −}, to an estimated accuracy of 0.1 mE{sub h} of the exact result in this basis.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2015

Towards understanding intrinsic degradation and breakdown mechanisms in SiOCH low-k dielectrics

The degradation and breakdown mechanisms of a SiOCH low-k material with k = 2.3 (25% porosity) and thicknesses ranging from 90 nm to 20 nm were investigated. By combining the time dependent dielectric breakdown data at positive/negative bias stress with the thickness scaling results, dielectric failure is proven to be intrinsic and not influenced by copper drift or metal barrier deposition induced dielectric damage. It is shown that stress induced leakage current (SILC) can be used as a measure of dielectric degradation. Therefore, low field lifetimes can be safely estimated using SILC extrapolation. Based on our results, both the impact damage model and the power law model have a good accuracy for low field lifetime prediction. Recovery and anneal experiments are used to study the physical nature causing the observed negative flatband voltage shifts in our metal-insulator-semiconductor planar capacitor structures, where hydrogen induced unstable fast and slow donor type interface states are hypothesized to be the root cause of the observed shifts. We suggest that atomic hydrogen is released from the dielectric during electron injection and migrates to the interfacial region. Our model is further supported by an observed irreversible SILC change during the recovery and anneal studies. The degradation mechanism proposed in this work, supported by the low field lifetime data, provides a feasible explanation for intrinsic low-k dielectric failure.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2015

Basis convergence of range-separated density-functional theory

Range-separated density-functional theory (DFT) is an alternative approach to Kohn-Sham density-functional theory. The strategy of range-separated density-functional theory consists in separating the Coulomb electron-electron interaction into long-range and short-range components and treating the long-range part by an explicit many-body wave-function method and the short-range part by a density-functional approximation. Among the advantages of using many-body methods for the long-range part of the electron-electron interaction is that they are much less sensitive to the one-electron atomic basis compared to the case of the standard Coulomb interaction. Here, we provide a detailed study of the basis convergence of range-separated density-functional theory. We study the convergence of the partial-wave expansion of the long-range wave function near the electron-electron coalescence. We show that the rate of convergence is exponential with respect to the maximal angular momentum L for the long-range wave function, whereas it is polynomial for the case of the Coulomb interaction. We also study the convergence of the long-range second-order Møller-Plesset correlation energy of four systems (He, Ne, N{sub 2}, and H{sub 2}O) with cardinal number X of the Dunning basis sets cc − p(C)V XZ and find that the error in the correlation energy is best fitted by an exponential in X. This leads us to propose a three-point complete-basis-set extrapolation scheme for range-separated density-functional theory based on an exponential formula.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2015

Prestack reverse time migration for 3D marine reflection seismic data

Prestack reverse time migration (RTM) is a method for imaging the subsurface using the inner product of wavefield extrapolation in shot domain and in receiver domain. It is well known that RTM is better for preserving amplitudes and phases than other prestack migrations. Since 3D seismic data is huge data volume and it needs heavy computing works, it requires parallel computing in order to have a meaningful depth image of the 3D subsurface. We implemented a parallelized version of 3D RTM for prestack depth migration. The results of numerical example for 3D SEG/EAGE salt model showed good agreement with the original geological model. We applied RTM to offshore 3D seismic reflection data. The study area is 12 × 25 km with 120 survey lines. Shot and receiver spacing is 25 m and 12.5 m. The line spacing is 100 m. Shot gathers were preprocessed to enhance signal to noise ratio and velocity model was calculated from conventional stack velocity. Both of them were used to obtain 3D image using RTM. The results show reasonable subsurface image.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2015

A preliminary prediction of the (ν{sub 2} + ν{sub 4}) combination band of {sup 34}SF{sub 6}

Sulfur hexafluoride appears to be a species of growing importance in the field of atmospheric physics and chemistry as a pollutant that can contribute to the greenhouse effect. The infrared spectroscopy of SF6 is essential for quantitative measurements in the atmosphere, is still insufficiently known. In this work, we present a new contribution to this topic concerning (ν{sub 2}+ν{sub 4}) combination band. The prediction of this band is based on the tensorial formalism and vibrational extrapolation methods developed in Dijon. The effective Hamiltonian contains contribution from the ground state level GS, ν{sub 2} and ν{sub 4} witch are known from previous studies. The analysis has been performed using nonlinear least squares fit procedure included in the XTDS program to Jup=95.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2015

Benchmark of numerical tools simulating beam propagation and secondary particles in ITER NBI

Injection of high energy beams of neutral particles is a method for plasma heating in fusion devices. The ITER injector, and its prototype MITICA (Megavolt ITER Injector and Concept Advancement), are large extrapolations from existing devices: therefore numerical modeling is needed to set thermo-mechanical requirements for all beam-facing components. As the power and charge deposition originates from several sources (primary beam, co-accelerated electrons, and secondary production by beam-gas, beam-surface, and electron-surface interaction), the beam propagation along the beam line is simulated by comprehensive 3D models. This paper presents a comparative study between two codes: BTR has been used for several years in the design of the ITER HNB/DNB components; SAMANTHA code was independently developed and includes additional phenomena, such as secondary particles generated by collision of beam particles with the background gas. The code comparison is valuable in the perspective of the upcoming experimental operations, in order to prepare a reliable numerical support to the interpretation of experimental measurements in the beam test facilities. The power density map calculated on the Electrostatic Residual Ion Dump (ERID) is the chosen benchmark, as it depends on the electric and magnetic fields as well as on the evolution of the beam species via interaction with the gas. Finally the paper shows additional results provided by SAMANTHA, like the secondary electrons produced by volume processes accelerated by the ERID fringe-field towards the Cryopumps.

43 PARTICLE ACCELERATORS↗

DOE OSTI2015

Electromagnetic effects on turbulent transport in high-performance ASDEX Upgrade discharges

Modern tokamak H-mode discharges routinely operate at high plasma beta. Dedicated experiments performed on multiple machines measure contradicting dependence of the plasma confinement on this important parameter. In view of designing high-performance scenarios for next-generation devices like ITER, a fundamental understanding of the involved physics is crucial. Theoretical results—most of which have been obtained for simplified setups—indicate that increased beta does not only modify the characteristics of microturbulence but also potentially introduces fundamentally new physics. Empowered by highly accurate measurements at ASDEX Upgrade, the GENE turbulence code is used to perform a comprehensive gyrokinetic study of dedicated H-Mode plasmas. We find the stabilization of ion-temperature-gradient driven turbulence to be the most pronounced beta effect in these experimentally relevant cases. The resulting beta-improved core confinement should thus be considered for extrapolations to future machines.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2015

Temperature dependence of the electronic structure of semiconductors and insulators

The renormalization of electronic eigenenergies due to electron-phonon coupling (temperature dependence and zero-point motion effect) is sizable in many materials with light atoms. This effect, often neglected in ab initio calculations, can be computed using the perturbation-based Allen-Heine-Cardona theory in the adiabatic or non-adiabatic harmonic approximation. After a short description of the recent progresses in this field and a brief overview of the theory, we focus on the issue of phonon wavevector sampling convergence, until now poorly understood. Indeed, the renormalization is obtained numerically through a slowly converging q-point integration. For non-zero Born effective charges, we show that a divergence appears in the electron-phonon matrix elements at q → Γ, leading to a divergence of the adiabatic renormalization at band extrema. This problem is exacerbated by the slow convergence of Born effective charges with electronic wavevector sampling, which leaves residual Born effective charges in ab initio calculations on materials that are physically devoid of such charges. Here, we propose a solution that improves this convergence. However, for materials where Born effective charges are physically non-zero, the divergence of the renormalization indicates a breakdown of the adiabatic harmonic approximation, which we assess here by switching to the non-adiabatic harmonic approximation. Also, we study the convergence behavior of the renormalization and develop reliable extrapolation schemes to obtain the converged results. Finally, the adiabatic and non-adiabatic theories, with corrections for the slow Born effective charge convergence problem (and the associated divergence) are applied to the study of five semiconductors and insulators: α-AlN, β-AlN, BN, diamond, and silicon. For these five materials, we present the zero-point renormalization, temperature dependence, phonon-induced lifetime broadening, and the renormalized electronic band structure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2015

Effect of Sn doping on nonlinear optical properties of quaternary Se-Sn-(Bi,Te) chalcogenide thin films

The aim of this work is to report the effect of Sn doping on the third order nonlinear optical properties of chalcogenide Se{sub 84-x}Te{sub 15}Bi{sub 1.0}Sn{sub x} thin films. Melt quenching technique has been used for the preparation of bulk chalcogenide glasses. Thin films of the studied composition are deposited on cleaned glass substrate by thermal evaporation technique. Optical band gap (E{sub g}) is calculated by using Tauc extrapolation method and is found to increase from 1.27 eV to 1.64 eV with the incorporation of Sn content. Stryland approach is utilized for the calculation of two photon absorption coefficient (β{sub 2}). The nonlinear refractive index (n{sub 2}) and third order susceptibility (χ{sup (3}) are calculated using Tichy and Ticha approach. The result shows that nonlinear refractive index (n{sub 2}) follows the same trend as that of linear refractive index (n). The values of n{sub 2} of studied composition as compared to pure silica are 1000-5000 times higher.

36 MATERIALS SCIENCE↗

DOE OSTI2015

Accurate and precise determination of critical properties from Gibbs ensemble Monte Carlo simulations

Since the seminal paper by Panagiotopoulos [Mol. Phys. 61, 813 (1997)], the Gibbs ensemble Monte Carlo (GEMC) method has been the most popular particle-based simulation approach for the computation of vapor–liquid phase equilibria. However, the validity of GEMC simulations in the near-critical region has been questioned because rigorous finite-size scaling approaches cannot be applied to simulations with fluctuating volume. Valleau [Mol. Simul. 29, 627 (2003)] has argued that GEMC simulations would lead to a spurious overestimation of the critical temperature. More recently, Patel et al. [J. Chem. Phys. 134, 024101 (2011)] opined that the use of analytical tail corrections would be problematic in the near-critical region. To address these issues, we perform extensive GEMC simulations for Lennard-Jones particles in the near-critical region varying the system size, the overall system density, and the cutoff distance. For a system with N = 5500 particles, potential truncation at 8σ and analytical tail corrections, an extrapolation of GEMC simulation data at temperatures in the range from 1.27 to 1.305 yields T{sub c} = 1.3128 ± 0.0016, ρ{sub c} = 0.316 ± 0.004, and p{sub c} = 0.1274 ± 0.0013 in excellent agreement with the thermodynamic limit determined by Potoff and Panagiotopoulos [J. Chem. Phys. 109, 10914 (1998)] using grand canonical Monte Carlo simulations and finite-size scaling. Critical properties estimated using GEMC simulations with different overall system densities (0.296 ≤ ρ{sub t} ≤ 0.336) agree to within the statistical uncertainties. For simulations with tail corrections, data obtained using r{sub cut} = 3.5σ yield T{sub c} and p{sub c} that are higher by 0.2% and 1.4% than simulations with r{sub cut} = 5 and 8σ but still with overlapping 95% confidence intervals. In contrast, GEMC simulations with a truncated and shifted potential show that r{sub cut} = 8σ is insufficient to obtain accurate results. Additional GEMC simulations for hard-core square-well particles with various ranges of the attractive well and for n-decane molecules represented by the TraPPE force field yield data that support the trends observed for Lennard-Jones particles. The finite-size dependence of the critical properties obtained from GEMC simulations is significantly smaller than those from grand-canonical ensemble simulations. Thus, when resources are not available for a rigorous finite-size scaling study, GEMC simulations provide a straightforward route to determine fairly accurate critical properties using relatively small system sizes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2015

The benzene⋯naphthalene complex: A more challenging system than the benzene dimer for newly developed computational methods

High-level coupled cluster singles, doubles, and perturbative triples [CCSD(T)] computations with up to the aug-cc-pVQZ basis set (1924 basis functions) and various extrapolations toward the complete basis set (CBS) limit are presented for the sandwich, T-shaped, and parallel-displaced benzene⋯naphthalene complex. Using the CCSD(T)/CBS interaction energies as a benchmark, the performance of some newly developed wave function and density functional theory methods has been evaluated. The best performing methods were found to be the dispersion-corrected PBE0 functional (PBE0-D3) and spin-component scaled zeroth-order symmetry-adapted perturbation theory (SCS-SAPT0). The success of SCS-SAPT0 is very encouraging because it provides one method for energy component analysis of π-stacked complexes with 200 atoms or more. Most newly developed methods do, however, overestimate the interaction energies. The results of energy component analysis show that interaction energies are overestimated mainly due to the overestimation of dispersion energy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2015

Resonance strength measurement at astrophysical energies: The {sup 17}O(p,α){sup 14}N reaction studied via Trojan Horse Method

In recent years, the Trojan Horse Method (THM) has been used to investigate the low-energy cross sections of proton-induced reactions on {sup 17}O nuclei, overcoming extrapolation procedures and enhancement effects due to electron screening. We will report on the indirect study of the {sup 17}O(p,α){sup 14}N reaction via the THM by applying the approach developed for extracting the resonance strength of narrow resonance in the ultralow energy region. Two measurements will be described and the experimental THM cross sections will be shown for both experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2015

Nanocluster ionization energies and work function of aluminum, and their temperature dependence

Ionization threshold energies of Al{sub n} (n = 32-95) nanoclusters are determined by laser ionization of free neutral metal clusters thermalized to several temperatures in the range from 65 K to 230 K. The photoion yield curves of cold clusters follow a quadratic energy dependence above threshold, in agreement with the Fowler law of surface photoemission. Accurate data collection and analysis procedures make it possible to resolve very small (few parts in a thousand) temperature-induced shifts in the ionization energies. Extrapolation of the data to the bulk limit enables a determination of the thermal shift of the polycrystalline metal work function, found to be in excellent agreement with theoretical prediction based on the influence of thermal expansion. Small clusters display somewhat larger thermal shifts, reflecting their greater susceptibility to thermal expansion. Ionization studies of free size-resolved nanoclusters facilitate understanding of the interplay of surface, electronic, and lattice properties under contamination-free conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2015

On a simple way to calculate electronic resonances for polyatomic molecules

We propose a simple method for calculation of low-lying shape electronic resonances of polyatomic molecules. The method introduces a perturbation potential and requires only routine bound-state type calculations in the real domain of energies. Such a calculation is accessible by most of the free or commercial quantum chemistry software. The presented method is based on the analytical continuation in a coupling constant model, but unlike its previous variants, we experience a very stable and robust behavior for higher-order extrapolation functions. Moreover, the present approach is independent of the correlation treatment used in quantum many-electron computations and therefore we are able to apply Coupled Clusters (CCSD-T) level of the correlation model. We demonstrate these properties on determination of the resonance position and width of the {sup 2}Π{sub u} temporary negative ion state of diacetylene using CCSD-T level of theory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗