Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Kernel”

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,729 records · Page 96
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Anderson metal-insulator transitions with classical magnetic impurities

We study the effects of classical magnetic impurities on the Anderson metal-insulator transition (AMIT) numerically. In particular we find that while a finite concentration of Ising impurities lowers the critical value of the site-diagonal disorder amplitude W{sub c}, in the presence of Heisenberg impurities, W{sub c} is first increased with increasing exchange coupling strength J due to time-reversal symmetry breaking. The resulting scaling with J is compared to analytical predictions by Wegner [1]. The results are obtained numerically, based on a finite-size scaling procedure for the typical density of states [2], which is the geometric average of the local density of states. The latter can efficiently be calculated using the kernel polynomial method [3]. Although still suffering from methodical shortcomings, our method proves to deliver results close to established results for the orthogonal symmetry class [4]. We extend previous approaches [5] by combining the KPM with a finite-size scaling analysis. We also discuss the relevance of our findings for systems like phosphor-doped silicon (Si:P), which are known to exhibit a quantum phase transition from metal to insulator driven by the interplay of both interaction and disorder, accompanied by the presence of a finite concentration of magnetic moments [6].

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Excitonic effects in oxyhalide scintillating host compounds

Ab-initio calculations based on density functional theory have been performed to study the electronic, optical, mechanical, and vibrational properties of scintillator host compounds YOX (X = F, Cl, Br, and I). Semiempirical dispersion correction schemes are used to find the effect of van der Waals forces on these layered compounds and we found this effect to be negligible except for YOBr. Calculations of phonons and elastic constants showed that all the compounds studied here are both dynamically and mechanically stable. YOF and YOI are found to be indirect band gap insulators while YOCl and YOBr are direct band gap insulators. The band gap is found to decrease as we move from fluorine to iodine, while the calculated refractive index shows the opposite trend. As the band gap decreases on going down the periodic table from YOF to YOI, the luminescence increases. The excitonic binding energy calculated, within the effective mass approximation, is found to be more for YOF than the remaining compounds, suggesting that the excitonic effect to be more in YOF than the other compounds. The optical properties are calculated within the Time-Dependent Density Functional Theory (TDDFT) and compared with results obtained within the random phase approximation. The TDDFT calculations, using the newly developed bootstrap exchange-correlation kernel, showed significant excitonic effects in all the compounds studied here.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Computational study of collisions between O({sup 3}P) and NO({sup 2}Π) at temperatures relevant to the hypersonic flight regime

Reactions involving N and O atoms dominate the energetics of the reactive air flow around spacecraft when reentering the atmosphere in the hypersonic flight regime. For this reason, the thermal rate coefficients for reactive processes involving O({sup 3}P) and NO({sup 2}Π) are relevant over a wide range of temperatures. For this purpose, a potential energy surface (PES) for the ground state of the NO{sub 2} molecule is constructed based on high-level ab initio calculations. These ab initio energies are represented using the reproducible kernel Hilbert space method and Legendre polynomials. The global PES of NO{sub 2} in the ground state is constructed by smoothly connecting the surfaces of the grids of various channels around the equilibrium NO{sub 2} geometry by a distance-dependent weighting function. The rate coefficients were calculated using Monte Carlo integration. The results indicate that at high temperatures only the lowest A-symmetry PES is relevant. At the highest temperatures investigated (20 000 K), the rate coefficient for the “O1O2+N” channel becomes comparable (to within a factor of around three) to the rate coefficient of the oxygen exchange reaction. A state resolved analysis shows that the smaller the vibrational quantum number of NO in the reactants, the higher the relative translational energy required to open it and conversely with higher vibrational quantum number, less translational energy is required. This is in accordance with Polanyi's rules. However, the oxygen exchange channel (NO2+O1) is accessible at any collision energy. Finally, this work introduces an efficient computational protocol for the investigation of three-atom collisions in general.

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

Hysteresis modeling in graphene field effect transistors

Graphene field effect transistors with an Al{sub 2}O{sub 3} gate dielectric are fabricated on H-intercalated bilayer graphene grown on semi-insulating 4H-SiC by chemical vapour deposition. DC measurements of the gate voltage v{sub g} versus the drain current i{sub d} reveal a severe hysteresis of clockwise orientation. A capacitive model is used to derive the relationship between the applied gate voltage and the Fermi energy. The electron transport equations are then used to calculate the drain current for a given applied gate voltage. The hysteresis in measured data is then modeled via a modified Preisach kernel.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Strain softening mechanism at meso scale during micro-compression in an ultrafine-grained pure copper

Strain softening behavior has been found at meso scale using micro-compression testing in an ultrafine-grained (UFG) pure copper by comparison with the typical strain hardening in conventional coarse-grained (CG) material. Microstructural observations show that grain size remains nearly the same including the fraction of high-angle grain boundaries during micro-compression in UFG pure copper. The Kernel average misorientation(KAM) distribution measured by electron backscatter diffraction (EBSD), as a statistical method, is applied to qualitatively evaluate dislocation density in the interior of the grains. It is suggested that the deformation mechanisms are dominated by grain boundary sliding and grain rotation accompanied by dislocation slip in UFG pure copper, which demonstrates that the strain softening behavior is primarily caused by dislocation annihilation during micro-compression.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Many-body Green’s function theory for electron-phonon interactions: The Kadanoff-Baym approach to spectral properties of the Holstein dimer

We present a Kadanoff-Baym formalism to study time-dependent phenomena for systems of interacting electrons and phonons in the framework of many-body perturbation theory. The formalism takes correctly into account effects of the initial preparation of an equilibrium state and allows for an explicit time-dependence of both the electronic and phononic degrees of freedom. The method is applied to investigate the charge neutral and non-neutral excitation spectra of a homogeneous, two-site, two-electron Holstein model. This is an extension of a previous study of the ground state properties in the Hartree (H), partially self-consistent Born (Gd) and fully self-consistent Born (GD) approximations published in Säkkinen et al. [J. Chem. Phys. 143, 234101 (2015)]. Here, the homogeneous ground state solution is shown to become unstable for a sufficiently strong interaction while a symmetry-broken ground state solution is shown to be stable in the Hartree approximation. Signatures of this instability are observed for the partially self-consistent Born approximation but are not found for the fully self-consistent Born approximation. By understanding the stability properties, we are able to study the linear response regime by calculating the density-density response function by time-propagation. This amounts to a solution of the Bethe-Salpeter equation with a sophisticated kernel. The results indicate that none of the approximations is able to describe the response function during or beyond the bipolaronic crossover for the parameters investigated. Overall, we provide an extensive discussion on when the approximations are valid and how they fail to describe the studied exact properties of the chosen model system.

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

Simulation of macromolecule self-assembly in solution: A multiscale approach

One of the most common processes to produce polymer nanoparticles is to induce self-assembly by using the solvent-displacement method, in which the polymer is dissolved in a “good” solvent and the solution is then mixed with an “anti-solvent”. The polymer ability to self-assemble in solution is therefore determined by its structural and transport properties in solutions of the pure solvents and at the intermediate compositions. In this work, we focus on poly-ε-caprolactone (PCL) which is a biocompatible polymer that finds widespread application in the pharmaceutical and biomedical fields, performing simulation at three different scales using three different computational tools: full atomistic molecular dynamics (MD), population balance modeling (PBM) and computational fluid dynamics (CFD). Simulations consider PCL chains of different molecular weight in solution of pure acetone (good solvent), of pure water (anti-solvent) and their mixtures, and mixing at different rates and initial concentrations in a confined impinging jets mixer (CIJM). Our MD simulations reveal that the nano-structuring of one of the solvents in the mixture leads to an unexpected identical polymer structure irrespectively of the concentration of the two solvents. In particular, although in pure solvents the behavior of the polymer is, as expected, very different, at intermediate compositions, the PCL chain shows properties very similar to those found in pure acetone as a result of the clustering of the acetone molecules in the vicinity of the polymer chain. We derive an analytical expression to predict the polymer structural properties in solution at different solvent compositions and use it to formulate an aggregation kernel to describe the self-assembly in the CIJM via PBM and CFD. Simulations are eventually validated against experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Inter-phase heat transfer and energy coupling in turbulent dispersed multiphase flows

Here, the present paper addresses important fundamental issues of inter-phase heat transfer and energy coupling in turbulent dispersed multiphase flows through scaling analysis. In typical point-particle or two-fluid approaches, the fluid motion and convective heat transfer at the particle scale are not resolved and the momentum and energy coupling between fluid and particles are provided by proper closure models. By examining the kinetic energy transfer due to the coupling forces from the macroscale to microscale fluid motion, closure models are obtained for the contributions of the coupling forces to the energy coupling. Due to the inviscid origin of the added-mass force, its contribution to the microscale kinetic energy does not contribute to dissipative transfer to fluid internal energy as was done by the quasi-steady force. Time scale analysis shows that when the particle is larger than a critical diameter, the diffusive-unsteady kernel decays at a time scale that is smaller than the Kolmogorov time scale. As a result, the computationally costly Basset-like integral form of diffusive-unsteady heat transfer can be simplified to a non-integral form. Conventionally, the fluid-to-particle volumetric heat capacity ratio is used to evaluate the relative importance of the unsteady heat transfer to the energy balance of the particles. Therefore, for gas-particle flows, where the fluid-to-particle volumetric heat capacity ratio is small, unsteady heat transfer is usually ignored. However, the present scaling analysis shows that for small fluid-to-particle volumetric heat capacity ratio, the importance of the unsteady heat transfer actually depends on the ratio between the particle size and the Kolmogorov scale. Furthermore, the particle mass loading multiplied by the heat capacity ratio is usually used to estimate the importance of the thermal two-way coupling effect. Through scaling argument, improved estimates are established for the energy coupling parameters of each energy exchange mechanism between the phases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Exactly solvable problems in the momentum space with a minimum uncertainty in position

A new approach in solution of simple quantum mechanical problems in deformed space with minimal length is presented. We propose the generalization of Schrödinger equation in momentum representation on the case of deformed Heisenberg algebra with minimal length. Assuming that the kernel of potential energy operator does not change in the case of deformation, we obtain exact solution of eigenproblem of a particle in delta potential as well as double delta potential. Particle in Coulomb like potential is revisited and the problem of inversibility and hermicity of inverse coordinate operator is solved.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Relativistic dynamics, Green function and pseudodifferential operators

The central role played by pseudodifferential operators in relativistic dynamics is known very well. In this work, operators like the Schrodinger one (e.g., square root) are treated from the point of view of the non-local pseudodifferential Green functions. Starting from the explicit construction of the Green (semigroup) theoretical kernel, a theorem linking the integrability conditions and their dependence on the spacetime dimensions is given. Relativistic wave equations with arbitrary spin and the causality problem are discussed with the algebraic interpretation of the radical operator and their relation with coherent and squeezed states. Also we perform by means of pure theoretical procedures (based in physical concepts and symmetry) the relativistic position operator which satisfies the conditions of integrability: it is a non-local, Lorentz invariant and does not have the same problems as the “local”position operator proposed by Newton and Wigner. Physical examples, as zitterbewegung and rogue waves, are presented and deeply analyzed in this theoretical framework.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Coherent states, quantum gravity, and the Born- Oppenheimer approximation. II. Compact Lie groups

In this article, the second of three, we discuss and develop the basis of a Weyl quantisation for compact Lie groups aiming at loop quantum gravity-type models. This Weyl quantisation may serve as the main mathematical tool to implement the program of space adiabatic perturbation theory in such models. As we already argued in our first article, space adiabatic perturbation theory offers an ideal framework to overcome the obstacles that hinder the direct implementation of the conventional Born-Oppenheimer approach in the canonical formulation of loop quantum gravity. Additionally, we conjecture the existence of a new form of the Segal-Bargmann-Hall “coherent state” transform for compact Lie groups G, which we prove for G = U(1){sup n} and support by numerical evidence for G = SU(2). The reason for conjoining this conjecture with the main topic of this article originates in the observation that the coherent state transform can be used as a basic building block of a coherent state quantisation (Berezin quantisation) for compact Lie groups G. But, as Weyl and Berezin quantisation for ℝ{sup 2d} are intimately related by heat kernel evolution, it is natural to ask whether a similar connection exists for compact Lie groups as well. Moreover, since the formulation of space adiabatic perturbation theory requires a (deformation) quantisation as minimal input, we analyse the question to what extent the coherent state quantisation, defined by the Segal-Bargmann-Hall transform, can serve as basis of the former.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

The fermionic projector in a time-dependent external potential: Mass oscillation property and Hadamard states

We give a non-perturbative construction of the fermionic projector in Minkowski space coupled to a time-dependent external potential which is smooth and decays faster than quadratically for large times. The weak and strong mass oscillation properties are proven. We show that the integral kernel of the fermionic projector is of the Hadamard form, provided that the time integral of the spatial sup-norm of the potential satisfies a suitable bound. This gives rise to an algebraic quantum field theory of Dirac fields in an external potential with a distinguished pure quasi-free Hadamard state.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Ion parallel closures

Ion parallel closures are obtained for arbitrary atomic weights and charge numbers. For arbitrary collisionality, the heat flow and viscosity are expressed as kernel-weighted integrals of the temperature and flow-velocity gradients.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Charge transfer excitations from particle-particle random phase approximation—Opportunities and challenges arising from two-electron deficient systems

The particle-particle random phase approximation (pp-RPA) is a promising method for studying charge transfer (CT) excitations. Through a detailed analysis on two-electron deficient systems, we show that the pp-RPA is always able to recover the long-distance asymptotic –1/R trend for CT excitations as a result of the concerted effect between orbital energies and the pp-RPA kernel. We also provide quantitative results for systems with relatively short donor-acceptor distances. With conventional hybrid or range-separated functionals, the pp-RPA performs much better than time-dependent density functional theory (TDDFT), although it still gives underestimated results which are not as good as TDDFT with system-dependent tuned functionals. For pp-RPA, there remain three great challenges in dealing with CT excitations. First, the delocalized frontier orbitals in strongly correlated systems often lead to difficulty with self-consistent field convergence as well as an incorrect picture with about half an electron transferred. Second, the commonly used density functionals often underestimate the energy gap between the highest occupied molecular orbital and the lowest unoccupied molecular orbital (LUMO) for the two-electron deficient species, resulting in systems with delocalized orbitals. Third, the performance of pp-RPA greatly depends on the energy difference between the LUMO and a higher virtual orbital. However, the meaning of the orbital energies for higher virtual orbitals is still not clear. We also discuss the performance of an approximate pp-RPA scheme that uses density functional tight binding (pp-DFTB) as reference and demonstrate that the aforementioned challenges can be overcome by adopting suitable range-separated hybrid functionals. Here, the pp-RPA and pp-DFTB are thus promising general approaches for describing charge transfer excitations.

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

Excitation energies from particle-particle random phase approximation with accurate optimized effective potentials

The optimized effective potential (OEP) that gives accurate Kohn-Sham (KS) orbitals and orbital energies can be obtained from a given reference electron density. These OEP-KS orbitals and orbital energies are used here for calculating electronic excited states with the particle-particle random phase approximation (pp-RPA). Our calculations allow the examination of pp-RPA excitation energies with the exact KS density functional theory (DFT). Various input densities are investigated. Specifically, the excitation energies using the OEP with the electron densities from the coupled-cluster singles and doubles method display the lowest mean absolute error from the reference data for the low-lying excited states. This study probes into the theoretical limit of the pp-RPA excitation energies with the exact KS-DFT orbitals and orbital energies. We believe that higher-order correlation contributions beyond the pp-RPA bare Coulomb kernel are needed in order to achieve even higher accuracy in excitation energy calculations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2018

Electron parallel closures for the 3 + 1 fluid model

Linear closures are obtained for arbitrary collisionality for the 3 + 1 fluid model which includes the evolution of density, flow velocity, and pressure both parallel and perpendicular to a preferred direction, usually a magnetic field. A large set of 6400 moment equations is solved to provide closures that are accurate in the collisional regime and well into the collisionless regime. The closures in the collisionless limit are determined by solving the kinetic equation with a model collision operator. Simple fits for the kernel functions that define the closures are obtained for arbitrary collisionality in wave number space. The results are linearly accurate to within 3% across the entire range of collisionality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Accurate double excitations from ensemble density functional calculations

The recent use of a new ensemble in density functional theory (DFT) to produce direct corrections to the Kohn-Sham transitions yields the elusive double excitations that are missed by time-dependent DFT (TDDFT) with the standard adiabatic approximation. But accuracies are lower than for single excitations, and formal arguments about TDDFT suggest that a correction kernel is needed. In principle, ensemble DFT with direct corrections at the exchange level must yield accurate doubles in the weakly correlated limit. We illustrate with exact calculations and analytic results on the Hubbard dimer. Lastly, we also explain the error in formal arguments in TDDFT.

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

Modeling L 2,3 -edge X-ray absorption spectroscopy with linear response exact two-component relativistic time-dependent density functional theory

X-ray absorption spectroscopy (XAS) is a powerful tool that can provide physical insights into element-specific chemical processes and reactivities. Although relativistic time-dependent density functional theory (TDDFT) has been previously applied to model the L-edge region in XAS, there has not been a more comprehensive study of the choices of basis sets and density functional kernels available for variational relativistic excited state methods. In this work, we introduce the implementation of the generalized preconditioned locally harmonic residual algorithm to solve the complex-valued relativistic TDDFT for modeling the L-edge X-ray absorption spectra. Further, we investigate the L 2,3 -edge spectra of a series of molecular complexes using relativistic linear response TDDFT with a hybrid iterative diagonalization algorithm. A systematic error analysis was carried out with a focus on the energetics, intensities, and magnitude of L 2 –L 3 splitting compared to experiments. Additionally, the results from relativistic TDDFT calculations are compared to those computed using other theoretical methods, and the multideterminantal effects on the L-edge XAS were investigated.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗