SEARCH · Search NASA
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.
We have carried out extensive ab initio calculations of the electronic structure of methane, and these results are used to compute vibrational energy levels. We include basis set extrapolations, core-valence correlation, relativistic effects, and Born-Oppenheimer breakdown terms in our calculations. Our ab initio predictions of the lowest lying levels are superb.
Current interest in clusters stems from their role as novel materials as well as a possible extension of cluster results to bulk systems. Experimental investigations on clusters have been carried out using laser spectroscopy, microwave spectroscopy, heavy-particle collisions, as well as electron collisions with earlier experimental work on electron attachment and ionization having been reviewed previously. Recently, Mark and coworkers studied the decay channels of cluster ions following electron impact ionization. Rauth et al. reported the formation of the superhalogen ion SF7(-) and other nonstoichiometric cluster ions in their study of electron attachment to SF6 clusters. Kresin et al. measured the absolute electron-impact depletion cross section of metal clusters Na8, Na(20), and Na(40). They found that the inelastic scattering cross section increased with cluster size and was considerably greater than the hard sphere collision cross sections. They hypothesized that electron attachment and collision-induced fragmentation were the dominant physical processes responsible for this effect. For the two smaller clusters, they also found a sharp increase in the cross section near threshold. Most theoretical studies of clusters have been devoted to their electronic structures, vibrational relaxation, and predissociation while investigations of electron scattering from clusters has been lacking. In view of this, we recently undertook an ab initio study of electron scattering from small Be clusters and BeCO. Beryllium was chosen because it is readily amenable to ab t'nitio calculations. Moreover, the electronic structure of Be clusters has been studied extensively, showing that the Be-Be bond is relatively weak in comparison with a normal chemical bond. Our investigation focuses on how the cross sections change with cluster size and geometry. The range of energy studied, 0.05 - 5.0 eV, is chosen because of the ubiquitous resonance in the low-energy scattering of Be. Hence it can be determined if it is possible, as a consequence of the weaker bonds in the Be clusters, to identify the atomic origin of the cluster resonance.
The goal of this project was to develop computational tools for the calculation of the electronic structure of molecules containing heavy atoms, and to use these tools in the study of catalytic processes, with the overall objective of gaining an understanding of the catalytic process which could be used to design more efficient catalysts. The main catalytic system of interest was the combustion of hydrogen on platinum surfaces. Under this project, a flexible Dirac-Hartree-Fock (DHF) program has been developed, a code to calculate DHF correlation energies at the second-order Moeller-Plesset perturbation (PP2) level is almost complete, and code to include correlation at the MCSCF and MCSCF/MP2 level is planned. The tools so far developed have been validated and used to calibrate some more approximate methods, and applied to investigate the importance of relativistic effects in the bonding of hydrogen to platinum.
Calculations of fine structure and transition probabilities, using the superstructure program developed at University College London, have been carried out for the 2p(k) and 2p(k-1)3l configurations of ions with six, seven, or eight electrons. The predicted line patterns are compared with existing observations for silicon, sulfur, argon, and calcium to yield some new identifications in laboratory and solar spectra and to indicate some transitions, as yet unobserved, that may be of observational importance.
Pressure-dependent product yields have been experimentally determined for the cross-radical reaction C2H5 + C2H3. These results have been extended by calculations. It is shown that the chemically activated combination adduct, 1-C4H8*, is either stabilized by bimolecular collisions or subject to a variety of unimolecular reactions including cyclizations and decompositions. Therefore the "apparent" combination/disproportionation ratio exhibits a complex pressure dependence. The experimental studies were performed at 298 K and at selected pressures between about 4 Torr (0.5 kPa) and 760 Torr (101 kPa). Ethyl and vinyl radicals were simultaneously produced by 193 nm excimer laser photolysis of C2H5COC2H3 or photolysis of C2H3Br and C2H5COC2H5. Gas chromatograph/mass spectrometry/flame ionization detection (GC/MS/FID) were used to identify and quantify the final reaction products. The major combination reactions at pressures between 500 (66.5 kPa) and 760 Torr are (1c) C2H5 + C2H3 yields 1-butene, (2c) C2H5 + C2H5 yields n-butane, and (3c) C2H3 + C2H3 yields 1,3-butadiene. The major products of the disproportionation reactions are ethane, ethylene, and acetylene. At moderate and lower pressures, secondary products, including propene, propane, isobutene, 2-butene (cis and trans), 1-pentene, 1,4-pentadiene, and 1,5-hexadiene are also observed. Two isomers of C4H6, cyclobutene and/or 1,2-butadiene, were also among the likely products. The pressure-dependent yield of the cross-combination product, 1-butene, was compared to the yield of n-butane, the combination product of reaction (2c), which was found to be independent of pressure over the range of this study. The [ 1-C4H8]/[C4H10] ratio was reduced from approx.1.2 at 760 Torr (101 kPa) to approx.0.5 at 100 Torr (13.3 kPa) and approx.0.1 at pressures lower than about 5 Torr (approx.0.7 kPa). Electronic structure and RRKM calculations were used to simulate both unimolecular and bimolecular processes. The relative importance of C-C and C-H bond ruptures, cyclization, decyclization, and complex decompositions are discussed in terms of energetics and structural properties. The pressure dependence of the product yields were computed and dominant reaction paths in this chemically activated system were determined. Both modeling and experiment suggest that the observed pressure dependence of [1-C4H8]/[C4H10] is due to decomposition of the chemically activated combination adduct 1-C4H8* in which the weaker allylic C-C bond is broken: H2C=CHCH2CH3 yields C3H5 + CH3. This reaction occurs even at moderate pressures of approx.200 Torr (26 kPa) and becomes more significant at lower pressures. The additional products detected at lower pressures are formed from secondary radical-radical reactions involving allyl, methyl, ethyl, and vinyl radicals. The modeling studies have extended the predictions of product distributions to different temperatures (200-700 K) and a wider range of pressures (10(exp -3) - 10(exp 5) Torr). These calculations indicate that the high-pressure [1-C4H8]/[C4H10] yield ratio is 1.3 +/- 0.1.
The current capabilities of relativistic electronic structure codes for calculations on molecules such as TlF and YbF, which are of interest in the study of parity non-conservation, will be presented. The issues of basis set size and correlation methods will be discussed, and planned developments of the codes which would extend the range of their capabilities will be outlined. Issues of the sensitivity of the desired properties to the form of the nuclear potential and the flexibility of the wave function near the nucleus will be raised to stimulate discussion.
Polycyclic aromatic hydrocarbons and related species have been suggested to play a key role in the astrochemical evolution of the interstellar medium, but the formation mechanism of even their simplest building block-the aromatic benzene molecule-has remained elusive for decades. Here we demonstrate in crossed molecular beam experiments combined with electronic structure and statistical calculations that benzene (C6H6) can be synthesized via the barrierless, exoergic reaction of the ethynyl radical and 1,3- butadiene, C2H + H2CCHCHCH2 --> C6H6, + H, under single collision conditions. This reaction portrays the simplest representative of a reaction class in which aromatic molecules with a benzene core can be formed from acyclic precursors via barrierless reactions of ethynyl radicals with substituted 1,3-butadlene molecules. Unique gas-grain astrochemical models imply that this low-temperature route controls the synthesis of the very first aromatic ring from acyclic precursors in cold molecular clouds, such as in the Taurus Molecular Cloud. Rapid, subsequent barrierless reactions of benzene with ethynyl radicals can lead to naphthalene-like structures thus effectively propagating the ethynyl-radical mediated formation of aromatic molecules in the interstellar medium.
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.
For several years, scientists at NASA Ames have been studying the properties of carbon nanotubes using various experimental and computational methods. In this talk, I will compare different strategies for using quantum chemistry calculations to describe the electronic structure, deformation and chemical functionalization of single wall carbon nanotubes (SWNT) and the physisorption of small molecules on nanotube surfaces. The SWNT can be treated as an infinite (periodic) or finite length carbon cylinder or as a polycyclic aromatic hydrocarbon (PAH) molecule with an imposed curvature maintained by external constraints (as if it were cut out of the SWNT surface). Calculations are carried out using DFT and MP2 methods and a variety of atomic orbital basis sets from minimal (STO-3G) to valence triple zeta. The optimal approach is based on the particular SWNT property of interest. Examples to be discussed include: nanotube fluorination and other functionalization reactions; coating of nanotubes by water vapor and low-molecular weight organic molecules; and the nature of the interface between SWNT and liquids such as water and amines. In many cases, the quantum chemistry calculations are used to parameterize or validate force fields for molecular dynamics simulations. The results of these calculations have helped explain experimental data and contributed to the design of novel materials and sensors based on carbon nanotubes. Some of this research is described in the following papers:
The density of the corona over the active regions generating the type III exciters are examined, taking into account radio data together with daily white light coronal observations obtained with a K-coronameter. The intensity profiles and the synoptic maps derived from the K-corona data make it possible to distinguish between dense structures and low intensity regions. A deconvolution technique developed by Leblanc et al. (1970) is applied to calculate the electron density of these structures. For the low intensity regions van de Hulst's (1950) method has been used to compute the electron density of models.
Stable midlatitude red arc observation by Alouette 2 satellite for electron temperature and density structure measurement, calculating intensity and extent
Automatic calculations of strengths of force transmitting components of thin-walled aircraft structures
To provide theoretical insights into the stability and dynamics of the new rare gas compounds HArF and HKrF, reaction paths for decomposition processes HRgF to Rg + HF and HRgF to H + Rg + F (Rg = Ar, Kr) are calculated using ab initio electronic structure methods. The bending channels, HRgF to Rg + HF, are described by single-configurational MP2 and CCSD(T) electronic structure methods, while the linear decomposition paths, HRgF to H + Rg + F, require the use of multi-configurational wave functions that include dynamic correlation and are size extensive. HArF and HKrF molecules are found to be energetically stable with respect to atomic dissociation products (H + Rg + F) and separated by substantial energy barriers from Rg + HF products, which ensure their kinetic stability. The results are compatible with experimental data on these systems.
Satellite observations of electron temperature and the topside electron density structure are used to calculate the red arc according to the thermal conduction model. In this model, energy from the magnetosphere flows in the electron gas along geomagnetic field lines into the ionosphere. This energy heats the ambient F-region electrons sufficiently to excite the oxygen atoms to the O I(1D) level by collisional impact giving rise to the 6300 A emission characteristic of the arc. For the 31 October/1 November 1968 red arc, the calculated emission rate, geographical position, and horizontal extent of the red arc are in agreement with the photometric data obtained by airglow observatories. Almost all of the energy conducted into the red arc is ultimately transferred to the neutral gas through elastic and inelastic collisions. This energy drives a large thermal cell and the circulation extends the influence of the arc thousands of kilometers beyond the region of direct heat input. The calculated neutral gas temperature response to electron heating within the arc is small.
Anharmonic vibrational frequencies and intensities are computed for hydrogen fluoride clusters (HF)n with n=3,4 and mixed clusters of hydrogen fluoride with water (HF)n(H2O)n where n=1,2. For the (HF)4(H2O)4 complex, the vibrational spectra are calculated at the harmonic level, and anharmonic effects are estimated. Potential energy surfaces for these systems are obtained at the MP2/TZP level of electronic structure theory. Vibrational states are calculated from the potential surface points using the correlation-corrected vibrational self-consistent field (CC-VSCF) method. The method accounts for the anharmonicities and couplings between all vibrational modes and provides fairly accurate anharmonic vibrational spectra that can be directly compared with experimental results without a need for empirical scaling. For (HF)n, good agreement is found with experimental data. This agreement shows that the MP2 potential surfaces for these systems are reasonably reliable. The accuracy is best for the stiff intramolecular modes, which indicates the validity of MP2 in describing coupling between intramolecular and intermolecular degrees of freedom. For (HF)n(H2O)n experimental results are unavailable. The computed intramolecular frequencies show a strong dependence on cluster size. Intensity features are predicted for future experiments.
Anharmonic vibrational frequencies and intensities are computed for hydrogen fluoride clusters (HF)n, with n = 3, 4 and mixed clusters of hydrogen fluoride with water (HF)n(H2O)n where n = 1, 2. For the (HF)4(H2O)4 complex, the vibrational spectra are calculated at the harmonic level, and anharmonic effects are estimated. Potential energy surfaces for these systems are obtained at the MP2/TZP level of electronic structure theory. Vibrational states are calculated from the potential surface points using the correlation-corrected vibrational self-consistent field method. The method accounts for the anharmonicities and couplings between all vibrational modes and provides fairly accurate anharmonic vibrational spectra that can be directly compared with experimental results without a need for empirical scaling. For (HF)n, good agreement is found with experimental data. This agreement shows that the Moller-Plesset (MP2) potential surfaces for these systems are reasonably reliable. The accuracy is best for the stiff intramolecular modes, which indicates the validity of MP2 in describing coupling between intramolecular and intermolecular degrees of freedom. For (HF)n(H2O)n experimental results are unavailable. The computed intramolecular frequencies show a strong dependence on cluster size. Intensity features are predicted for future experiments.
DNA is damaged both by the direct and indirect effects of radiation. In the direct effect, the DNA itself is ionized, whereas the indirect effect involves the radiolysis of the water molecules surrounding the DNA and the subsequent reaction of the DNA with radical products. While this problem has been studied for many years, many unknowns still exist. To study this problem, we have developed the computer code RITRACKS [1], which simulates the radiation track structure for heavy ions and electrons, calculating all energy deposition events and the coordinates of all species produced by the water radiolysis. In this work, we plan to simulate DNA damage by using the crystal structure of a nucleosome and calculations performed by RITRACKS. The energy deposition events are used to calculate the dose deposited in nanovolumes [2] and therefore can be used to simulate the direct effect of the radiation. Using the positions of the radiolytic species with a radiation chemistry code [3] it will be possible to simulate DNA damage by indirect effect. The simulation results can be compared with results from previous calculations such as the frequencies of simple and complex strand breaks [4] and with newer experimental data using surrogate markers of DNA double ]strand breaks such as . ]H2AX foci [5].
The role of Quantum Theoretical Methods is both predictive and supportive of experimental results in Chemistry. Present day methods are able to calculate vibrational spectra and stereochemical interactions for molecules of moderate size (up to 20 atoms). As for the predictive side, the electronic structure of molecules and polymers can be calculated in order to narrow down the field of many potential candidates, which would have the novel properties looked for. The following has been accomplished at the Rutgers Camden Chemistry Department as results of calculations on molecular and polymeric systems of interest to the Polymers Branch of the NASA Lewis Research Center, under Grant NAG3-956.