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 atomic structure in a 110 screw dislocation core for aluminum is obtained by computer simulation. The lattice statics technique is employed since it entails no artificially imposed elastic boundary around the defect. The interatomic potential has no adjustable parameters and was derived from pseudopotential theory. The resulting atomic displacements were allowed to relax in all three dimensions.
The prolongation structure of a closed ideal of exterior differential forms is further discussed, and its use illustrated by application to an ideal (in six dimensions) representing the cubically nonlinear Schroedinger equation. The prolongation structure in this case is explicitly given, and recurrence relations derived which support the conjecture that the structure is open - i.e., does not terminate as a set of structure relations of a finite-dimensional Lie group. We introduce the use of multiple pseudopotentials to generate multiple Baecklund transformation, and derive the double Baecklund transformation. This symmetric transformation concisely expresses the (usually conjectured) theorem of permutability, which must consequently apply to all solutions irrespective of asymptotic constraints.
The atomic structure of a screw dislocation core of the 110 line type in aluminum is calculated by the modified lattice-statics method developed in the preceding paper. The method includes anharmonic as well as harmonic forces and permits relaxation of the atoms in all three dimensions. All forces used in the present calculations were derived from a first-principles interatomic pair potential obtained via pseudopotential theory. Several significant differences from the ordinary lattice statics results are noted, including the displacement field, Peierl's energy barrier, and the equilibrium core-center location.
The criterion for the existence of vortex-like ion phase-space configurations, as obtained by a standard pseudopotential method, is found to coincide with the criterion for the linear instability for two (cold) counterstreaming ion beams. A nonlinear equation is derived, which demonstrates that this instability actually evolves into such phase-space configurations. A small, but nonzero, ion temperature turns out to be essential for the saturation into stationary structures
Quasi-geostrophic dynamics, as usually formulated, centers on the conservation of pseudopotential vorticity (PV) on isobaric surfaces. The present study has the objective to investigate some simple constraints which conservation of PV imposes on internal Rossby waves. Attention is given to some of the limitations that conservation of potential vorticity imposes on the dynamics of the stratosphere. Another possible constraint imposed by hydrodynamic stability is also discussed. The quantitative limitations these constraints put on the amplitudes of internal Rossby waves are assessed and compared with observations of the lower stratosphere.
The high pressure structural phase transition in Ge has been studied using the energy dispersive X-ray diffraction technique and a synchrotron radiation source. Ge was observed to transform to the beta-Sn tetragonal structure in agreement with the earlier results of Jamieson, but the phase transition began at 80 + or - 5 kilobars, a somewhat lower value than generally reported. These experimental diffraction results are compared with the recent self-consistent pseudopotential calculations of Yin and Cohen (1981) and with the observed transition pressure for shock wave loaded Ge.
Stevens et al. (1981) considered CsH as a two electron problem, but simulated the Cs core electrons by an empirical pseudopotential following the work by Bardsley (1970). These potentials, since they are derived empirically, include the relativistic effects and the atomic core-valence interactions. However, molecular core-core interactions are not accounted for. Stevens et al. obtained an R(e) value which is too small. This result was attributed to uncorrected core (proton)-core interactions in the molecule. The present investigation is concerned with a resolution of the discrepancies between the results of earlier studies. The X 1 Sigma + ground state of CsH is recalculated. The calculation employs a nine valence electron relativistic effective core potential (RECP) for Cs. The bonding in CsH is found to involve a Cs(6s)-H(1s) bond but with a significant ionic (Cs+H-) component.
A detailed calculation of the index refraction of various GaAs-AlAs superlattices is presented for the first time. The calculation is performed by using a hybrid approach which combines the k-p method with the pseudopotential technique. Appropriate quantization conditions account for the influence of the superstructures on the electronic properties of the systems. The results of the model are in very good agreement with the experimental data. In comparison with the index of refraction of the corresponding AlGaAs alloy, characterized by the same average mole fraction of Al, the results indicate that the superlattice index of refraction values attain maxima at the various quantized transition energies. For certain structures the difference can be as large as 2 percent. These results suggest that the waveguiding and dispersion relation properties of optoelectronic devices can be tailored to design for specific optical application by an appropriate choice of the superlattice structure parameters.
A major goal of this project is to understand the surface tension and other thermophysical properties of liquid metals and alloys from a fundamental viewpoint. The approach is to calculate these quantities by a first principles technique which combines the statistical-mechanical theory of the liquid state with an electronic pseudopotential theory of electrons in metals. The inhomogeneity of the surface is treated using an ionic-density-functional formalism developed with the support of NASA. Of particular interest are the variation of surface tension with temperature and impurity concentration: such variations strongly influence the types of convection which make take place in a low-gravity environment. Some progress has already been achieved in computing the reduction of surface tension due to the presence of low-surface-tension impurities, and the corresponding surface segregation of such impurities. In the coming year, it is planned to concentrate on the surface properties of materials of particular interest to the MSA program: Si, Ga and GaSn alloys. An additional goal is to gain some theoretical understanding of the high temperature thermophysical properties of liquid metals, particularly high melting point materials which have not been studied extensively from a theoretical viewpoint.
The optical properties of GaAs-Al sub x Ga sub 1-xAs superlattices are calculated as a function of the frequency and superlattice structure. The comutations are performed using a partition method which combines the vectors k.p method with the pseudopotential technique. The influence of the super-structure on the electronic properties of the systems is accounted for by appropriate quantization conditions. The anisotropy and structure dependence of the dielectric constant result mainly from the contribution of the gamma region while the contributions of the other regions of the Brillouin zone are rather insensitive to the superlattice structure. The superlattice index of refraction values are shown to attain maxima at the various quantized transition energies, where for certain structures, the difference between the refractive indices of the superlattices and its corresponding Al sub x Ga sub 1-xAs alloy can be as large as 2%. In general results are in good agreement with the experimental data.
The optical properties of III to V binary and ternary compounds and GaAs-Al(x)Ga(1-x)As superlattices are determined by calculating the real and imaginary parts of the transverse dielectric constant. Emphasis is given to determining the influence of different material and superlattice parameters on the values of the index of refraction and absorption coefficient. In order to calculate the optical properties of a material, it is necessary to compute its electronic band structure. This was accomplished by introducing a partition band structure approach based on a combination of the vector k x vector p and nonlocal pseudopotential techniques. The advantages of this approach are that it is accurate, computationally fast, analytical, and flexible. These last two properties enable incorporation of additional effects into the model, such as disorder scattering, which occurs for alloy materials and excitons. Furthermore, the model is easily extended to more complex structures, for example multiple quantum wells and superlattices. The results for the transverse dielectric constant and absorption coefficient of bulk III to V compounds compare well with other one-electron band structure models and the calculations show that for small frequencies, the index of refraction is determined mainly by the contibution of the outer regions of the Brillouin zone.
A general model of the transverse dielectric constant of III-V compounds is developed using a hybrid method which combines the kp method with a nonlocal pseudopotential calculation. In this method the Brillouin zone is partitioned into three regions by expanding the energy bands and matrix elements about the F, X, and L symmetry points. The real and imaginary parts of the dielectric constant are calculated as a sum of the individual contributions of each region. By using this partition method, it is possible to get good insight into the dependence of the dielectric constant on the shape of the band structure.
Stationary solutions of the derivative nonlinear Schroedinger equation are discussed and classified by using a pseudopotential formulation. The solutions consist of a rich family of nonlinear Alfven waves and solitons with parallel and oblique propagation directions. Expressions for the envelope and the phase of nonlinear waves with periodic envelope modulation, and 'hyperbolic' and 'algebraic' solitons are given. The propagation angle for the slightly modulated elliptic, periodic waves and for oblique solitons is evaluated.
A composite velocity procedure for the three-dimensional reduced Navier-Stokes equations is developed. In the spirit of matched asymptotic expansions, the velocity components are written as a combination multiplicative and additive composite of viscouslike velocities and pseudopotential or inviscid velocities. The solution procedure is then consistent with both asymptotic inviscid flow and boundary layer theory. For transonic flow cases, the Enquist-Osher flux biasing scheme developed for the full potential equation is used. A quasi-conservation form of the governing equations is used in the shock region to capture the correct rotational shock with the standard nonconservation form of the equations used in nonshock regions. The consistent coupled strongly implicit procedure coupled with a plane relaxation procedure is used to solve the discretized equations.
The flow over afterbody geometries was investigated using the reduced Navier-Stokes (RNS) approximation. Both pressure velocity flux-split and composites velocity primitive variable formulations were considered. Pressure or pseudopotential relaxation procedures are combined with sparse matrix or coupled strongly implicit algorithms to form a three-dimensional solver for general non-orthogonal coordinates. Three-dimensional subsonic and transonic viscous/inviscid interacting flows were evaluated. Solutions with and without regions of recirculation were obtained.
A reduced form of the Navier-Stokes equations, defined by a single composite of the Euler, boundary layer, and triple deck approximations, is considered for the computation of viscous interacting flows. Global pressure or pseudopotential relaxation methods are integrated with coupled sparse matrix direct solvers or coupled strongly implicit ILU algorithms to efficiently capture sharp shocks and regions of recirculation. Solutions are obtained for a variety of 2D and 3D geometries and for Mach numbers (M) spanning the range from incompressible (M = 0) to supersonic (M = 6) flow.
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.
Since the original chapter, we briefly review two calculational methods used for basic applications in autoionization of few body systems: (a) complex rotation and (b) a pseudopotential method. Complex rotation, briefly mentioned above, has been extensively applied with great accuracy. Two additional basic systems to be mentioned here are H(sup -) and Ps(sup -) (Ps=positronium). In complex rotation the particle distances are multiplied by the common phase factor given.