Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Function approximation”

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,243 records · Page 69
All sourcesNASA NTRSDOE OSTI
DOE OSTI2012

Spring constant calibration of atomic force microscope cantilevers of arbitrary shape

The spring constant of an atomic force microscope cantilever is often needed for quantitative measurements. The calibration method of Sader et al. [Rev. Sci. Instrum. 70, 3967 (1999)] for a rectangular cantilever requires measurement of the resonant frequency and quality factor in fluid (typically air), and knowledge of its plan view dimensions. This intrinsically uses the hydrodynamic function for a cantilever of rectangular plan view geometry. Here, we present hydrodynamic functions for a series of irregular and non-rectangular atomic force microscope cantilevers that are commonly used in practice. Cantilever geometries of arrow shape, small aspect ratio rectangular, quasi-rectangular, irregular rectangular, non-ideal trapezoidal cross sections, and V-shape are all studied. This enables the spring constants of all these cantilevers to be accurately and routinely determined through measurement of their resonant frequency and quality factor in fluid (such as air). An approximate formulation of the hydrodynamic function for microcantilevers of arbitrary geometry is also proposed. Implementation of the method and its performance in the presence of uncertainties and non-idealities is discussed, together with conversion factors for the static and dynamic spring constants of these cantilevers. These results are expected to be of particular value to the design and application of micro- and nanomechanical systems in general.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2015

First-order derivative couplings between excited states from adiabatic TDDFT response theory

We present a complete derivation of derivative couplings between excited states in the framework of adiabatic time-dependent density functional response theory. Explicit working equations are given and the resulting derivative couplings are compared with derivative couplings from a pseudo-wavefunction ansatz. For degenerate excited states, i.e., close to a conical intersection (CI), the two approaches are identical apart from an antisymmetric overlap term. However, if the difference between two excitation energies equals another excitation energy, the couplings from response theory exhibit an unphysical divergence. This spurious behavior is a result of the adiabatic or static kernel approximation of time-dependent density functional theory leading to an incorrect analytical structure of the quadratic response function. Numerical examples for couplings close to a CI and for well-separated electronic states are given.

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

First-Principles Investigation of Magnetic Properties and Faraday Rotation of Co Doping CdTe

In this paper, we aimed to study the effect of doping of the compound CdTe with the cobalt impurity (Co), as well as the vacancy defects in Cd sites. On one hand, this leads to the investigation of the magnetic properties and the Faraday rotation effect for the studied alloy, doped with different concentrations of cobalt (0.01, 0.05, 0.1, 0.15, 0.2, 0.25). On the other hand, we have created 0.01 of vacancy defects in Cd sites. Then, we raised the vacancy defect concentration to 0.05, keeping the same concentrations of cobalt. As a result, we have investigated that there is a magnetism appearing with the cobalt doping, while the vacancy defects in Cd sites affect the stability of the magnetic states. In general, it improves the ferromagnetic state that will be well explained in the discussion. Our calculations were performed using the KKR-CPA method within the spin-polarized density functional theory (DFT) with the local density approximation (LDA). The curves of the density of states (DOS) illustrate the results of this study which has been discussed, analyzed, and explained below. In addition, the energy of each case was calculated and given in the tables below for both the DLM (half of the cobalt spins are up while the other half are down) and the ferromagnetic state in order to confirm which one of them is stable.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Structural and electronic properties of TiX (X=N, As) in rock salt and zinc blende phase: A DFT study

Quantum mechanical first principle calculations have been performed to study the electronic and structural properties of TiN and TiAs in zinc blende (ZB) and rock salt (RS) structures. The full-potential linearized augmented plane wave (FP-LAPW) method has been used within the framework of density functional theory (DFT). The exchange correlation functional has been solved employing generalized gradient approximation (GGA). Our predicted results for lattice constants are in good agreement with the earlier findings. The electronic band structures of TiX are metallic in both the phases.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Ab initio theory of the equation of state for compressed rare gas crystals

Nonempirical equations of state of compressed rare gas crystals Ne, Ar, Kr, and Xe are studied on the basis of the earlier-obtained ab initio adiabatic potential. The paired and three-body short-range repulsive potentials are calculated by the Hartree–Fock method in the basis of localized functions with their exact mutual orthogonalization and do not contain experimentally determined parameters. The theory is compared with the experiment and results of calculations by other authors. Analysis of the proposed equations of state for large compressions has shown the importance of taking into account the three-body interaction and the terms of the higher order in the overlap integral in compressed neon and the sufficiency of the quadratic approximation in the orthogonalization of functions in heavy rare gas crystals.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2011

Extremal trajectories in a nilpotent sub-Riemannian problem on the Engel group

On the Engel group a nilpotent sub-Riemannian problem is considered, a 4-dimensional optimal control problem with a 2-dimensional linear control and an integral cost functional. It arises as a nilpotent approximation to nonholonomic systems with 2-dimensional control in a 4-dimensional space (for example, a system describing the navigation of a mobile robot with trailer). A parametrization of extremal trajectories by Jacobi functions is obtained. A discrete symmetry group and its fixed points, which are Maxwell points, are described. An estimate for the cut time (the time of the loss of optimality) on extremal trajectories is derived on this basis. Bibliography: 25 titles.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Electronic and optical properties of NbO{sub 2}

In the present study, we combine theoretical and experimental approaches in order to gain insight into the electronic properties of both the high-temperature, rutile (metallic) and low-temperature, body-centered tetragonal (insulating) phase of niobium dioxide (NbO{sub 2}) as well as the optical properties of the low-temperature phase. Theoretical calculations performed at the level of the local density approximation, Hubbard U correction, and hybrid functional are complemented with the spectroscopic ellipsometry (SE) of epitaxial films grown by molecular beam epitaxy. For the rutile phase, the local density approximation (LDA) gives the best description and predicts Fermi surface nesting consistent with wave vectors that lead to niobium-niobium dimerization during the phase transition. For the insulating phase, LDA provides a good quantitative description of the lattice, but only a qualitative description for the band gap. Including a Hubbard U correction opens the band gap at the expense of correctly describing the valence band and lattice of both phases. The hybrid functional slightly overestimates the band gap. Ellipsometric measurement is consistent with insulating behavior with a 1.0 eV band gap. Comparison with the theoretical dielectric functions, obtained utilizing a scissors operator to adjust the LDA band gap to reproduce the ellipsometry data, allows for identification of the SE peak features.

36 MATERIALS SCIENCE↗
DOE OSTI2025

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↗
DOE OSTI2024

Local convergence analysis of an inexact trust-region method for nonsmooth optimization

In Baraldi, we introduced an inexact trust-region algorithm for minimizing the sum of a smooth nonconvex function and a nonsmooth convex function in Hilbert space—a class of problems that is ubiquitous in data science, learning, optimal control, and inverse problems. Furthermore, this algorithm has demonstrated excellent performance and scalability with problem size. In this paper, we enrich the convergence analysis for this algorithm, proving strong convergence of the iterates with guaranteed rates. In particular, we demonstrate that the trust-region algorithm recovers superlinear, even quadratic, convergence rates when using a second-order Taylor approximation of the smooth objective function term.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Enabling attractive-repulsive potentials in binary-collision-approximation monte-carlo codes for ion-surface interactions

Abstract Binary Collision Approximation (BCA) codes for ion-material interactions, such as SRIM, Tridyn, F-TRIDYN, and SDtrimSP, have historically been limited to screened Coulomb potentials even at low energies due to the difficulty in numerically solving the Distance of Closest Approach (DOCA) problem for attractive-repulsive potentials. Techniques such as direct n-body simulation or modifications to Newton’s method are either prohibitively costly or not guaranteed to work for all potentials. Advanced rootfinding techniques, such as companion matrix solvers, offer a solution. For many attractive-repulsive potentials, however, a companion matrix cannot be used directly, because there is no way to put the associated functions into a monomial basis form. A complementary technique is proxy rootfinding—by finding the best-fit polynomial approximant of a function, the zeros of the approximant can be guaranteed to be close to the zeros of the function. Using the Chebyshev basis and grid offers additional guarantees with regards to the quality of the approximation, the speed of convergence, and the avoidance of Runge’s phenomenon. By finding Chebyshev interpolants and using the Chebyshev-Frobenius companion matrix, the zeros of any real function on a bounded domain can be found. Here we show that using an Adaptive Chebyshev Proxy Rootfinder with Automatic Subdivision (ACPRAS) with appropriate scaling functions, numerical issues presented by attractive-repulsive potentials, including those of scale, can be handled. Using these techniques, we show that it is possible to include any physically reasonable interatomic potential in a BCA code, and to guarantee correctness of the resulting scattering angle calculations.

Materials Science↗
DOE OSTI2025

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↗
DOE OSTI2016

Low rank approximation in G 0 W 0 calculations

The single particle energies obtained in a Kohn-Sham density functional theory (DFT) calculation are generally known to be poor approximations to electron excitation energies that are measured in tr ansport, tunneling and spectroscopic experiments such as photo-emission spectroscopy. The correction to these energies can be obtained from the poles of a single particle Green’s function derived from a many-body perturbation theory. From a computational perspective, the accuracy and efficiency of such an approach depends on how a self energy term that properly accounts for dynamic screening of electrons is approximated. The G 0 W 0 approximation is a widely used technique in which the self energy is expressed as the convolution of a noninteracting Green’s function (G 0 ) and a screened Coulomb interaction (W 0 ) in the frequency domain. The computational cost associated with such a convolution is high due to the high complexity of evaluating W 0 at multiple frequencies. In this paper, we discuss how the cost of G 0 W 0 calculation can be reduced by constructing a low rank approximation to the frequency dependent part of W 0 . In particular, we examine the effect of such a low rank approximation on the accuracy of the G 0 W 0 approximation. We also discuss how the numerical convolution of G 0 and W 0 can be evaluated efficiently and accurately by using a contour deformation technique with an appropriate choice of the contour.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

NbTi_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 niobium-titanium (Nb-Ti). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP), https://github.com/ORNL/ScalableWorkflow_VASP_Calculations. 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 Nb 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↗
DOE OSTI2026

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↗
DOE OSTI2014

Density functional theory for the description of spherical non-associating monomers in confined media using the SAFT-VR equation of state and weighted density approximations

As a first step of an ongoing study of thermodynamic properties and adsorption of complex fluids in confined media, we present a new theoretical description for spherical monomers using the Statistical Associating Fluid Theory for potential of Variable Range (SAFT-VR) and a Non-Local Density Functional Theory (NLDFT) with Weighted Density Approximations (WDA). The well-known Modified Fundamental Measure Theory is used to describe the inhomogeneous hard-sphere contribution as a reference for the monomer and two WDA approaches are developed for the dispersive terms from the high-temperature Barker and Henderson perturbation expansion. The first approach extends the dispersive contributions using the scalar and vector weighted densities introduced in the Fundamental Measure Theory (FMT) and the second one uses a coarse-grained (CG) approach with a unique weighted density. To test the accuracy of this new NLDFT/SAFT-VR coupling, the two versions of the theoretical model are compared with Grand Canonical Monte Carlo (GCMC) molecular simulations using the same molecular model. Only the version with the “CG” approach for the dispersive terms provides results in excellent agreement with GCMC calculations in a wide range of conditions while the “FMT” extension version gives a good representation solely at low pressures. Hence, the “CG” version of the theoretical model is used to reproduce methane adsorption isotherms in a Carbon Molecular Sieve and compared with experimental data after a characterization of the material. The whole results show an excellent agreement between modeling and experiments. Thus, through a complete and consistent comparison both with molecular simulations and with experimental data, the NLDFT/SAFT-VR theory has been validated for the description of monomers.

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

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↗
DOE OSTI2015

Palladium dimers adsorbed on graphene: A DFT study

The 2D structure of graphene shows a great promise for enhanced catalytic activity when adsorbed with palladium. We performed a systematic density functional theory (DFT) study of the adsorption of palladium dimer (Pd{sub 2}) on graphene using SIESTA package, in the generalized gradient approximation (GGA). The adsorption energy, geometry, and charge transfer of Pd{sub 2}-graphene system are calculated. Both horizontal and vertical orientations of Pd{sub 2} on graphene are studied. Our calculations revealed that the minimum energy configuration for Pd dimer is parallel to the graphene sheet with its two atoms occupying centre of adjacent hexagonal rings of graphene sheet. Magnetic moment is induced for Pd dimer adsorbed on graphene in vertical orientation while horizontal orientation of Pd dimer on graphene do not exhibit magnetism. Insignificant energy differences among adsorption sites means that dimer mobility on the graphene sheet is high. There is imperceptible distortion of graphene sheet perpendicular to its plane. However, some lateral displacements are seen.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

A Variational Nodal Approach to 2D/1D Pin Resolved Neutron Transport for Pressurized Water Reactors

A two-dimensional/one-dimensional (2D/1D) variational nodal approach is presented for pressurized water reactor core calculations without fuel-moderator homogenization. A 2D/1D approximation to the within-group neutron transport equation is derived and converted to an even-parity form. The corresponding nodal functional is presented and discretized to obtain response matrix equations. Within the nodes, finite elements in the x-y plane and orthogonal functions in z are used to approximate the spatial flux distribution. On the radial interfaces, orthogonal polynomials are employed; on the axial interfaces, piecewise constants corresponding to the finite elements eliminate the interface homogenization that has been a challenge for method of characteristics (MOC)-based 2D/1D approximations. The angular discretization utilizes an even-parity integral method within the nodes, and low-order spherical harmonics (P N ) on the axial interfaces. The x-y surfaces are treated with high-order P N combined with quasi-reflected interface conditions. Furthermore, the method is applied to the C5G7 benchmark problems and compared to Monte Carlo reference calculations.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗