Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “basis sets”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS1996

The Calculation of Accurate Harmonic Frequencies of Large Molecules: The Polycyclic Aromatic Hydrocarbons, a Case Study

The vibrational frequencies and infrared intensities of naphthalene neutral and cation are studied at the self-consistent-field (SCF), second-order Moller-Plesset (MP2), and density functional theory (DFT) levels using a variety of one-particle basis sets. Very accurate frequencies can be obtained at the DFT level in conjunction with large basis sets if they are scaled with two factors, one for the C-H stretches and a second for all other modes. We also find remarkably good agreement at the B3LYP/4-31G level using only one scale factor. Unlike the neutral PAHs where all methods do reasonably well for the intensities, only the DFT results are accurate for the PAH cations. The failure of the SCF and MP2 methods is caused by symmetry breaking and an inability to describe charge delocalization. We present several interesting cases of symmetry breaking in this study. An assessment is made as to whether an ensemble of PAH neutrals or cations could account for the unidentified infrared bands observed in many astronomical sources.

Bauschlicher, Charles W., Jr.↗
NASA NTRS2010

Nonlinear Dynamics of Cantilever-Sample Interactions in Atomic Force Microscopy

The interaction of the cantilever tip of an atomic force microscope (AFM) with the sample surface is obtained by treating the cantilever and sample as independent systems coupled by a nonlinear force acting between the cantilever tip and a volume element of the sample surface. The volume element is subjected to a restoring force from the remainder of the sample that provides dynamical equilibrium for the combined systems. The model accounts for the positions on the cantilever of the cantilever tip, laser probe, and excitation force (if any) via a basis set of set of orthogonal functions that may be generalized to account for arbitrary cantilever shapes. The basis set is extended to include nonlinear cantilever modes. The model leads to a pair of coupled nonlinear differential equations that are solved analytically using a matrix iteration procedure. The effects of oscillatory excitation forces applied either to the cantilever or to the sample surface (or to both) are obtained from the solution set and applied to the to the assessment of phase and amplitude signals generated by various acoustic-atomic force microscope (A-AFM) modalities. The influence of bistable cantilever modes of on AFM signal generation is discussed. The effects on the cantilever-sample surface dynamics of subsurface features embedded in the sample that are perturbed by surface-generated oscillatory excitation forces and carried to the cantilever via wave propagation are accounted by the Bolef-Miller propagating wave model. Expressions pertaining to signal generation and image contrast in A-AFM are obtained and applied to amplitude modulation (intermittent contact) atomic force microscopy and resonant difference-frequency atomic force ultrasonic microscopy (RDF-AFUM). The influence of phase accumulation in A-AFM on image contrast is discussed, as is the effect of hard contact and maximum nonlinearity regimes of A-AFM operation.

Cantrell, John H.↗
NASA NTRS1992

The vibrational frequencies of difluoroethyne

Ab initio coupled-cluster calculations with single and double excitations and with a perturbational treatment of connected triple excitations are reported for difluoroethyne using large basis sets. The results for the transbending mode nu-4 are extremely sensitive to electron correlation and basis set effects. The best theoretical and experimental estimates for the fundamental vibrational frequencies are in excellent agreement.

Breidung, Juergen↗
NASA NTRS1983

On the role of the seasonal cycle

The definition of the seasonal cycle as the projection of an atmospheric time series onto a suitably defined subset of orthogonal basis functions guarantees that any atmospheric covariance may be expressed as the sum of a seasonal cycle part and a transient part (where transient refers to departures from the seasonal cycle rather than the time mean). Results are presented for the seasonal cycle contribution to the zonally averaged fluxes of momentum and heat, and for the zonally averaged height variance, in data for: the winter season, the winter season with the Fourier basis set appropriate for the seven-year time series, and the seven-year mean. The Fourier basis set calculation indicates that the interannual variability of the momentum flux is dominated by the interactions between the seasonal cycle and the meteorologically low frequency flow.

Straus, D. M.↗
NASA NTRS1988

Computed potential energy surfaces for chemical reactions

The minimum energy path for the addition of a hydrogen atom to N2 is characterized in CASSCF/CCI calculations using the (4s3p2d1f/3s2p1d) basis set, with additional single point calculations at the stationary points of the potential energy surface using the (5s4p3d2f/4s3p2d) basis set. These calculations represent the most extensive set of ab initio calculations completed to date, yielding a zero point corrected barrier for HN2 dissociation of approx. 8.5 kcal mol/1. The lifetime of the HN2 species is estimated from the calculated geometries and energetics using both conventional Transition State Theory and a method which utilizes an Eckart barrier to compute one dimensional quantum mechanical tunneling effects. It is concluded that the lifetime of the HN2 species is very short, greatly limiting its role in both termolecular recombination reactions and combustion processes.

Walch, Stephen P.↗
NASA NTRS1987

An initio calculations on C2, Si2, and SiC

Full-configuration-interaction (FCI) calculations have been performed for the Z 1Sigma(g)(+), a 3Pi(u), and b 3Sigma(g)(-) states of C2 and the X 3Sigma(g)(-) and A 3Pi(u) states of Si2 in valence double-zeta plus polarization one-particle basis sets. By carrying out the CASSCF multireference singles plus doubles CI (MRCI) procedure in a very large Gaussian basis set and correcting for errors with the FCI, the triplet-triplet separation in C2 is reproduced to within 47/cm. Analogous calculations on Si2 definitively determine the ground state to be 3Sigma(g)(-); the 3Pi(u) state is only 440 + or - 100/cm higher in energy. For the SiC molecule, which has never been observed spectroscopically but is believed to be an important component in carbon stars and to be present in interstellar space, the A 3Sigma(-) state is nearly 4000/cm above the X 3Pi ground state.

Bauschlicher, Charles W., Jr.↗
NASA NTRS1987

Theoretical studies of the potential surface for the F - H2 greater than HF + H reaction

The F + H2 yields HF + H potential energy hypersurface was studied in the saddle point and entrance channel regions. Using a large (5s 5p 3d 2f 1g/4s 3p 2d) atomic natural orbital basis set, a classical barrier height of 1.86 kcal/mole was obtained at the CASSCF/multireference CI level (MRCI) after correcting for basis set superposition error and including a Davidson correction (+Q) for higher excitations. Based upon an analysis of the computed results, the true classical barrier is estimated to be about 1.4 kcal/mole. The location of the bottleneck on the lowest vibrationally adiabatic potential curve was also computed and the translational energy threshold determined from a one-dimensional tunneling calculation. Using the difference between the calculated and experimental threshold to adjust the classical barrier height on the computed surface yields a classical barrier in the range of 1.0 to 1.5 kcal/mole. Combining the results of the direct estimates of the classical barrier height with the empirical values obtained from the approximation calculations of the dynamical threshold, it is predicted that the true classical barrier height is 1.4 + or - 0.4 kcal/mole. Arguments are presented in favor of including the relatively large +Q correction obtained when nine electrons are correlated at the CASSCF/MRCI level.

Bauschlicher, Charles W., Jr.↗
NASA NTRS1997

Accurate Heats of Formation for SiF(n) and SiF(n+), for N=1-4

Accurate heats of formation are computed for SiFn and SiFn+, for n=1-4. The vibrational frequencies are determined at the B3LYP level of theory. The energetics are determined at the CCSD(T) level of theory. Basis set limit values are obtained by extrapolation. In those cases where the CCSD(T) calculations become prohibitively large, the basis set extrapolation is performed at the MP2 level. The temperature dependence of the heat of formation, heat capacity, and entropy are computed for the temperature range 300 to 4000 K and fit to a polynomial. The CCSD(T) bond energies are compared with those obtained at the B3LYP, MP2, G2, and G2MP2 levels of theory.

Bauschlicher, Charles W., Jr.↗
NASA NTRS1998

Heats of Formation for CF(sub n) (n = 1 - 4), CF(sup +, sub n) (n = 1 - 4), and CF(sup -, sub n) (n = 1 - 3)

Accurate heats of formation are computed for CF(sub n) (n = 1 - 4), CF(sup +, sub n) (n = 1 - 4), and CF(sup -, sub n) (n = 1 - 3). The geometries and vibrational frequencies are determined at the B3LYP level of theory. The energetics are determined at the CCSD(T) level of theory. Basis set limit values are obtained by extrapolation. In those cases where the CCSD(T) calculations become prohibitively large, the basis set extrapolation is performed at the MP2 level. The temperature dependence of the heat of formation, heat capacity, and entropy are computed for the temperature range 300 to 4000 K and fit to a polynomial.

Ricca, Alessandra↗
NASA NTRS1998

Heats of NF(sub n) (n= 1-3) and NF(sub n)(+)(n = 1-3)

Accurate heats of formation are computed for NF(sub n) and NF(sub n)(+), for n = 1-3. The geometries and the vibrational frequencies are determined at the B3LYP level of theory. The energetics are determined at the CCSD(T) level of theory. Basis set limit values are obtained by extrapolation. In those cases where the CCSD(T) calculations become prohibitively large, the basis set extrapolation is performed at the MP2 level. The temperature dependence of the heat of formation, heat capacity, and entropy are computed for the temperature range 300 to 4000 K and fit to a polynomial.

Ricca, Alessandra↗
NASA NTRS1987

Mixed Cu-simple metal dimers and trimers - CuLi, CuLi2, CuNa, CuK, CuBe, CuBe2, Cu2Be, CuAl, and CuAl2

Theoretical studies of selected diatomic and triatomic molecules containing copper and the simple metals Li, Na, K, Be, and Al are presented, with emphasis on elucidating the nature of the bonding in mixed transition metal-simple metal systems. Large Gaussian basis sets are used in the diatomic calculations, and are used to calibrate the triatomic calculations, in which somewhat smaller Gaussian basis sets are employed. Electron correlation is incorporated using both the single-reference singles plus doubles configuration interaction and coupled pair functional methods. It is found that alkali atoms form very polar sigma bonds with copper, and that the ionicity increases with the inclusion of higher excitations because they improve the electron affinity of copper, which in turn allows a larger negative charge on copper. Aluminum is found to form stronger bonds than beryllium, since it does not have to undergo sp hybridization. Some of the trimers bond by forming three-center three-electron bonds. These multicenter bonds are quite strong even when compared to the two-electron bonds in the dimers or to other bonding mechanisms in the trimers.

Bauschlicher, Charles W., Jr.↗
NASA NTRS1987

Accurate ab initio calculations for the ground states of N2, O2 and F2

Spectroscopic constants and dissociation energies for the ground states of N2, O2, and F2 determined at the CAS SCF MRCI correlation level are in excellent agreement with experiment when very large primitive valence and polarization one-particle Gaussian basis sets are employed. The dissociation energy (De) for N2 is larger than experiment unless the 2s electrons are correlated. The basis set requirement for an accurate determination of De is found to increase with the degree of multiple-bond character in the molecule.

Langhoff, Stephen R.↗
NASA NTRS1994

Computed Minimum Energy Pathway for Isomerization in Ketene

A minimum energy pathway for interchange of the CH groups in ketene via a C2v structure has been obtained using complete active space self consistent field (CASSCF) derivative methods with a polarized valence double zeta basis set to define the reaction pathway followed by multi-reference internally contracted configuration interaction (ICCI) calculations with a [3s3p2d/3s2p] basis set to determine the energetics. Qualitatively, the C2v structure is found to be a shallow minimum on the potential energy surface separated from ketene by a small barrier (0.2 kcal/mol), a second minimum, and a larger barrier (3.0 kcal/mol). The minimum energy pathway leading from the C2v minimum to ketene starts by simultaneous rotation of the farther CH group out of the plane and away from the oxygen followed by increase of the CCO angle and subsequent 1,2-migration of the H of the nearer CH group toward the carbon of the farther CH group.

Walch, Stephen P.↗
NASA NTRS1986

Accurate ab initio calculations which demonstrate a 3 Pi u ground state for Al2

The spectroscopic parameters and separations between the three low-lying X 3 Pi u, A 3 Sigma g -, and a 1 Sigma g + states of Al2 are studied as a function of both the one-particle and n-particle basis set. Approximate correlation treatments are calibrated against full Cl calculations correlating the six valence electrons in a double-zeta plus two d-function basis set. Since the CASSCF/MRCI 3 Pi u to 3 Sigma g - separation is in excellent agreement wtih the FCI value, the MRCI calculations were carried out in an extended (20s13p6d4f)/(6s5p3d2f) gaussian basis. Including a small correction for relativistic effects, the best estimate is that 3 Sigma g - state lies 174/cm above the 3 Pi u ground state. The 1 Sigma g + state lies at least 2000/cm higher in energy. At the CPF level, inclusion of 2s and 2p correlation has little effect on D sub e, reduces T sub e by only 26/cm, and shortens the bond lengths by about 0.02 a sub o. Further strong support for a 3 Pi u ground state comes from the experimental absorption spectra, since both observed transitions can be convincingly assigned as 3 Pi u yields 3 Pi g. The (2) 3 Pi g state is observed to be sensitive to the level of correlation treatment, and to have its minimum shifted to shorter rho values, such that the strongest experimental absorption peak probably corresponds to the 0 yields 2 transition.

Bauschlicher, Charles W., Jr.↗
NASA NTRS1986

On the bonding in Fe2(CO)9

The present paper is concerned with an analysis of the bonding in Fe2(CO)9, taking into account the results of self-consistent-field (SCF) calculations. The basis sets used are the same as those employed in the study of Fe(CO)5 conducted by Bauschlicher and Bagus (1985). The Fe basis set starts with the 14s9p5d primitive set optimized by Wachters (1970), which is contracted to 8s4p5d. An SCF wave function, consisting of only the core electrons, was determined, and then electrons were slowly added to the virtual orbitals with the lowest eigenvalues. The conducted calculations show no direct Fe-Fe bonding, but rather bonding resulting from the bridging CO's. The bonding is found to involve both the 5 sigma and 2 pi(asterisk) orbitals of the bridging CO's.

Bauschlicher, C. W., Jr.↗
NASA NTRS2012

Accurate ab initio Quartic Force Fields of Cyclic and Bent HC2N Isomers

Highly correlated ab initio quartic force field (QFFs) are used to calculate the equilibrium structures and predict the spectroscopic parameters of three HC2N isomers. Specifically, the ground state quasilinear triplet and the lowest cyclic and bent singlet isomers are included in the present study. Extensive treatment of correlation effects were included using the singles and doubles coupled-cluster method that includes a perturbational estimate of the effects of connected triple excitations, denoted CCSD(T). Dunning s correlation-consistent basis sets cc-pVXZ, X=3,4,5, were used, and a three-point formula for extrapolation to the one-particle basis set limit was used. Core-correlation and scalar relativistic corrections were also included to yield highly accurate QFFs. The QFFs were used together with second-order perturbation theory (with proper treatment of Fermi resonances) and variational methods to solve the nuclear Schr dinger equation. The quasilinear nature of the triplet isomer is problematic, and it is concluded that a QFF is not adequate to describe properly all of the fundamental vibrational frequencies and spectroscopic constants (though some constants not dependent on the bending motion are well reproduced by perturbation theory). On the other hand, this procedure (a QFF together with either perturbation theory or variational methods) leads to highly accurate fundamental vibrational frequencies and spectroscopic constants for the cyclic and bent singlet isomers of HC2N. All three isomers possess significant dipole moments, 3.05D, 3.06D, and 1.71D, for the quasilinear triplet, the cyclic singlet, and the bent singlet isomers, respectively. It is concluded that the spectroscopic constants determined for the cyclic and bent singlet isomers are the most accurate available, and it is hoped that these will be useful in the interpretation of high-resolution astronomical observations or laboratory experiments.

Inostroza, Natalia↗
NASA NTRS1985

Theoretical transition probabilities between the lowest 2S, 2P and 2D states of Na, K, Rb and Cs

Theoretical transition probabilities between the lowest 2S, 2P and 2D states of the alkali atoms Na through Cs have been computed using near Hartree-Fock quality Slater basis sets. The important core-valence correlation effects are incorporated explicitly by a configuration-interaction procedure. For Cs, the calculations were repeated using a Gaussian basis set so that relativistic effects could be incorporated through an effective core potential procedure. The best calculated electric quadrupole Einstein coefficients are Na(196.3/s), K(103.6/s), Rb(72.4/s) and Cs(19.7/s). Core-valence effects become increasingly important down the column, and reduce the quadrupole transition strengths to about the same degree as for the 2P-2S and 2D-2P dipole-allowed transitions. Relativistic effects increase the quadrupole moment of Cs, but less so than in Ba, presumably because the alkali 2D states are more diffuse.

Langhoff, S. R.↗
NASA NTRS1986

A full CI treatment of Ne atom - A benchmark calculation performed on the NAS CRAY 2

Full CI calculations are performed for Ne atom using Gaussian basis sets of up to triple-zeta plus double polarization quality. The total valence correlation energy through double, triple, quadruple and octuple excitations is compared for eight different basis sets. These results are expected to be an important benchmark for calibrating methods for estimating the importance of higher excitations.

Bauschlicher, C. W., Jr.↗