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,531 records · Page 85

Periodicity property of relativistic Thomson scattering with application to exact calculations of angular and spectral distributions of the scattered field

We prove that the analytical expression of the intensity of the relativistic Thomson scattered field for a system composed of an electron interacting with a plane electromagnetic field can be written in the form of a composite periodic function of only one variable, that is, the phase of the incident field. This property is proved without using any approximation in the most general case in which the field is elliptically polarized, the initial phase of the incident field and the initial velocity of the electron are taken into consideration, and the direction in which the radiation is scattered is arbitrary. This property leads to an exact method for calculating the angular and spectral distributions of the scattered field, which reveals a series of physical details of these distributions, such as their dependence on the components of the initial electron velocity. Since the phase of the field is a relativistic invariant, it follows that the periodicity property is also valid when the analysis is made in the inertial system in which the initial velocity of the electron is zero in the case of interactions between very intense electromagnetic fields and relativistic electrons. Consequently, the calculation method can be used for the evaluation of properties of backscattered hard radiations generated by this type of interaction. The theoretical evaluations presented in this paper are in good agreement with the experimental data from literature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Analytic continuation of scattering data to the region of negative energies for systems that have one and two bound states

An exactly solvable potential model is used to study the possibility of deducing information about the features of bound states for the system under consideration (binding energies and asymptotic normalization coefficients) on the basis of data on continuum states. The present analysis is based on an analytic approximation and on the subsequent continuation of a partial-wave scattering function from the region of positive energies to the region of negative energies. Cases where the system has one or two bound states are studied. The α+d and α+{sup 12}C systems are taken as physical examples. In the case of one bound state, the scattering function is a smooth function of energy, and the procedure of its analytic continuation for different polynomial approximations leads to close results, which are nearly coincident with exact values. In the case of two bound states, the scattering function has two poles—one in the region of positive energies and the other in the region of negative energies between the energies corresponding to the two bound states in question. Padéapproximants are used to reproduce these poles. The inclusion of these poles proves to be necessary for correctly describing the properties of the bound states.

74 ATOMIC AND MOLECULAR PHYSICS↗

Electronic Structure and Related Band Parameters of Hexagonal Wurtzite Zn{sub 1−x}Mn{sub x}O Magnetic Semiconducting Alloys

The electronic band structure and related band parameters such as, lattice parameters, energy band gaps, density of states and magnetic moment of Zn{sub 1−x}Mn{sub x}O ternary system in the hypothetical hexagonal wurtzite structure have been investigated. The calculations are performed using a full-potential linearized augmented plane wave method based on density functional theory in the generalized gradient approximation of Perdew-Burke-Ernzerhof. The evolution of the features of interest with increasing the Mn doping concentration from 0 to 1 in the ZnO structure has been examined and discussed. Our findings show that the behaviour of all studied band parameters with respect to the alloy composition x is consistent with those reported experimentally in the literature. The present study may be useful for spintronics applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Plasmoid instability in high-Lundquist-number magnetic reconnection

Our understanding of magnetic reconnection in resistive magnetohydrodynamics has gone through a fundamental change in recent years. The conventional wisdom is that magnetic reconnection mediated by resistivity is slow in laminar high Lundquist (S) plasmas, constrained by the scaling of the reconnection rate predicted by Sweet-Parker theory. However, recent studies have shown that when S exceeds a critical value ∼10{sup 4}, the Sweet-Parker current sheet is unstable to a super-Alfvénic plasmoid instability, with a linear growth rate that scales as S{sup 1/4}. In the fully developed statistical steady state of two-dimensional resistive magnetohydrodynamic simulations, the normalized average reconnection rate is approximately 0.01, nearly independent of S, and the distribution function f(ψ) of plasmoid magnetic flux ψ follows a power law f(ψ)∼ψ{sup −1}. When Hall effects are included, the plasmoid instability may trigger onset of Hall reconnection even when the conventional criterion for onset is not satisfied. The rich variety of possible reconnection dynamics is organized in the framework of a phase diagram.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Temperature dependence of nuclear quantum effects on liquid water via artificial neural network model based on SCAN meta-GGA functional

We investigate the temperature dependence of nuclear quantum effects (NQEs) on structural and dynamic properties of liquid water by training a neural network force field using first-principles molecular dynamics (FPMD) based on the strongly constrained and appropriately normed meta-generalized gradient approximation exchange-correlation approximation. The FPMD simulation based on density functional theory has become a powerful computational approach for studying a wide range of condensed phase systems. However, its large computational cost makes it difficult to incorporate NQEs in the simulation and investigate temperature dependence of various properties. To circumvent this difficulty, we use an artificial neural network model and employ the thermostatted ring polymer MD approach for studying the temperature dependence of NQEs on various properties. The NQEs generally bring the radial distribution functions closer to the experimental measurements. Translational diffusivity and rotational dynamics of water molecules are both slowed down by the NQEs. The competing inter-molecular and intra-molecular quantum effects on hydrogen bonds, as discussed by Habershon, Markland, and Manolopoulos [J. Chem. Phys. 131(2), 024501 (2019)], can explain the observed temperature dependence of the NQEs on the dynamical properties in our simulation.

Yao, Yi (ORCID:0000000299576979)↗

Range-separated thermal hybrid exchange-correlation density functional for accurate band-gap calculations of warm dense matter

Here, in this study, we present a range-separated thermal hybrid exchange-correlation (XC) functional, termed RS-KDT0, which provides highly accurate band-gap predictions at elevated temperature while reducing the computational cost relative to thermal global hybrid functionals. At zero temperature, RS-KDT0 recovers the widely used ground-state range-separated hybrid functional HSE06. To assess its performance across temperature regimes, we conducted static calculations of electronic pressure, band gaps, and band structures. RS-KDT0 exhibits the same level of accuracy as HSE06 in the low-temperature regime and demonstrates improved agreement with finite-temperature GW calculations at elevated temperatures when compared to other thermal XC functionals such as the generalized gradient approximation (GGA)-level KDT16 and the global hybrid KDT0.

density functional theory↗

Double beta decays of {sup 106}Cd

The two-neutrino (2{nu}2{beta}) and neutrinoless (0{nu}2{beta}) double beta decays of {sup 106}Cd are studied for the transitions to the ground state 0{sub gs}{sup +} and 0{sup +} and 2{sup +} excited states in {sup 106}Pd by using realistic many-body wave functions calculated in the framework of the quasiparticle random-phase approximation. Effective, G-matrix-derived nuclear forces are used in realistic single-particle model spaces. All the possible channels, {beta}{sup +}{beta}{sup +}, {beta}{sup +}EC, and ECEC, are discussed for both the 2{nu}2{beta} and 0{nu}2{beta} decays. The associated half-lives are computed and particular attention is devoted to the study of the detectability of the resonant neutrinoless double electron capture (R0{nu}ECEC) process in {sup 106}Cd. The calculations of the present article constitute the thus far most complete and up-to-date investigation of the double-beta-decay properties of {sup 106}Cd.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Self-interaction effects on charge-transfer collisions

In this article, we investigate the role of the self-interaction error in the simulation of collisions using time-dependent density functional theory (TDDFT) and Ehrenfest dynamics. In addition, we compare many different approximations of the exchange and correlation potential, using as a test system the collision of H + + CH 4 at 30 eV. We find that semilocal approximations, like the Perdew-Burke- Ernzerhof (PBE), and even hybrid functionals, such as the Becke, 3-parameter, Lee-Yang-Parr (B3LYP), produce qualitatively incorrect predictions for the scattering of the proton. This discrepancy appears because the self-interaction error allows the electrons to jump too easily to the proton, leading to radically different forces with respect to the non-self-interacting case. Lastly, from our results, we conclude that using a functional that is self-interaction free is essential to properly describing charge-transfer collisions between ions and molecules in TDDFT.

74 ATOMIC AND MOLECULAR PHYSICS↗

Energy band alignment and electronic states of amorphous carbon surfaces in vacuo and in aqueous environment

In this paper, we obtain the energy band positions of amorphous carbon (a–C) surfaces in vacuum and in aqueous environment. The calculations are performed using a combination of (i) classical molecular dynamics (MD), (ii) Kohn-Sham density functional theory with the Perdew-Burke-Ernzerhof (PBE) exchange-correlation functional, and (iii) the screened-exchange hybrid functional of Heyd, Scuseria, and Ernzerhof (HSE). PBE allows an accurate generation of a-C and the evaluation of the local electrostatic potential in the a-C/water system, HSE yields an improved description of energetic positions which is critical in this case, and classical MD enables a computationally affordable description of water. Our explicit calculation shows that, both in vacuo and in aqueous environment, the a-C electronic states available in the region comprised between the H{sub 2}/H{sub 2}O and O{sub 2}/H{sub 2}O levels of water correspond to both occupied and unoccupied states within the a-C pseudogap region. These are localized states associated to sp{sup 2} sites in a-C. The band realignment induces a shift of approximately 300 meV of the a-C energy band positions with respect to the redox levels of water.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Investigation of Transport Parameters of Graphene-Based Nanostructures

The paper presents results of computer simulation of the main transport parameters of nanostructures obtained through the row-by-row removal of carbon atoms from graphene ribbon. Research into the electrical parameters is carried out within the density functional theory using the non-equilibrium Green functions in the local-density approximation. Virtual NanoLab based on Atomistix ToolKit is used to construct structures and analyze simulation results. Current-voltage characteristics, differential conductivity and transmittance spectra of nanostructures are calculated at different values of bias voltage. It is found that there is a large region of negative differential resistance in current-voltage characteristics of nanostructures caused by resonant tunneling of quasi-particles. Differential (dI/dV) characteristic also has similar changes. The obtained results can be useful for building novel electronic devices in the field of nanoelectronics.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

A phase-based interior penalty discontinuous Galerkin method for the Helmholtz equation with spatially varying wavenumber

This paper is concerned with an interior penalty discontinuous Galerkin (IPDG) method based on a flexible type of non-polynomial local approximation space for the Helmholtz equation with varying wavenumber. The local approximation space consists of multiple polynomial-modulated phase functions which can be chosen according to the phase information of the solution. We obtain some approximation properties for this space and a priori L 2 error estimates for the h-convergence of the IPDG method using duality argument. We also provide ample numerical examples to show that, building phase information into the local spaces often gives more accurate results comparing to using the standard polynomial spaces.

42 ENGINEERING↗

Structural properties of small rhodium clusters

We report a systematic study of the structural properties of rhodium clusters at the atomistic level. A novel global-minimum search algorithm, known as parallel tempering multicanonical basin hopping plus genetic algorithm (PTMBHGA), is used to obtain the geometrical structures with lowest minima at the semi-empirical level where Gupta potential is used to describe the atomic interaction among the rhodium atoms. These structures are then re-optimized at the density functional theory (DFT) level with exchange-correlation energy approximated by Perdew-Burke-Ernzerhof (PBE) generalized gradient approximation (GGA). The structures are optimized for different spin multiplicities. The ones with lowest energies will be taken as ground-state structures. In most cases, we observe only minor changes in the geometry and bond length of the clusters as a result of DFT-level re-optimization. Only in some limited cases, the initial geometries obtained from the PTMBHGA are modified by the re-optimization. The variation of structural properties, such as ground-state geometry, symmetry and binding energy, with respect to the cluster size is studied and agreed well with other results available in the literature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The snapshot Hubble U-band cluster survey (SHUCS). II. The star cluster population of NGC 2997

We study the star cluster population of NGC 2997, a giant spiral galaxy located at 9.5 Mpc and targeted by the Snapshot Hubble U-band Cluster Survey (SHUCS). Combining our U-band imaging from SHUCS with archival BVI imaging from the Hubble Space Telescope, we select a high confidence sample of clusters in the circumnuclear ring and disk through a combination of automatic detection procedures and visual inspection. The cluster luminosity functions in all four filters can be approximated by power laws with indices of –1.7 to –2.3. Some deviations from pure power-law shape are observed, hinting at the presence of a high-mass truncation in the cluster mass function. However, upon inspection of the cluster mass function, we find it is consistent with a pure power law of index –2.2 ± 0.2 despite a slight bend at ∼2.5 × 10{sup 4} M {sub ☉}. No statistically significant truncation is observed. From the cluster age distributions, we find a low rate of disruption (ζ ∼ –0.1) in both the disk and circumnuclear ring. Finally, we estimate the cluster formation efficiency (Γ) over the last 100 Myr in each region, finding 7% ± 2% for the disk, 12% ± 4% for the circumnuclear ring, and 10% ± 3% for the entire UBVI footprint. This study highlights the need for wide-field UBVI coverage of galaxies to study cluster populations in detail, though a small sample of clusters can provide significant insight into the characteristics of the population.

79 ASTRONOMY AND ASTROPHYSICS↗

Hydrogen Effect on Electron-Phonon Interactions in L1{sub 0} FePd

The effect of H position in L1{sub 0}FePd magnetic equiatomic layered structure on electron-phonon interactions is studied via ab initio pseudopotentials density functional approach in the general gradient approximation. The structures investigated are L1{sub 0}FePd (parent structure), and the specific hydrides 2(FePd)1H (H@Fe layer) and 2(FePd)2H (H@Pd layer). At ground state (0 K), electronic band structure results demonstrate the emergence of Fermi surfaces (FS) in topologies of small isolated multipockets in the hydrides. Using the harmonic approximation I find the phonon energy dispersion of parent and H@Fe simple and classical but for H@Pd complex. I linked this latter structure to its topology of FS, known as Kohn anomaly. Using McMillan-Eliashberg model, I get superconducting transition temperature (T{sub c}) values ≈ 0.00 K, 3.48 K, and 19.75 K for the parent, H@Fe and H@Pd respectively. Therefore, these hypothetical hydrides prove that H position in the structure has a direct influence on T{sub c} values and consequently the suppression of magnetism.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Quasielastic electron-deuteron scattering in the weak binding approximation

We perform a global analysis of all available electron-deuteron quasielastic scattering data using Q^2-dependent smearing functions that describe inclusive inelastic e-d scattering within the weak binding approximation. We study the dependence of the cross sections on the deuteron wave function and the off-shell extrapolation of the elastic electron-nucleon cross section, which show particular sensitivity at x >> 1. The excellent overall agreement with data over a large range of Q^2 and x suggest a limited need for effects beyond the impulse approximation, with the exception of the very high-x or very low-Q^2 regions, where short-distance effects in the deuteron become more relevant.

Ethier, Jacob J.↗

Electronic structure calculations of PbS quantum rods and tubes

We study absorption spectra, optical and HOMO-LUMO gaps, and the density of states for PbS quantum rods (QRs) and tubes (QTs). We find some similarities and also differences in QR and QT properties. For both QRs and QTs, the optical and HOMO-LUMO gaps reach the plateaus for small lengths. We find that tubes are as stable as rods. The optical spectra exhibit a peak that can be due to the electron-hole interaction or be a prototype of an S{sub e}–S{sub h} transition in the effective mass approximation. We also calculate the density of states by the density functional theory (DFT) and time-dependent density functional theory (TDDFT) methods. The TDDFT density of states function is shifted towards the red side by 0.5 eV indicating the strong e-h interaction.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Approximating Photo- z PDFs for Large Surveys

Modern galaxy surveys produce redshift probability density functions (PDFs) in addition to traditional photometric redshift (photo-z) point estimates. However, the storage of photo-z PDFs may present a challenge with increasingly large catalogs, as we face a trade-off between the accuracy of subsequent science measurements and the limitation of finite storage resources. This work presents qp, a Python package for manipulating parameterizations of one-dimensional PDFs, as suitable for photo-z PDF compression. We use qp to investigate the performance of three simple PDF storage formats (quantiles, samples, and step functions) as a function of the number of stored parameters on two realistic mock data sets, representative of upcoming surveys with different data qualities. We propose some best practices for choosing a photo-z PDF approximation scheme and demonstrate the approach on a science case using performance metrics on both ensembles of individual photo-z PDFs and an estimator of the overall redshift distribution function. We show that both the properties of the set of PDFs we wish to approximate and the fidelity metric(s) chosen affect the optimal parameterization. Additionally, we find that quantiles and samples outperform step functions, and we encourage further consideration of these formats for PDF approximation.

79 ASTRONOMY AND ASTROPHYSICS↗

Revisiting the Fundamental Analytical Solutions of Heat and Mass Transfer: The Kernel of Multirate and Multidimensional Diffusion

There are two types of analytical solutions of temperature/concentration in and heat/mass transfer through boundaries of regularly shaped 1D, 2D, and 3D blocks. These infinite-series solutions with either error functions or exponentials exhibit highly irregular but complementary convergence at different dimensionless times, t d0 . In this paper, approximate solutions were developed by combining the error-function-series solutions for early times and the exponential-series solutions for late times and by using time partitioning at the switchover time, t d0 . The combined solutions contain either the leading term of both series for normal-accuracy approximations (with less than 0.003 relative error) or the first two terms for high-accuracy approximations (with less than 10-7 relative error) for 1D isotropic (spheres, cylinders, slabs) and 2D/3D rectangular blocks (squares, cubes, rectangles, and rectangular parallelepipeds). This rapid and uniform convergence for rectangular blocks was achieved by employing the same time partitioning with individual dimensionless times for different directions and the product of their combined 1D slab solutions. The switchover dimensionless time was determined to minimize the maximum approximation errors. Furthermore, the analytical solutions of first-order heat/mass flux for 2D/3D rectangular blocks were derived for normal-accuracy approximations. These flux equations contain the early-time solution with a three-term polynomial in √td and the late-time solution with the limited-term exponentials for rectangular blocks. The heat/mass flux equations and the combined temperature/concentration solutions form the ultimate kernel for fast simulations of multirate and multidimensional heat/mass transfer in porous/fractured media with millions of low-permeability blocks of varying shapes and sizes.

42 ENGINEERING↗