Search NASA⌕ Search

SEARCH · Search NASA

Results for “quantum chemical calculations”

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 19 records

Accurate quantum chemical calculations

An important goal of quantum chemical calculations is to provide an understanding of chemical bonding and molecular electronic structure. A second goal, the prediction of energy differences to chemical accuracy, has been much harder to attain. First, the computational resources required to achieve such accuracy are very large, and second, it is not straightforward to demonstrate that an apparently accurate result, in terms of agreement with experiment, does not result from a cancellation of errors. Recent advances in electronic structure methodology, coupled with the power of vector supercomputers, have made it possible to solve a number of electronic structure problems exactly using the full configuration interaction (FCI) method within a subspace of the complete Hilbert space. These exact results can be used to benchmark approximate techniques that are applicable to a wider range of chemical and physical problems. The methodology of many-electron quantum chemistry is reviewed. Methods are considered in detail for performing FCI calculations. The application of FCI methods to several three-electron problems in molecular physics are discussed. A number of benchmark applications of FCI wave functions are described. Atomic basis sets and the development of improved methods for handling very large basis sets are discussed: these are then applied to a number of chemical and spectroscopic problems; to transition metals; and to problems involving potential energy surfaces. Although the experiences described give considerable grounds for optimism about the general ability to perform accurate calculations, there are several problems that have proved less tractable, at least with current computer resources, and these and possible solutions are discussed.

Bauschlicher, Charles W., Jr.↗

Quantum Chemical Calculations

The current methods of quantum chemical calculations will be reviewed. The accent will be on the accuracy that can be achieved with these methods. The basis set requirements and computer resources for the various methods will be discussed. The utility of the methods will be illustrated with some examples, which include the calculation of accurate bond energies for SiF$_n$ and SiF$_n^+$ and the modeling of chemical data storage.

Bauschlicher, Charles W.↗

Quantum chemical calculation of the equilibrium structures of small metal atom clusters

The application of ab initio quantum mechanical approaches in the study of metal atom clusters requires simplifying techniques that do not compromise the reliability of the calculations. Various aspects of the implementation of the effective core potential (ECP) technique for the removal of the metal atom core electrons from the calculation were examined. The ECP molecular integral formulae were modified to bring out the shell characteristics as a first step towards fulfilling the increasing need to speed up the computation of the ECP integrals. Work on the relationships among the derivatives of the molecular integrals that extends some of the techniques pioneered by Komornicki for the calculation of the gradients of the electronic energy was completed and a formulation of the ECP approach that quite naturally unifies the various state-of-the-art "shape- and Hamiltonian-consistent" techniques was discovered.

Kahn, L. R.↗

Quantum chemical calculation of the equilibrium structures of small metal atom clusters

Metal atom clusters are studied based on the application of ab initio quantum mechanical approaches. Because these large 'molecular' systems pose special practical computational problems in the application of the quantum mechanical methods, there is a special need to find simplifying techniques that do not compromise the reliability of the calculations. Research is therefore directed towards various aspects of the implementation of the effective core potential technique for the removal of the metal atom core electrons from the calculations.

Kahn, L. R.↗

Quantum chemical calculations of interatomic potentials for computer simulation of solids

A comprehensive mathematical model by which the collective behavior of a very large number of atoms within a metal or alloy can accurately be simulated was developed. Work was done in order to predict and modify the strength of materials to suit our technological needs. The method developed is useful in studying atomic interactions related to dislocation motion and crack extension.

Source record↗

Quantum chemical calculation of the equilibrium structures of small metal atom clusters

A decomposition of the molecular energy is presented that is motivated by the atom superposition and electron delocalization physical model of chemical binding. The energy appears in physically transparent form consisting of a classical electrostatic interaction, a zero order two electron exchange interaction, a relaxation energy, and the atomic energies. Detailed formulae are derived in zero and first order of approximation. The formulation extends beyond first order to any chosen level of approximation leading, in principle, to the exact energy. The structure of this energy decomposition lends itself to the fullest utilization of the solutions to the atomic sub problems to simplify the calculation of the molecular energy. If nonlinear relaxation effects remain minor, the molecular energy calculation requires at most the calculation of two center, two electron integrals. This scheme thus affords the prospects of substantially reducing the computational effort required for the calculation of molecular energies.

Kahn, L. R.↗

Quantum chemical calculations for polymers and organic compounds

The relativistic effects of the orbiting electrons on a model compound were calculated. The computational method used was based on 'Modified Neglect of Differential Overlap' (MNDO). The compound tetracyanoplatinate was used since empirical measurement and calculations along "classical" lines had yielded many known properties. The purpose was to show that for large molecules relativity effects could not be ignored and that these effects could be calculated and yield data in closer agreement to empirical measurements. Both the energy band structure and molecular orbitals are depicted.

Lopez, J.↗

Metastable isomers - A new class of interstellar molecules

The abundances of a variety of metastable isomers of small organic molecules, analogous to HNC/HCN, in dense interstellar clouds are considered. These metastable species, some of which are thought to exist as intermediates in laboratory organic chemical reactions, are of considerable interest to chemists. Current ideas of gas-phase, ion-molecule chemistry are utilized to demonstrate that such metastable species should often be present in dense clouds in sufficient abundance to be observed. Unfortunately, the spectral constants of metastable isomers have rarely been determined in the laboratory, and quantum chemical calculations of a varying degree of accuracy must be utilized; results are included of some new quantum chemical calculations. The interstellar chemistry and expected microwave spectra of a representative sample of possibly important interstellar metastable isomers are discussed.

Green, S.↗

Applicability of PM3 to transphosphorylation reaction path: Toward designing a minimal ribozyme

A growing body of evidence shows that RNA can catalyze many of the reactions necessary both for replication of genetic material and the possible transition into the modern protein-based world. However, contemporary ribozymes are too large to have self-assembled from a prebiotic oligonucleotide pool. Still, it is likely that the major features of the earliest ribozymes have been preserved as molecular fossils in the catalytic RNA of today. Therefore, the search for a minimal ribozyme has been aimed at finding the necessary structural features of a modern ribozyme (Beaudry and Joyce, 1990). Both a three-dimensional model and quantum chemical calculations are required to quantitatively determine the effects of structural features of the ribozyme on the reaction it catalyzes. Using this model, quantum chemical calculations must be performed to determine quantitatively the effects of structural features on catalysis. Previous studies of the reaction path have been conducted at the ab initio level, but these methods are limited to small models due to enormous computational requirements. Semiempirical methods have been applied to large systems in the past; however, the accuracy of these methods depends largely on a simple model of the ribozyme-catalyzed reaction, or hydrolysis of phosphoric acid. We find that the results are qualitatively similar to ab initio results using large basis sets. Therefore, PM3 is suitable for studying the reaction path of the ribozyme-catalyzed reaction.

Manchester, John I.↗

Calculations on the rate of the ion-molecule reaction between NH3(+) and H2

The rate coefficient for the ion-molecule reaction NH3(+) + H2 yields NH4(+) + H has been calculated as a function of temperature with the use of the statistical phase space approach. The potential surface and reaction complex and transition state parameters used in the calculation have been taken from ab initio quantum chemical calculations. The calculated rate coefficient has been found to mimic the unusual temperature dependence measured in the laboratory, in which the rate coefficient decreases with decreasing temperature until 50-100 K and then increases at still lower temperatures. Quantitative agreement between experimental and theoretical rate coefficients is satisfactory given the uncertainties in the ab initio results and in the dynamics calculations. The rate coefficient for the unusual three-body process NH3(+) + H2 + He yields NH4(+) + H + He has also been calculated as a function of temperature and the result found to agree well with a previous laboratory determination.

Herbst, Eric↗

New calculations on the ion-molecule processes C2H2(+) + H2 --> C2H3(+) + H and C2H2(+) + H2 --> C2H4+

New high-level quantum chemical calculations have been undertaken to understand the rates and mechanisms of the reactive and associative channels for the reactants C2H2(+) + H2. The reactive channel, which produces C2H3(+) + H, has been shown to be slightly endothermic, confirming earlier calculations at a somewhat lower level and in agreement with some recent experimental work. The associative channel, leading to C2H4+, has been shown to proceed via a transition state with negative energy relative to the reactants, so that association is predicted to be efficient. This result is in conflict with an earlier theoretical study but in agreement with low-temperature experimental measurements.

NASA Discipline Exobiology↗

Ab initio studies of dissociative recombination

Quantum chemical calculations of the dissociative recombination of O2(+) and N2(+) are reported. An approach for calculating autoionization widths from high-principal-quantum-number Rydberg states is summarized, and an example is presented for the lowest dissociative state of O2. For O2(+), the 1Sigma(+)u state is the sole source of O(1S) from the lowest 10 vibrational levels of the ion. Rate coefficients for generating O(1S) and O(1D) at ionospheric temperatures are reported.

Guberman, Steven L.↗

Mechanism for the Green Glow of the Upper Ionosphere

The generation of the green line of atomic oxygen by dissociative recombination of 02 plus occurs by the capture of an electron into a repulsive state of 02 followed by dissociation along another state of a different electronic symmetry. The two states are coupled together by mixed symmetry Rydberg states. Quantum chemical calculations give a rate coefficient at room temperature of (0.39 (+ 0.31 or -0.19)) x 10 exp -8 cubic centimeters per second. The quantum yield of excited oxygen is within the range deduced from ground, rocket, and satellite observations. The rate coefficients and yields are needed in models of the optical emission, chemistry, and energy balance of planetary ionospheres.

Guberman, Steven L.↗

Quantum mechanical calculations to chemical accuracy

The accuracy of current molecular-structure calculations is illustrated with examples of quantum mechanical solutions for chemical problems. Two approaches are considered: (1) the coupled-cluster singles and doubles (CCSD) with a perturbational estimate of the contribution of connected triple excitations, or CCDS(T); and (2) the multireference configuration-interaction (MRCI) approach to the correlation problem. The MRCI approach gains greater applicability by means of size-extensive modifications such as the averaged-coupled pair functional approach. The examples of solutions to chemical problems include those for C-H bond energies, the vibrational frequencies of O3, identifying the ground state of Al2 and Si2, and the Lewis-Rayleigh afterglow and the Hermann IR system of N2. Accurate molecular-wave functions can be derived from a combination of basis-set saturation studies and full configuration-interaction calculations.

Bauschlicher, Charles W., Jr.↗

Cavity Ring-Down Spectroscopy of Anthracene, 9-Methylanthracene, and 2 Methylanthracene in Supersonic Expansion

The S0–S1 absorption spectra of anthracene (C14H10), 9-methylanthracene (C15H12), and 2-methylanthracene (C15H12) are measured in the ultraviolet region between 330 and 375 nm (26,666 to 30,303 cm–1) with cavity ring-down spectroscopy in supersonic free-jet expansions of argon. The associated vibronic band systems and their spectroscopic assignments are discussed and compared to previous studies performed using fluorescence excitation and dispersed fluorescence techniques. Density functional theory (DFT) calculations were carried out to study the structures and evaluate the vibrational transitions of the ground and excited states. Time-dependent DFT calculations of the first electronic excited states and Franck–Condon factor calculations were carried out to assist in the assignment of the experimentally measured vibronic bands. The vibronic spectra obtained in absorption agree well with fluorescence excitation spectra in terms of peak positions but exhibit different relative band intensities. We find a very good match between experimentally obtained vibronic line positions and the peak positions of the quantum chemically calculated Franck–Condon excitation lines.

Salma Bejaoui↗

Theoretical Studies of Important Processes in Planetary and Comet Atmospheres

Using theoretical quantum chemical calculations, I have successfully described the dissociative recombination (DR) of O2(+) leading to the excited S-1 state of atomic oxygen, the upper state of the well known green line emission. The process is described by O2(+) + e(-) yields O(S-1)) + O(D-1) (1) where e(-) is an electron and the product oxygen atoms are both excited. This process is important in the atmospheres of Venus, Mars and Earth. I have shown in prior calculations that only one repulsive potential curve of O2, f(sup l)Sigma(sub u, sup +) can generate O(S-1) from DR of the lowest vibrational levels of O2(+). However, in the prior results, the calculated quantum yield (i.e. the number of O(S-1) atoms produced for every two product atoms) from the O2(+) v = 0 level was smaller than the laboratory and atmospheric measurements by more than an order of magnitude. Including only direct recombination, the calculated quantum yield for O(S-1) is only 0.0016. In a calculation that accounts for both direct and indirect recombination, the quantum yield is 0.0012. The range of experimentally determined quantum yields is between 0.01 and 0.23. Because of this large difference between the theoretical and experimental quantum yields, it was thought for some time that the ionospheric and laboratory O2(+) must be vibrationally excited since for excited levels, theory gave quantum yields that are similar to experimental yields. It was also suggested that some other process was generating O(S-1) but none could be identified. Under current NASA support, I have found that reaction (1) proceeds via an unusual mechanism. The f(sup 1)Sigma(sub u, sup,+) state does not cross the ion between the turning points of the v = 0 level of the O2(+) ground state. The lack of a favorable crossing leads to a very small calculated DR rate coefficient. However, this mechanism assumes that initial electron capture must occur into the repulsive state that leads to )(S-1). Instead, I have found that initial electron capture occurs mostly into the B(sup 3)Sigma(sup u, sub -) state which crosses the ion between the turning points of the v = 0 ion level and has a large DR rate coefficient. The B state dissociates to O(D-1) and )(P-3). After capture, some of the flux is transferred to the f(sup 1)Sigma(sup u, sub +) state via symmetry mixed intermediate Rydberg states. The neutral Rydberg states are a mixture of Sigma(sub u, sup +)-1 and Sigma(sub u, sup -)-3 symmetry.

Guberman, Steven L.↗

Computational Study of Electron-Molecule Collisions Related to Low-Temperature Plasmas

Computational study of electron-molecule collisions not only complements experimental measurements, but can also be used to investigate processes not readily accessible experimentally. A number of ab initio computational methods are available for these types of calculations. Here we describe a recently developed technique, the finite element Z-matrix method, Analogous to the R-matrix, method, it partitions the space into regions and employs real matrix elements. However, unlike the implementation of the R-matrix method commonly used in atomic and molecular physics, the Z-matrix method is fully variational. In the present implementation, a mixed basis of finite elements and Gaussians is used to represent the continuum electron, thus offering full flexibility without imposing fixed boundary conditions. Numerical examples include the electron-impact dissociation of N2 via the metastable A3Su+ state, a process which may be important in the lower thermosphere, and the dissociation of the CF radical, a process of interest to plasma etching. To understand the dissociation pathways, large scale quantum chemical calculations have been carried out for all target states which dissociate to the lowest five limits in the case of N2, and to the lowest two limits in the case of CF. For N2, the structural calculations clearly show the preference for predissociation if the initial state is the ground X1 Sg+ state, but direct dissociation appears to be preferable if the initial state is the A3Su+ state. Multi-configuration SCF target functions are used in the collisional calculation.

Huo, Winifred M.↗

Comment on the relation between the nonadiabatic coupling and the complex intersection of potential energy curves

Simple relations are discussed that provide a correspondence between the complex intersection of two potential surfaces and the nonadiabatic coupling matrix element between those surfaces. These are key quantities in semiclassical and quantum mechanical theories of collision induced electronic transitions. Within the two state approximation, the complex intersection is shown to be directly related to the location and magnitude of the peak in the nonadiabatic coupling. Two cases are considered: the avoided crossing between two potential surfaces; and the spin orbit interaction due to a P-2 halogen atom. Comparisons are made between the results of the two-state model and the results of ab initio quantum chemical calculations.

Jaffe, R. L.↗