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.
Cation and anion vacancy formation energies are calculated for HgTe, ZnTe, CdTe, and their dilute alloys. Harrison's tight-bonding theory is extended to a cluster embedded in an extended crystal. Only neutral vacancies have been considered, and two final states for the removed atom have been addressed: a free atom in vacuum and an atom on an ideal (111) surface.
Nonpairwise van der Waal interactions effects on vacancy formation energies in monatomic solids, using Lorentz oscillator and two body models for comparison
Isotropic elastic continuum model application to calculate energy and entropy of vacancy formation in metal crystals
The alloy variation of the band gap and the electron and hole effective masses have been calculated for HgCdTe and HgZnTe. Band-gap bowing is larger in HgZnTe than in HgCdTe because of the larger bond length mismatch of HgTe and ZnTe; electron and hole effective masses are found to be comparable for the two alloys for a given band gap. The electron mobility was calculated in both alloys with contributions from phonon, impurity, and alloy scattering. Contributions to the E1 line width due to alloy and impurity scattering in Hg(0.7)Cd(03)Te have been calculated. Results of calculations of the vacancy formation energies in HgTe, ZnTe, and CdTe are discussed.
A model for calculating the extraction energies and vacancy-formation energies in diamond-cubic and zinc-blende semiconductors is developed on the basis of Harrison's (1980, 1983) tight-binding theory. The extraction energies provide a reference from which other final states of the removed atoms can be calculated. The results of calculations show that, in a given compound, the calculated extraction energies are larger for the anion than for the cation, with the difference between the cation and the anion being larger in the II-VI than in the III-V compounds. This is in agreement with experimental results.
Structural and energetics properties of atoms near a grain boundary are of great importance from theoretical and experimental standpoints. From various experimental work it is concluded that diffusion at low temperatures at polycrystalline materials take place near grain boundary. Experimental and theoretical results also indicate changes of up to 70 percent in physical properties near a grain boundary. The Embedded Atom Method (EAM) calculations on structural properties of Au twist grain boundaries are in quite good agreement with their experimental counterparts. The EAM is believed to predict reliable values for the single vacancy formation energy as well as migration energy. However, it is not clear whether the EAM functions which are fitted to the bulk properties of a perfect crystalline solid can produce reliable results on grain boundaries. One of the objectives of this work is to construct the EAM functions for Cu and use them in conjunction with the molecular static simulation to study structures and energetics of atoms near twist grain boundaries in Cu. This provides tests of the EAM functions near a grain boundary. In particular, we determine structure, single vacancy formation energy, migration energy, single vacancy activation energy, and interlayer spacing as a function of distance from grain boundary. Our results are compared with the available experimental and theoretical results from grain boundaries and bulk.
Ytterbium disilicate is of interest as a potential environmental barrier coating for aerospace applications, notably for use in next generation jet turbine engines. In such applications, the transport of oxygen and water vapor through these coatings to the ceramic substrate is undesirable if high temperature oxidation is to be avoided. In an effort to understand the diffusion process in these materials, we have performed kinetic Monte Carlo simulations of vacancy-mediated and interstitial oxygen diffusion in Ytterbium disilicate. Oxygen vacancy and interstitial site energies, vacancy and interstitial formation energies, and migration barrier energies were computed using Density Functional Theory. We have found that, in the case of vacancy-mediated diffusion, many potential diffusion paths involve large barrier energies, but some paths have barrier energies smaller than one electron volt. However, computed vacancy formation energies suggest that the intrinsic vacancy concentration is small. In the case of interstitial diffusion, migration barrier energies are typically around one electron volt, but the interstitial defect formation energies are positive, with the result that the disilicate is unlikely to exhibit experience significant oxygen permeability except at very high temperature.
We have derived the embedding energy functional and two-body potential of the embedded atom method using decreasing exponentials for both the electron density and the two body potential. The embedding function was obtained from the equation of state given by Rose et al. (1984). Because of the form of the embedding function, the equilibrium lattice constant, cohesive energy, and bulk modulus are automatically satisfied. The two parameters phi(e) and gamma of the two-body potential were determined by fitting to shear modulus and the single vacancy formation energy. Contributions of up to the third nearest neighbors were included in the evaluation of the charge density rho and the two-body potential phi. The stability and anisotropy of each structure were estimated and compared with the available experimental data.
Density functional theory was used to screen eleven refractory materials – two pure metals, six nitrides, and three carbides–as high-temperature hydrogen permeation barriers to prevent hydrogen embrittlement. Activation energies were calculated for atomic hydrogen (H) diffusion into the first subsurface layer from the lowest energy surface of the high-temperature phase of each candidate material. The candidate barrier materials with the highest activation energies are h-BN, c-BN, HfN, and ZrN with predicted barriers of 3.25 eV, 3.23 eV, 3.14 eV, and 2.76 eV, respectively. Strain energies, Bader charges, and density of states were calculated for the diffusing H at the relaxed initial state and the transition state to provide insight into contributing factors to high energy barriers. The diffusing H atom in materials with the highest predicted barriers are protic. In addition, interstitial H atoms induce mid-gap states in the density of states of both BN polymorphs. The nitrogen retention of each nitride material at high temperatures was predicted using nitrogen vacancy formation energies with respect to gaseous nitrogen. Experimental evaluation of nitrogen retention in h-BN, ZrN, and TiN confirmed their resistance to nitrogen loss at 1773 K. However, of these nitrides, TiN is predicted to be the least stable. This work identifies multiple promising materials that are predicted to be effective hydrogen barriers at high temperatures and that are stable at temperatures above 2700 K, with BN predicted to perform best.
The initial goal was to grow and characterize HgZnTe and determine if it indeed had the advantageous properties that were predicted. Researchers grew both bulk and liquid phase epitaxial HgZnTe. It was determined that HgZnTe had the following properties: (1) microhardness at least 50 percent greater than HgCdTe of equivalent bandgap; (2) Hg annealing rates of at least 2 to 4 times longer than HgCdTe; and (3) higher Hg vacancy formation energies. This early work did not focus on one specific composition (x-value) of HgZnTe since NASA was interested in HgZnTe's potential for a variety of applications. Since the beginning of 1989, researchers have been concentrating, however, on the liquid phase growth of very long wavelength infrared (VLWIR) HgZnTe (cutoff approx. equals 17 microns at 65K) to address the requirements of the Earth Observing System (EOS). Since there are no device models to predict the advantages in reliability one can gain with increased microhardness, surface stability, etc., one must fabricate HgZnTe detectors and assess their relative bake stability (accelerated life test behavior) compared with HgCdTe devices fabricated in the same manner. Researchers chose to fabricate HIT detectors as a development vehicle for this program because high performance in the VLWIR has been demonstrated with HgCdTe HIT detectors and the HgCdTe HIT process should be applicable to HgZnTe. HIT detectors have a significant advantage for satellite applications since these devices dissipate much less power than conventional photoconductors to achieve the same responsivity.
Models of encapsulated 1/2 nuclear spin H-1 and P-31 atoms in fullerene and diamond nanocrystallite, respectively, are proposed and examined with ab-initio local density functional method for possible applications as single quantum bits (qubits) in solid-state quantum computers. A H-1 atom encapsulated in a fully deuterated fullerene, C(sub 20)D(sub 20), forms the first model system and ab-initio calculation shows that H-1 atom is stable in atomic state at the center of the fullerene with a barrier of about 1 eV to escape. A P-31 atom positioned at the center of a diamond nanocrystallite is the second model system, and 3 1P atom is found to be stable at the substitutional site relative to interstitial sites by 15 eV, Vacancy formation energy is 6 eV in diamond so that substitutional P-31 atom will be stable against diffusion during the formation mechanisms within the nanocrystallite. The coupling between the nuclear spin and weakly bound (valance) donor electron coupling in both systems is found to be suitable for single qubit applications, where as the spatial distributions of (valance) donor electron wave functions are found to be preferentially spread along certain lattice directions facilitating two or more qubit applications. The feasibility of the fabrication pathways for both model solid-state qubit systems within practical quantum computers is discussed with in the context of our proposed solid-state qubits.
We have used a Green's function technique to calculate the energy levels and formation energy of deep defects in the narrow gap semiconductors mercury cadmium telluride (MCT), mercury zinc telluride (MZT) and mercury zinc selenide (MZS). The formation energy is calculated from the difference between the total energy with an impurity cluster and the total energy for the perfect crystal. Substitutional (including antisite), interstitial (self and foreign), and vacancy deep defects are considered. Relaxation effects are calculated (with molecular dynamics). By use of a pseudopotential, we generalize the ideal vacancy model so as to be able to consider relaxation for vacancies. Different charge states are considered and the charged state energy shift (as computed by a modified Haldane-Anderson model) can be twice that due to relaxation. Different charged states for vacancies were not calculated to have much effect on the formation energy. For all cases we find deep defects in the energy gap only for cation site s-like orbitals or anion site p-like orbitals, and for the substitutional case only the latter are appreciably effected by relaxation. For most cases for MCT, MZT, MZS, we consider x (the concentration of Cd or Zn) in the range appropriate for a band gap of 0.1 eV. For defect energy levels, the absolute accuracy of our results is limited, but the precision is good, and hence chemical trends are accurately predicted. For the same reason, defect formation energies are more accurately predicted than energy level position. We attempt, in Appendix B, to calculate vacancy formation energies using relatively simple chemical bonding ideas due to Harrison. However, these results are only marginally accurate for estimating vacancy binding energies. Appendix C lists all written reports and publications produced for the grant. We include abstracts and a complete paper that summarizes our work which is not yet available.
Refractory materials such as metallic borides, often considered as ultra high temperature ceramics (UHTC), are characterized by high melting point, high hardness, and good chemical inertness. These materials have many applications which require high temperature materials that can operate with no or limited oxidation. Ab initio, first principles methods are the most accurate modeling approaches available and represent a parameter free description of the material based on the quantum mechanical equations. Using these methods, many of the intrinsic properties of these material can be obtained. We performed ab initio calculations based on density functional theory for the UHTC materials ZrB2 and HfB2. Computational results are presented for structural information (lattice constants, bond lengths, etc), electronic structure (bonding motifs, densities of states, band structure, etc), thermal quantities (phonon spectra, phonon densities of states, specific heat), as well as information about point defects such as vacancy and antisite formation energies.
Understanding surface diffusion is essential in understanding surface phenomena, such as crystal growth, thin film growth, corrosion, physisorption, and chemisorption. Because of its importance, various experimental and theoretical efforts have been directed to understand this phenomena. The Field Ion Microscope (FIM) has been the major experimental tool for studying surface diffusion. FIM have been employed by various research groups to study surface diffusion of adatoms. Because of limitations of the FIM, such studies are only limited to a few surfaces: nickel, platinum, aluminum, iridium, tungsten, and rhodium. From the theoretical standpoint, various atomistic simulations are performed to study surface diffusion. In most of these calculations the Embedded Atom Method (EAM) along with the molecular static (MS) simulation are utilized. The EAM is a semi-empirical approach for modeling the interatomic interactions. The MS simulation is a technique for minimizing the total energy of a system of particles with respect to the positions of its particles. One of the objectives of this work is to develop the EAM functions for Cu and use them in conjunction with the molecular static (MS) simulation to study diffusion of a Cu atom on a perfect as well as stepped Cu(100) surfaces. This will provide a test of the validity of the EAM functions on Cu(100) surface and near the stepped environments. In particular, we construct a terrace-ledge-kink (TLK) model and calculate the migration energies of an atom on a terrace, near a ledge site, near a kink site, and going over a descending step. We have also calculated formation energies of an atom on the bare surface, a vacancy in the surface, a stepped surface, and a stepped-kink surface. Our results are compared with the available experimental and theoretical results.
As crewed travel becomes more common to the Moon and beyond, space nuclear propulsion-based craft offer reduced travel times and double the specific impulse compared to conventional chemical rockets. The nuclear core and its integral hydrogen propellant channels will require advanced refractory carbide coatings over a carbon/carbon substrate for hydrogen resistance. Using ab initio and other computational techniques, we studied select (Zr/Ta/Nb)C mixed carbide systems for erosion resistance and hydrogen diffusion using realistic levels of carbon vacancies and experimentally obtainable stoichiometries. Various contributions to the formation and hydrogen erosion free energies were calculated, including configurational, electronic, and vibrational terms. Our calculations indicate that ZrC is less likely to lose carbon than NbC or TaC, providing an explanation for heritage NERVA data on ZrC and NbC mass loss performance. Additionally, binary mixtures show compositions that are more stable than their end-member counterparts. Hydrogen diffusion in stoichiometric ZrC was found to be slightly higher than previously theorized, with the caveat that hydrogen saturation of carbon vacancies will hinder further hydrogen diffusion through the bulk, while hydrogen diffusion in ZrNbC is slightly enhanced over ZrC. And finally, we show that as carbon is depleted, the deformation mechanism of ZrC changes from brittle to ductile. Though many questions still remain as to the bulk performance of these materials, these trends and results are important in terms of materials selection efforts for hydrogen-resistant coatings on carbon/carbon substrates.
As crewed travel becomes more common to the Moon and beyond, space nuclear propulsion-based craft offer reduced travel times and double the specific impulse compared to conventional chemical rockets. The nuclear core and its integral hydrogen propellant channels will require advanced refractory carbide coatings over a carbon/carbon substrate for hydrogen resistance. Using \textit{ab initio} and other computational techniques, we studied select (Zr/Ta/Nb)C mixed carbide systems for erosion resistance and hydrogen diffusion using realistic levels of carbon vacancies and experimentally obtainable stoichiometries. Various contributions to the formation and hydrogen erosion free energies were calculated, including configurational, electronic, and vibrational terms. Our calculations indicate that ZrC is less likely to lose carbon than NbC or TaC, providing an explanation for heritage NERVA data on ZrC and NbC mass loss performance. Additionally, binary mixtures show compositions that are more stable than their end-member counterparts. Hydrogen diffusion in stoichiometric ZrC was found to be slightly higher than previously theorized, with the caveat that hydrogen saturation of carbon vacancies will hinder further hydrogen diffusion through the bulk, while hydrogen diffusion in ZrNbC is slightly enhanced over ZrC. And finally, we show that as carbon is depleted, the deformation mechanism of ZrC changes from brittle to ductile. Though many questions still remain as to the bulk performance of these materials, these trends and results are important in terms of materials selection efforts for hydrogen-resistant coatings on carbon/carbon substrates.
As crewed travel becomes more common to the Moon and beyond, space nuclear propulsion-based craft offer reduced travel times and double the specific impulse compared to conventional chemical rockets. The nuclear core and its integral hydrogen propellant channels will require advanced refractory carbide coatings over a carbon/carbon substrate for hydrogen resistance. Using ab initio and other computational techniques, we studied select (Zr/Ta/Nb)C mixed carbide systems for erosion resistance and hydrogen diffusion using realistic levels of carbon vacancies and experimentally obtainable stoichiometries. Various contributions to the formation and hydrogen erosion free energies were calculated, including configurational, electronic, and vibrational terms. Our calculations indicate that ZrC is less likely to lose carbon than NbC or TaC, providing an explanation for heritage NERVA data on ZrC and NbC mass loss performance. Additionally, binary mixtures show compositions that are more stable than their end-member counterparts. Hydrogen diffusion in stoichiometric ZrC was found to be slightly higher than previously theorized, with the caveat that hydrogen saturation of carbon vacancies will hinder further hydrogen diffusion through the bulk, while hydrogen diffusion in ZrNbC is slightly enhanced over ZrC. And finally, we show that as carbon is depleted, the deformation mechanism of ZrC changes from brittle to ductile. Though many questions still remain as to the bulk performance of these materials, these trends and results are important in terms of materials selection efforts for hydrogen-resistant coatings on carbon/carbon substrates.
CdTe-HgTe heterostructures preparation with minimum energy gaps, showing photovoltaic effect caused by undesirable p-n junction due to Cd vacancy formation