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,711 records · Page 95

Stochastic modeling of sunshine number data

In this paper, we will present a unified statistical modeling framework for estimation and forecasting sunshine number (SSN) data. Sunshine number has been proposed earlier to describe sunshine time series in qualitative terms (Theor Appl Climatol 72 (2002) 127-136) and since then, it was shown to be useful not only for theoretical purposes but also for practical considerations, e.g. those related to the development of photovoltaic energy production. Statistical modeling and prediction of SSN as a binary time series has been challenging problem, however. Our statistical model for SSN time series is based on an underlying stochastic process formulation of Markov chain type. We will show how its transition probabilities can be efficiently estimated within logistic regression framework. In fact, our logistic Markovian model can be relatively easily fitted via maximum likelihood approach. This is optimal in many respects and it also enables us to use formalized statistical inference theory to obtain not only the point estimates of transition probabilities and their functions of interest, but also related uncertainties, as well as to test of various hypotheses of practical interest, etc. It is straightforward to deal with non-homogeneous transition probabilities in this framework. Very importantly from both physical and practical points of view, logistic Markov model class allows us to test hypotheses about how SSN dependents on various external covariates (e.g. elevation angle, solar time, etc.) and about details of the dynamic model (order and functional shape of the Markov kernel, etc.). Therefore, using generalized additive model approach (GAM), we can fit and compare models of various complexity which insist on keeping physical interpretation of the statistical model and its parts. After introducing the Markovian model and general approach for identification of its parameters, we will illustrate its use and performance on high resolution SSN data from the Solar Radiation Monitoring Station of the West University of Timisoara.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Analytic cubic and quartic force fields using density-functional theory

We present the first analytic implementation of cubic and quartic force constants at the level of Kohn–Sham density-functional theory. The implementation is based on an open-ended formalism for the evaluation of energy derivatives in an atomic-orbital basis. The implementation relies on the availability of open-ended codes for evaluation of one- and two-electron integrals differentiated with respect to nuclear displacements as well as automatic differentiation of the exchange–correlation kernels. We use generalized second-order vibrational perturbation theory to calculate the fundamental frequencies of methane, ethane, benzene, and aniline, comparing B3LYP, BLYP, and Hartree–Fock results. The Hartree–Fock anharmonic corrections agree well with the B3LYP corrections when calculated at the B3LYP geometry and from B3LYP normal coordinates, suggesting that the inclusion of electron correlation is not essential for the reliable calculation of cubic and quartic force constants.

74 ATOMIC AND MOLECULAR PHYSICS↗

Stochastic, real-space, imaginary-time evaluation of third-order Feynman–Goldstone diagrams

A new, alternative set of interpretation rules of Feynman–Goldstone diagrams for many-body perturbation theory is proposed, which translates diagrams into algebraic expressions suitable for direct Monte Carlo integrations. A vertex of a diagram is associated with a Coulomb interaction (rather than a two-electron integral) and an edge with the trace of a Green's function in real space and imaginary time. With these, 12 diagrams of third-order many-body perturbation (MP3) theory are converted into 20-dimensional integrals, which are then evaluated by a Monte Carlo method. It uses redundant walkers for convergence acceleration and a weight function for importance sampling in conjunction with the Metropolis algorithm. The resulting Monte Carlo MP3 method has low-rank polynomial size dependence of the operation cost, a negligible memory cost, and a naturally parallel computational kernel, while reproducing the correct correlation energies of small molecules within a few mE{sub h} after 10{sup 6} Monte Carlo steps.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

2D stochastic-integral models for characterizing random grain noise in titanium alloys

We extend our previous work, in which we applied high-dimensional model representation (HDMR) and analysis of variance (ANOVA) concepts to the characterization of a metallic surface that has undergone a shot-peening treatment to reduce residual stresses, and has, therefore, become a random conductivity field. That example was treated as a onedimensional problem, because those were the only data available. In this study, we develop a more rigorous two-dimensional model for characterizing random, anisotropic grain noise in titanium alloys. Such a model is necessary if we are to accurately capture the 'clumping' of crystallites into long chains that appear during the processing of the metal into a finished product. The mathematical model starts with an application of the Karhunen-Loève (K-L) expansion for the random Euler angles, θ and φ, that characterize the orientation of each crystallite in the sample. The random orientation of each crystallite then defines the stochastic nature of the electrical conductivity tensor of the metal. We study two possible covariances, Gaussian and double-exponential, which are the kernel of the K-L integral equation, and find that the double-exponential appears to satisfy measurements more closely of the two. Results based on data from a Ti-7Al sample will be given, and further applications of HDMR and ANOVA will be discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Using indium tin oxide material to implement the imaging of microwave plasma ignition process

In this paper, a method is introduced to get global observation of microwave plasma ignition process at high pressure. A microwave resonator was designed with an indium tin oxide coated glass at bottom. Microwave plasma ignition was implemented in methane and air mixture at 10 bars by a 2 ms-3 kW-2.45 GHz microwave pulse, and the high speed images of the ignition process were obtained. The images visually proved that microwave plasma ignition could lead to a multi-point ignition. The system may also be applied to obtain Schlieren images, which is commonly used to observe the development of flame kernel in an ignition process.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Assessing electronic structure approaches for gas-ligand interactions in metal-organic frameworks: The CO{sub 2}-benzene complex

Adsorption of gas molecules in metal-organic frameworks is governed by many factors, the most dominant of which are the interaction of the gas with open metal sites, and the interaction of the gas with the ligands. Herein, we examine the latter class of interaction in the context of CO{sub 2} binding to benzene. We begin by clarifying the geometry of the CO{sub 2}–benzene complex. We then generate a benchmark binding curve using a coupled-cluster approach with single, double, and perturbative triple excitations [CCSD(T)] at the complete basis set (CBS) limit. Against this ΔCCSD(T)/CBS standard, we evaluate a plethora of electronic structure approximations: Hartree-Fock, second-order Møller-Plesset perturbation theory (MP2) with the resolution-of-the-identity approximation, attenuated MP2, and a number of density functionals with and without different empirical and nonempirical van der Waals corrections. We find that finite-basis MP2 significantly overbinds the complex. On the other hand, even the simplest empirical correction to standard density functionals is sufficient to bring the binding energies to well within 1 kJ/mol of the benchmark, corresponding to an error of less than 10%; PBE-D in particular performs well. Methods that explicitly include nonlocal correlation kernels, such as VV10, vdW-DF2, and ωB97X-V, perform with similar accuracy for this system, as do ωB97X and M06-L.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ultra-nonlocality in density functional theory for photo-emission spectroscopy

We derive an exact expression for the photocurrent of photo-emission spectroscopy using time-dependent current density functional theory (TDCDFT). This expression is given as an integral over the Kohn-Sham spectral function renormalized by effective potentials that depend on the exchange-correlation kernel of current density functional theory. We analyze in detail the physical content of this expression by making a connection between the density-functional expression and the diagrammatic expansion of the photocurrent within many-body perturbation theory. We further demonstrate that the density functional expression does not provide us with information on the kinetic energy distribution of the photo-electrons. Such information can, in principle, be obtained from TDCDFT by exactly modeling the experiment in which the photocurrent is split into energy contributions by means of an external electromagnetic field outside the sample, as is done in standard detectors. We find, however, that this procedure produces very nonlocal correlations between the exchange-correlation fields in the sample and the detector.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The rigorous stochastic matrix multiplication scheme for the calculations of reduced equilibrium density matrices of open multilevel quantum systems

Understanding the roles of the temporary and spatial structures of quantum functional noise in open multilevel quantum molecular systems attracts a lot of theoretical interests. I want to establish a rigorous and general framework for functional quantum noises from the constructive and computational perspectives, i.e., how to generate the random trajectories to reproduce the kernel and path ordering of the influence functional with effective Monte Carlo methods for arbitrary spectral densities. This construction approach aims to unify the existing stochastic models to rigorously describe the temporary and spatial structure of Gaussian quantum noises. In this paper, I review the Euclidean imaginary time influence functional and propose the stochastic matrix multiplication scheme to calculate reduced equilibrium density matrices (REDM). In addition, I review and discuss the Feynman-Vernon influence functional according to the Gaussian quadratic integral, particularly its imaginary part which is critical to the rigorous description of the quantum detailed balance. As a result, I establish the conditions under which the influence functional can be interpreted as the average of exponential functional operator over real-valued Gaussian processes for open multilevel quantum systems. I also show the difference between the local and nonlocal phonons within this framework. With the stochastic matrix multiplication scheme, I compare the normalized REDM with the Boltzmann equilibrium distribution for open multilevel quantum systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Two-level system in spin baths: Non-adiabatic dynamics and heat transport

We study the non-adiabatic dynamics of a two-state subsystem in a bath of independent spins using the non-interacting blip approximation, and derive an exact analytic expression for the relevant memory kernel. We show that in the thermodynamic limit, when the subsystem-bath coupling is diluted (uniformly) over many (infinite) degrees of freedom, our expression reduces to known results, corresponding to the harmonic bath with an effective, temperature-dependent, spectral density function. We then proceed and study the heat current characteristics in the out-of-equilibrium spin-spin-bath model, with a two-state subsystem bridging two thermal spin-baths of different temperatures. We compare the behavior of this model to the case of a spin connecting boson baths, and demonstrate pronounced qualitative differences between the two models. Specifically, we focus on the development of the thermal diode effect, and show that the spin-spin-bath model cannot support it at weak (subsystem-bath) coupling, while in the intermediate-strong coupling regime its rectifying performance outplays the spin-boson model.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Electronic nose with a new feature reduction method and a multi-linear classifier for Chinese liquor classification

An electronic nose (e-nose) was designed to classify Chinese liquors of the same aroma style. A new method of feature reduction which combined feature selection with feature extraction was proposed. Feature selection method used 8 feature-selection algorithms based on information theory and reduced the dimension of the feature space to 41. Kernel entropy component analysis was introduced into the e-nose system as a feature extraction method and the dimension of feature space was reduced to 12. Classification of Chinese liquors was performed by using back propagation artificial neural network (BP-ANN), linear discrimination analysis (LDA), and a multi-linear classifier. The classification rate of the multi-linear classifier was 97.22%, which was higher than LDA and BP-ANN. Finally the classification of Chinese liquors according to their raw materials and geographical origins was performed using the proposed multi-linear classifier and classification rate was 98.75% and 100%, respectively.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A fast non-Fourier method for Landau-fluid operators

An efficient and versatile non-Fourier method for the computation of Landau-fluid (LF) closure operators [Hammett and Perkins, Phys. Rev. Lett. 64, 3019 (1990)] is presented, based on an approximation by a sum of modified-Helmholtz-equation solves (SMHS) in configuration space. This method can yield fast-Fourier-like scaling of the computational time requirements and also provides a very compact data representation of these operators, even for plasmas with large spatial nonuniformity. As a result, the method can give significant savings compared with direct application of “delocalization kernels” [e.g., Schurtz et al., Phys. Plasmas 7, 4238 (2000)], both in terms of computational cost and memory requirements. The method is of interest for the implementation of Landau-fluid models in situations where the spatial nonuniformity, particular geometry, or boundary conditions render a Fourier implementation difficult or impossible. Systematic procedures have been developed to optimize the resulting operators for accuracy and computational cost. The four-moment Landau-fluid model of Hammett and Perkins has been implemented in the BOUT++ code using the SMHS method for LF closure. Excellent agreement has been obtained for the one-dimensional plasma density response function between driven initial-value calculations using this BOUT++ implementation and matrix eigenvalue calculations using both Fourier and SMHS non-Fourier implementations of the LF closures. The SMHS method also forms the basis for the implementation, which has been carried out in the BOUT++ code, of the parallel and toroidal drift-resonance LF closures. The method is a key enabling tool for the extension of gyro-Landau-fluid models [e.g., Beer and Hammett, Phys. Plasmas 3, 4046 (1996)] to codes that treat regions with strong profile variation, such as the tokamak edge and scrapeoff-layer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Stochastic models for surface diffusion of molecules

We derive a stochastic model for the surface diffusion of molecules, starting from the classical equations of motion for an N-atom molecule on a surface. The equation of motion becomes a generalized Langevin equation for the center of mass of the molecule, with a non-Markovian friction kernel. In the Markov approximation, a standard Langevin equation is recovered, and the effect of the molecular vibrations on the diffusion is seen to lead to an increase in the friction for center of mass motion. This effective friction has a simple form that depends on the curvature of the lowest energy diffusion path in the 3N-dimensional coordinate space. We also find that so long as the intramolecular forces are sufficiently strong, memory effects are usually not significant and the Markov approximation can be employed, resulting in a simple one-dimensional model that can account for the effect of the dynamics of the molecular vibrations on the diffusive motion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

On the oscillator realization of conformal U(2, 2) quantum particles and their particle-hole coherent states

We revise the unireps. of U(2, 2) describing conformal particles with continuous mass spectrum from a many-body perspective, which shows massive conformal particles as compounds of two correlated massless particles. The statistics of the compound (boson/fermion) depends on the helicity h of the massless components (integer/half-integer). Coherent states (CS) of particle-hole pairs (“excitons”) are also explicitly constructed as the exponential action of exciton (non-canonical) creation operators on the ground state of unpaired particles. These CS are labeled by points Z (2×2 complex matrices) on the Cartan-Bergman domain D₄=U(2,2)/U(2)², and constitute a generalized (matrix) version of Perelomov U(1, 1) coherent states labeled by points z on the unit disk D₁=U(1,1)/U(1)². First, we follow a geometric approach to the construction of CS, orthonormal basis, U(2, 2) generators and their matrix elements and symbols in the reproducing kernel Hilbert space H{sub λ}(D₄) of analytic square-integrable holomorphic functions on D₄, which carries a unitary irreducible representation of U(2, 2) with index λϵN (the conformal or scale dimension). Then we introduce a many-body representation of the previous construction through an oscillator realization of the U(2, 2) Lie algebra generators in terms of eight boson operators with constraints. This particle picture allows us for a physical interpretation of our abstract mathematical construction in the many-body jargon. In particular, the index λ is related to the number 2(λ – 2) of unpaired quanta and to the helicity h = (λ – 2)/2 of each massless particle forming the massive compound.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

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↗

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↗

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↗

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↗

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↗