Search NASA⌕ Search

SEARCH · Search NASA

Results for “quartic force field”

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 37 records · Page 2

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.↗

The Anharmonic Force Field of Ethylene, C2H4, by Means of Accurate Ab Initio Calculations

The quartic force field of ethylene, C2H4, has been calculated ab initio using augmented coupled cluster, CCSD(T), methods and correlation consistent basis sets of spdf quality. For the C-12 isotopomers C2H4, C2H3D, H2CCD2, cis-C2H2D2, trans-C2H2D2, C2HD3, and C2D4, all fundamentals could be reproduced to better than 10 per centimeter, except for three cases of severe Fermi type 1 resonance. The problem with these three bands is identified as a systematic overestimate of the Kiij Fermi resonance constants by a factor of two or more; if this is corrected for, the predicted fundamentals come into excellent agreement with experiment. No such systematic overestimate is seen for Fermi type 2 resonances. Our computed harmonic frequencies suggest a thorough revision of the accepted experimentally derived values. Our computed and empirically corrected re geometry differs substantially from experimentally derived values: both the predicted rz geometry and the ground-state rotational constants are, however, in excellent agreement with experiment, suggesting revision of the older values. Anharmonicity constants agree well with experiment for stretches, but differ substantially for stretch-bend interaction constants, due to equality constraints in the experimental analysis that do not hold. Improved criteria for detecting Fermi and Coriolis resonances are proposed and found to work well, contrary to the established method based on harmonic frequency differences that fails to detect several important resonances for C2H4 and its isotopomers. Surprisingly good results are obtained with a small spd basis at the CCSD(T) level. The well-documented strong basis set effect on the v8 out-of-plane motion is present to a much lesser extent when correlation-optimized polarization functions are used. Complete sets of anharmonic, rovibrational coupling, and centrifugal distortion constants for the isotopomers are available as supplementary material to the paper.

Martin, Jan M. 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↗

Excited Vibrational Level Rotational Constants for SiC2: A Sensitive Molecular Diagnostic for Astrophysical Conditions

Silacyclopropynylidene, SiC2, is a known and highly abundant circumstellar molecule. Its spectrum has been established as a major component of lines observed toward the carbon-rich star IRC +10216 (CW Leonis). It has been detected in its low-lying v(sub 3) = 1 and 2 vibrational states as well as in various isotopic compositions. Increasing sensitivity and spatial resolution will enable many more emission or absorption lines to be detected. In order to detect new molecular species, unassigned lines of known species must be identified. This work uses established ab initio quartic force fields to produce data necessary for this classification of lines related to SiC2. Agreement between the theoretical vibrational frequencies and known rotational and spectroscopic constants is quite good, as good as 5 cm(exp -1) and 3 MHz, respectively in some cases. Additionally, experimentally unknown vibrational frequencies and rotational constants are provided for the first overtones and combination bands in addition to 3(sub v3), the second overtone of the low-lying antisymmetric stretch/carbide rotation mode. Frequencies of v(sub 3) = 3 low-J rotational transitions of the main isotopic species are also estimated from published data for v(sub 3) ≤ 2. Further, we determine rotational and centrifugal distortion parameters for which in most cases vibrational effects due to the v(sub 3) mode were reduced to first, and in several cases also to second order. These values may approximate equilibrium values better than the ground state values. The data produced herein will aid in the experimental and observational characterization of this known astromolecule in order to identify some of the unassigned lines for a known entity.

Astrochemistry↗

Helium droplet infrared spectroscopy of the butyl radicals

Butyl radicals (n-, s-, i-, and tert-butyl) are formed from the pyrolysis of stable precursors (1-pentyl nitrite, 2-methyl-1-butyl nitrite, isopentyl nitrite, and azo-tert-butane, respectively). The radicals are doped into a beam of liquid helium droplets and probed with infrared action spectroscopy from 2700 to 3125 cm -1 , allowing for a low temperature measurement of the CH stretching region. The presence of anharmonic resonance polyads in the 2800–3000 cm -1 region complicates its interpretation. To facilitate spectral assignment, the anharmonic resonances are modeled with two model Hamiltonian approaches that explicitly couple CH stretch fundamentals to HCH bend overtones and combinations: a VPT2+K normal mode model based on coupled-cluster with single, double, and perturbative triple excitations [CCSD(T)] quartic force fields and a semi-empirical local mode model. Both of these computational methods provide generally good agreement with the experimental spectra.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Finite-temperature vibrational full configuration interaction

Thermodynamic functions of an ideal molecular gas due to its anharmonic vibrations are evaluated in a wide range of temperature (T) by the vibrational full-configuration-interaction (FCI) method using a quartic force field and a finite number (N) of harmonic-oscillator basis functions along each normal mode. The thermodynamic functions considered are the grand potential (Ω), internal energy (U), and entropy (S). They are compared with those obtained from the Bose–Einstein theory with or without truncation of the harmonic-oscillator basis functions after quantum number N–1. The comparison reveals that the finite-basis-set errors in Ω and U are, respectively, k B Tln(k B T/Nℏω) and k B T per mode in the high-T limit, obscuring anharmonic effects when k B T > ℏω, where ω is the lowest mode frequency. Here, the benchmark data for several low-order perturbation corrections to Ω, U, and S are also obtained as the numerical derivatives of their FCI values with respect to dimensionless perturbation strength, and the domain of T and N in which these data are reliable (for the N → ∞ limits) is discussed.

74 ATOMIC AND MOLECULAR PHYSICS↗

Variational calculations of vibrational properties of ozone

A variational method is used to obtain vibrational-rotational properties for ozone from an experimental quartic force field. Band positions, average structures, matrix elements for calculating infrared intensities, and effective rotational constants are reported for (O-16)3 and its O-18 isotopic species. Also, the degree to which the vibrational energies and properties are converged is investigated as a function of the basis set parameters and basis set size, and of the method of obtaining the Hamiltonian matrix elements. A convenient procedure for assigning the vibrational states is developed for variational vibrational wavefunction expansions.

Carney, G. D.↗

Vibrations in small Mg clusters

The equilibrium geometries, binding energies, and harmonic frequencies of Mg3 and Mg4 have been determined using large atomic natural orbital basis sets in conjunction with high levels of electron correlation. The correlation treatments comprise multireference configuration interaction (MRCI), singles and doubles coupled-cluster (CCSD) theory and the CCSD(T) extension that includes a perturbational estimate of connected triple excitations. As with previous studies of small Be clusters, the CCSD(T) method is found to reproduce the MRCI results with a remarkable degree of accuracy and at a fraction of the computational cost. Using the CCSD(T) method, full quartic force fields for Mg3 and Mg4 have been determined and anharmonic analyses have been performed using second-order perturbation theory. Vibration-rotation interaction constants and centrifugal distortion constants have also been determined.

Lee, Timothy J.↗

The structures, binding energies and vibrational frequencies of Ca3 and Ca4: An application of the CCSD(T) method

The Ca3 and Ca4 metallic clusters have been investigated using state-of-the-art ab initio quantum mechanical methods. Large atomic natural orbital basis sets have been used in conjunction with the singles and doubles coupled-cluster (CCSD) method, a coupled-cluster method that includes a perturbational estimate of connected triple excitations, denoted CCSD(T), and the multireference configuration interaction (MRCI) method. The equilibrium geometries, binding energies and harmonic vibrational frequencies have been determined with each of the methods so that the accuracy of the coupled-cluster methods may be assessed. Since the CCSD(T) method reproduces the MRCI results very well, cubic and quartic force fields of Ca3 and Ca4 have been determined using this approach and used to evaluate the fundamental vibrational frequencies. The infrared intensities of both the e' mode of Ca3 and the t2 mode of Ca4 are found to be small. The results obtained in this study are compared and contrasted with those from our earlier studies on small Be and Mg clusters.

Lee, Timothy J.↗

The Heat of Formation of HNO

The HNO molecule is of interest in both combustion and atmospheric chemistry. For example, Guadagnini et al. have recently presented ab initio potential energy surfaces for the three lowest lying electronic states of HNO and then used these in examining several chemical reactions that take place in the combustion of nitrogen containing fuels and in the oxidation of atmospheric nitrogen. We have previously studied the ground state potential energy surface (i.e., stationary points along the HNO rev. reaction HON path), vibrational spectrum (using an accurate quartic force field), zero-point energy, and bonding of HNO using coupled-cluster ab initio methods. HNO is also very interesting because of the unique nature of its bonding characteristics. That is, the potential energy surface is very flat along the H-N bonding coordinate thereby giving unusual harmonic and fundamental vibrational frequencies, and the H-N bond energy is rather weak in comparison to other H-N bond energies. In fact, using experimental heats of formation for HO, H, and NO, the H- bond energy is computed to be only 49.9 kcal/ mol (298 K). However, ab initio calculations of isodesmic reaction energies involving HNO, FNO, ClNO, and several other molecules have shown that there is an inconsistency in the experimental heats of formation of the XNO (X=H, F, and Cl) species. Hence the motivation for this study was to determine a very accurate(DELTA)H(sup o)(sub f) value for HNO using state of-the-art ab initio methods. Based on many recent studies it is evident that the singles and doubles coupled-cluster method that includes a perturbational estimate of the effects of connected triple excitations, denoted CCSD(T), in conjunction with large one-particle basis sets should be reliable to better than +0.8 kcal/mol for this quantity. The computational methodology is described in the next section followed by our results and discussion. Conclusions are presented in the final section.

Lee, Timothy J.↗

The Heat of Formation of HNO

The HNO molecule is of interest in both combustion and atmospheric chemistry. For example, Guadagnini et al. have recently presented ab initio potential energy surfaces for the three lowest lying electronic states of HNO and then used these in examining several chemical reactions that take place in the combustion of nitrogen containing fuels and in the oxidation of atmospheric nitrogen. We have previously studied the ground state potential energy surface (i.e., stationary points along the HNO reversible reaction HON path), vibrational spectrum (using an accurate quartic force field), zero-point energy, and bonding of HNO using coupled-cluster ab initio methods. HNO is also very interesting because of the unique nature of its bonding characteristics. That is, the potential energy surface is very flat along the H-N bonding coordinate thereby giving unusual harmonic and fundamental vibrational frequencies, and the H-N bond energy is rather weak in comparison to other H-N bond energies. In fact, using experimental heats of formation for HNO, H, and NO, the H- bond energy is computed to be only 49.9 kcal/ mol (298 K). However, ab initio calculations of isodesmic reaction energies involving HNO, FNO, ClNO, and several other molecules have shown that there is an inconsistency in the experimental heats of formation of the XNO (X double bond H, F, and Cl) species. Hence the motivation for this study was to determine a very accurate (Delta)H(sub f, sup o) value for HNO using state-of-the-art ab initio methods. Based on many recent studies it is evident that the singles and doubles coupled-cluster method that includes a perturbational estimate of the effects of connected triple excitations, denoted CCSD(T), in conjunction with large one-particle basis sets should be reliable to better than +/- 0.8 kcal/mol for this quantity. The computational methodology is described in the next section followed by our results and discussion. Conclusions are presented in the final section.

Lee, Timothy J.↗

Ab Initio Potential Energy Surfaces and the Calculation of Accurate Vibrational Frequencies

Due to advances in quantum mechanical methods over the last few years, it is now possible to determine ab initio potential energy surfaces in which fundamental vibrational frequencies are accurate to within plus or minus 8 cm(exp -1) on average, and molecular bond distances are accurate to within plus or minus 0.001-0.003 Angstroms, depending on the nature of the bond. That is, the potential energy surfaces have not been scaled or empirically adjusted in any way, showing that theoretical methods have progressed to the point of being useful in analyzing spectra that are not from a tightly controlled laboratory environment, such as vibrational spectra from the interstellar medium. Some recent examples demonstrating this accuracy will be presented and discussed. These include the HNO, CH4, C2H4, and ClCN molecules. The HNO molecule is interesting due to the very large H-N anharmonicity, while ClCN has a very large Fermi resonance. The ab initio studies for the CH4 and C2H4 molecules present the first accurate full quartic force fields of any kind (i.e., whether theoretical or empirical) for a five-atom and six-atom system, respectively.

Lee, Timothy J.↗

Ab Initio Calculation of Accurate Vibrational Frequencies for Molecules of Interest in Atmospheric Chemistry

Due to advances in quantum mechanical methods over the last few years, it is now possible to determine ab initio potential energy surfaces in which fundamental vibrational frequencies are accurate to within +/- 8 cm(sup -1) on average, and molecular bond distances are accurate to within +/- 0.001-0.003 A, depending on the nature of the bond. That is, the potential energy surfaces have not been scaled or empirically adjusted in any way, showing that theoretical methods have progressed to the point of being useful in analyzing spectra that are not from a tightly controlled laboratory environment, such as rovibrational spectra from the interstellar medium. Some recent examples demonstrating this accuracy win be presented and discussed. These include the HNO, CH4, C2H4, and ClCN molecules. The HNO molecule is interesting due to the very large H-N anharmonicity, while ClCN has a very large Fermi resonance. The ab initio studies for the CH4 and C2H4 molecules present the first accurate full quartic force fields of any kind (i.e., whether theoretical or empirical) for a five-atom and six-atom system, respectively.

Lee, Timothy J.↗

The Effect of Approximating Some Molecular Integrals in Coupled-Cluster Calculations: Fundamental Frequencies and Rovibrational Spectroscopic Constants of Cyclopropenylidene

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

Lee, Timothy J.↗

Reliable Spectroscopic Constants for CCH-, NH2- and Their Isotopomers from an Accurate Potential Energy Function

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.

Lee, Timothy J.↗

Rovibrational Spectroscopic Constants and Fundamental Vibrational Frequencies for Isotopologues of Cyclic and Bent Singlet HC2N isomers

Through established, highly-accurate ab initio quartic force fields (QFFs), a complete set of fundamental vibrational frequencies, rotational constants, and rovibrational coupling and centrifugal distortion constants have been determined for both the cyclic 1(sup 1) 1A' and bent 2(sup 1)A' DCCN, H(C13)CCN, HC(C-13)N, and HCC(N-15) isotopologues of HCCN. Spectroscopic constants are computed for all isotopologues using second-order vibrational perturbation theory (VPT2), and the fundamental vibrational frequencies are computed with VPT2 and vibrational configuration interaction (VCI) theory. Agreement between VPT2 and VCI results is quite good with the fundamental vibrational frequencies of the bent isomer isotopologues in accord to within a 0.1 to 3.2 / cm range. Similar accuracies are present for the cyclic isomer isotopologues. The data generated here serve as a reference for astronomical observations of these closed-shell, highly-dipolar molecules using new, high-resolution telescopes and as reference for laboratory studies where isotopic labeling may lead to elucidation of the formation mechanism for the known interstellar molecule: X 3A0 HCCN.

molecular data↗

The Structure, Anharmonic Vibrational Frequencies, and Intensities of NNHNN+

A semi-global potential energy surface (PES) and quartic force field (QFF) based on fitting high-level electronic structure energies are presented to describe the structures and spectroscopic properties of NNHNN+. The equilibrium structure of NNHNN+ is linear with the proton equidistant between the two nitrogen groups and thus of D(sub ∞h) symmetry. Vibrational second-order perturbation theory (VPT2) calculations based on the QFF fails to describe the proton “rattle” motion, i.e., the antisymmetric proton stretch, due to the very flat nature of PES around the global minimum, but performs properly for other modes with sharper potential wells. Vibrational self-consistent field/virtual state configuration interaction (VSCF/VCI) calculations using a version of MULTIMODE without angular momentum terms successfully describe this motion and predict the fundamental to be at 759 cm(exp -1). This is in good agreement with the value of 746 cm(exp -1) from a fixed-node diffusion Monte Carlo calculation and the experimental Ar-tagged result of 743 cm(exp -1). Other VSCF/VCI energies are in good agreement with other experimentally reported ones. Both double-harmonic intensity and rigorous MULTIMODE intensity calculations show the proton transfer fundamental has a very strong intensity.

Astrochemistry↗