Coupled multiconfigurational self-consistent-field method of atomic dipole polarizabilities. I - Theory and application to carbon.
Explore the source record for details and available documents.
Billingsley, F. P., II↗
SEARCH · Engineering Papers
Results for “multiconfigurational methods”
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.
Plasma-environment effects on K lines of astrophysical interest III. IPs, K thresholds, radiative rates, and Auger widths in Fe ix – Fe xvi
Aims. In the context of black-hole accretion disks, we aim to compute the plasma-environment e ects on the atomic parameters used to model the decay of K-vacancy states in moderately charged iron ions, namely Fe ix – Fe xvi. Methods. We used the fully relativistic multiconfiguration Dirac–Fock method approximating the plasma electron–nucleus and electron–electron screenings with a time-averaged Debye–Hückel potential. Results. We report modified ionization potentials, K-threshold energies, wavelengths, radiative emission rates, and Auger widths for plasmas characterized by electron temperatures and densities in the ranges 10^(5)-10^(7) K and 10^(18)-10^(22)/cu. cm. Conclusions. This study confirms that the high-resolution X-ray spectrometers onboard the future XRISM and Athena space missions will be capable of detecting the lowering of the K edges of these ions due to the extreme plasma conditions occurring in accretion disks around compact objects.
J. Deprince↗
Plasma environment effects on K lines of astrophysical interest IV. IPs, K thresholds, radiative rates, and Auger widths in Fe ii–Fe viii
Aims. Within the framework of compact-object accretion disks, we calculate plasma environment e ects on the atomic structure and decay parameters used in the modeling of K lines in lowly charged iron ions, namely Fe ii–Fe viii. Methods. For this study, we used the fully relativistic multiconfiguration Dirac–Fock method approximating the plasma electron–nucleus and electron-electron screenings with a time-averaged Debye-Hückel potential. Results. We report modified ionization potentials, K-threshold energies, wavelengths, radiative emission rates, and Auger widths for plasmas characterized by electron temperatures and densities in the ranges 10(exp 5)-10(exp 7) K and 10(exp 18)-10(exp 22) per cu.cm. In addition, we propose two universal fitting formulae to predict the IP and K-threshold lowerings in any elemental ion. Conclusions. We conclude that the high-resolution X-ray spectrometers onboard the future XRISM and ATHENA space missions will be able to detect the lowering of the K edges of these Fe ions due to the extreme plasma conditions occurring in the accretion disks around compact objects.
J. Deprince↗
The Fraternal Twins of Quartet O4(+)
Eleven stationary geometries of quartet O4(+) have been studied by ab initio methods. The geometries were optimized at the Complete Active Space Self-Consistent Field (CASSCF) level of theory and the energies were calculated by the multiconfigurational second order pertubation method (CASPT2), using Double-Zeta Plus polarization (DZP), Triple-Zeta Plus Double Polarization (TZ2P), average Atomic Natural Orbital (ANO) (5s4p2d) and average ANO (6s5p3d2f) basis sets. The rectangular and trans-planar structures are found to be the most stable, with an energy barrier to conversion between the two at the threshold of dissociation. Both have a delocalized hole and are stable relative to separated 02 and 02(+)by 11.0 and 11.5 kcal/mol for the rectangular and the 2 trans-planar structure, respectively, compared with the experimentally deduced energy in the range of 9.2 to 10.8 kcal/mol. The adiabatic ionization potentials of 04 and 02 are computed to be 11.67 and 12.21 eV, while experimental values are 11.66 and 12.07 eV, respectively. The vibrational frequencies have been computed for all degrees of freedom at the CASSCF level of theory. Symmetry breaking is found to be a particular problem in the computation of the antisymmetric stretch frequency for the delocalized structures at the CASSCF level of theory. Attempts to rectify these problems using the restricted active space self-consistent field (RASSCF) method leads to additional difficulties, but further analysis yields insight into the symmetry breaking and problems with earlier calculations. Finally, a nonorthogonal configuration interaction (CI) calculation based on the interaction of localized CASSCF wave functions using the Complete Active Space State Interation (CASSI) method leads to a balanced treatment of the antisymmetric stretch which is free from symmetry breaking. The study explains the four most prominent absorption frequencies observed in the partially unassigned IR spectrum of O4(+) isolated in solid neon as the antisymmetric OO stretch, and the combination band of the symmetric and antisymmetric 00 stretch of both the rectangular and trans-planar structures.
Lindh, Roland↗
Plasma environment effects on K lines of astrophysical interest I. Atomic structure, radiative rates, and Auger widths of oxygen ions
Aims. In the context of black-hole accretion disks, the main goal of the present study is to estimate the plasma environment effects on the atomic structure and radiative parameters associated with the K-vacancy states in ions of the oxygen isonuclear sequence. Methods. We used a time-averaged Debye–Hückel potential for both the electron–nucleus and the electron–electron interactions implemented in the fully relativistic multiconfiguration Dirac–Fock (MCDF) method. Results. Modified ionization potentials, K thresholds, Auger widths, and radiative transition wavelengths and rates are reported for O I–O VII in plasma environments with electron temperature and density ranges 10(exp 5)−10(exp 7) K and 10(exp 18)−10(exp 22)/cu.cm.
J. Deprince↗
Plasma environment effects on K lines of astrophysical interest II. Ionization potentials, K thresholds, radiative rates, and Auger widths in Ne- through He-like iron ions (Fe xvii – Fe xxv)
Aims. In the context of accretion disks around black holes, we estimate plasma-environment effects on the atomic parameters associated with the decay of K-vacancy states in highly charged iron ions, namely Fe XVII – Fe XXV. Methods. Within the relativistic multiconfiguration Dirac–Fock (MCDF) framework, the electron–nucleus and electron–electron plasma screenings were approximated with a time-averaged Debye–Hückel potential. Results. Modified ionization potentials, K thresholds, wavelengths, radiative emission rates, and Auger widths are reported for astrophysical plasmas characterized by electron temperatures and densities in the ranges 10^5 − 10^7 K and 10^18 − 10^22 per cu.cm, respectively. Conclusions. We conclude that the high-resolution microcalorimeters on board future X-ray missions such as XRISM and ATHENA are expected to be sensitive to the lowering of the iron K edge due to the extreme plasma conditions occurring in accretion disks around compact objects.
J. Deprince↗
Multiconfiguration self-consistent-field calculation of the dipole moment function of CO/X 1 sigma +/
Using the optimized valence configurations (OVC) multiconfiguration self-consistent-field (MCSCF) method, the dipole moment function for the ground state of CO in the vicinity of the equilibrium internuclear distance has been calculated. The OVC MCSCF calculation results are compared with existing Hartree-Fock and configuration interaction treatments of this molecule at single points and also the dipole moment function deduced from experimental infrared intensities. A general prescription for constructing OVC wavefunctions for diatomic molecules is also presented.
Billingsley, F. P., II↗
Dissociative Ionization and Product Distributions of Benzene and Pyridine by Electron Impact
We report a theoretical study of the dissociative ionization (DI) and product distributions of benzene (C6H6) and pyridine (C5H5N) from their low-lying ionization channels. Our approach makes use of the fact that electronic motion is much faster than nuclear motion allowing DI to be treated as a two-step process. The first step is the electron-impact ionization resulting in an ion with the same nuclear geometry as the neutral molecule. In the second step, the nuclei relax from the initial geometry and undergo unimolecular dissociation. For the ionization process we use the improved binary-encounter dipole (iBED) model [W.M. Huo, Phys. Rev. A64,042719-I (2001)]. For the unimolecular dissociation, we use multiconfigurational self-consistent field (MCSCF) methods to determine the steepest descent pathways to the possible product channels. More accurate methods are then used to obtain better energetics of the paths which are used to determine unimolecular dissociation probabilities and product distributions. Our analysis of the dissociation products and the thresholds of their productions for benzene are compared with the recent dissociative photoionization meausurements of benzene by Feng et al. [R. Feng, G. Cooper, C.E. Brion, J. Electron Spectrosc. Relat. Phenom. 123,211 (2002)] and the dissociative photoionization measurements of pyridine by Tixier et al. [S. Tixier, G. Cooper, R. Feng, C.E. Brion, J. Electron Spectrosc. Relat. Phenom. 123,185 (2002)] using dipole (e,e+ion) coincidence spectroscopy.
Dateo, Christopher E.↗
Design Of Combined Stochastic Feedforward/Feedback Control
Methodology accommodates variety of control structures and design techniques. In methodology for combined stochastic feedforward/feedback control, main objectives of feedforward and feedback control laws seen clearly. Inclusion of error-integral feedback, dynamic compensation, rate-command control structure, and like integral element of methodology. Another advantage of methodology flexibility to develop variety of techniques for design of feedback control with arbitrary structures to obtain feedback controller: includes stochastic output feedback, multiconfiguration control, decentralized control, or frequency and classical control methods. Control modes of system include capture and tracking of localizer and glideslope, crab, decrab, and flare. By use of recommended incremental implementation, control laws simulated on digital computer and connected with nonlinear digital simulation of aircraft and its systems.
Halyo, Nesim↗
Orbital configuration assignment of 4S/e/ doubly excited states of lithium isoelectronic sequence
Results of Hartree-Fock and multiconfiguration Hartree-Fock calculations using only two configurations (degenerate in zero order within the context of the nuclear-charge-expansion perturbation method) for 4S(e) doubly excited states of lithium isoelectronic sequence are used to propose an orbital configuration classification of these states extending the similar classification by Holoien and Geltman for 4S(e) states of lithium only.
Ali, M. A.↗
Integral processing in beyond-Hartree-Fock calculations
The increasing rate at which improvements in processing capacity outstrip improvements in input/output performance of large computers has led to recent attempts to bypass generation of a disk-based integral file. The direct self-consistent field (SCF) method of Almlof and co-workers represents a very successful implementation of this approach. This paper is concerned with the extension of this general approach to configuration interaction (CI) and multiconfiguration-self-consistent field (MCSCF) calculations. After a discussion of the particular types of molecular orbital (MO) integrals for which -- at least for most current generation machines -- disk-based storage seems unavoidable, it is shown how all the necessary integrals can be obtained as matrix elements of Coulomb and exchange operators that can be calculated using a direct approach. Computational implementations of such a scheme are discussed.
Taylor, P. R.↗
Radiative lifetimes and transition probabilities for electric-dipole delta n equals zero transitions in highly stripped sulfur ions
The beam-foil time-of-flight method has been used to investigate radiative lifetimes and transition rates involving allowed intrashell transitions within the L shell of highly ionized sulfur. The results for these transitions, which can be particularly correlation-sensitive, are compared with current calculations based upon multiconfigurational models.
Pegg, D. J.↗
FUTURE LEVELS OF AEROSOLS & RADIATIVE FORCING OVER THE ARCTIC
The Arctic Monitoring and Assessment Programme (AMAP) is currently assessing the impacts of Short-Lived Climate Forcers (SLCF) on Arctic climate and air quality. In support of the assessment, we used the NASA Goddard Institute of Space Sciences (GISS) Earth System Model (modelE2.1), with prescribed sea surface temperature and sea-ice fraction, to simulate SLCF concentrations globally between 1995 and 2015. Two simulations were conducted, using the One-Moment Aerosol (OMA) and the Multiconfiguration Aerosol TRacker of mIXing state (MATRIX) aerosol modules. OMA is a mass-based scheme in which aerosols are assumed to remain externally mixed and have a prescribed and constant size distribution, while MATRIX is an aerosol microphysics scheme based on the quadrature method of moments, which is able to explicitly simulate the mixing state of aerosols. Anthropogenic emissions from the ECLIPSE v6b emissions database were used, along with emissions from aircrafts and open biomass burning from the Coupled Model Intercomparison Project Phase 6 (CMIP6), while the natural emissions of sea salt, DMS, isoprene and dust are calculated interactively. The simulated monthly surface concentrations of sulfate (SO4), black carbon (BC), organic carbon (OA), and ozone (O3) are compared with observations from a set of Arctic stations, extracted from the EBAS and IMPROVE databases, as well as a few additional locations. Simulated aerosol optical depths (AOD) are also compared with Advanced Very-High Resolution Radiometer (AVHRR). The study will present the evaluation of the modelE2.1 in simulating SLCF levels over the Arctic using different aerosol schemes, along with observed and simulated trends of SLCFs over the Arctic between 1995 and 2015.
Aersosols↗
Potential surfaces for O atom-polymer reactions
Ab initio quantum chemistry methods are used to study the energetics of interactions of O atoms with organic compounds. Polyethylene (CH2)n has been chosen as the model system to study the interactions of O(3P) and O(1D) atoms with polymers. In particular, H abstraction is investigated and polyethylene is represented by a C3 (propane) oligomeric model. The gradient method, as implemented in the GRADSCF package of programs, is used to determine the geometries and energies of products and reactants. The saddle point, barrier geometry is determined by minimizing the squares of the gradients of the potential with respect to the internal coordinates. To correctly describe the change in bonding during the reaction at least a two configuration MCSCF (multiconfiguration self consistent field) or GVB (generalized valence bond) wave function has to be used. Basis sets include standard Pople and Dunning sets, however, increased with polarization functions and diffuse p functions on both the C and O atoms. The latter is important due to the O(-) character of the wave function at the saddle point and products. Normal modes and vibrational energy levels are given for the reactants, saddle points and products. Finally, quantitative energetics are obtained by implementing a small CAS (complete active space) approach followed by limited configuration interaction (CI) calculations. Comparisons are made with available experimental data.
Laskowski, B. C.↗