Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “band structure methods”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Electron–phonon physics from first principles using the EPW code

EPW is an open-source software for ab initio calculations of electron–phonon interactions and related materials properties. The code combines density functional perturbation theory and maximally localized Wannier functions to efficiently compute electron–phonon coupling matrix elements, and to perform predictive calculations of temperature-dependent properties and phonon-assisted quantum processes in bulk solids and low-dimensional materials. Here, we report on significant developments in the code since 2016, namely: a transport module for the calculation of charge carrier mobility under electric and magnetic fields using the Boltzmann transport equation; a superconductivity module for calculations of phonon-mediated superconductors using the anisotropic multi-band Eliashberg theory; an optics module for calculations of phonon-assisted indirect transitions; a module for the calculation of small and large polarons without supercells; and a module for calculating band structure renormalization and temperature-dependent optical spectra using the special displacement method. For each capability, we outline the methodology and implementation and provide example calculations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Probing topological phase transitions using high-harmonic generation

We report the prediction and realization of topological insulators have sparked great interest in experimental approaches to the classification of materials. The phase transition between non-trivial and trivial topological states is important, not only for basic materials science but also for next-generation technology, such as dissipation-free electronics. It is therefore crucial to develop advanced probes that are suitable for a wide range of samples and environments. Here we demonstrate that circularly polarized laser-field-driven high-harmonic generation is distinctly sensitive to the non-trivial and trivial topological phases in the prototypical three-dimensional topological insulator bismuth selenide. The phase transition is chemically initiated by reducing the spin–orbit interaction strength through the substitution of bismuth with indium atoms. We find strikingly different high-harmonic responses of trivial and non-trivial topological surface states that manifest themselves as a conversion efficiency and elliptical dichroism that depend both on the driving laser ellipticity and the crystal orientation. The origins of the anomalous high-harmonic response are corroborated by calculations using the semiconductor optical Bloch equations with pairs of surface and bulk bands. As a purely optical approach, this method offers sensitivity to the electronic structure of the material, including its nonlinear response, and is compatible with a wide range of samples and sample environments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Effects of spin-orbit interaction and electron correlations in strontium titanate

Here, we show that the Bloch states in the conduction band of SrTiO 3 arise from the interplay between highly anisotropic hopping in subbands derived from the Ti t2g orbitals and spin-orbit coupling that mixes these orbitals. Because of the nearly flat-band characteristics for one of the principal axes, at sufficiently high doping these Bloch states become unstable with respect to electron interactions, resulting in Mott-like singlet correlations. These findings may be relevant to the anomalous electronic properties of SrTiO 3 , including its unusual superconductivity.

36 MATERIALS SCIENCE↗
NASA NTRS1991

Comparing direct and iterative equation solvers in a large structural analysis software system

Two direct Choleski equation solvers and two iterative preconditioned conjugate gradient (PCG) equation solvers used in a large structural analysis software system are described. The two direct solvers are implementations of the Choleski method for variable-band matrix storage and sparse matrix storage. The two iterative PCG solvers include the Jacobi conjugate gradient method and an incomplete Choleski conjugate gradient method. The performance of the direct and iterative solvers is compared by solving several representative structural analysis problems. Some key factors affecting the performance of the iterative solvers relative to the direct solvers are identified.

Poole, E. L.↗
NASA NTRS1997

Doping Scheme in Atomic Chain Electronics

Due to the dramatic reduction in MOS size, there appear many unwanted effects. In these small devices, the number of dopant atoms in the channel is not macroscopic and electrons may suffer significantly different scattering from device to device since the spatial distribution of dopant atoms is no longer regarded as continuous. This prohibits integration, while it is impossible to control such dopant positions within atomic scale. A fundamental solution is to create electronics with simple but atomically precise structures, which could be fabricated with recent atom manipulation technology. All the constituent atoms are placed as planned, and then the device characteristics are deviation-free, which is mandatory for integration. Atomic chain electronics belongs to this category. Foreign atom chains or arrays form devices, and they are placed on the atomically flat substrate surface. We can design the band structure and the resultant Fermi energy of these structures by manipulating the lattice constant. Using the tight-binding theory with universal parameters, it has been predicted that isolated Si chains and arrays are metallic, Mg chains are insulating, and Mg arrays have metallic and insulating phases [1]. The transport properties along a metallic chain have been studied, emphasizing the role of the contact to electrodes [2]. For electronic applications, it is essential to establish a method to dope a semiconducting chain, which is to control the Fermi energy position without altering the original band structure. If we replace some of the chain atoms with dopant atoms randomly, the electrons will see random potential along die chain and will be localized strongly in space (Anderson localization). However, if we replace periodically, although the electrons can spread over the chain, there will generally appear new bands and band gaps reflecting the new periodicity of dopant atoms. This will change the original band structure significantly. In order to overcome this dilemma, we may place a dopant atom beside the chain at every N lattice periods (N > 1). Because of the periodic arrangement of pant atoms, we can avoid the unwanted Anderson localization. Moreover, since the dopant atoms do not constitute the chain, the overlap interaction between them is minimized, and the band structure modification can be made smallest. Some tight-binding results will be discussed to demonstrate the present idea.

Toshishige, Yamada↗
NASA NTRS2010

Estimation of Signal Coherence Threshold and Concealed Spectral Lines Applied to Detection of Turbofan Engine Combustion Noise

Combustion noise from turbofan engines has become important, as the noise from sources like the fan and jet are reduced. An aligned and un-aligned coherence technique has been developed to determine a threshold level for the coherence and thereby help to separate the coherent combustion noise source from other noise sources measured with far-field microphones. This method is compared with a statistics based coherence threshold estimation method. In addition, the un-aligned coherence procedure at the same time also reveals periodicities, spectral lines, and undamped sinusoids hidden by broadband turbofan engine noise. In calculating the coherence threshold using a statistical method, one may use either the number of independent records or a larger number corresponding to the number of overlapped records used to create the average. Using data from a turbofan engine and a simulation this paper shows that applying the Fisher z-transform to the un-aligned coherence can aid in making the proper selection of samples and produce a reasonable statistics based coherence threshold. Examples are presented showing that the underlying tonal and coherent broad band structure which is buried under random broadband noise and jet noise can be determined. The method also shows the possible presence of indirect combustion noise. Copyright 2011 Acoustical Society of America. This article may be downloaded for personal use only. Any other use requires prior permission of the author and the Acoustical Society of America.

Miles, Jeffrey Hilton↗
DOE OSTI2021

Monte Carlo transport analysis to assess intensity dependent response of a carbon-doped GaN photoconductor

In this report evaluation of the photoresponse in wurtzite GaN photoconductive switches is presented based on kinetic Monte Carlo simulations. The focus is on electron transport physics and assessment of high frequency operation. The roles of GaN band structure, Pauli exclusion, and treatment of internal fields based on the fast multipole method are all comprehensively included. The implementation was validated through comparisons of velocity-field characteristics for GaN with computational results in the literature. Photocurrent widths of less than ~7 ps for the 1 μm device can be expected, which translates into a 100 GHz upper bound. Photocurrent pulse compression below the laser full width at half maxima at high applied fields are predicted based on the interplay of space-charge effects and the negative differential velocity characteristics of GaN.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Robust analytic continuation of Green's functions via projection, pole estimation, and semidefinite relaxation

Green's functions of fermions are described by matrix-valued Herglotz-Nevanlinna functions. Since analytic continuation is fundamentally an ill-posed problem, the causal space described by the matrix-valued Herglotz-Nevanlinna structure can be instrumental in improving the accuracy and in enhancing the robustness with respect to noise. We demonstrate a three-pronged procedure for robust analytic continuation called PES: (1) projection of data to the causal space; (2) estimation of pole locations; and (3) semidefinite relaxation within the causal space. We compare the performance of PES with the recently developed Nevanlinna and Carathéodory continuation methods and find that PES is more robust in the presence of noise and does not require the usage of extended precision arithmetics. We also demonstrate that a causal projection improves the performance of the Nevanlinna and Carathéodory methods. The PES method is generalized to bosonic response functions, for which the Nevanlinna and Carathéodory continuation methods have not yet been developed. Furthermore, it is particularly useful for studying spectra with sharp features, as they occur in the study of molecules and band structures in solids.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Monte Carlo Explicitly Correlated Second-Order Many-Body Green’s Function Calculations of Semiconductor Band Gaps

A systematically converging series of ab initio, post-density-functional, size-consistent, electron-correlated approximations is desired for predictive computing of felectronic band structures of insulating, semiconducting, and metallic solids. A series that meets all of these desiderata (except the applicability to metals) is ab initio many-body Green's function theory based on Gaussian-type-orbital (GTO) basis sets. Here, its leading-order approximation, the second-order Green's function (GF2) method in the diagonal and frequency-independent approximations with the aug-cc-pVDZ basis set, is applied to the fundamental band gaps of three semiconductors (diamond, silicon, and silicon carbide in the zincblende structure) using cluster models. Corrections are made to the basis-set-incompleteness errors by the explicit-correlation (F12) ansatz (GF2-F12) for the valence band edges. The crystals are modeled as surface-passivated clusters of increasing sizes, whose wave functions are expanded by up to 2709 GTO basis functions. Immense computational costs of these calculations are overcome by the highly scalable stochastic algorithm of the Monte Carlo GF2-F12 method, whose operation cost per state increases only as a cubic power of system size, which has a tiny memory footprint and easily achieves near-perfect parallel efficiency on thousands of CPUs or on hundreds of GPUs. The correlated, F12-corrected highest-occupied and lowest-unoccupied molecular-orbital energy (HOMO-LUMO) gap is 5.78 ± 0.07 eV for C 87 H 76 as compared with the experimental value of the fundamental (indirect) band gap of bulk diamond at 5.48 eV. The correlated, F12-corrected HOMO-LUMO gaps for Si 75 H 76 and Si 32 C 43 H 76 are 2.56 ± 0.15 eV and 3.50 ± 0.12 eV, respectively, which are expected to decrease further with increasing cluster sizes. As a result, the experimental fundamental (indirect) band gaps of bulk silicon and silicon carbide are 1.17 eV and 2.42 eV, respectively.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Ultrafast high-harmonic spectroscopy of solids

High-harmonic spectroscopy, an ultrafast all-optical technique initially conceptualized in atomic and molecular systems, has now emerged as a powerful platform for studying the structure and dynamics of condensed matter. Unlike that in the gas phase, solid-state high-harmonic generation relies on the fundamental response from high atomic density and periodicity, leading to interband transitions and coherent driving of electrons and holes in their respective bands. These mechanisms make high-harmonic spectroscopy particularly sensitive to the electronic band structure, topological properties and many-body correlations in condensed media. An advantage of high-harmonic spectroscopy over other spectroscopic methods is its ability to probe ultrafast phenomena, capturing femto- to attosecond dynamics of multi-band and strongly correlated electron interactions in solids. Furthermore, in this Review, we discuss the latest experimental and theoretical advances in ultrafast high-harmonic spectroscopy of solids and provide perspectives for future research in this field.

High-harmonic generation↗
DOE OSTI2024

Comparison of interband related optical transitions and excitons in ZnGeN 2 and GaN

The optical dielectric function of ZnGeN 2 is calculated from the interband transitions using the energy bands calculated in the quasiparticle self-consistent (QS)GW^ method using two different levels of approximation: the independent particle approximation (IPA) and the Bethe-Salpeter Equation (BSE) approach. The first allows us to relate peaks in ε 2 to specific bands and k-points but does not include electron-hole interaction effects. The second includes electron-hole interaction or excitonic effects. The corresponding changes in the shape of ε 2 (ω) are found to be similar to those in GaN. The screened Coulomb interaction W^ is here calculated including electron-hole interactions in the polarization function and gives a band structure already going beyond the random phase approximation. The static dielectric constants including only electronic screening, commonly referred to as ε∞ , were calculated separately by extrapolating the wave vector dependent macroscopic dielectric function in the long-wavelength limit. Below the quasiparticle gap, we find three bound excitons optically active for different polarization. The convergence of these bound excitons with respect to the density of the k-mesh used in the BSE is studied and found to require a fine mesh. It is also found that these bound excitons originate from only the lowest conduction band and the top three valence bands. To incorporate the lattice screening, we include a scaling factor (ε∞/ε 0 ) 2 , which allows us to obtain exciton binding energies of the correct order of magnitude similar to those in GaN. Here, the excitons are related to each of the three fold split valence bands and the splittings of the latter are also studied as function of strain. Finally, a relation between the anisotropic effective masses and the valence band splitting is pointed out and explained.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Spin switchable optical phenomena in Rashba band structures through intersystem crossing in momentum space in solution-processing 2D-superlattice perovskite film

Spin-switchable phenomena are a critical element for the development of spintronic and chiroptic devices. Herein we combine a 2D-superlattice perovskite (4,4-DFPD2PbI4) film with a ferromagnetic cobalt (Co) layer to form a multiferroic perovskite/Co interface, and demonstrate spin-switchable circularly polarized luminescence (CPL) between right-handed σ+ and left-handed σ− polarizations. When the ferromagnetic spins of Co at the Co/perovskite interface are altered between positive and negative magnetic field directions, the CPL from the 2D-superlattice perovskite switches from σ+ to σ− polarization. The magnetic field effects present a unique method to confirm that CPL is generated by the circular-orbital momentum of light-emitting excitons within Rashba band structures, eliminating artifacts involving structural birefringence. Our polarized neutron reflectometry measurements confirm a super long-range spin-orbit interaction occurring in the 2D-superlattice perovskite films. The temperature dependence of spin-switchable phenomenon indicates an extraordinarily long orbital polarization lifetime, reaching microseconds at room temperature and milliseconds at 5 K.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Quasiparticle self-consistent GW band structures and phase transitions of LiAlO 2 in tetrahedrally and octahedrally coordinated structures

In this work, a first-principles computational study is presented of various phases of LiAlO 2 . The relative total energies and equations of state of the α, β, and γ phases are determined after structural relaxation of each phase. The β and γ tetrahedral phases are found to be very close in energy and lattice volume with the γ phase having the lowest energy. The octahedral α phase is a high-pressure phase and the transition pressure from the γ and β phases to α is determined to be about 1 GPa. The electronic band structures of each phase at their own equilibrium volume are determined using the quasiparticle self-consistent (QS) GW method as well as using the 0.8 Σ approach in which the QS GW self-energy is reduced by a factor of 0.8 to correct for the underscreening of W in QS GW. The effective masses of the band edges and the nature of the band gaps are presented. The lowest energy γ phase is found to have a pseudodirect gap of 7.69 eV. The gap is direct at Γ but corresponds to a dipole forbidden transition. The imaginary part of the dielectric function and the absorption coefficient are calculated in the long-wavelength limit and the random phase approximation, without local field or electron-hole interaction effects for each phase, and their anisotropies are discussed. Si doping on the Al site is investigated as a possible n-type dopant in γ-LiAlO 2 using a 128-atom supercell corresponding to 3.125% Si on the Al sublattice in the generalized gradient approximation and a smaller 16-atom cell with 25% Si in the QS GW approximation. The Si is found to significantly perturb the conduction band and lower the gap but a clearly separated deep donor defect level is not found. However, the donor binding energy is still expected to be relatively deep, on the order of a few tenths eV in the hydrogenic effective mass approximation.

36 MATERIALS SCIENCE↗
NASA NTRS1986

Solving large-scale dynamic systems using band Lanczos method in Rockwell NASTRAN on CRAY X-MP

The improved cost effectiveness using better models, more accurate and faster algorithms and large scale computing offers more representative dynamic analyses. The band Lanczos eigen-solution method was implemented in Rockwell's version of 1984 COSMIC-released NASTRAN finite element structural analysis computer program to effectively solve for structural vibration modes including those of large complex systems exceeding 10,000 degrees of freedom. The Lanczos vectors were re-orthogonalized locally using the Lanczos Method and globally using the modified Gram-Schmidt method for sweeping rigid-body modes and previously generated modes and Lanczos vectors. The truncated band matrix was solved for vibration frequencies and mode shapes using Givens rotations. Numerical examples are included to demonstrate the cost effectiveness and accuracy of the method as implemented in ROCKWELL NASTRAN. The CRAY version is based on RPK's COSMIC/NASTRAN. The band Lanczos method was more reliable and accurate and converged faster than the single vector Lanczos Method. The band Lanczos method was comparable to the subspace iteration method which was a block version of the inverse power method. However, the subspace matrix tended to be fully populated in the case of subspace iteration and not as sparse as a band matrix.

Gupta, V. K.↗
DOE OSTI2020

Many-body electronic structure of LaScO 3 by real-space quantum Monte Carlo

We present real space quantum Monte Carlo (QMC) calculations of the scandate LaSc0 3 that proved to be challenging for traditional electronic structure approaches due to strong correlation effects resulting in inaccurate band gaps from DFT and GW methods when compared with existing experimental data. Besides calculating an accurate QMC band gap corrected for supercell size biases and in agreement with numerous experiments, we also predict the cohesive energy of the crystal using the standard fixed-node QMC without any empirical or non-variational parameters. We show that promotion (optical) gap and fundamental gap agree with each other illustrating a clear absence of significant excitonic effects in the ideal crystal. We obtained these results in perfect consistency in two independent tracks that employ different basis sets (plane wave vs. localized gaussians), different codes for generating orbitals (QUANTUM ESPRESSO vs. CRYSTAL), different QMC codes (QMCPACK vs. QWALK) and different high-accuracy pseudopotentials (ccECPs vs. Troullier-Martins) presenting the maturity and consistency of QMC methodology and tools for studies of strongly correlated problems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1977

Variational calculations of vibrational properties of ozone

A variational method is used to obtain vibrational-rotational properties for ozone from an experimental quartic force field. Band positions, average structures, matrix elements for calculating infrared intensities, and effective rotational constants are reported for (O-16)3 and its O-18 isotopic species. Also, the degree to which the vibrational energies and properties are converged is investigated as a function of the basis set parameters and basis set size, and of the method of obtaining the Hamiltonian matrix elements. A convenient procedure for assigning the vibrational states is developed for variational vibrational wavefunction expansions.

Carney, G. D.↗
DOE OSTI2020

Monitoring of dielectric permittivity in accelerated alkali-silica reaction concrete with microwave backscattering

Deterioration of concrete due to the alkali–silica reaction (ASR) involves a reaction between alkaline ions in the cement pore solution and non-crystalline silica found in many aggregates. The product of reaction is a porous hydrophilic ASR gel, which causes expansion and cracking of concrete structure. Currently, there is no reliable method for non-destructive evaluation of ASR. We have investigated the feasibility of using relative dielectric permittivity obtained from X-band microwave backscattering of concrete. Advantages of this method of concrete evaluation are that measurements are non-contact and one-sided, and the approach is scalable to arbitrary-size concrete structures. The exploratory study was performed using an unrestrained set of accelerated ASR concrete prism specimens developed according to the ASTM1293 standard. One set of specimens contained no entrained air, while another set was contained entrained air. The specimens were removed from the environmental chamber at different times during a year-long study to create a set of different ASR maturity specimens. Strain measurements were performed on all specimens. Strain data both air entrained and non-air entrained sets were shown to be in close agreement with Larive model of isothermal unrestrained ASR expansion with similar fitting parameters. Furthermore, relative dielectric permittivity of concrete specimens was correlated with strain data. Both sets show correlation between dielectric permittivity and strain, with stronger linear correlation observed for the air-entrained specimens. Development of such correlation would enable estimating expansion of actual concrete structures from microwave backscattering measurements, which could be used for risk stratification to guide and minimize coring.

36 MATERIALS SCIENCE↗