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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3

All sources NASA NTRS DOE OSTI

NASA NTRS1984

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

NASA NTRS1986

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

NASA NTRS1987

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

NASA NTRS1987

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

NASA NTRS1987

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

NASA NTRS1990

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

NASA NTRS1991

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

NASA NTRS1992

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

NASA NTRS1992

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

NASA NTRS1990

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

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

DOE OSTI2024

ASMS 2024 Investigation of Uranyl Perchlorate Anion Complexes in the Gas Phase via Infrared Multiphoton Dissociation and Collision Induced Dissociation

Investigation of Uranyl Perchlorate Anion Complexes in the Gas Phase via Infrared Multiphoton Dissociation and Collision Induced Dissociation Brittany D. M. Hodges, Christopher A. Zarzana, JungSoo Kim, Jonathan Martens, and W. C. M. Berden Introduction (120 words max) Effects of electronic structure on chemical bonding and reactivity play critical roles shaping the chemical bonding and reactivity behaviors of heavy elements. Understanding the role of f electrons in bond formation between the actinide-series elements like uranium with other ligands is critical for solving technical challenges associated with these heavy elements, important to nuclear fuel cycles, efficient separations of rare earth metals, and understanding the chemistry of stored nuclear fuels and related environmental management sites. In this study, we further examine the interactions between uranyl and the perchlorate ion in order to understand the structures of these ions through the use of IRMPD. Here we report the IRMPD spectra of [UO2(ClO4)3]-, [UO3(ClO4)2]-, and a proposed transition state. Methods (120 word max) IRMPD spectra and CID spectra were acquired using a Bruker amaZon QIT/MS instrument at the Free-Electron Lasers for Infrared eXperiments (FELIX) laboratory at Radboud University. The FELIX QIT/MS is modified to allow for the high-intensity tunable IR beam from FELIX to be directed into the ion packet, resulting in multiphoton dissociation that is measured only when the IR frequency is in resonance with an adequately high absorption vibrational mode of the mass-selected complex. DFT geometry optimizations and frequency calculations using the Gaussian suite of programs were performed using B3LYP, TPSSh, and PBE0 level of theory with 6-31+G(d) basis for the O, C, H, and N atoms and the SDD basis set for U. The SDD basis set employs the Stuttgart/Dresden effective core potential. Preliminary Data or Plenary Speakers Abstract (300 words max) Metal ion clusters of uranyl perchlorate were formed via direct electrospray ionization. For each metal ligand complex of interest, the parent ion was isolated and collision induced dissociation fragmentation and Infrared Multiphoton Dissociation (IRMPD) fragmentation spectra were acquired. Results presented here are the first look at the IRMPD spectra of [UO3(ClO4)2]-, [UO2(ClO4)3]-. Structures were examined using Gaussian at different levels of theory B3LYP level of theory, TPPSh and PBE0 levels, to reflect the behaviors of uranium metal ligand complexes most accurately. In these structures, we identified an overlap between each of these uranyl stretches resulting in their largely being obscured by a perchlorate mode. The CID product spectra agree with similar structures reported by Groenewold for uranyl nitrate in 2006 (10.1021/ja058106n).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Effect of Ionic Core on the Properties of Non-Isothermal Plasmas

This study investigates the effects of ionic cores on non-isothermal plasmas using a novel ion-ion interaction potential that incorporates screening effects from both ion cores and exchange-correlation interactions. Our findings indicate that with increasing distance, the effective potential approaches a Yukawa-like screening potential, while at shorter distances, strong electron binding weakens screening. The different values of the cutoff radius ' cut and core edge steepness α significantly influence the potential behavior and radial distribution functions (RDFs). Higher coupling parameters (Γ i ) strengthen the electron-ion interactions, leading to deeper potential wells and more pronounced non-ideality corrections. Increasing Θ decreases the absolute values of non-ideality corrections, indicating fewer interactions in the system. A larger cutoff radius r cut due to weaker screening effects. As Γ i at a fixed parameter also reduces corrections increases, non-ideality corrections grow, reflecting stronger coupling. The results show the importance of taking into account the ion core effects in dense non-isothermal plasma research.

dense plasmas↗

DOE OSTI2022

Quasi-linear toroidal simulations of resonant magnetic perturbations in eight ITER H-mode scenarios

Abstract Both linear and quasi-linear aspects of the plasma response to the resonant magnetic perturbation (RMP) field are numerically investigated for various H-mode scenarios in ITER, covering the pre-fusion power operation and the fusion power operation phases. Linear response computations for eight ITER scenarios, with varying plasma current and toroidal magnetic field, reveal that the best coil current phasing for controlling the type-I edge localized modes (ELMs) scales roughly linearly with the edge safety factor. The coil phasing is defined as the relative toroidal phase of the coil currents between different rows, for a given toroidal harmonic. Quasi-linear initial value simulation, which is the focus of the present study, shows that application of the n = 3 ( n is the toroidal mode number) RMP field has a minimum side effect on the plasma core momentum confinement but potentially a large effect on the global particle transport. Generally, the RMP field with the best (worst) coil phasing for ELM control produces the strongest (weakest) effect on the plasma edge flow and the overall density. This robustly holds for all eight ITER scenarios. Consequently, in order to minimize the RMP induced side effects while achieving ELM control (suppression) in ITER, a compromise is necessary in choosing the coil current configuration.

Physics↗

NASA NTRS1979

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

NASA NTRS1991

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↗

DOE OSTI2024

Ab initio leading order effective potential for elastic proton scattering based on the symmetry-adapted no-core shell model

Calculating microscopic optical potentials for elastic scattering at intermediate energies from light nuclei in an ab initio fashion within the Watson expansion has been established within the last few years. Based on the Watson expansion of the multiple scattering series, we employ a nonlocal translationally invariant nuclear density derived within the symmetry-adapted no-core shell model (SA-NCSM) framework from a chiral next-to-next-to-leading order (NNLO) nucleon-nucleon interaction and the very same interaction for a consistent full-folding calculation of the effective (optical) potential for nucleon-nucleus scattering for medium-heavy nuclei. The leading order effective (optical) folding potential is computed by integrating over a translationally invariant SA-NCSM one-body scalar density, spin-projected momentum distribution, and the Wolfenstein amplitudes 𝐴, 𝐶, and 𝑀. The resulting nonlocal potentials serve as input for a momentum space Lippmann-Schwinger equation. In the SA-NCSM, the model space is systematically up-selected using Sp⁡(3,ℝ) symmetry considerations. For the light nucleus of 6 He, we establish a systematic selection scheme in the SA-NCSM for scattering observables. Then, we apply this scheme to calculations of scattering observables, such as differential cross sections, analyzing powers, and spin rotation functions for elastic proton scattering from 20 Ne and 40 Ca in the energy regime between 65 and 200 MeV, and compare to available data. Furthermore, our calculations show that the leading order effective nucleon-nucleus potential in the Watson expansion of multiple scattering theory obtained from an up-selected SA-NCSM model space describes 40 Ca elastic scattering observables reasonably well to about 60 degrees in the center-of-mass frame, which coincides roughly with the validity of the NNLO chiral interaction used to calculate both the nucleon-nucleon amplitudes and the one-body scalar and spin nuclear densities.

Ab initio calculations↗

NASA NTRS1994

Vertical velocity in oceanic convection off tropical Australia

Time series of 1-Hz vertical velocity data collected during aircraft penetrations of oceanic cumulonimbus clouds over the western Pacific warm pool as part of the Equatorial Mesoscale Experiment (EMEX) are analyzed for updraft and downdraft events called cores. An updraft core is defined as occurring whenever the vertical velocity exceeds 1 m/sec for at least 500 m. A downdraft core is defined analogously. Over 19,000 km of straight and level flight legs are used in the analysis. Five hundred eleven updraft cores and 253 downdraft cores are included in the dataset. Core properties are summarized as distributions of average and maximum vertical velocity, diameter, and mass flux in four altitude intervals between 0.2 and 5.8 km. Distributions are approximately lognormal at all levels. Examination of the variation of the statistics with height suggests a maximum in vertical velocity between 2 and 3 km; slightly lower or equal vertical velocity is indicated at 5 km. Near the freezing level, virtual temperature deviations are found to be slightly positive for both updraft and downdraft cores. The excess in updraft cores is much smaller than that predicted by parcel theory. Comparisons with other studies that use the same analysis technique reveal that EMEX cores have approximately the same strength as cores of other oceanic areas, despite warmer sea surface temperatures. Diameter and mass flux are greater than those in the Global Atmospheric Research Program (GATE) but smaller than those in hurricane rainbands. Oceanic cores are much weaker and appear to be slightly smaller than those observed over land during the Thunderstorm Project. The markedly weaker oceanic vertical velocities below 5.8 km (compared to the continental cores) cannot be attributed to smaller total convective available potential energy or to very high water loading. Rather, it is suggested that water loading, although less than adiabatic, is more effective in reducing buoyancy of oceanic cores because of the smaller potential buoyancy below 5.8 km. Entrainment appears to be more effective in reducing buoyancy to well below adiabatic values in oceanic cores, a result consistent with the smaller oceanic core diameters in the lower cloud layer. It is speculated further that core diameters are related to boundary layer depth, which is clearly smaller over the oceans.

Lucas, Christopher↗