Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “effective core potentials”

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

Oscillator strengths of some Ba lines - A treatment including core-valence correlation and relativistic effects

Theoretical calculations of selected excitation energies and oscillator strengths for Ba are presented that overcome the difficulties of previous theoretical treatments. A relativistic effective-core potential treatment is used to account for the relativistic core contraction, but the outermost ten electrons are treated explicitly. Core-valence correlation can be included in this procedure in a rigorous and systematic way through a configuration-interaction calculation. Insight is gained into the importance of relativistic effects by repeating many of the calculations using an all-electron nonrelativistic treatment employing an extended Slater basis set. It is found that the intensity of the intercombination line 3P1-1S0 is accurately determined by accounting for the deviation from LS coupling through spin-orbit mixing with the 1P1 state, and that deviations from the Lande interval rule provide an accurate measure of the degree of mixing.

Bauschlicher, C. W., Jr.↗

Correlation Consistent Basis Sets for In

Correlation consistent basis sets are developed for a 13 electron correlation treatment of In atom. The basis sets are used in all-electron (AE) and effective core potential (ECP) calculations. The atomization energies of InCl and InCl3 are used as a test of the basis set. The results are extrapolated to the complete basis set (CBS) limit. The AE and ECP CCSD(T) CBS values are in very good mutal when the AE values are corrected for scalar relativistic effects. The best computed values are in good agreement with experiment.

Bauschlicher, Charles W., Jr.↗

Theoretical calculation of low-lying states of NaAr and NaXe

Potential curves as well as dipole moments and linking transition moments are calculated for the ground X 2 Sigma + and low lying excited A 2 Pi, B 2 Sigma +, C 2 Sigma +, (4) 2 Sigma +, (2) 2 Pi and (1) 2 Delta states of NaAr and NaXe. Calculations are performed using a self-consistent field plus configuration-interaction procedure with the core electrons replaced by an ab initio effective core potential. The potential curves obtained are found to be considerably less repulsive than the semiempirical curves of Pascale and Vandeplanque (1974) and to agree well with existing experimental data, although the binding energies of those states having potential minima due to van der Waals interactions are underestimated. Emission bands are also calculated for the X 2 Sigma + - C 2 Sigma + excimer transitions of NaAr and NaXe using the calculated transition moments and potential curves, and shown to agree well with experiment on the short-wavelength side of the maximum.

Laskowski, B. C.↗

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

A final report on the extension of the MCSCF/CI method

For the period of July 1982 through January 1983, the research carried out focused on the nature of metal adsorbate bonding. Two main areas were investigated in detail: (1) the adsorption of CO on Al, Na, Mg and Ni; and (2) bonding in transition metal oxides, CuO, CrO and MoO. These theoretical studies have lead to a new understanding of CO-metal bonding as well as to a better understanding of the d involvement in the bonding in transition metal oxides. Four areas associated with cluster model studies are discussed in detail: hybridization and changing in cluster model studies, the adsorption of CO on metals (Al, Na, Mg), the adsorption of CO on Ni, and finally the use of effective core potentials (ECP) in cluster models.

Nelin, C. J.↗

The structure of small metal clusters

One metal atom surrounded by its 12 nearest neighbors is considered for both D(3d) (face-centered cubic-like) and D(3h) (hexagonal close-packed-like) geometries. For Al and Be, the neutral cluster and the positive and negative ions are considered for idealized (all bonds equal) and distorted geometries. The D(3d) geometry is found to be the lowest for Be13, while the D(3h) geometry is lower for Al13. This is the reverse of what is expected based upon the bulk metal structures, Be(hcp) and Al(fcc). Al13 is found to have only small distortions, while Be13 shows large distortions for both the D(3d) and D(3h) geometries. The ions have geometries which are similar to those found for the neutral systems. Both all-electron and effective core potential calculations were carried out on the X13 clusters; the agreement is very good.

Bauschlicher, C. W., Jr.↗

Theoretical spectroscopic parameters for the low-lying states of the second-row transition metal hydrides

A systematic analysis of the low-lying states of all of the second-row transition metal (TM) hydrides except CdH is reported. The calculations included the dominant relativistic contributions through the use of the relativistic effective core potentials of Hay and Wadt (1985). Electron correlation was incorporated, using single-plus-double configuration interaction, the coupled pair functional (CPF) formalism of Ahlrichs et al. (1985), and the Chong and Langhoff (1986) modified version of the CPF method. The spectroscopic parameters D(e), r(e), and mu(e) determined for the low-lying states are compared with the available experimental data and previous theoretical results. In contrast to the first-row TM hydrides studied earlier (Chong et al., 1986), the spectroscopic constants for the second-row TM hydrides were found to be much less sensitive to the level of correlation treatment.

Langhoff, Stephen R.↗

Positive ions of the first- and second-row transition metal hydrides

Theoretical dissociation energies for the first- and second-row transition metal hydride positive ions are critically compared against recent experimental values obtained from ion beam reactive scattering methods. Theoretical spectroscopic parameters and dipole moments are presented for the ground and several low-lying excited states. The calculations employ large Gaussian basis sets and account for electron correlation using the single-reference single- and double-excitation configuration interaction and coupled-pair-functional methods. The Darwin and mass-velocity contributions to the relativistic energy are included in the all-electron calculations on the first-row systems using first-order perturbation theory, and in the second-row systems using the Hay and Wadt relativistic effective core potentials. The theoretical D(0) values for the second-row transition metal hydride positive ions should provide a critical measure of the experimental values, which are not as refined as many of those in the first transition row.

Pettersson, Lars G. M.↗

Small Al clusters. I - The effect of basis set and correlation on the geometry of small Al clusters

A detailed study is presented of the basis set requirements and effects of correlation on the geometry and structure of small Al(n) clusters n = 2, 4, and 13. An effective core potential (ECP) is developed from the Huzinaga basis which accurately reproduces the all-electron results. It is found that depolarization functions are very important in computing the bond length, and that the bond contraction obtained is about as large for a 13-atom cluster as for the Al(4) rhombus structure. With d functions on each center in Al(13), a bond distance shorter than the bulk is obtained, as expected. The inclusion of correlation in Al(4) is found to be less important than the addition of polarization functions for determining the bond length, but increases substantially the calculated atomization energy. These calibration calculations allow an accurate treatment for Al(n) clusters to be designed.

Bauschlicher, Charles W., Jr.↗

All-electron molecular Dirac-Hartree-Fock calculations: The group 4 tetrahydrides CH4, SiH4, GeH4, SnH4 and PbH4

A basis-set-expansion Dirac-Hartree-Fock program for molecules is described. Bond lengths and harmonic frequencies are presented for the ground states of the group 4 tetrahydrides, CH4, SiH4, GeH4, SnH4, and PbH4. The results are compared with relativistic effective core potential (RECP) calculations, first-order perturbation theory (PT) calculations and with experimental data. The bond lengths are well predicted by first-order perturbation theory for all molecules, but non of the RECP's considered provides a consistent prediction. Perturbation theory overestimates the relativistic correction to the harmonic frequencies; the RECP calculations underestimate the correction.

Dyall, Kenneth G.↗

All-electron molecular Dirac-Hartree-Fock calculations - The group IV tetrahydrides CH4, SiH4, GeH4, SnH4, and PbH4

A basis-set-expansion Dirac-Hartree-Fock program for molecules is described. Bond lengths and harmonic frequencies are presented for the ground states of the group 4 tetrahydrides, CH4, SiH4, GeH4, SnH4, and PbH4. The results are compared with relativistic effective core potential (RECP) calculations, first-order perturbation theory (PT) calculations and with experimental data. The bond lengths are well predicted by first-order perturbation theory for all molecules, but none of the RECP's considered provides a consistent prediction. Perturbation theory overestimates the relativistic correction to the harmonic frequencies; the RECP calculations underestimate the correction.

Dyall, Kenneth G.↗

Theoretical investigation of gas-surface interactions

The investigation into the appearance of intruder states from the negative continuum when some of the two-electron integrals were omitted was completed. The work shows that, provided all integrals involving core contracted functions in an atomic general contraction are included, or that the core functions are radially localized, meaningful results are obtained and intruder states do not appear. In the area of program development, the Dirac-Hartree-Fock (DHF) program for closed-shell polyatomic molecules was extended to permit Kramers-restricted open-shell DHF calculations with one electron in an open shell or one hole in a closed shell, or state-averaged DHF calculations over several particle or hole doublet states. One application of the open-shell code was to the KO molecule. Another major area of program development is the transformation of integrals from the scalar basis in which they are generated to the 2-spinor basis employed in parts of the DHF program, and hence to supermatrix form. Particularly concerning the omission of small component integrals, and with increase in availability of disk space, it is now possible to consider transforming the integrals. The use of ordered integrals, either in the scalar basis or in the 2-spinor basis, would considerably speed up the construction of the Fock matrix, and even more so if supermatrices were constructed. A considerable amount of effort was spent on analyzing the integral ordering and tranformation for the DHF program. The work of assessing the reliability of the relativistic effective core potentials (RECPs) was continued with calculation of the group IV monoxides. The perturbation of the metal atom provided by oxygen is expected to be larger than that provided by hydrogen and thus provide a better test of the qualification of the RECPs. Calculations on some platinum hydrides were carried out at nonrelativistic (NR), perturbation theory (PT) and DHF levels. Reprints of four papers describing this work are included.

Dyall, Kenneth G.↗

All-electron molecular Dirac-Hartree-Fock calculations - Properties of the XH4 and XH2 molecules and the reaction energy XH4 yields XH2 + H2, X = Si, Ge, Sn, Pb

Relativistic corrections to a number of properties of the Group IV hydrides are calculated using the Dirac-Hartree-Fock method. The use of first-order perturbation theory is sufficient to obtain relativistic corrections for Ge, but the effects of spin-orbit interaction and other higher-order effects begin to show for Sn and become important for Pb. The energy of the reaction XH4 yields XH2 + H2 (X = Si, Ge, Sn, and Pb) is also calculated. The results are compared with relativistic effective core potential calculations, first-order perturbation theory calculations, and limited experimental data.

Dyall, Kenneth G.↗

An Ab Initio Study of the Low-Lying Doublet States of AgO and AgS

Spectroscopic constants (D(sub o), r(sub e), mu(sub e), T(sub e)) are determined for the doublet states of AgO and AgS below approx. = 30000/cm. Large valence basis sets are employed in conjunction with relativistic effective core potentials (RECPs). Electron correlation is included using the modified coupled-pair functional (MCPF) and multireference configuration interaction (MRCI) methods. The A(sup 2)Sigma(sup +) - X(sup 2)Pi band system is found to occur in the near infrared (approx. = 9000/cm) and to be relatively weak with a radiative lifetime of 900 microns for A(sup 2)Sigma(sup +) (upsilon = 0). The weakly bound C(sup 2)Pi state (our notation), the upper state of the blue system, is found to require high levels of theoretical treatment to determine a quantitatively accurate potential. The red system is assigned as a transition from the C(sup 2)Pi state to the previously unobserved A(sup 2)Sigma(sup +) state. Several additional transitions are identified that should be detectable experimentally. A more limited study is performed for the vertical excitation spectrum of AgS. In addition, a detailed all-electron study of the X(sup 2)Pi and A(sup 2)Sigma(sup +) states of AgO is carried out using large atomic natural orbital (ANO) basis sets. Our best calculated D(sub o) value for AgO is significantly less than the experimental value, which suggests that there may be some systematic error in the experimental determination.

Bauschlicher, Charles W., Jr.↗

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

A new theory of simple classical fluids

The paper presents a unified structural and thermodynamic theory of simple classical fluids in which the interactions between the particles can be represented by spherically symmetric pairwise potentials. Both the excess entropy and the gross features of the radial distribution function are determined mainly by excluded volume effects, which are in turn governed by a hard-core property intrinsic to any pair potential. The potential beyond this effective hard core is considered relatively weak and can be treated as a perturbation. It is also considered essential to sum all subclasses of diagrams to infinite order. The basic form of a diagrammatic scheme which allows both summation and the determination of the structure and excess entropy is presented, and a statement of universality is derived. The statement has been confirmed within the accuracy of present-day computer simulations, and a possible procedure for calculating both the structure and thermodynamics of every physically conceivable pair potential is presented.

Rosenfeld, Y.↗

Gravitational lensing frequencies - Galaxy cross-sections and selection effects

Four issues - (1) the best currently available data on the galaxy velocity-dispersion distribution, (2) the effects of finite core radii potential ellipticity on lensing cross sections, (3) the predicted distribution of lens image separations compared to observational angular resolutions, and (4) the preferential inclusion of lens systems in flux limited samples - are considered in order to facilitate more realistic predictions of multiple image galaxy-quasar lensing frequencies. It is found that (1) the SIS lensing parameter F equals 0.047 +/-0.019 with almost 90 percent contributed by E and S0 galaxies, (2) observed E and S0 core radii are remarkably small, yielding a factor of less than about 2 reduction in total lensing cross sections, (3) 50 percent of galaxy-quasar lenses have image separations greater than about 1.3 arcsec, and (4) amplification bias factors are large and must be carefully taken into account. It is concluded that flat universe models excessively dominated by the cosmological constant are not favored by the small observed galaxy-quasar lensing rate.

Fukugita, Masataka↗