Search NASA⌕ Search

SEARCH · Search NASA

Results for “density functional calculation”

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 91 records · Page 5

Stratospheric eddy diffusion coefficients from tracer data

Global distributions of nitrous oxide, methane, ozone, and carbon 14 are used to estimate four sets of stratospheric eddy diffusion coefficients. A photochemical equilibrium model calculates O(3P), O(1D), H, HO2, OH, H2O2, NO, and NO2 densities, as a function of altitude, latitude, and time. The calculated O(1D), OH, and observed Cl densities are used to obtain the eddy profiles associated with the methane and nitrous oxide distributions, for altitudes between 10 and 40 km. Application of a constant flux condition to the seasonally averaged ozone data yields eddy values below 20 km. Time-dependent carbon 14 calculations produce eddy coefficients between 13 and 27 km. A composite profile is obtained by comparing the four sets of coefficients. Further, carbon 14 computations are used to test these profiles as well as those recommended in reports issued by the National Academy of Sciences in 1976 and 1979. The composite eddy profile produces the best agreement.

Massie, S. T.↗

Radiation pressure-driven shocks in winds from hot stars

The dynamical theory of radiation-pressure-driven shocks in winds from early-type stars is outlined. After the principal parameters needed for a description of these shocks are isolated (speeds relative to incoming material and their column densities), the shocks' principal properties as functions of these parameters are estimated. The special case of time-steady shocks is examined, and calculations are given for density, ionization state, temperature, and flow velocity as functions of position; also calculated are X-ray emissivity, UV line emissivity, and total radiation force as functions of shock velocity and postshock column density.

Krolik, J. H.↗

Electromagnetic processes in the atmosphere of pulsars

The work consists of two parts. The first deals with the fine structure of radio pulses. Based on kinetic theory, processes occurring in the plasma shell of a pulsar when external electromagnetic radiation is present are investigated. It is shown that electromagnetic waves cause electrons to drift relative to ions, and initiate longitudinal oscillations. A dispersion equation describing the longitudinal oscillations in magnetized plasma is derived. Conditions for excitation of oscillations are found. Correlation functions of electron density are calculated, along with the coefficients of electromagnetic wave scattering. It is shown that variations in the amplitude of pulsar pulses are associated with scintillations caused by fluctuations in the plasma electron density. The second part of the study presents a mechanism for the radio emission of pulsars. The model of a rotating and a pulsating star, a neutron star with dipolar or more complex magnetic field, is examined.

Yukhimuk, A. K.↗

Infrared forbidden lines in H II regions and planetary nebulae

Emissivity coefficients are calculated as functions of electron density and temperature for the IR forbidden lines between 2 and 300 microns of the elements and ionization stages most abundant in H II regions and planetary nebulae. The effect of self-absorption is investigated, and a method is presented for estimating when self-absorption may be important. The IR line spectrum of the Orion Nebula is predicted using the emissivity coefficients, and ionic abundances are derived for the Ar III, S IV, and Ne II forbidden lines in a number of planetary nebulae. The sulfur abundance and sulfur/oxygen ratio are discussed for H II regions and planetary nebulae. The average sulfur/oxygen ratio is found to be about 0.05 for gaseous nebulae.

Simpson, J. P.↗

Estimation of minority carrier diffusion lengths in InP/GaAs solar cells

Minority carrier diffusion length is one of the most important parameters affecting the solar cell performance. An attempt is made to estimate the minority carrier diffusion lengths is the emitter and base of InP/GaAs heteroepitaxial solar cells. The PC-1D computer model was used to simulate the experimental cell results measured at NASA Lewis under AMO (air mass zero) spectrum at 25 C. A 16 nm hole diffusion length in the emitter and a 0.42 micron electron diffusion length in the base gave very good agreement with the I-V curve. The effect of varying minority carrier diffusion lengths on cell short current, open circuit voltage, and efficiency was studied. It is also observed that the front surface recombination velocity has very little influence on the cell performance. The poor output of heteroepitaxial cells is caused primarily by the large number of dislocations generated at the interfaces that propagate through the bulk indium phosphide layers. Cell efficiency as a function of dislocation density was calculated and the effect of improved emitter bulk properties on cell efficiency is presented. It is found that cells with over 16 percent efficiencies should be possible, provided the dislocation density is below 10(exp 6)/sq cm.

Jain, R. K.↗

What is the Structure of FeC5H6(+)?

Calculations performed using both density functional theory (DFT) and the modified coupled-pair functional (MCPF) approach show that FeC5H6(+) is more stable than HFeC5H5(+) by about 10 kcal/mol. The ground state of FeC5H6(+) is a quartet state derived from the 3d(sup 7) occupation of Fe+. For HFeC5H6(+), the MCPF 5 approach yields a sextet ground state while DFT yields a quartet; however, these two states are close in energy at both levels of theory.

Bauschlicher, Charles W., Jr.↗

Theoretical Characterization of the NNO...O3 and NNO...O Complexes

The bonding mechanisms and bond strengths of complexes between two molecular and/or atomic species will be discussed. The implications of complex formation for stratospheric chemistry will be presented. Results of theoretical calculations on NNO ... O3 and NNO ...O will be used as specific examples to illustrate the principles of complex formation. The equilibrium structures and binding energies of these complexes will be presented as well as shifts in other molecular properties due to complexation, such as vibrational frequencies. Thermal dissociation rates will be presented. The theoretical methodology used in this study consists of high-level coupled-cluster calculations and more approximate density functional theory approaches, together with large one-particle basis sets.

Lee, Timothy J.↗

Infrared Spectroscopy of Matrix-Isolated Polycyclic Aromatic Compounds and their Ions: Polycyclic Aromatic Nitrogen Heterocycles - 6

The matrix-isolation technique has been employed to measure the mid-infrared spectra of several polycyclic aromatic nitrogen heterocycles in both neutral and cationic forms. The species studied include: 7,8 benzoquinoline (C13H9N); 2-azapyrene (C15H9N); 1- and 2-azabenz(a)anthracene (C17H11N); and 1-, 2-, and 4-azachrysene (also C17H11N). The experimentally measured band frequencies and intensities for each molecule are tabulated and compared with their theoretically calculated values computed using density functional theory at the B3LYP/4-31G level. The overall agreement between experiment and theory is quite good, in keeping with previous investigations involving the parent aromatic hydrocarbons. Several interesting spectroscopic trends are found to accompany nitrogen substitution into the aromatic framework of these compounds. First, for the neutral species, the nitrogen atom produces a significant increase in the total integrated infrared intensity across the 1600 - 1100/cm region and plays an essential role in the molecular vibration that underlies an uncharacteristically intense, discrete feature that is observed near 1400/cm in the spectra of 7,8 benzoquinoline, 1-azabenz(a)anthracene, and 4-azachrysene. The origin of this enhanced infrared activity and the nature of the new 1400/cm vibrational mode are explored. Finally, in contrast to the parent hydrocarbon species, these aromatic nitrogen heterocycles possess a significant permanent dipole moment. Consequently, these dipole moments and the rotational constants are reported for these species in their neutral and ionized forms.

Mattioda, A. L.↗

Collaborative Simulation Grid: Multiscale Quantum-Mechanical/Classical Atomistic Simulations on Distributed PC Clusters in the US and Japan

A multidisciplinary, collaborative simulation has been performed on a Grid of geographically distributed PC clusters. The multiscale simulation approach seamlessly combines i) atomistic simulation backed on the molecular dynamics (MD) method and ii) quantum mechanical (QM) calculation based on the density functional theory (DFT), so that accurate but less scalable computations are performed only where they are needed. The multiscale MD/QM simulation code has been Grid-enabled using i) a modular, additive hybridization scheme, ii) multiple QM clustering, and iii) computation/communication overlapping. The Gridified MD/QM simulation code has been used to study environmental effects of water molecules on fracture in silicon. A preliminary run of the code has achieved a parallel efficiency of 94% on 25 PCs distributed over 3 PC clusters in the US and Japan, and a larger test involving 154 processors on 5 distributed PC clusters is in progress.

Kikuchi, Hideaki↗

Scalability of a Low-Cost Multi-Teraflop Linux Cluster for High-End Classical Atomistic and Quantum Mechanical Simulations

Scalability of a low-cost, Intel Xeon-based, multi-Teraflop Linux cluster is tested for two high-end scientific applications: Classical atomistic simulation based on the molecular dynamics method and quantum mechanical calculation based on the density functional theory. These scalable parallel applications use space-time multiresolution algorithms and feature computational-space decomposition, wavelet-based adaptive load balancing, and spacefilling-curve-based data compression for scalable I/O. Comparative performance tests are performed on a 1,024-processor Linux cluster and a conventional higher-end parallel supercomputer, 1,184-processor IBM SP4. The results show that the performance of the Linux cluster is comparable to that of the SP4. We also study various effects, such as the sharing of memory and L2 cache among processors, on the performance.

Kikuchi, Hideaki↗

Negative Ion Chemistry in the Coma of Comet 1P/Halley

Negative ions (anions) were identified in the coma of comet 1P/Halley from in-situ measurements performed by the Giotto spacecraft in 1986. These anions were detected with masses in the range 7-110 amu, but with insufficient mass resolution to permit unambiguous identification. We present details of a new chemical-hydrodynamic model for the coma of comet Halley that includes - for the first time - atomic and molecular anions, in addition to a comprehensive hydrocarbon chemistry. Anion number densities arc calculated as a function of radius in the coma, and compared with the Giotto results. Important anion production mechanisms arc found to include radiative electron attachment, polar photodissociation, dissociative electron attachment, and proton transfer. The polyyne anions C4H(-) and C6H(-) arc found to be likely candidates to explain the Giotto anion mass spectrum in the range 49-73 amu. Thc CN(-) anion probably makes a significant contribution to the mass spectrum at 26 amu. Larger carbon-chain anions such as C8H(1) can explain the peak near 100 amu provided there is a source of large carbon-chain-bearing molecules from the cometary nucleus.

Cordiner, M. A.↗

Application of Ab Initio Methods in the Development of Advanced Technical Ceramics

A significant need exists to develop materials not only capable of providing desired electronic and mechanical properties but also survival at extreme temperatures during service and device fabrication, such as co-firing. Computational methods offers an efficient and systematic manner to design new materials and guide their development. As an example computational -based material approaches can be used to determine the suitability of a given materials as a practical thermoelectric for energy harvesting. In this presentation examples are given of applications to relevant technical ceramics such as thermoelectrics, dielectrics, and magnetic systems. The calculation were based on density functional theory and carried out with norm conserving and projector augmented wave (PAW) methods using commercial codes Materials Studio ( Biovia, Inc) with the Cambridge Serial Total Energy Package (CASTEP) and MedeA (Materials Design Inc.) utilizing the Vienna Ab-initio Simulation Package (VASP) as the respective computational engines. This study makes predictions of relevant technical properties of ceramic materials.

thermoelectric↗

Theoretical Infrared Spectra for Polycyclic Aromatic Hydrocarbon Neutrals, Cations and Anions

Calculations are carried out using density functional theory (DFT) to determine the harmonic frequencies and intensities of the neutrals and cations of thirteen polycyclic aromatic hydrocarbons (PAHs) up to the size of ovalene. Calculations are also carried out for a few PAH anions. The DFT harmonic frequencies, when uniformly scaled by the factor of 0.958 to account primarily for anharmonicity, agree with the matrix isolation fundamentals to within an average error of about 10 per centimeter. Electron correlation is found to significantly reduce the intensities of many of the cation harmonics, bringing them into much better agreement with the available experimental data. While the theoretical infrared spectra agree well with the experimental data for the neutral systems and for many of the cations, there are notable discrepancies with the experimental matrix isolation data for some PAH cations that are difficult to explain in terms of limitations in the calculations. In agreement with previous theoretical work, the present calculations show that the relative intensities for the astronomical unidentified infrared (UIR) bands agree reasonably well with those for a distribution of polycyclic aromatic hydrocarbon (PAH) cations, but not with a distribution of PAH neutrals. We also observe that the infrared spectra of highly symmetrical cations such as coronene agree much better with astronomical observations than do those of, for example, the polyacenes such as tetracene and pentacene. The total integrated intensities for the neutral species are found to increase linearly with size, while the total integrated intensities are much larger for the cations and scale more nearly quadratically with size. We conclude that emission from moderate-sized highly symmetric PAH cations such as coronene and larger could account for the UIR bands.

Langhoff, Stephen R.↗

Simulations of Spray Reacting Flows in a Single Element LDI Injector With and Without Invoking an Eulerian Scalar PDF Method

This paper presents the numerical simulations of the Jet-A spray reacting flow in a single element lean direct injection (LDI) injector by using the National Combustion Code (NCC) with and without invoking the Eulerian scalar probability density function (PDF) method. The flow field is calculated by using the Reynolds averaged Navier-Stokes equations (RANS and URANS) with nonlinear turbulence models, and when the scalar PDF method is invoked, the energy and compositions or species mass fractions are calculated by solving the equation of an ensemble averaged density-weighted fine-grained probability density function that is referred to here as the averaged probability density function (APDF). A nonlinear model for closing the convection term of the scalar APDF equation is used in the presented simulations and will be briefly described. Detailed comparisons between the results and available experimental data are carried out. Some positive findings of invoking the Eulerian scalar PDF method in both improving the simulation quality and reducing the computing cost are observed.

Shih, Tsan-Hsing↗

Short wavelength laser calculations for electron pumping in neon-like krypton (Kr XXVII)

Calculations of electron impact collision strengths and spontaneous radiative decay rates are made for neon-like krypton (Kr XXVII) for the 2s2 2p6, 2s2 2p5 3s, 2s2 2p5 3p, and 2s2 2p5 3d configurations. From these atomic data, the level populations as a function of the electron density are calculated at two temperatures, 1 x 10 to the 7th K and 3 x 10 to the 7th K. An analysis of level populations reveals that a volume of krypton in which a significant number of the ions are in the Kr XXVII degree of ionization can produce a significant gain in transition between the 2s2 2p5 3s and 2s2 2p5 3p configurations. At an electron density of 1 x 10 to the 19th/cu cm the plasma length has to be of the order of 1 m; at a density of 1 x 10 to the 21st/cu cm the length is reduced to approximately 0.5 cm; and at an electron density of 1 x 10 to the 22nd/cu cm the length of the plasma is further reduced to approximately 1 mm.

Feldman, U.↗

Theoretical study of the freezing of polystyrene sphere suspensions

A theoretical study of melting in aqueous suspensions of polystyrene spheres is presented. The Helmholtz free energies of the liquid and solid phase of the suspensions are calculated as functions of sphere number density and electrolyte concentrations. The results tend to show that the freezing curve of such suspensions is that of a conventional classical liquid with repulsive short-range interactions. The interactions can be treated satisfactorily within a slightly modified Debye-Hueckel approximation. The modifications include size correction and a correction for nonlinear screening. The results are confirmed by analogous calculations for the solid phase, and for the line along which liquid and solid free energies are equal.

Shih, W.-H.↗

Test particle calculations of pick-up ions in the vicinity of Comet Giacobini-Zinner

Ionization of cometary neutral molecules produces ions which are picked-up by the solar wind. The cometary ion pick-up process for Comet Giacobini-Zinner is studied using a test particle method in which trajectories are numerically calculated for several thousand ions whose initial locations were chosen randomly with probability proportional to the neutral density. The cometary ion distribution function was calculated at several locations upstream of the bow shock. Several types of magnetic fluctuations were included. Circularly and linearly polarized monochromatic waves result in less pitch-angle scattering of pick-up ions than turbulent fluctuations which have about the same amplitude. Some anisotropy remains in the distribution function just upstream of the shock, even for the turbulent fluctuations.

Cravens, T. E.↗