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,441 records · Page 80

Investigation on phonon scattering in a GaAs nanowire field effect transistor using the non-equilibrium Green's function formalism

Using quantum transport simulations, the impact of electron-phonon scattering on the transfer characteristic of a gate-all-around nanowire (GaAs) field effect transistor (NWFET) has been thoroughly investigated. The Non-Equilibrium Green's Function formalism in the effective mass approximation using a decoupled mode decomposition has been deployed. NWFETs of different dimensions have been considered, and scattering mechanisms including acoustic, optical and polar optical phonons have been included. The effective masses were extracted from tight binding simulations. High and low drain bias have been considered. We found substantial source to drain tunnelling current and significant impact of phonon scattering on the performance of the NWFET. At low drain bias, for a 2.2 × 2.2 nm{sup 2} cross-section transistor, scattering caused a 72%, 77%, and 81% decrease in the on-current for a 6 nm, 10 nm, and 20 nm channel length, respectively. This reduction in the current due to scattering is influenced by the increase in the tunnelling current. We include the percentage tunnelling for each valley at low and high drain bias. It was also found that the strong quantisation caused the relative position of the valleys to vary with the cross-section. This had a large effect on the overall tunnelling current. The phonon-limited mobility was also calculated, finding a mobility of 950 cm{sup 2}/V s at an inversion charge density of 10{sup 12 }cm{sup −2} for a 4.2 × 4.2 nm{sup 2} cross-section device.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Non-intrusive nonlinear model reduction via machine learning approximations to low-dimensional operators

Abstract Although projection-based reduced-order models (ROMs) for parameterized nonlinear dynamical systems have demonstrated exciting results across a range of applications, their broad adoption has been limited by their intrusivity: implementing such a reduced-order model typically requires significant modifications to the underlying simulation code. To address this, we propose a method that enables traditionally intrusive reduced-order models to be accurately approximated in a non-intrusive manner. Specifically, the approach approximates the low-dimensional operators associated with projection-based reduced-order models (ROMs) using modern machine-learning regression techniques. The only requirement of the simulation code is the ability to export the velocity given the state and parameters; this functionality is used to train the approximated low-dimensional operators. In addition to enabling nonintrusivity, we demonstrate that the approach also leads to very low computational complexity, achieving up to $$10^3{\times }$$ 10 3 × in run time. We demonstrate the effectiveness of the proposed technique on two types of PDEs. The domain of applications include both parabolic and hyperbolic PDEs, regardless of the dimension of full-order models (FOMs).

42 ENGINEERING↗

Twisted bilayer graphene. I. Matrix elements, approximations, perturbation theory, and a k · p two-band model

We investigate the twisted bilayer graphene (TBG) model of Bistritzer and MacDonald (BM) [Bistritzer and MacDonald, Proc. Natl. Acad. Sci. 108, 12233 (2011)] to obtain an analytic understanding of its energetics and wave functions needed for many-body calculations. We provide an approximation scheme for the wave functions of the BM model, which first elucidates why the BM K M -point centered original calculation containing only four plane waves provides a good analytical value for the first magic angle (θ M ≈ 1°). The approximation scheme also elucidates why most of the many-body matrix elements in the Coulomb Hamiltonian projected to the active bands can be neglected. By applying our approximation scheme at the first magic angle to a Γ M -point centered model of six plane waves, we analytically understand the reason for the small Γ M -point gap between the active and passive bands in the isotropic limit w 0 = w 1 . Furthermore, we analytically calculate the group velocities of the passive bands in the isotropic limit, and show that they are almost doubly degenerate, even away from the Γ M point, where no symmetry forces them to be. Furthermore, moving away from the Γ M and K M points, we provide an explicit analytical perturbative understanding as to why the TBG bands are flat at the first magic angle, despite the first magic angle is defined by only requiring a vanishing K M -point Dirac velocity. We derive analytically a connected “magic manifold” w 1 = $2\sqrt{1 + w^{2}_{0}}$ $-\sqrt{2 + 3w^2_0}$, on which the bands remain extremely flat as w 0 is tuned between the isotropic (w 0 = w 1 ) and chiral (w 0 = 0) limits. We analytically show why going away from the isotropic limit by making w 0 less (but not larger) than w 1 increases the Γ M -point gap between the active and the passive bands. Finally, by perturbation theory, we provide an analytic Γ M point k ∙ p two-band model that reproduces the TBG band structure and eigenstates within a certain w 0 , w 1 parameter range. Further refinement of this model are discussed, which suggest a possible faithful representation of the TBG bands by a two-band Γ M point k ∙ p model in the full w 0 , w 1 parameter range.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Modeling of non-equilibrium phenomena in expanding flows by means of a collisional-radiative model

The effects of non-equilibrium in a quasi-one-dimensional nozzle flow are investigated by means of a collisional-radiative model. The gas undergoing the expansion is an air plasma and consists of atoms, molecules, and free electrons. In the present analysis, the electronic excited states of atomic and molecular species are treated as separate pseudo-species. Rotational and vibrational energy modes are assumed to be populated according to Boltzmann distributions. The coupling between radiation and gas dynamics is accounted for, in simplified manner, by using escape factors. The flow governing equations for the steady quasi-one-dimensional flow are written in conservative form and discretized in space by means of a finite volume method. Steady-state solutions are obtained by using a fully implicit time integration scheme. The analysis of the evolution of the electronic distribution functions reveals a substantial over-population of the high-lying excited levels of atoms and molecules in correspondence of the nozzle exit. The influence of optical thickness is also studied. The results clearly demonstrate that the radiative transitions, within the optically thin approximation, drastically reduce the over-population of high-lying electronic levels.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Rational RBF-based partition of unity method for efficiently and accurately approximating 3D objects

We consider the problem of reconstructing 3D objects via meshfree interpolation methods. In this framework, we usually deal with large data sets and thus develop an efficient local scheme via the well-known partition of unity method. The main contribution in this paper consists in constructing the local interpolants for the implicit interpolation by means of rational radial basis functions. Numerical experiments, devoted to test the accuracy of the scheme, confirm that the proposed method is particularly performing when 3D objects, or more in general implicit functions defined by scattered points, need to be approximated.

97 MATHEMATICS AND COMPUTING↗

A course-grained model for polyethylene glycol polymer

A coarse-grained (CG) model of polyethylene glycol (PEG) was developed and implemented in CG molecular dynamics (MD) simulations of PEG chains with degree of polymerization (DP) 20 and 40. In the model, two repeat units of PEG are grouped as one CG bead. Atomistic MD simulation of PEG chains with DP = 20 was first conducted to obtain the bonded structural probability distribution functions (PDFs) and nonbonded pair correlation function (PCF) of the CG beads. The bonded CG potentials are obtained by simple inversion of the corresponding PDFs. The CG nonbonded potential is parameterized to the PCF using both an inversion procedure based on the Ornstein-Zernike equation with the Percus-Yevick approximation (OZPY{sup -1}) and a combination of OZPY{sup -1} with the iterative Boltzmann inversion (IBI) method (OZPY{sup -1}+IBI). As a simple one step method, the OZPY{sup -1} method possesses an advantage in computational efficiency. Using the potential from OZPY{sup -1} as an initial guess, the IBI method shows fast convergence. The coarse-grained molecular dynamics (CGMD) simulations of PEG chains with DP = 20 using potentials from both methods satisfactorily reproduce the structural properties from atomistic MD simulation of the same systems. The OZPY{sup -1}+IBI method yields better agreement than the OZPY{sup -1} method alone. The new CG model and CG potentials from OZPY{sup -1}+IBI method was further tested through CGMD simulation of PEG with DP = 40 system. No significant changes are observed in the comparison of PCFs from CGMD simulations of PEG with DP = 20 and 40 systems indicating that the potential is independent of chain length.

36 MATERIALS SCIENCE↗

Ultra-fast computation of electronic spectra for large systems by tight-binding based simplified Tamm-Dancoff approximation (sTDA-xTB)

The computational bottleneck of the extremely fast simplified Tamm-Dancoff approximated (sTDA) time-dependent density functional theory procedure [S. Grimme, J. Chem. Phys. 138, 244104 (2013)] for the computation of electronic spectra for large systems is the determination of the ground state Kohn-Sham orbitals and eigenvalues. This limits such treatments to single structures with a few hundred atoms and hence, e.g., sampling along molecular dynamics trajectories for flexible systems or the calculation of chromophore aggregates is often not possible. The aim of this work is to solve this problem by a specifically designed semi-empirical tight binding (TB) procedure similar to the well established self-consistent-charge density functional TB scheme. The new special purpose method provides orbitals and orbital energies of hybrid density functional character for a subsequent and basically unmodified sTDA procedure. Compared to many previous semi-empirical excited state methods, an advantage of the ansatz is that a general eigenvalue problem in a non-orthogonal, extended atomic orbital basis is solved and therefore correct occupied/virtual orbital energy splittings as well as Rydberg levels are obtained. A key idea for the success of the new model is that the determination of atomic charges (describing an effective electron-electron interaction) and the one-particle spectrum is decoupled and treated by two differently parametrized Hamiltonians/basis sets. The three-diagonalization-step composite procedure can routinely compute broad range electronic spectra (0-8 eV) within minutes of computation time for systems composed of 500-1000 atoms with an accuracy typical of standard time-dependent density functional theory (0.3-0.5 eV average error). An easily extendable parametrization based on coupled-cluster and density functional computed reference data for the elements H–Zn including transition metals is described. The accuracy of the method termed sTDA-xTB is first benchmarked for vertical excitation energies of open- and closed-shell systems in comparison to other semi-empirical methods and applied to exemplary problems in electronic spectroscopy. As side products of the development, a robust and efficient valence electron TB method for the accurate determination of atomic charges as well as a more accurate calculation scheme of dipole rotatory strengths within the Tamm-Dancoff approximation is proposed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Conductance fluctuation of edge-disordered graphene nanoribbons: Crossover from diffusive transport to Anderson localization

Conductance fluctuation of edge-disordered graphene nanoribbons (ED-GNRs) is examined using the non-equilibrium Green's function technique combined with the extended Hückel approximation. The mean free path λ and the localization length ξ of the ED-GNRs are determined to classify the quantum transport regimes. In the diffusive regime where the length L{sub c} of the ED-GNRs is much longer than λ and much shorter than ξ, the conductance histogram is given by a Gaussian distribution function with universal conductance fluctuation. In the localization regime where L{sub c}≫ξ, the histogram is no longer the universal Gaussian distribution but a lognormal distribution that characterizes Anderson localization.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

KICKSTARTING REIONIZATION WITH THE FIRST BLACK HOLES: THE EFFECTS OF SECOND-ORDER PERTURBATION THEORY IN PRE-REIONIZATION VOLUMES

We explore structure formation in the dark ages (z {approx} 30-6) using two well-known methods for initializing cosmological N-body simulations. Overall, both the Zel'dovich approximation and second-order Lagrangian perturbation theory (2LPT) are known to produce accurate present-day dark matter halo mass functions. However, since the 2LPT method drives more rapid evolution of dense regions, it increases the occurrence of rare massive objects-an effect that is most pronounced at high redshift. We find that 2LPT produces more halos that could harbor Population III stars and their black hole remnants, and they produce them earlier. Although the differences between the 2LPT and Zel'dovich approximation mass functions are nearly erased by z = 6, this small boost to the number and mass of black holes more than doubles the reionized volume of the early universe. We discuss the implications for reionization and massive black hole growth.

79 ASTRONOMY AND ASTROPHYSICS↗

Two charges on plane in a magnetic field I. “Quasi-equal” charges and neutral quantum system at rest cases

Low-lying bound states for the problem of two Coulomb charges of finite masses on a plane subject to a constant magnetic field B perpendicular to the plane are considered. Major emphasis is given to two systems: two charges with the equal charge-to-mass ratio (quasi-equal charges) and neutral systems with concrete results for the hydrogen atom and two electrons (quantum dot). It is shown that for these two cases, when a neutral system is at rest (the center-of-mass momentum is zero), some outstanding properties occur: in double polar coordinates in CMS (R,ϕ) and relative (ρ,φ) coordinate systems (i) the eigenfunctions are factorizable, all factors except for ρ-dependent are found analytically, they have definite relative angular momentum, (ii) dynamics in ρ-direction is the same for both systems being described by a funnel-type potential; (iii) at some discrete values of dimensionless magnetic fields b≤1 the system becomes quasi-exactly-solvable and a finite number of eigenfunctions in ρ are polynomials. The variational method is employed. Trial functions are based on combining for the phase of a wavefunction (a) the WKB expansion at large distances, (b) the perturbation theory at small distances (c) with a form of the known analytically (quasi-exactly-solvable) eigenfunctions. Such a form of trial function appears as a compact uniform approximation for lowest eigenfunctions. For the lowest states with relative magnetic quantum numbers s=0,1,2 this approximation gives not less than 7 s.d., 8 s.d., 9 s.d., respectively, for the total energy E(B) for magnetic fields 0.049a.u.<B<2000a.u. (hydrogen atom) and 0.025a.u.⋜B⋜1000a.u. (two electrons). The evolution of nodes of excited states with the magnetic field change is indicated. In the framework of convergent perturbation theory the corrections to proposed approximations are evaluated. -- Highlights: •Approximate solution for the low-lying eigenstates for two Coulomb charges given. •Factorization of eigenfunctions in double polar center-of-mass coordinates uncovered. •A magnetic field range varies from weak to ultra-strong, 1000–2000 a.u.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Topological phase transition in LaAs under pressure

Using density functional theory (DFT) calculations with the screened hybrid functional of Heyd, Scuseria, and Ernzerhof (HSE06), we study the effects of hydrostatic pressure on the electronic properties of LaAs. We focus on the band crossing near the X point that can make LaAs a topological semimetal, discussing results of both DFT within the generalized gradient approximation (GGA) and the HSE06 hybrid functional. We find that in DFT-GGA, under the calculated equilibrium lattice parameter, LaAs displays a crossing between the highest As p band and the lowest La d band near the X point due to the overestimated p–d band overlap. Such crossing does not occur when the band overlap is corrected in the HSE06 calculation. However, we find that the p–d crossing can be induced in LaAs under hydrostatic pressure, showing a topological phase transition at ~7GPa. Furthermore, the rocksalt crystal structure of LaAs is predicted to be stable under applied pressure up to 20 GPa, in good agreement with experimental observations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Benchmark tests and spin adaptation for the particle-particle random phase approximation

The particle-particle random phase approximation (pp-RPA) provides an approximation to the correlation energy in density functional theory via the adiabatic connection [H. van Aggelen, Y. Yang, and W. Yang, Phys. Rev. A 88, 030501 (2013)]. It has virtually no delocalization error nor static correlation error for single-bond systems. However, with its formal O(N{sup 6}) scaling, the pp-RPA is computationally expensive. In this paper, we implement a spin-separated and spin-adapted pp-RPA algorithm, which reduces the computational cost by a substantial factor. We then perform benchmark tests on the G2/97 enthalpies of formation database, DBH24 reaction barrier database, and four test sets for non-bonded interactions (HB6/04, CT7/04, DI6/04, and WI9/04). For the G2/97 database, the pp-RPA gives a significantly smaller mean absolute error (8.3 kcal/mol) than the direct particle-hole RPA (ph-RPA) (22.7 kcal/mol). Furthermore, the error in the pp-RPA is nearly constant with the number of atoms in a molecule, while the error in the ph-RPA increases. For chemical reactions involving typical organic closed-shell molecules, pp- and ph-RPA both give accurate reaction energies. Similarly, both RPAs perform well for reaction barriers and nonbonded interactions. These results suggest that the pp-RPA gives reliable energies in chemical applications. The adiabatic connection formalism based on pairing matrix fluctuation is therefore expected to lead to widely applicable and accurate density functionals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Nuclear theory and event generators for charge-changing neutrino reactions

Semi-inclusive $CC\nu$ cross sections based on factorized cross sections are studied for a selection of spectral function models with the objective of facilitating the choice of models for use as input into event generators. The basic formalism for such cross sections is presented along with an introduction to constructing spectral functions for simple models based on the independent particle shell model (IPSM), the relativistic Fermi gas model (RFG) and a local density approximation (LDA) based on the RFG. Spectral functions for these models are shown for $^{16}$O along with a more sophisticated model which includes nucleon-nucleon interactions.. Inclusive and semi-inclusive cross sections are calculated for these models. Although the inclusive cross sections are all of similar size and shape, the semi-inclusive cross sections are subtantially different depending upon whether the spectral functions contain some features associated with the nuclear shell model or are based on the RFG and LDA models. Here, calculations of average values and standard deviations of the initial neutrino energy using the semi-inclusive cross section for the various models are presented and indicate that there may be simple kinematical descriptions of the average neutrino energy which is common to all of these models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Bayesian Monte Carlo Evaluation Framework for Imperfect Data and Models [Abstract]

Nuclear data evaluation methods conventionally make the following assumptions: prior and posterior probability distribution functions (PDFs) of all model parameters and data are normal (Gaussian); the linear approximation is sufficiently accurate for minimization of a cost function (even for non-linear models); and that both the model (of, e.g., neutron cross section) and experimental data (including their covariance data) are perfect. These assumptions are inherent to the well-known generalized linear least squares (GLLS) minimization method commonly used for evaluations of resolved resonance region (RRR) neutron cross sections. However, these assumptions are often not justified due to the presence of non-normal PDFs, non-linear models (e.g. R -matrix formalism), and inherent imperfections in data and models (e.g. discrepant data sets, discrepancies between the previous evaluation and newly measured data, or imperfect covariance data). We remove the said assumptions in a mathematical framework of Bayes’ theorem, and implement it using the Metropolis-Hastings Monte Carlo method. Parameters of a new kind are introduced to parameterize inherent imperfections, e.g. , any discrepancies between the theoretical model and measured data. These new parameters enable evaluators to quantify their expert judgement about any discrepancies or imperfections in a reproducible manner. We demonstrate the framework with an ongoing evaluation of 233 U in the eV region using the ENDF-B/VIII library and transmission data measured by Guber, et al. , and compare the posterior parameters to those obtained by conventional evaluation methods.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Electronic and magnetic properties of Mo doped graphene; full potential approach

The electronic and magnetic properties of Pristine and Mo doped Graphene have been calculated using WIEN2k implementation of full potential linearized augmented plane wave (FPLAPW) method based on Density Functional Theory (DFT). The exchange and correlation (XC) effects were taken into account by generalized gradient approximation (GGA). The calculated results show that Mo doping creates magnetism in Graphene by shifting the energy levels at E{sub F} and opens up a channel for Graphene to be used in real nanoscale device applications. The unpaired d-electrons of Mo atom are responsible for induced magnetism in Graphene. Magnetic ordering created in Graphene in this way makes it suitable for recording media, magnetic sensors, magnetic inks and spintronic devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Optimal Bandwidth for Multitaper Spectrum Estimation

A systematic method for bandwidth parameter selection is desired for Thomson multitaper spectrum estimation. We give a method for determining the optimal bandwidth based on a mean squared error (MSE) criterion. When the true spectrum has a second-order Taylor series expansion, one can express quadratic local bias as a function of the curvature of the spectrum, which can be estimated by using a simple spline approximation. This is combined with a variance estimate, obtained by jackknifing over individual spectrum estimates, to produce an estimated MSE for the log spectrum estimate for each choice of time-bandwidth product. The bandwidth that minimizes the estimated MSE then gives the desired spectrum estimate. Additionally, the bandwidth obtained using our method is also optimal for cepstrum estimates. We give an example of a damped oscillatory (Lorentzian) process in which the approximate optimal bandwidth can be written as a function of the damping parameter. Furthermore, the true optimal bandwidth agrees well with that given by minimizing estimated the MSE in these examples.

97 MATHEMATICS AND COMPUTING↗

Using spatial information about recurrence risk for robust optimization of dose-painting prescription functions

Purpose: To develop a robust method for deriving dose-painting prescription functions using spatial information about the risk for disease recurrence. Methods: Spatial distributions of radiobiological model parameters are derived from distributions of recurrence risk after uniform irradiation. These model parameters are then used to derive optimal dose-painting prescription functions given a constant mean biologically effective dose. Results: An estimate for the optimal dose distribution can be derived based on spatial information about recurrence risk. Dose painting based on imaging markers that are moderately or poorly correlated with recurrence risk are predicted to potentially result in inferior disease control when compared the same mean biologically effective dose delivered uniformly. A robust optimization approach may partially mitigate this issue. Conclusions: The methods described here can be used to derive an estimate for a robust, patient-specific prescription function for use in dose painting. Two approximate scaling relationships were observed: First, the optimal choice for the maximum dose differential when using either a linear or two-compartment prescription function is proportional to R, where R is the Pearson correlation coefficient between a given imaging marker and recurrence risk after uniform irradiation. Second, the predicted maximum possible gain in tumor control probability for any robust optimization technique is nearly proportional to the square of R.

61 RADIATION PROTECTION AND DOSIMETRY↗