Search

High levels of theory have been used to compute quartic force fields (QFFs) for the cyclic and linear forms of the C H + molecular cation, referred to as c-C H + and I-C H +. Specifically the 33 3333 singles and doubles coupled-cluster method that includes a perturbational estimate of connected triple excitations, CCSD(T), has been used in conjunction with extrapolation to the one-particle basis set limit and corrections for scalar relativity and core correlation have been included. The QFFs have been used to compute highly accurate fundamental vibrational frequencies and other spectroscopic constants using both vibrational 2nd-order perturbation theory and variational methods to solve the nuclear Schroedinger equation. Agreement between our best computed fundamental vibrational frequencies and recent infrared photodissociation experiments is reasonable for most bands, but there are a few exceptions. Possible sources for the discrepancies are discussed. We determine the energy difference between the cyclic and linear forms of C H +, 33 obtaining 27.9 kcal/mol at 0 K, which should be the most reliable available. It is expected that the fundamental vibrational frequencies and spectroscopic constants presented here for c-C H + 33 and I-C H + are the most reliable available for the free gas-phase species and it is hoped that 33 these will be useful in the assignment of future high-resolution laboratory experiments or astronomical observations.

A feasible initiation step for particle growth in the interstellar medium (ISM) is simulated by means of ab quantum chemistry methods. The systems studied are dimer ions formed by pairing nitrogen containing small molecules known to exist in the ISM with ions of unsaturated hydrocarbons or vice versa. Complexation energies, structures of ensuing complexes and electronic excitation spectra of the encounter complexes are estimated using various quantum chemistry methods. Moller-Plesset perturbation theory (MP2, Z-averaged perturbation theory (ZAP2), coupled cluster singles and doubles with perturbative triples corrections (CCSD(T)), and density functional theory (DFT) methods (B3LYP) were employed along with the correlation consistent cc-pVTZ and aug-cc-pVTZ basis sets. Two types of complexes are predicted. One type of complex has electrostatic binding with moderate (7-20 kcal per mol) binding energies, that are nonetheless significantly stronger than typical van der Waals interactions between molecules of this size. The other type of complex develops strong covalent bonds between the fragments. Cyclic isomers of the nitrogen containing complexes are produced very easily by ion-molecule reactions. Some of these complexes show intense ultraviolet visible spectra for electronic transitions with large oscillator strengths at the B3LYP, omegaB97, and equations of motion coupled cluster (EOM-CCSD) levels. The open shell nitrogen containing carbonaceous complexes especially exhibit a large oscillator strength electronic transition in the visible region of the electromagnetic spectrum.

We have investigated the physical characteristics of greenhouse gases (GHGs) to assess which properties are most important in determining the efficiency of a GHG. Chlorofluorcarbons (CFCs), hydrofluorocarbons (HFCs), perfluorocarbons (PFCs), nitrogen fluorides, and various other known atmospheric trace molecules have been included in this study. Compounds containing the halogens F or Cl have in common very polar X-F or X-Cl bonds, particularly the X-F bonds. It is shown that as more F atoms bond to the same central atom, the bond dipoles become larger as a result of the central atom becoming more positive. This leads to a linear increase in the total or integrated XF bond dipole derivatives for the molecule, which leads to a non-linear (quadratic) increase in infrared (IR) intensity. Moreover, virtually all of the X-F bond stretches occur in the atmospheric IR window as opposed to X-H stretches, which do not occur in the atmospheric window. It is concluded that molecules possessing several F atoms will always have a large radiative forcing parameter in the calculation of their global warming potential. Some of the implications for global warming and climate change are discussed.

This report will discuss the use of advanced simulation techniques to optimize the performance of the proposed Orion Crew Module airbag landing system design. The Boeing Company and the National Aeronautic and Space Administration s Langley Research Center collaborated in the analysis of the proposed airbag landing system for the next generation space shuttle replacement, the Orion spacecraft. Using LS-DYNA to simulate the Crew Module landing impacts, two main objectives were established and achieved: the investigation of potential methods of optimizing the airbag performance in order to reduce rebound on the anti-bottoming bags, lower overall landing loads, and increase overall Crew Module stability; and the determination of the Crew Module stability and load boundaries using the optimized airbag design, based on the potential Crew Module landing pitch angles and ground slopes in both the center of gravity forward and aft configurations. This paper describes the optimization and stability and load boundary studies and presents a summary of the results obtained and key lessons learned from this analysis.

A paper [Chem. Phys. Lett. 345, 295 (2001)] describes theoretical studies of excited electronic states of nitrogen molecules, with a view toward utilizing those states in synthesizing tetrahedral N4, or Td N4 a metastable substance under consideration as a high-energy-density rocket fuel. Several ab initio theoretical approaches were followed in these studies, including complete active space self-consistent field (CASSCF), state-averaged CASSCF (SA-CASSCF), singles configuration interaction (CIS), CIS with second-order and third-order correlation corrections [CIS(D) and CIS(3)], and linear response singles and doubles coupled-cluster (LRCCSD). Standard double zeta polarized and triple zeta double polarized one-particle basis sets were used. The CASSCF calculations overestimated the excitation energies, while SACASSCF calculations partly corrected these overestimates. The accuracy of the CIS calculations varied, depending on the particular state, while the CIS(D), CIS(3), and LRCCSD results were in generally good agreement. The energies of the lowest six excited singlet states of Td N4 as calculated by the LRCCSD were compared with the energies of possible excited states of N2 + N2 fragments, leading to the conclusion that the most likely route for synthesis of Td N4 would involve a combination of two bound quintet states of N2.

The singles and doubles coupled-cluster method that includes a perturbational estimate of connected triple excitations, denoted CCSD(T), has been used, in conjunction with approximate integral techniques, to compute highly accurate rovibrational spectroscopic constants of cyclopropenylidene, C3H2. The approximate integral technique was proposed in 1994 by Rendell and Lee in order to avoid disk storage and input/output bottlenecks, and today it will also significantly aid in the development of algorithms for distributed memory, massively parallel computer architectures. It is shown in this study that use of approximate integrals does not impact the accuracy of CCSD(T) calculations. In addition, the most accurate spectroscopic data yet for C3H2 is presented based on a CCSD(T)/cc-pVQZ quartic force field that is modified to include the effects of core-valence electron correlation. Cyclopropenylidene is of great astronomical and astrobiological interest because it is the smallest aromatic ringed compound to be positively identified in the interstellar medium, and is thus involved in the prebiotic processing of carbon and hydrogen. The singles and doubles coupled-cluster method that includes a perturbational estimate of

Accurate quartic force fields have been determined for the CCH- and NH2- molecular anions using the singles and doubles coupled-cluster method that includes a perturbational estimate of the effects of connected triple excitations, CCSD(T). Very large one-particle basis sets have been used including diffuse functions and up through g-type functions. Correlation of the nitrogen and carbon core electrons has been included, as well as other "small" effects, such as the diagonal Born-Oppenheimer correction, and basis set extrapolation, and corrections for higher-order correlation effects and scalar relativistic effects. Fundamental vibrational frequencies have been computed using standard second-order perturbation theory as well as variational methods. Comparison with the available experimental data is presented and discussed. The implications of our research for the astronomical observation of molecular anions will be discussed.

Previous work has shown that pentacoordinated bromine compounds have their lowest excited electronic states shifted to the blue relative to monocoordinated bromine molecules, and that this shift may be large enough to render them photostable in the lower stratosphere. Our earlier work has also shown that certain pentacoordinated bromine compounds are thermodynamically stable relative to their mono- or tricoordinated isomers, suggesting that if a bromine stratospheric reservoir species exists, then it is most likely a pentacoordinated compound. In this study we have examined the singlet excited electronic states of several bromine compounds in order to assess their photostability excited states in mono-, tri-, and pentacoordinated bromine molecules. Due to the strong spin-orbit mixing in bromine, we have also examined the lowest triplet excited state.

The anharmonic force field of BeH2 has been calculated near the basis set and n-particle space limits. The computed antisymmetric stretch frequencies of BeH2 and BeD2 are in excellent agreement with recent high-resolution gas-phase measurements. The agreement between theory and experiment for the other spectroscopic constants is also excellent, except for omega(sub 3) and X(sub 33) for BeH2 and G(sub 22) for BeD2. It is concluded that further experimental work is needed in order to resolve these discrepancies.

Recent observations from the Galileo satellite indicate that three of the Jupiter moons, Europa, Ganymede, and Callisto, may have subsurface oceans. Possible existence of such ocean and the nature of its composition are of great interest to astrobiologists. Data from Galileo's NIMS spectrometer indicate the possibility of hydrated salts on Europa's surface. To aid in the design of future missions, we investigated infrared spectra of MgSO4-nH20, n=1-3 using ab initio calculations. Geometry, energetics, dipole moments, vibrational frequencies and infrared intensities of pure and hydrated MgSO4 salts were determined. Significant differences are found between vibrational spectra of water molecules in complexes with MgSO4 and pure water. Some of the O-H stretching frequencies in the complexes are shifted to the red by up to 1,500 - 2,000 per cm. In addition, the SO2 stretching vibrations are found at lower frequency regions than the water vibrations. The calculated bands of water and SO2 fragments can serve as markers for the existence of the salt-water complexes on the surface of Jupiter's moon.

Harmonic and anharmonic vibrational frequencies, as well as infrared and Raman intensities, are calculated for MgSO4.nH20 (n=1-3). Electronic structure theory at the second order Moller-Plesset perturbation theory (MP2) level with a triple-zeta + polarization (TZP) basis set is used to determine the geometry, properties, and vibrational spectra of pure and hydrated MgSO4 salts. The direct vibrational self-consistent field (VSCF) method and its correlation corrected (CC-VSCF) extension are used to determine anharmonic corrections to vibrational frequencies and intensities for the pure MgSO4 and its complex with one water molecule. Very significant differences are found between vibrational of water molecules in complexes with MgSO4 and pure water. Some of the O-H stretching frequencies are shifted to the red very significantly (by up to 1500-2000/cm) upon complexation with magnesium sulfate. They should be observed between 1700 and 3000/cm in a region very different from the corresponding O-H stretch frequency region of pure water (3700-3800/cm). In addition, the SO2 stretching vibrations are found at lower frequency regions than the water vibrations. They can serve as unique identifiers for the presence of sulfate salts. The predicted infrared and Raman spectra should be of valuable help in the design of future missions and analysis of observed data from the ice surface of Jupiter's moon Europa that possibly contains hydrated MgSO4 salts.

It is shown that the T(sub 1) operator used in a previous study to define the open-shell D(sub 1) diagnostic is invalid, and leads to an arbitrary definition of the open-shell D(sub 1) diagnostic. A new definition is proposed that eliminates this ambiguity and approximately restores the mathematical relationship previously noted between the closed-shell D(sub 1) and T(sub 1) diagnostics. Statistical comparison of the T(sub 1) and D(sub 1) diagnostics shows a very high degree of correlation between them, although it is argued that both diagnostics used together can provide more information than either can separately.

We investigate the electronic absorption spectra of several maximally pericondensed polycyclic aromatic hydrocarbon radical cations with time dependent density functional theory calculations. We find interesting trends in the vertical excitation energies and oscillator strengths for this series containing pyrene through circumcoronene, the largest species containing more than 50 carbon atoms. We discuss the implications of these new results for the size and structure distribution of the diffuse interstellar band carriers.

We present a full experimental and theoretical study of an interesting series of polycyclic aromatic hydrocarbons, the oligorylenes. The absorption spectra of perylene, terrylene and quaterrylene in neutral, cationic and anionic charge states are obtained by matrix-isolation spectroscopy in Ne. The experimental spectra are dominated by a bright state that red shifts with growing molecular size. Excitation energies and state symmetry assignments are supported by calculations using time dependent density functional theory methods. These calculations also provide new insight into the observed trends in oscillator strength and excitation energy for the bright states: the oscillator strength per unit mass of carbon increases along the series.

The CCSD(T) method has been used to compute a highly accurate quartic force field and fundamental frequencies for all N-14 and N-15 isotopomers of the high energy density material T(sub d)N(sub 4). The computed fundamental frequencies show beyond doubt that the bands observed in a matrix isolation experiment by Radziszewski and coworkers are not due to different isotopomers of T(sub d)N(sub 4). The most sophisticated thermochemical calculations to date yield a N(sub 4) -> 2N(sub 2) heat of reaction of 182.22 +/- 0.5 kcal/mol at 0 K (180.64 +/- 0.5 at 298 K). It is hoped that the data reported herein will aid in the ultimate detection of T(sub d)N(sub 4).

Vertical electronic excitation energies for single states have been computed for the high energy density material (HEDM) Td N4 in order to assess possible synthetic routes that originate from excited electronic states of N2 molecules. Several ab initio theoretical approaches have been used, including complete active space self-consistent field (CASSCF), state averaged CASSCF (SA-CASSCF), singles configuration interaction (CIS), CIS with second-order and third-order correlation corrections [CIS(D)) and CIS(3)], and linear response singles and doubles coupled-cluster (LRCCSD), which is the highest level of theory employed. Standard double zeta polarized (DZP) and triple zeta double polarized (TZ2P) one-particle basis sets were used. The CASSCF calculations are found to overestimate the excitation energies, while the SA-CASSCF approach rectifies this error to some extent, but not completely. The accuracy of the CIS calculations varied depending on the particular state, while the CIS(D), CIS(3), and LRCCSD results are in generally good agreement. Based on the LRCCSD calculations, the lowest six excited singlet states are 9.35(l(sup)T1), 10.01(l(sup)T2), 10.04(1(sup)A2), 10.07(1(sup)E), 10.12(2(sup)T1), and 10.42(2(sup)T2) eV above the ground state, respectively. Comparison of these excited state energies with the energies of possible excited states of N2+N2 fragments, leads us to propose that the most likely synthetic route for Td N4 involving this mechanism arises from combination of two bound quintet states of N2.