Search NASA⌕ Search

SEARCH · Search NASA

Results for “ab initio theory”

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 73 records · Page 4

Positive Anharmonicities: The Oxonide Anion as an Example

An accurate ab initio quartic force field for the ozonide anion has been determined at the singles and doubles coupled-cluster level of theory that includes a perturbational estimate of the effects of connected triple excitations, denoted CCSD(T), using the augmented valence triple-zeta correlation consistent one-particle basis set. Convergence of the harmonic frequencies with respect to the one-particle basis set has been demonstrated by computing quadratic force fields at the CCSD(T) level using augmented valence double-zeta and quadruple-zeta basis sets. Fundamental vibrational frequencies have been determined via second-order ro-vibrational perturbation theory and also using exact variational methods. Agreement is very good, and they both show that the antisymmetric O-O stretch, v 3, possesses a positive anharmonic correction (that is, the fundamental frequency is predicted to be higher in energy than the harmonic frequency). Comparison of the O_3 and O3 quartic force fields shows that the positive anharmonic correction is the result of a particularly large F3333 symmetry internal coordinate force constant. The reasonableness of this force constant has been tested by computing both F33 and F3333 at the CCSD, CCSD(T), and CASPT2 levels of theory. A discussion of known positive anharmonicities for stretching modes is presented.

Lee, Timothy J.↗

Infrared laser spectroscopy of jet-cooled carbon clusters: the nu 5 band of linear C9

The nu 5 antisymmetric stretching vibration of 1 sigma+g C9 has been observed using direct infrared diode laser absorption spectroscopy of a pulsed supersonic cluster beam. Twenty-eight rovibrational transitions measured in the region of 2079-2081 cm-1 were assigned to this band. A combined least squares fit of these transitions with previously reported nu 6 transitions yielded the following molecular constants for the nu 5 band: nu 0 = 2 079.673 58(17) cm-1, B"= 0.014 321 4(10) cm-1, and B'=0.014 288 9(10) cm-1. The IR intensity of the nu 5 band relative to nu 6 was found to be 0.108 +/- 0.006. Theoretical predictions for the relative intensities vary widely depending upon the level of theory employed, and the experimental value reported here is in reasonable agreement only with the result obtained from the most sophisticated ab initio calculation considered (CCSD).

Non-NASA Center↗

A reinvestigation of the rate of the C/+/ + H2 radiative association reaction

It is noted that new experimental results and statistical theories have prompted a reinvestigation of the rate coefficient of the interstellar reaction in which C(+) + H2 yields CH2(+) + h(nu) in the 10-100 K temperature range. The results presented here indicate a rate coefficient between 10 to the -16 and 10 to the -15 cu cm/s at all temperatures studied. In applying the modified thermal and phase space theories, it is expected that they will be as accurate as for radiative association as for three-body association, provided the ab initio value for the radiative decay rate is correct. It is expected that the calculated values of the rate coefficient will be accurate to within an order of magnitude and will have the correct temperature dependence.

Herbst, E.↗

AB Initio Study of the Structure and Spectroscopic Properties of Halogenated Thioperoxy Radicals

Thioperoxy (XSO or XOS) radicals exist in a variety of chemical environments, and they have as a consequence drawn some interest. HSO, an important species in the chemistry of the troposphere, has been examined both experimentally. The halogenated (X = F, Cl or Br) peroxy species and isovalent thioperoxy species have been studied less, but they too are potentially interesting because oxidized sulfur species and halogen sources are present in the atmosphere. Learning the fate of XSO and XOS radicals is important to understanding the atmospheric oxidation chemistry of sulfur compounds. Of these, FSO and ClSO are particularly interesting because they have been directly detected spectroscopically. Recent studies in our laboratory on the photochemistry of thionyl halides (X2SO; where X = F or Cl) have suggested new ways to generate XSO species. The laser-induced photodissociation of thionyl fluoride, F2SO, at 193 nm and thionyl chloride, ClSO, at 248 nm is characterized by a radical mechanism, X2SO -> XSO + X. The structure of FSO has been characterized experimentally by Endo et cd. employing microwave spectroscopy. Using the unrestricted Hartree-Fock (UHF) self-consistent field (SCF) method, Sakai and Morokuma computed the electronic structure of the ground (sup 2)A" and the first excited (sup 2)A' states of FSO. Electron correlation was not taken into account in their study. In a laser photodissociation experiment, Huber et al. identified ClSO mass spectromctrically. ClSO has also been detected in low temperature matrices by EPR and in the gas phase by far IR laser magnetic resonance. Although the structure of FSO is known in detail, the only study, experimental or theoretical, of CISO has been an ab initio HFSCF study by Hinchliffe. Electron correlation corrections were also excluded from this study. In order to better understand the isomerization and dissociation dynamics of the radical species, we have performed ab initio correlated studies of the potential energy surfaces (PES) of ClSO and its isomer ClOS at the QCISD(T)/6-31 G* level of theory. For FSO and FOS, more extensive QCISD/6-31 1G(2df) calculations have been possible, and the results are summarized here.

Munoz, Luis A.↗

Formation of Carbon Nanotube Based Gears: Quantum Chemistry and Molecular Mechanics Study of the Electrophilic Addition of o-Benzyne to Fullerenes, Graphene, and Nanotubes

Considerable progress has been made in recent years in chemical functionalization of fullerene molecules. In some cases, the predominant reaction products are different from those obtained (using the same reactants) from polycyclic aromatic hydrocarbons (PAHs). One such example is the cycloaddition of o-benzyne to C60. It is well established that benzyne adds across one of the rings in naphthalene, anthracene and other PAHs forming the [2+4] cycloaddition product (benzobicyclo[2.2.2.]-octatriene with naphthalene and triptycene with anthracene). However, Hoke et al demonstrated that the only reaction path for o-benzyne with C60 leads to the [2+2] cycloaddition product in which benzyne adds across one of the interpentagonal bonds (forming a cyclobutene ring in the process). Either reaction product results in a loss of aromaticity and distortion of the PAH or fullerene substrate, and in a loss of strain in the benzyne. It is not clear, however, why different products are preferred in these cases. In the current paper, we consider the stability of benzyne-nanotube adducts and the ability of Brenner's potential energy model to describe the structure and stability of these adducts. The Brenner potential has been widely used for describing diamondoid and graphitic carbon. Recently it has also been used for molecular mechanics and molecular dynamics simulations of fullerenes and nanotubes. However, it has not been tested for the case of functionalized fullerenes (especially with highly strained geometries). We use the Brenner potential for our companion nanogear simulations and believe that it should be calibrated to insure that those simulations are physically reasonable. In the present work, Density Functional theory (DFT) calculations are used to determine the preferred geometric structures and energetics for this calibration. The DFT method is a kind of ab initio quantum chemistry method for determining the electronic structure of molecules. For a given basis set expansion, it is comparable in accuracy to the MP2 method (better than Hartree Fock, but less accurate than more extensive electron correlation methods such as MP4 or CCSD). However, for systems with large numbers of basis functions it more efficient than any other methods that include electron correlation effects. In this presentation we show the results of DFT calculations for the reaction of benzyne with naphthalene, C60, and nanotube models. We compare energies for [2+2] and [2+4] cycloaddition products. The preferred products for the naphthalene and C60 reactions have been determined by experiment and, thus, these cases serve as a validation of our quantum chemical approach. We also compare the DFT and Brenner potential results. Finally we can predict the likelihood of reaction between benzyne and nanotubes.

Jaffe, Richard↗

Equilibrium properties of transition-metal ion-argon clusters via simulated annealing

The geometrical structures of M(+) (Ar)n ions, with n = 1-14, have been studied by the minimization of a many-body potential surface with a simulated annealing procedure. The minimization method is justified for finite systems through the use of an information theory approach. It is carried out for eight potential-energy surfaces constructed with two- and three-body terms parametrized from experimental data and ab initio results. The potentials should be representative of clusters of argon atoms with first-row transition-metal monocations of varying size. The calculated geometries for M(+) = Co(+) and V(+) possess radial shells with small (ca. 4-8) first-shell coordination number. The inclusion of an ion-induced-dipole-ion-induced-dipole interaction between argon atoms raises the energy and generally lowers the symmetry of the cluster by promoting incomplete shell closure. Rotational constants as well as electric dipole and quadrupole moments are quoted for the Co(+) (Ar)n and V(+) (Ar)n predicted structures.

Asher, Robert L.↗

Transport properties in the Jovian atmosphere

Transport properties in a Jupiter-like atmosphere (89 mol % hydrogen and 11 mol % helium) are obtained by using the method of the kinetic theory of gases. The transport collision integrals are calculated by fitting various two-body semiempirical interaction potentials for which the collision integrals are tabulated to ab initio quantum mechanical calculations of the two-body interactions. The collision integrals are used to calculate the binary diffusion coefficients, viscosity, and 'total' thermal conductivity of the pure gases and the gas mixtures at 1-atm pressure from 1000 K to 25,000 K.

Biolsi, L.↗

Anharmonic Rovibrational Calculations of Singlet Cyclic C4 Using a New Ab Initio Potential and a Quartic Force

We report a CCSD(T)/cc-pCV5Z quartic force field (QFF) and a semi-global CCSD(T)-F12b/aug-cc-pVTZ potential energy surface (PES) for singlet, cyclic C4. Vibrational fundamentals, combinations and overtones are obtained using vibrational second-order perturbation theory (VPT2) and the vibrational configurationinteraction (VCI) approach. Agreement is within 10 cm(exp −1) between the VCI calculated fundamentals on the QFF and PES using the MULTIMODE (MM) program, and VPT2 and VCI results agree for the fundamentals. The agreement between VPT2- QFF and MM-QFF results is also good for the C4 combinations and overtones. The J = 1 and J = 2 rovibrational energies are reported from both VCI (MM) on the PES and VPT2 on the QFF calculations. The spectroscopic constants of (12)C4 and two C(sub 2v)-symmetry, single (13)C-substituted isotopologues are presented, which may help identification of cyclic C4 in future experimental analyses or astronomical observations.

quartic force field↗

Towards Accurate Predictions of Martensitic Transition Temperatures for Shape Memory Alloys from Ab Initio Simulations

Experimentally, NiTi undergoes a single martensitic phase transition around 341 K from the lowtemperature (T) monoclinic B19’ phase (P21/m) to the high-temperature cubic B2 phase (Pm3m). Theoretically, an orthorhombic B33 (Cmcm) has also been proposed as the T=0 ground state structure, although this phase has never been observed. Accurate predictions of martensitic transition temperatures (MTT) have remained elusive in part due to several well-known theoretical complexities of these systems including low temperature instabilities of the B2 phase. Recently, we proposed a rigorous thermodynamic integration approach based on ab initio simulations to resolve many of these difficulties [1,2]. However, an unsatisfying overprediction of the MTT relative to experiment (by ~100 K) means a fully quantitative theory is still lacking. In this work, we report several new developments to our method that bring first principles theory and experiment much closer into agreement. Our calculations indicate that phonon free energies at low temperature stabilizes B19’ over B33, rationalizing B19’ as the ground state down to T=0. We also find that accurate computations of the electronic free energy, i.e. the change in energy and the appearance of electronic configurational entropy due to finite temperature, is crucial to obtain accurate MTT. Incorporating these corrections results in an MTT prediction of 365 K for binary NiTi, which is in very close agreement with experiment. Our theoretical approach is expected to be a broadly applicable and predictive theory for MTT of SMAs.

Zhigang Wu↗

Ab initio ground states of strongly-correlated materials on quantum computers

The accurate first-principles description of strongly-correlated materials is an important and challenging problem in condensed matter physics. Ab initio downfolding has emerged as a way of deriving accurate many-body Hamiltonians including strong correlations, representing a subspace of interest of a material, using density functional theory calculations as a starting point. However, the solution of these material-specific models can scale exponentially on classical computers, constituting a challenge. Here we propose that utilizing quantum computers for obtaining the properties of downfolded Hamiltonians yields an accurate description of the ground state properties of strongly-correlated systems, while circumventing the exponential scaling problem. We benchmark the solution of Hubbard-like models obtained through downfolding by utilizing a classical tensor network implementation of variational quantum eigensolvers (VQE), and we reveal a strategy for driving the optimization through a hybrid minimization of the energy and maximization of the overlap with an approximate solution obtained through low-cost computational methods. This results in a reduction of the energy error by orders of magnitude compared to conventional VQE approaches, and allows us to reproduce long-range correlations for the first time. We demonstrate our first-principles approach for diverse strongly-correlated materials, correctly predicting the antiferromagnetic state of one-dimensional cuprate Ca 2 CuO 3 , the excitonic ground state of monolayer WTe2, and the charge-ordered state of correlated metal SrVO 3 . Our efficient computational implementation allows us to simulate large systems with up to 54 qubits and encompassing up to four correlated bands, which is indicative of the complexity that our framework can address.

Antonios M Alvertis↗

Theory of collision-induced translation-rotation spectra: H2-He

An adiabatic quantal theory of spectral line shapes in collision-induced absorption and emission is presented which incorporates the induced translation-rotation and translation-vibration spectra. The generalization to account for the anisotropy of the scattering potential is given. Calculations are carried out of the collision-induced absorption spectra of He in collisions with H2 with ab initio electric dipole functions and realistic potentials. The anisotropy of the interaction potential is small and is not included in the calculations. The predicted spectra are in satisfactory agreement with experimental data though some deviations occur which may be significant. The rotational line shapes have exponential wings and are not Lorentzian. The connection between the quantal and classical theories is written out explicitly for the isotropic overlap induction.

Birnbaum, G.↗

Analysis of the far infrared H2-He spectrum

Previous measurements of the far infrared absorption due to H2-He collisions at the temperatures of 77, 195, and 292 K are analyzed. The spectra are fitted by a semiempirical line shape representing the isotropic induced overlap component and combined anisotropic quadrupolar and overlap components. The experimental spectral moments are evaluated and compared with theory for several induced-dipole and potential models. From the isotropic contribution, the range and strength of the induced dipole is evaluated and compared with the results of ab initio calculations. Fitting parameters are obtained with physically plausible temperature dependences which allow simple and accurate representation of the spectra and of their moments at temperatures different from those of the measurements.

Cohen, E. R.↗

Calculations of rate constants for the three-body recombination of H2 in the presence of H2

A new global potential energy hypersurface for H2 + H2 is constructed and quasiclassical trajectory calculations performed using the resonance complex theory and energy transfer mechanism to estimate the rate of three body recombination over the temperature range 100 to 5000 K. The new potential is a faithful representation of ab initio electron structure calculations, is unchanged under the operation of exchanging H atoms, and reproduces the accurate H3 potential as one H atom is pulled away. Included in the fitting procedure are geometries expected to be important when one H2 is near or above the dissociation limit. The dynamics calculations explicitly include the motion of all four atoms and are performed efficiently using a vectorized variable-stepsize integrator. The predicted rate constants are approximately a factor of two smaller than experimental estimates over a broad temperature range.

Schwenke, David W.↗

Calculations of rate constants for the three-body recombination of H2 in the presence of H2

A new global potential energy hypersurface for H2 + H2 is constructed and quasiclassical trajectory calculations performed using the resonance complex theory and energy transfer mechanism to estimate the rate of three body recombination over the temperature range 100 to 5000 K. The new potential is a faithful representation of ab initio electron structure calculations, is unchanged under the operation of exchanging H atoms, and reproduces the accurate H3 potential as one H atom is pulled away. Included in the fitting procedure are geometries expected to be important when one H2 is near or above the dissociation limit. The dynamics calculations explicitly include the motion of all four atoms and are performed efficiently using a vectorized variable-stepsize integrator. The predicted rate constants are approximately a factor of two smaller than experimental estimates over a broad temperature range.

Schwenke, David W.↗

Rototranslational absorption spectra of H2-H2 pairs in the far infrared

On the basis of ab initio induced dipole components and a well-tested isotropic potential, the collision-induced absorption rototranslational profiles of hydrogen have been determined using an exact wave-mechanical theory. The main induction mechanisms, particularly the quadrupole-induced 0223,2023 components, are found to be consistent with observations. It is demonstrated that the present theory can be used to predict spectra at temperatures where no measurements exist, as is needed for the modeling of planetary and stellar atmospheres.

Meyer, Wilfried↗

Accurate Quartic Force Fields and Vibrational Frequencies for HCN and HNC

The quartic force fields of HCN and HNC are determined using atomic natural orbital one-particle basis sets of spdf/spd and spdfg/spdf quality in conjunction with the CCSD(T) electron correlation method (singles and doubles coupled-cluster theory plus a perturbational estimate of the effects of connected triple excitations). The HCN force field is in good agreement with a recent experimentally derived force field and also with the force field recently computed by Wong and Bacskay. On the basis of the good agreement obtained for HCN, it is argued that the ab initio quartic force field for HNC is superior to a prior force field derived from experiment. The harmonic frequencies of HNC are predicted to be 3822 +/- 10,472 +/- 5, and 2051 +/-10/cm for omega1, omega2, and omega3, respectively; the experimentally derived values are above these values and fall outside the estimated uncertainties. Using the quartic force field, spectroscopic constants are predicted for HNC based on a vibrational second-order perturbation theory analysis. It is also asserted that the gas-phase fundamental v(sub 3) for HNC is slightly lower than the matrix isolation value. The range of validity of the quartic force fields is investigated by comparison of variational vibrational energies computed with the quartic force fields to those obtained from our recently reported global HCN/HNC potential energy surface and also to experimental data.

Lee, Timothy J.↗

Molecular orbital predictions of the vibrational frequencies of some molecular ions

The initial detections of IR vibration-rotation bands in polyatomic molecular ions by recent spectroscopic advances were guided by ab initio prediction of vibrational frequencies. The present calculations predict the vibrational frequencies of additional ions which are candidates for laboratory analysis. Neutral molecule vibrational frequencies were computed at three levels of theory and then compared with experimental data; the effect of scaling was also investigated, in order to determine how accurately vibrational frequencies could be predicted. For 92 percent of the frequencies examined, the relatively simple HF/6-31G theory's vibrational frequencies were within 100/cm of experimental values, with a mean absolute error of 49/cm. On this basis, the frequencies of 30 molecular ions (many possessing astrophysical significance) were computed.

Defrees, D. J.↗

Theoretical determination of chemical rate constants using novel time-dependent methods

The work completed within the grant period 10/1/91 through 12/31/93 falls primarily in the area of reaction dynamics using both quantum and classical mechanical methodologies. Essentially four projects have been completed and have been or are in preparation of being published. The majority of time was spent in the determination of reaction rate coefficients in the area of hydrocarbon fuel combustion reactions which are relevant to NASA's High Speed Research Program (HSRP). These reaction coefficients are important in the design of novel jet engines with low NOx emissions, which through a series of catalytic reactions contribute to the deterioration of the earth's ozone layer. A second area of research studied concerned the control of chemical reactivity using ultrashort (femtosecond) laser pulses. Recent advances in pulsed-laser technologies have opened up a vast new field to be investigated both experimentally and theoretically. The photodissociation of molecules adsorbed on surfaces using novel time-independent quantum mechanical methods was a third project. And finally, using state-of-the-art, high level ab initio electronic structure methods in conjunction with accurate quantum dynamical methods, the rovibrational energy levels of a triatomic molecule with two nonhydrogen atoms (HCN) were calculated to unprecedented levels of agreement between theory and experiment.

Dateo, Christopher E.↗