Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Moller-Plesset perturbation 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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Second-Order Moller-Plesset Perturbation Theory for Molecular Dirac-Hartree-Fock Wave Functions

Moller-Plesset perturbation theory is developed to second order for a selection of Kramers restricted Dirac-Hartree-Fock closed and open-shell reference wave functions. The open-shell wave functions considered are limited to those with no more than two electrons in open shells, but include the case of a two-configuration SCF reference. Denominator shifts are included in the style of Davidson's OPT2 method. An implementation which uses unordered integrals with labels is presented, and results are given for a few test cases.

Dyall, Kenneth G.↗
NASA NTRS2002

Theoretical Study of Infrared and Raman Spectra of Hydrated Magnesium Sulfate Salts

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.

Chaban, Galina M.↗
NASA NTRS1998

C2 Fragmentation Energy of C60 Revisited: Theory Disagrees with Most Experiments

Following our earlier work on the subject, we have carried out density functional theory (DFT) and second-order Moller-Plesset perturbation theory (MP2) calculations of the dissociation energy of the reaction C60 yields C58 + C2 using polarized basis sets and geometries optimized with DFT methods. The present theoretical results support an electronic fragmentation energy D(sub e) around 10-11 eV in disagreement with most experimental results that place the dissociation energy D(sub o) (including zero point energy) around 7-8 eV. The plausible errors remaining in the theoretical calculations are unlikely to account for this big difference (2-4 eV).

Boese, A. Daniel↗
NASA NTRS1995

A Comparison of the Accuracy of Different Functionals

The errors in the computed geometries, zero-point energies, and atomization energies of molecules containing only first and second row atoms are compared for several levels of theory, including Hartree-Fock, second-order Moller-Plesset perturbation theory (MP2), and density functional theory (DFT) using five different functionals, including two hybrid functionals. The 6-31G* and 6-311+G(3df,2p) basis set are used. Overall, the B3LYP hybrid functional yields the best results.

Bauschlicher, Charles W., Jr.↗
NASA NTRS1997

The Vibrational Frequencies of CaO2, ScO2, and TiO2: A Comparison of Theoretical Methods

The vibrational frequencies of several states of CaO2, ScO2, and TiO2 are computed at using density functional theory (DFT), the Hatree-Fock approach, second order Moller-Plesset perturbation theory (MP2), and the complete-active-space self-consistent-field theory. Three different functionals are used in the DFT calculations, including two hybrid functionals. The coupled cluster singles and doubles approach including the effect of unlinked triples, determined using perturbation theory, is applied to selected states. The Becke-Perdew 86 functional appears to be the cost effective method of choice, although even this functional does not perform well for one state of CaO2. The MP2 approach is significantly inferior to the DFT approaches.

Rosi, Marzio↗
NASA NTRS1991

A parallel vectorized implementation of triple excitations in CCSD(T) - Application to the binding energies of the AlH3, AlH2F, AlHF2 and AlF3 dimers

An efficient method for various noniterative estimates of connected triple excitations in coupled-cluster theory is outlined and related to a similar expression occurring in Moller-Plesset perturbation theory. The method is highly vectorized and capable of utilizing multiple processors on a shared-memory machine, leading to computational rates in excess of one billion floating-point operations per second on four processors of a CRAY Y-MP. Using the new procedure, the binding energies of the D(2h) diborane-type dimers of AlH3, AlH2F, AlHF2, and AlF3 have been determined to be 32, 40, 20, and 47 kcal/mol, respectively. For Al2F6, the correlation procedure includes 232 molecular orbitals and over 1.5 x 10 to the 6th single and double coupled-cluster amplitudes, effectively accounting for over 2 x 10 to the 9th connected triple excitations.

Rendell, Alistair P.↗
NASA NTRS1996

An ab initio benchmark study of the H + CO --> HCO reaction

The H + CO --> HCO reaction has been characterized with correlation consistent basis sets at five levels of theory in order to benchmark the sensitivities of the barrier height and reaction ergicity to the one-electron and n-electron expansions of the electronic wave function. Single and multireference methods are compared and contrasted. The coupled cluster method RCCSD(T) was found to be in very good agreement with Davidson-corrected internally-contracted multireference configuration interaction (MRCI+Q). Second-order Moller-Plesset perturbation theory (MP2) was also employed. The estimated complete basis set (CBS) limits for the barrier height (in kcal/mol) for the five methods, including harmonic zero-point energy corrections, are MP2, 4.66; RCCSD, 4.78; RCCSD(T), 4.15; MRCI, 5.10; and MRCI+Q, 4.07. Similarly, the estimated CBS limits for the ergicity of the reaction are: MP2, -17.99; RCCSD, -13.34; RCCSD(T), -13.79; MRCI, -11.46; and MRCI+Q, -13.70. Additional basis set explorations for the RCCSD(T) method demonstrate that aug-cc-pVTZ sets, even with some functions removed, are sufficient to reproduce the CBS limits to within 0.1-0.3 kcal/mol.

NASA Discipline Exobiology↗
NASA NTRS1989

A high level Ab initio study of the anionic hydrogen-bonded complexes FH-CN-, FH-NC-, H2O-CN- and H2O-NC-

HF, H2O, CN- and their hydrogen-bonded complexes were studied using state-of-the-art ab initio quantum mechanical methods. A large Gaussian one particle basis set consisting of triple zeta plus double polarization plus diffuse s and p functions (TZ2P + diffuse) was used. The theoretical methods employed include self consistent field, second order Moller-Plesset perturbation theory, singles and doubles configuration interaction theory and the singles and doubles coupled cluster approach. The FH-CN- and FH-NC- and H2O-CN-, H2O-NC- pairs of complexes are found to be essentially isoenergetic. The first pair of complexes are predicted to be bound by approx. 24 kcal/mole and the latter pair bound by approximately 15 kcal/mole. The ab initio binding energies are in good agreement with the experimental values. The two being shorter than the analogous C-N hydrogen bond. The infrared (IR) spectra of the two pairs of complexes are also very similar, though a severe perturbation of the potential energy surface by proton exchange means that the accurate prediction of the band center of the most intense IR mode requires a high level of electronic structure theory as well as a complete treatment of anharmonic effects. The bonding of anionic hydrogen-bonded complexes is discussed and contrasted with that of neutral hydrogen-bonded complexes.

Lee, Timothy J.↗
NASA NTRS1989

A high-level ab initio study of the anionic hydrogen-bonded complexes FH-CN(-), FH-NC(-), H2O-CN(-), and H2O-NC(-)

HF, H2O, CN- and their hydrogen-bonded complexes were studied using state-of-the-art ab initio quantum mechanical methods. A large Gaussian one particle basis set consisting of triple zeta plus double polarization plus diffuse s and p functions (TZ2P + diffuse) was used. The theoretical methods employed include self consistent field, second order Moller-Plesset perturbation theory, singles and doubles configuration interaction theory and the singles and doubles coupled cluster approach. The FH-CN- and FH-NC- and H2O-CN-, H2O-NC- pairs of complexes are found to be essentially isoenergetic. The first pair of complexes are predicted to be bound by approx. 24 kcal/mole and the latter pair bound by approximately 15 kcal/mole. The ab initio binding energies are in good agreement with the experimental values. The two being shorter than the analogous C-N hydrogen bond. The infrared (IR) spectra of the two pairs of complexes are also very similar, though a severe perturbation of the potential energy surface by proton exchange means that the accurate prediction of the band center of the most intense IR mode requires a high level of electronic structure theory as well as a complete treatment of anharmonic effects. The bonding of anionic hydrogen-bonded complexes is discussed and contrasted with that of neutral hydrogen-bonded complexes.

Lee, Timothy J.↗
NASA NTRS1991

Analytic gradients for coupled-cluster energies that include noniterative connected triple excitations - Application to cis- and trans-HONO

An efficient formulation of the analytic energy gradient for the single and double excitation coupled-cluster method that includes a perturbational estimate of the effects of connected triple excitations is presented. The formulation has a small computational cost, and the algebraic manipulations may be applied generally to the analytic gradient of Moller-Plesset perturbation theory energies. The new formulation has been implemented in an efficient set of programs that utilize highly vectorized algorithms and has been used to investigate the equilibrium structures, harmonic vibrational frequencies, IR intensities, and energy separation of cis- and trans-HONO.

Lee, Timothy J.↗
NASA NTRS1992

Ab initio study of the molecular structure and vibrational spectrum of nitric acid and its protonated forms

The equilibrium structures, harmonic vibrational frequencies, IR intensities, and relative energetics of HNO3 and its protonated form H2NO3+ were investigated using double-zeta plus polarization and triple-zeta plus polarization basis sets in conjunction with high-level ab initio methods. The latter include second-order Moller-Plesset perturbation theory, the single and double excitation coupled cluster (CCSD) methods, a perturbational estimate of the effects of connected triple excitations (CCSD(T)), and the self-consistent field. To determine accurate energy differences CCSD(T) energies were computed using large atomic natural orbital basis sets. Four different isomers of H2NO3+ were considered. The lowest energy form of protonated nitric acid was found to correspond to a complex between H2O and NO2+, which is consistent with earlier theoretical and experimental studies.

Lee, Timothy J.↗
NASA NTRS1993

Theoretical Study of the H2-ML(+) Binding Energies

The cooperative ligand effects are studied in MLH2(+) and the results are compared to the recent experiments of Kemper et al. The bonding in these compounds is principally electrostatic in origin; however, ligand to metal and metal to ligand donations are important, especially for H2. We show that differences arise among the vanadium, cobalt, and copper complexes which are due to 3d donation to H2. Electron correlation is required to describe the dative interaction, and we find that second order Moller-Plesset perturbation theory (MP2) yields a good description of these systems compared with higher levels of correlation (such as the modified coupled pair functional and coupled cluster approaches) and experiment. However, obtaining quantitative results requires higher levels of theory than MP2.

Maitre, Philippe↗
NASA NTRS1979

Finite-field Calculations of Molecular Polarizabilities Using Field-Induced Polarization Functions: Second- and Third-order Perturbation Correlation Corrections to the Coupled Hartree-Fock Polarizability of H2O

Ordinary Rayleigh-Schroudinger perturbation theory with Moller-Plesset (RSMP) partitioning is used to calculate second- and third-order correlation corrections to the CHF polarizability and dipole moment of the water molecule by a finite-field procedure. Pade approximants are found to be useful in accelerating the convergence of the property perturbation expansions. Field-induced polarization functions suitable for polarizability calculations are determined. The average polarizability calculated, neglecting vibrational averaging, with Dunning's (9s5p/4s-4s2p/2s) contracted GTO basis set augmented by field-induced lslp2d/lp polarization functions is within 3 per cent of the experimental result. Correlation corrections to the dipole moment and polarizability of the water molecule calculated by the finite-field RSMP and single + double excitation CI(SDCI) methods for the same basis set are found to be in close agreement. The RSMP approach has the advantages of being size-consistent and of being capable of greater efficiency than the SCDI method. Comparative calculations carried out using Epstein-Nesbet partitioning show that through third order RSEN correlation perturbation expansions for the dipole moment and polarizability are less rapidly convergent than RSMP expansions. However, reasonable accord with RSMP results can be achieved by using Pade approximants to accelerate the convergence of RSEN energy perturbation expansions. The convergence of RSMP property correlation expansions based on the zeroth-order uncoupled-Hartree-Fock (UCHF) and coupled-Hartree-Fock (CHF) approximations are compared through third order. Whereas the CHF + RSMP expansions are for practical purposes fully converged, the UCHF + RSMP expansions are not adequately converged.

Langhoff, S. R.↗
DOE OSTI2024

A static quantum embedding scheme based on coupled cluster theory

Here, we develop a static quantum embedding scheme that utilizes different levels of approximations to coupled cluster (CC) theory for an active fragment region and its environment. To reduce the computational cost, we solve the local fragment problem using a high-level CC method and address the environment problem with a lower-level Møller–Plesset (MP) perturbative method. This embedding approach inherits many conceptual developments from the hybrid second-order Møller–Plesset (MP2) and CC works by Nooijen [J. Chem. Phys. 111, 10815 (1999)] and Bochevarov and Sherrill [J. Chem. Phys. 122, 234110 (2005)]. We go beyond those works here by primarily targeting a specific localized fragment of a molecule and also introducing an alternative mechanism to relax the environment within this framework. We will call this approach MP-CC. We demonstrate the effectiveness of MP-CC on several potential energy curves and a set of thermochemical reaction energies, using CC with singles and doubles as the fragment solver, and MP2-like treatments of the environment. The results are substantially improved by the inclusion of orbital relaxation in the environment. Using localized bonds as the active fragment, we also report results for N=N bond breaking in azomethane and for the central C–C bond torsion in butadiene. We find that when the fragment Hilbert space size remains fixed (e.g., when determined by an intrinsic atomic orbital approach), the method achieves comparable accuracy with both a small and a large basis set. Additionally, our results indicate that increasing the fragment Hilbert space size systematically enhances the accuracy of observables, approaching the precision of the full CC solver.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

The approximate second order coupled-cluster method based on a size-consistent Brillouin–Wigner partitioning

We present a variant of the approximate second order coupled-cluster method (CC2) with a two-parameter size-consistent Brillouin–Wigner (BW-s) partitioning instead of a Møller–Plesset (MP) partitioning for the unperturbed Hamiltonian, which we refer to as BWs-CC2. The computational complexity of this model scales identically to CC2 with molecular size. Conventional CC2 and its regularized BWs-CC2 variants, as well as conventional MP2 and two of its regularized BW-s2 variants, were assessed on a 535 element database spanning thermochemistry, non-covalent interactions, barrier heights, and isomerization energies. To ensure a well-defined model chemistry, the assessment was performed using internally stable spin-polarized Hartree–Fock (HF) orbitals in the finite aug-cc-pVQZ basis without counterpoise corrections. As a result of using stable orbitals, contrary to conventional wisdom, we find that CC2 substantially outperforms MP2 on molecules with significantly spin contaminated reference orbitals without a significant increase in error on systems with a spin-pure reference, showing the value of its single substitutions. While no single choice of regularization parameters can be optimal for all datasets, we find that BWs-CC2 generally outperforms both CC2 and BW-s2 with a single judicious parameter choice. Additional tests on dipole moments and bond lengths of diatomics provide further support for the utility of this choice. Furthermore, the main outliers and poorest performing cases are associated with large amounts of spin-contamination in the HF reference, which is indicative of systems with either strong correlation or extensive artificial symmetry breaking. Overall, these findings argue that the perception of the quality of the CC2 ground state should be reevaluated and that it can be further improved upon by the soundly based BWs-CC2 variant with the recommended parameter choice.

Correlation energy↗
NASA NTRS1994

Relativistic and Correlation Effects in CuH, AgH and AuH: Comparison of Various Relativistic Methods

The effects of relativity on the bond lengths, dissociation energies, and harmonic vibrational frequencies of the 1Epsilon(+) electronic ground states of the group IB hydrides CuH, AgH and AuH have been evaluated with a variety of ab initio methods. These properties were investigated with moderately-sized basis sets at the self-consistent field Hartree Fock (SCF HF) level and with second-order Moller-Plesset (MP2) perturbation theory for electron correlation. Comparisons were made between all-electron results using the nonrelativistic Hamiltonian, perturbation theory (PT) at first-order with only the one-electron non-fine structure terms of the Breit-Pauli Hamiltonian, the spin-free Douglas-Kroll (DK) transformed Dirac Hamiltonian and the untransformed Dirac Hamiltonian, and results using two sets of relativistic effective core potentials (RECPs). The expected trends of bond length decrease, dissociation energy increase and harmonic frequency increase with both relativity and correlation are found. Both sets of RECPs are shown to give good results, if accompanied by a reasonable basis set. The DK method is demonstrated to be an inexpensive, reliable approximation to the DHF method.

Collins, Charlene L.↗
NASA NTRS1994

Heat of formation determination of the ground and excited state of cyanomethylene (HCCN) radical

Ab initio electronic structure theory has been used to characterize the structure of the ground triplet and lowest singlet excited states of cyanomethylene. The geometries, vibrational frequencies, and heats of formation have been determined using second-order Moller-Plesset perturbation, single and double excitation configuration interaction, and quadratic configuration interaction theory. The heat of formation is predicted with isodesmic reaction and Gaussian-2 theory (G2) for the ground triplet and first excited singlet states of cyanomethylene. For the ground state Delta-H(sub 0)(sup f,0) is 114.8+/-2 kcal/mol while for the excited single state it is 126.5+/-2 kcal/mol.

Francisco, Joseph S.↗
NASA NTRS1994

On The Choice of a Zeroth-Order Hamiltonian for Second-Order Perturbation Theory with A CASSCF Reference Function

A new approach to perturbation theory based on a CASSCF reference function has been developed. The key to the approach is the definition of the zeroth order Hamiltonian, H(sub 0), which includes the full CI Hamiltonian for the active space. In the the inactive and secondary spaces, operators may be chosen which reduce to the usual Moller-Plesset or Epstein-Nesbet forms in the limit of a null active space. These operators are diagonal in the orbital indices and permit the block-diagonalization of H(sub 0). The reference is an eigenfunction of Ho without N-particle projection. H(sub 0) automatically incorporates denominator shifts in the style of those appearing in recent open-shell perturbation theories. Comparative results are presented for a few test cases.

Dyall, Kenneth G.↗