Search NASA⌕ Search

SEARCH · Search NASA

Results for “gradient calculation”

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,189 records · Page 66

Toward Ultraflat Surface Morphologies During Focused Electron Beam Induced Nanosynthesis: Disruption Origins and Compensation

Here we explore emerging applications for nanoscale materials demand precise deposit shape retention from design to deposition. This study investigates the effects that disrupt high-fidelity shapes during focused electron beam induced nanosynthesis. It is shown that process parameters, patterning strategies and deposit topography can impose lateral precursor coverage gradients during growth resulting in unwanted topographic artifacts. The study classifies the evolving surface shapes into four general types and explains the formation and transition from a fundamental point of view. Continuum model calculations and simulations expand the experimental results to provide a comprehensive insight into understand the disruption mechanism. The findings demonstrate that the well-established concept of growth regimes has to be expanded by its lateral gradients as they strongly influence final shape fidelities. Finally, the study is complemented by a compensation strategy that improves the edge fidelity on the lower nanoscale to further push this technique toward the intrinsic limitations.

36 MATERIALS SCIENCE↗

HfZr_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys hafnium-zirconium (Hf-Zr). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Hf and Zr. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions.

36 MATERIALS SCIENCE↗

HfTi_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys hafnium-titanium (Hf-Ti). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Hf and Ti. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions. The calculations have been collected on Air Force HPC11 cluster using the VASP 6.5.1. Additional methodology and file structure information is available in the dataset README.txt file.

36 MATERIALS SCIENCE↗

TiV_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys titanium-vanadium (Ti-V). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Ti and V. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions. The calculations have been collected on Air Force HPC11 cluster using the VASP 6.5.1. Additional methodology and file structure information is available in the dataset README.txt file.

36 MATERIALS SCIENCE↗

HfTa_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys hafnium-tantalum (Hf-Ta). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Hf and Ta. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect BCC lattice sites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,040 randomized atomic structures over 31 chemical compositions. Further methodological and structural information is contained in the dataset README.txt file.

36 MATERIALS SCIENCE↗

The role of cell size in density gradient electrophoretic separation of mouse leukemia cells according to position in the cell cycle

Cultured mouse leukemia cells line L5178Y were subjected to upward electrophoresis in a density gradient and the slower migrating cell populations were enriched in G2 cells. It is indicated that this cell line does not change electrophoretic mobility through the cell cycle. The possibility that increased sedimentation downward on the part of the larger G2 cells caused this separation was explored. Two different cell populations were investigated. The log phase population was found to migrate upward faster than the G2 population, and a similar difference between their velocities and calculated on the basis of a 1 um diameter difference between the two cell populations. The G2 and G1 enriched populations were isolated by Ficoll density gradient sedimentation. The bottom fraction was enriched in G2 cells and the top fraction was enriched with G1 cells, especially when compared with starting materials. The electrophoretic mobilities of these two cell populations did not differ significantly from one another. Cell diameter dependent migration curves were calculated and were found to be different. Families of migration curves that differ when cell size is considered as a parameter are predicted.

Plank, L. D.↗

Near total reflection x-ray photoelectron spectroscopy: quantifying chemistry at solid/liquid and solid/solid interfaces

Near total reflection regime has been widely used in x-ray science, specifically in grazing incidence small angle x-ray scattering and in hard x-ray photoelectron spectroscopy (XPS). In this work, we introduce some practical aspects of using near total reflection (NTR) in ambient pressure XPS and apply this technique to study chemical concentration gradients in a substrate/photoresist system. Experimental data are accompanied by x-ray optical and photoemission simulations to quantitatively probe the photoresist and the interface with the depth accuracy of ~1 nm. Together, our calculations and experiments confirm that NTR XPS is a suitable method to extract information from buried interfaces with highest depth-resolution, which can help address open research questions regarding our understanding of concentration profiles, electrical gradients, and charge transfer phenomena at such interfaces. Furthermore, the presented methodology is especially attractive for solid/liquid interface studies, since it provides all the strengths of a Bragg-reflection standing-wave spectroscopy without the need of an artificial multilayer mirror serving as a standing wave generator, thus dramatically simplifying the sample synthesis.

solid/solid interface↗

Structural, elastic, electronic and phonon properties of SnX{sub 2}O{sub 4} (X=Mg, Zn, Cd) spinel from density functional theory

First-principle calculations of structural, electronic, elastic and phonon properties of SnMg{sub 2}O{sub 4}, SnZn{sub 2}O{sub 4} and SnCd{sub 2}O{sub 4} compounds are presented, using the pseudo-potential plane waves approach based on density functional theory (DFT) within the generalized gradient approximation (GGA). The computed ground state structural parameters, i.e. lattice constants, internal free parameter and bulk modulus are in good agreement with the available theoretical results. Our calculated elastic constants are indicative of stability of SnX{sub 2}O{sub 4} (X=Mg, Zn, Cd) compounds in the spinel structure. The partial density of states (PDOS) of these compounds is in good agreement with the earlier ab-initio calculations. The phonon dispersion relations were calculated using the direct method. Phonon dispersion results indicate that SnZn{sub 2}O{sub 4} is dynamically stable, while SnMg{sub 2}O{sub 4} and SnCd{sub 2}O{sub 4} are unstable.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

SU-F-T-600: Influence of Acuros XB and AAA Dose Calculation Algorithms On Plan Quality Metrics and Normal Lung Doses in Lung SBRT

Purpose: To study the influence of superposition-beam model (AAA) and determinant-photon transport-solver (Acuros XB) dose calculation algorithms on the treatment plan quality metrics and on normal lung dose in Lung SBRT. Methods: Treatment plans of 10 Lung SBRT patients were randomly selected. Patients were prescribed to a total dose of 50-54Gy in 3–5 fractions (10?5 or 18?3). Doses were optimized accomplished with 6-MV using 2-arcs (VMAT). Doses were calculated using AAA algorithm with heterogeneity correction. For each plan, plan quality metrics in the categories- coverage, homogeneity, conformity and gradient were quantified. Repeat dosimetry for these AAA treatment plans was performed using AXB algorithm with heterogeneity correction for same beam and MU parameters. Plan quality metrics were again evaluated and compared with AAA plan metrics. For normal lung dose, V{sub 20} and V{sub 5} to (Total lung- GTV) were evaluated. Results: The results are summarized in Supplemental Table 1. PTV volume was mean 11.4 (±3.3) cm{sup 3}. Comparing RTOG 0813 protocol criteria for conformality, AXB plans yielded on average, similar PITV ratio (individual PITV ratio differences varied from −9 to +15%), reduced target coverage (−1.6%) and increased R50% (+2.6%). Comparing normal lung doses, the lung V{sub 20} (+3.1%) and V{sub 5} (+1.5%) were slightly higher for AXB plans compared to AAA plans. High-dose spillage ((V105%PD - PTV)/ PTV) was slightly lower for AXB plans but the % low dose spillage (D2cm) was similar between the two calculation algorithms. Conclusion: AAA algorithm overestimates lung target dose. Routinely adapting to AXB for dose calculations in Lung SBRT planning may improve dose calculation accuracy, as AXB based calculations have been shown to be closer to Monte Carlo based dose predictions in accuracy and with relatively faster computational time. For clinical practice, revisiting dose-fractionation in Lung SBRT to correct for dose overestimates attributable to algorithm may very well be warranted.

60 APPLIED LIFE SCIENCES↗

The lagRST Model: A Turbulence Model for Non-Equilibrium Flows

This study presents a new class of turbulence model designed for wall bounded, high Reynolds number flows with separation. The model addresses deficiencies seen in the modeling of nonequilibrium turbulent flows. These flows generally have variable adverse pressure gradients which cause the turbulent quantities to react at a finite rate to changes in the mean flow quantities. This "lag" in the response of the turbulent quantities can t be modeled by most standard turbulence models, which are designed to model equilibrium turbulent boundary layers. The model presented uses a standard 2-equation model as the baseline for turbulent equilibrium calculations, but adds transport equations to account directly for non-equilibrium effects in the Reynolds Stress Tensor (RST) that are seen in large pressure gradients involving shock waves and separation. Comparisons are made to several standard turbulence modeling validation cases, including an incompressible boundary layer (both neutral and adverse pressure gradients), an incompressible mixing layer and a transonic bump flow. In addition, a hypersonic Shock Wave Turbulent Boundary Layer Interaction with separation is assessed along with a transonic capsule flow. Results show a substantial improvement over the baseline models for transonic separated flows. The results are mixed for the SWTBLI flows assessed. Separation predictions are not as good as the baseline models, but the over prediction of the peak heat flux downstream of the reattachment shock that plagues many models is reduced.

Lillard, Randolph P.↗

Quasilinear Carbon Transport In An Impurity Hole Plasma In LHD

Comprehensive electrostatic gyrokinetic linear stability calculations for ion-scale microinstabilities in an LHD plasma with an ion-ITB and carbon "impurity hole" are used to make quasilinear estimates of particle flux to explore whether microturbulence can explain the observed outward carbon fluxes that flow "up" the impurity density gradient. The ion temperature is not stationary in the ion-ITB phase of the simulated discharge, during which the core carbon density decreases continuously. To fully sample these varying conditions the calculations are carried out at three radial locations and four times. The plasma parameter inputs are based on experimentally measured profiles of electron and ion temperature, as well as electron and carbon density. The spectroscopic line-average ratio of hydrogen and helium densities is used to set the density of these species. Three ion species (H,He,C) and the electrons are treated kinetically, including collisions. Electron instability drive does enhance the growth rate significantly, but the most unstable modes have characteristics of ion temperature gradient (ITG) modes in all cases. As the carbon density gradient is scanned between the measured value and zero, the quasilinear carbon flux is invariably inward when the carbon density profile is hollow, so turbulent transport due to the instabilities considered here does not explain the observed outward flux of impurities in impurity hole plasmas. The stiffness of the quasilinear ion heat flux is found to be 1.7-2.3, which is lower than several estimates in tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Adjoint-enabled multidimensional optimization of satellite electron/proton shields - 25452

The design of satellites usually includes the objective of minimizing mass due to high launch costs, which is complicated by the need to protect sensitive electronics from the space radiation environment. There is growing interest in automated design optimization techniques to help achieve that objective. Traditional optimization approaches that rely exclusively on response functions (e.g. dose calculations) can be quite expensive when applied to transport problems. Previously we showed how adjoint-based transport sensitivities used in conjunction with gradient-based optimization algorithms can be quite effective in designing mass-efficient electron/proton shields in one-dimensional slab geometries. In this paper we extend that work to two-dimensional Cartesian geometries. This consists primarily of deriving the sensitivities to geometric changes, given a particular prescription for parameterizing the shield geometry. We incorporate these sensitivities into our optimization process and demonstrate their effectiveness in such design calculations. (authors)

42 ENGINEERING↗

Achieving Hard X-ray Nanofocusing Using a Wedged Multilayer Laue Lens

Here, we report on the fabrication and the characterization of a wedged multilayer Laue lens for x-ray nanofocusing. The lens was fabricated using a sputtering deposition technique, in which a specially designed mask was employed to introduce a thickness gradient in the lateral direction of the multilayer. X-ray characterization shows an efficiency of 27% and a focus size of 26 nm at 14.6 keV, in a good agreement with theoretical calculations. Our results indicate that the desired wedging is achieved in the fabricated structure. Furthermore, we anticipate that continuous development on wedged MLLs will advance x-ray nanofocusing optics to new frontiers and enrich capabilities and opportunities for hard X-ray microscopy.

42 ENGINEERING↗

Finite difference solution for transient radiative cooling of a conducting semitransparent square region

Transient solutions were obtained for a square region of heat conducting semitransparent material cooling by thermal radiation. The region is in a vacuum environment, so energy is dissipated only by radiation from within the medium leaving through its boundaries. The effect of heat conduction during the transient is to partially equalize the internal temperature distribution. As the optical thickness of the region is increased, the temperature gradients increase near the boundaries and corners, unless heat conduction is large. The solution procedure must provide accurate temperature distributions in these regions to prevent error in the calculated radiation losses. Two-dimensional numerical Gaussian integration is used to obtain the local radiative source term. A finite difference procedure with variable space and time increments is used to solve the transient energy equation. Variable spacing was used to concentrate grid points in regions with large temperature gradients.

Siegel, R.↗

SU-E-J-249: Characterization of Gynecological Tumor Heterogeneity Using Texture Analysis in the Context of An 18F-FDG PET Adaptive Protocol

Purpose: We propose a method to examine gynecological tumor heterogeneity using texture analysis in the context of an adaptive PET protocol in order to establish if texture metrics from baseline PET-CT predict tumor response better than SUV metrics alone as well as determine texture features correlating with tumor response during radiation therapy. Methods: This IRB approved protocol included 29 women with node positive gynecological cancers visible on FDG-PET treated with EBRT to the PET positive nodes. A baseline and intra-treatment PET-CT was obtained. Tumor outcome was determined based on RECIST on posttreatment PET-CT. Primary GTVs were segmented using 40% threshold and a semi-automatic gradient-based contouring tool, PET Edge (MIM Software Inc., Cleveland, OH). SUV histogram features, Metabolic Volume (MV), and Total Lesion Glycolysis (TLG) were calculated. Four 3D texture matrices describing local and regional relationships between voxel intensities in the GTV were generated: co-occurrence, run length, size zone, and neighborhood difference. From these, 39 texture features were calculated. Prognostic power of baseline features derived from gradientbased and threshold GTVs were determined using the Wilcoxon rank-sum test. Receiver Operating Characteristics and logistic regression was performed using JMP (SAS Institute Inc., Cary, NC) to find probabilities of predicting response. Changes in features during treatment were determined using the Wilcoxon signed-rank test. Results: Of the 29 patients, there were 16 complete responders, 7 partial responders, and 6 non-responders. Comparing CR/PR vs. NR for gradient-based GTVs, 7 texture values, TLG, and SUV kurtosis had a p < 0.05. Threshold GTVs yielded 4 texture features and TLG with p < 0.05. From baseline to intra-treatment, 14 texture features, SUVmean, SUVmax, MV, and TLG changed with p < 0.05. Conclusion: Texture analysis of PET imaged gynecological tumors is an effective method for early prognosis and should be used complimentary to SUV metrics, especially when using gradient based segmentation.

60 APPLIED LIFE SCIENCES↗

Dissociative electron attachment to 5-bromo-uracil: non-adiabatic dynamics on complex-valued potential energy surfaces

Electron induced dissociation reactions are relevant to many fields, ranging from prebiotic chemistry to cancer treatments. However, the simulation of dissociation electron attachment (DEA) dynamics is very challenging because the auto-ionization widths of the transient negative ions must be accounted for. We propose an adaptation of the ab initio multiple spawning (AIMS) method for complex-valued potential energy surfaces, along the lines of recent developments based on surface hopping dynamics. Our approach combines models for the energy dependence of the auto-ionization widths, obtained from scattering calculations, with survival probabilities computed for the trajectory basis functions employed in the AIMS dynamics. The method is applied to simulate the DEA dynamics of 5-bromouracil in full dimensionality, i.e., taking all the vibrational modes into consideration. The propagation starts on the $π_{2}^{*}$ resonance state and describes the formation of Br - anions mediated by non-adiabatic couplings. The potential energies, gradients and non-adiabatic couplings were computed with the fractional-occupancy molecular orbital complete-active-space configuration-interaction method, and the calculated DEA cross section are consistent with the observed DEA intensities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Unraveling the intrinsic atomic physics behind x-ray absorption line shifts in warm dense silicon plasmas

Here, we present a novel free-energy density functional theory (DFT)-based methodology for optical property calculations of warm dense matter to cover a wide range of thermodynamic conditions and photon energies including the entire x-ray range. It uses Mermin–Kohn–Sham density functional theory with exchange-correlation (XC) thermal effects taken into account via a fully temperature dependent generalized gradient approximation XC functional. The methodology incorporates a combination of the ab initio molecular dynamics (AIMD) snapshotted Kubo–Greenwood optic data with a single-atom in simulation cell calculations to close the photon energy gap between the L and K edges and extend the K-edge tail toward many-keV photon energies. This gap arises in the standard scheme due to a prohibitively large number of bands required for the Kubo–Greenwood calculations with AIMD snapshots. Kubo–Greenwood data on snapshots provide an accurate description of optic properties at low photon frequencies slightly beyond the L-edge and x-ray absorption near edges structure (XANES) spectra, while data from periodic calculations with single atoms cover the tail regions beyond the edges. To demonstrate its applicability to mid-Z materials where the standard DFT-based approach is not computationally feasible, we have applied it to opacity calculations of warm dense silicon plasmas. These first-principles calculations revealed a very interesting phenomon of red-to-blue shift in K–L (1s → 2p) and K-edge absorptions along both isotherm and isochore, which are absent in most continuum-lowering models of traditional plasma physics. This new physics phenomenon can be attributed to the underlying competition between the screening of deeply bound core electrons and the screening of outer-shell electrons caused by warm-dense-plasma conditions. We further demonstrate that the ratio of 1s → 2p to the K-edge x-ray absorptions can be used to characterize warm-dense-plasma conditions. Eventually, based on our absorption calculations, we have established a first-principles opacity table (FPOT) for silicon in a wide range of material densities and temperatures.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Methodology to Generate Decay Gamma Sources for the Second Target Station Full Target Assembly

This report details the specifics of applying the Position-Averaged Method to the target wedges and staves. The methodology used to calculate the decay gamma sources for a single target wedge to the calculation of the Second Target Station (STS) target assembly activation and the resulting decay gamma source terms are detailed in this report. The methodology has been applied to the entire STS target assembly including all of the target wedges, staves, shaft, and drive. The target staves and shaft have been vertically segmented to ensure the decay gamma source gradient is sufficiently captured. MCNP ® Code Version 6.2.0 with the RNUCS patch coupled with CINDER2008 from the AARE V1.0 package are the main tools used to calculate the decay gamma source terms.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗