SEARCH · Engineering Papers
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.
The dynamics of a plasma crystal under the action of random external forces is considered. The pair correlation functions of the particle displacements are calculated in the harmonic approximation by using the Langevin equations. The case of a planar hexagonal lattice is analyzed in more detail. Analogues of the Van Hove singularities in the spectral densities are discovered.
Recent developments in theoretical techniques have significantly improved the predictive power of density-functional-based calculations. In this article, we discuss how such advancements have enabled improved understanding of native point defects in GaN. We review the methodologies for the calculation of point defects, and discuss how techniques for overcoming the band-gap problem of density functional theory affect native defect calculations. In particular, we examine to what extent calculations performed with semilocal functionals (such as the generalized gradient approximation), combined with correction schemes, can produce accurate results. The properties of vacancy, interstitial, and antisite defects in GaN are described, as well as their interaction with common impurities. We also connect the first-principles results to experimental observations, and discuss how native defects and their complexes impact the performance of nitride devices. Overall, we find that lower-cost functionals, such as the generalized gradient approximation, combined with band-edge correction schemes can produce results that are qualitatively correct. However, important physics may be missed in some important cases, particularly for optical transitions and when carrier localization occurs.
This paper presents an extension to the Wave Based Method to predict the absorption, reflection and transmission coefficients of a porous material with an embedded periodic set of inclusions. The porous unit cell is described using the Multi-Level methodology and by embedding Bloch–Floquet periodicity conditions in the weighted residual scheme. The dynamic pressure field in the semi-infinite acoustic domains is approximated using a novel wave function set that fulfils the Helmholtz equation, the Bloch–Floquet periodicity conditions and the Sommerfeld radiation condition. The method is meshless and computationally efficient, which makes it well suited for optimisation studies.
Ab-initio density functional theory based calculations of electronic properties of bulk and monolayer Molybdenum di-Sulfide (MoS{sub 2}) have been performed using all electron Full Potential Linearised Augmentad Plane Wave (FPLAPW) method using Elk code. We have used Generalised Gradient Approximation (GGA) for exchange and correlation functionals and performed calculaitons of Lattice parameters, Density Of States (DOS) and Band Structure (BS). Band structure calculations revealed that bulk MoS{sub 2} has indirect band gap of 0.97 eV and mono-layer MoS{sub 2} has direct band gap which has increased to 1.71 eV. These are in better agreement with experimental values as compared with the other calculations using pseudo-potential code. The effect of mechanical strain on the electronic properties of bulk MoS{sub 2} has also been studied. For the different values of compressive strain (varying from 2% to 8% in steps of 2%) along the c-axis, the corresponding DOS and BS are obtained. We observed that the band gap decreases by about 15% for every 2% increase in strain along the c-axis.
We perform first-principle calculations of electron-nucleus scattering on 3 He and 3 H using the Green's function Monte Carlo method and two approaches based on the factorization of the final hadronic state: the spectral-function formalism and the short-time approximation. These three methods are benchmarked among each other and compared to the experimental data for the longitudinal and transverse electromagnetic response functions of 3 He, and the inclusive cross sections of both 3 He and 3 H. Since these three approaches are based on the same description of nuclear dynamics of the initial target state, comparing their results enables a precise quantification of the uncertainties inherent to factorization schemes. At sufficiently large values of the momentum transfer, we find an excellent agreement of the Green's function Monte Carlo calculation with experimental data and with both the spectral-function formalism and the short-time approximation. Here, we also analyze the relevance of relativistic effects, whose inclusion becomes crucial to explain data at high momentum and energy transfer.
We perform first-principle calculations of electron-nucleus scattering on $^3$He and $^3$H using the Green's function Monte Carlo method and two approaches based on the factorization of the final hadronic state: the spectral-function formalism and the short-time approximation. These three methods are benchmarked among each other and compared to the experimental data for the longitudinal and transverse electromagnetic response functions of $^3$He, and the inclusive cross sections of both $^3$He and $^3$H. Since these three approaches are based on the same description of nuclear dynamics of the initial target state, comparing their results enables a precise quantification of the uncertainties inherent to factorization schemes. At sufficiently large values of the momentum transfer, we find an excellent agreement of the Green's function Monte Carlo calculation with experimental data and with both the spectral-function formalism and the short-time approximation. We also analyze the relevance of relativistic effects, whose inclusion becomes crucial to explain data at high momentum and energy transfer.
Recently, a series of double-perovskite halide compounds such as Cs 2 AgBiCl 6 and Cs 2 AgBiBr 6 have attracted intensive interest as promising alternatives to the solar absorber material CH 3 NH 3 PbI 3 because they are Pb-free and may exhibit enhanced stability. The thermodynamic stability of a number of double-perovskite halides has been predicted based on density functional theory (DFT) calculations of compound formation energies. In this paper, we found that the stability prediction can be dependent on the approximations used for the exchange-correlation functionals, e.g., the DFT calculations using the widely used Perdew, Burke, Ernzerhof (PBE) functional predict that Cs 2 AgBiBr 6 is thermodynamically unstable against phase-separation into the competing phases such as AgBr, Cs 2 AgBr 3 , Cs 3 Bi 2 Br 9 , etc., obviously inconsistent with the good stability observed experimentally. The incorrect prediction by the PBE calculation results from its failure to predict the correct ground-state structures of AgBr, AgCl, and CsCl. By contrast, the DFT calculations based on local density approximation, optB86b-vdW, and optB88-vdW functionals predict the ground-state structures of these binary halides correctly. Furthermore, the optB88-vdW functional is found to give the most accurate description of the lattice constants of the double-perovskite halides and their competing phases. Given these two aspects, we suggest that the optB88-vdW functional should be used for predicting thermodynamic stability in the future high-throughput computational material design or the construction of the Materials Genome database for new double-perovskite halides. As a result, using different exchange-correlation functionals has little influence on the dispersion of the conduction and the valence bands near the electronic bandgap; however, the calculated bandgap can be affected indirectly by the optimized lattice constant, which varies for different functionals.
Here, we explore non-adiabatic particle production in a de Sitter universe for a scalar spectator field, by allowing the effective mass m 2 (t) of this field and the cosmic time interval between non-adiabatic events to vary stochastically. Two main scenarios are considered depending on the (non-stochastic) mass M of the spectator field: the conformal case with M 2 =2H 2 , and the case of a massless field. We make use of the transfer matrix formalism to parametrize the evolution of the system in terms of the "occupation number", and two phases associated with the transfer matrix; these are used to construct the evolution of the spectator field. Assuming short-time interactions approximated by Dirac-delta functions, we numerically track the change of these parameters and the field in all regimes: sub- and super-horizon with weak and strong scattering. In all cases a log-normally distributed field amplitude is observed, and the logarithm of the field amplitude approximately satisfies the properties of a Wiener process outside the horizon. We derive a Fokker-Planck equation for the evolution of the transfer matrix parameters, which allows us to calculate analytically non-trivial distributions and moments in the weak-scattering limit.
The structural, electronic, magnetic, and mechanical properties of CoVTiAs and CoVTiSi quaternary Heusler alloys are studied using the full-potential linearized augmented plane-wave (FP-LAPW) method in the framework of the density functional theory (DFT). The generalized gradient approximation (GGA) is chosen for the exchange-correlation energy, whereas the modified Becke-Johnson (mBJ) formalism is applied for the electronic properties. Both CoVTiAs and CoVTiSi are stable in the type 2 (FM) structure. The mechanical properties of CoVTiX (X = As, Si) are predicted from the calculated elastic constants. Results identify that these compounds are mechanically stable. Our results with the mBJ approximation predict that CoVTiAs and CoVTiSi alloys are half-metallic ferromagnets with band gaps of 0.450 and 0.424 eV, respectively.
The F E integral for the normalized Boltzmann-weighted number of unimolecular states above the threshold energy is the key quantity for computing the collision efficiency in the pressure-dependent unimolecular rate theory developed by Troe, who calls this the energy dependence factor of the density of states. By using the Whitten–Rabinovitch approximation and assuming that the Whitten–Rabinovitch a(E) function is independent of energy, F E can be approximated by an analytical formula; this approximate formula is widely used because of its convenience and computational efficiency. Here we test its validity by comparing the rate constants computed by using the approximate F E to the ones determined by using the numerically integrated F E . For small-sized molecules and for reactions with high threshold energies E 0 , the differences are negligible at all temperatures, but in other cases, the approximate formula tends to underestimate F E and thus overestimates the collision efficiency, and this leads to smaller pressure falloff. When a(E) at high energies differs appreciably from a(E 0 ), we find that the underestimation of pressure-dependent rate constants by using the approximate formula can be greater than a factor of 5 at high temperatures. The physical insight we draw from this study is that, for reactions with threshold energies below about 30 kcal mol -1 , the rate of collisional energy transfer can be appreciably slowed down by the increase in the density of states at higher energies, and this increases the falloff effect by which finite-pressure rate constants become lower than the high-pressure limit, especially at higher temperatures.
The structural, electronic, thermal, and optical properties of cubic non magnetic CeX{sub 3}(X=Pd, Rh) compounds which crystallize in the Au{sub 3}Cu structure have been studied using the projected augmented wave (PAW) method within the density functional theory (DFT) with generalized gradient approximation (GGA) for exchange correlation potential. In this paper we have calculated the band structure which are interpreted using the density of states. The optical properties such as extinction coefficients clearly illustrate the changes in CeX{sub 3} due to intercalation of boron. Lattice instability is observed in CePd{sub 3}B from the calculated dynamical properties.
We introduce a novel black hole mass function that realistically models the physics of pair-instability supernovae with a minimal number of parameters. Applying this to all events in the LIGO-Virgo Gravitational-Wave Transient Catalog 2 (GWTC-2), we detect a peak at M BHMG = 46$^{+17}_{-6}$ M ⊙ . Repeating the analysis without the black holes from the event GW190521, we find this feature at M BHMG =54± 6 M ⊙ . These results establish the edge of the anticipated "black hole mass gap" at a value compatible with the expectation from standard stellar structure theory. The mass gap manifests itself as a discontinuity in the mass function and is populated by a distinct, less-abundant population of higher-mass black holes. We find that the primary black hole population scales with power-law index -1.95 ± 0.51 (-1.97 ± 0.44) with (without) GW190521, consistent with models of star formation. Using Bayesian techniques, we establish that our mass function fits a new catalog of black hole masses approximately as well as pre-existing phenomenological mass functions. We also remark on the implications of these results for constraining or discovering new phenomena in nuclear and particle physics.
In this study, a new computational approach is presented to solve the generalized pantograph-delay differential equations (PDDEs). The solutions obtained by our scheme represent by a linear combination of a special kind of basis functions, and can be deduced in a straightforward manner. Firstly, using the least squares approximation method and the Lagrange-multiplier method, the given PDDE is converted to a linear system of algebraic equations, and those unknown coefficients of the solution of the problem are determined by solving this linear system. Secondly, a PDDE related to the error function of the approximate solution is constructed based on the residual error function technique, and error estimation is presented for the suggested method. The convergence of the approximate solution is proved. Several numerical examples are given to demonstrate the accuracy and efficiency. Comparisons are made between the proposed method and other existing methods.
Two-neutrino double-beta (2νββ) decay has been used to constrain the neutron-proton part of effective interactions, which in turn is used to compute the nuclear matrix elements for neutrinoless double-beta decay, the observation of which would have important consequences for fundamental physics. We carefully examine 2νββ matrix elements within the proton-neutron quasiparticle random-phase approximation with nuclear energy density functionals. Here we work with functionals that are fit globally to single-beta-decay half-lives and charge-exchange giant-resonance energies, but not to 2νββ half-lives themselves, to evaluate the 2νββ nuclear matrix elements for all important nuclei, including those whose half-lives have not yet been measured. Such a comprehensive evaluation in large model spaces without configuration truncation requires an efficient computational scheme; we employ a double contour integration within the finite amplitude method. The results generally reproduce the nuclear matrix element extracted from half-lives well, without the use of any of those half-lives in the fitting procedure. We present predictions of the matrix elements in a total of 27 nuclei with half-lives that are still unmeasured.
In this work, we present the photoelectron momentum distributions (PMDs) of helium, neon, and argon atoms driven by a linearly polarized, visible (527-nm) or near-infrared (800-nm) laser pulse (20 optical cycles in duration) based on the time-dependent density-functional theory (TDDFT) under the local-density approximation with a self-interaction correction. A set of time-dependent Kohn-Sham equations for all electrons in an atom is numerically solved using the generalized pseudospectral method. An effect of the electron-electron interaction driven by a visible laser field is not recognizable in the helium and neon PMDs except for a reduction of the overall photoelectron yield, but there is a clear difference between the PMDs of an argon atom calculated with the frozen-core approximation and TDDFT, indicating an interference of its M-shell wave functions during the ionization. Furthermore, we find that the PMDs of degenerate p states are well separated in intensity when driven by a near-infrared laser field, so that the single-active-electron approximation can be adopted safely.
Highlights: • Structural stability of Na decorated single- and double- sided graphyne and its analogous. • Hydrogen storage capacity of decorated graphyne and its analogous, based on Density Functional Theory (DFT) by PBEsol, vdW-DF2-B86R and vdW-DF2 approximations. • Electrical properties like density of states and band gap after Na decoration and hydrogen adsorption. - Abstract: In this paper, we have addressed the structural stability and hydrogen storage capability of single side and double side Na-decorated γ-Graphyne, Graphyne-like Boron Nitride and CCBN-yne, using three methods of density functional theory approximation: PBEsol, vdW-DF2-B86R and vdW-DF2. We also investigated the electronic properties of these structures, decorated by Na atoms and adsorbing hydrogen, such as density of states and band structure. Our results showed that hydrogen can be adsorbed up to three molecules per Na atom. We found the optimal geometries of adsorbed hydrogen molecules on adsorption candidates. We also found that Direct band-gap at “M” point changes after decoration and adsorption in mentioned structures. Finally, we explored metallic structures for hydrogen storage which can be used in the industry of fuel cells.
Transition Edge Sensor microcalorimeters can measure X-ray and gamma-ray energies with very high energy resolution and high photon-collection efficiency. For this technology to reach its full potential in future X-ray observatories, each sensor must be able to measure hundreds or even thousands of photon energies per second. Current “optimal filtering” approaches to achieve the best possible energy resolution work only for photons that are well isolated in time, a requirement which is in direct conflict with the need for high-rate measurements. We describe a new analysis procedure to allow fitting for the pulse height of all photons even in the presence of heavy pulse pile-up. In the limit of isolated pulses, the technique reduces to standard optimal filtering with long records. We employ reasonable approximations to the noise covariance function in order to render this procedure computationally viable even for very long data records. The technique is employed to analyze X-ray emission spectra at 600 eV and 6 keV at rates up to 250 counts s{sup −1} in microcalorimeters having exponential signal decay times of approximately 1.2 ms.
In this paper, the electronic properties of TiB{sub 2} and Ti{sub 2}B are computed within the framework of Density Functional Theory (DFT) based on Linear Combination of Atomic Orbitals (LCAO) method. In the present calculations, the generalized gradient approximation (GGA) proposed by Perdew-Burke-Erenzerhof (PBE) and Becke’s scheme are considered to treat the correlation and exchange effects respectively. TiB{sub 2} crystallizes in the hexagonal AlB{sub 2} type structure which is designated as C32 with the space group P6/mmm. While the crystal structure of Ti{sub 2}B described in this paper is isomorphous to the body centered tetragonal (Al{sub 2}Cu - type) structure with space group I4/mcm.