Atomic structure of N-electron atoms is often determined using the Hartree-Fock method, which is an integro-differential equation. The exchange term of the Hartree-Fock equations is usually treated as an inhomogeneous term of a differential equation, or with a local density approximation. This work uses matrix methods to solve for the Hartree-Fock equations, rather than the more commonly-used shooting method to integrate an inhomogeneous differential equation. It is well known that a derivative operator can be expressed as a matrix made of finite-difference coefficients; energy eigenvalues and eigenvectors can be obtained by using computer linear-algebra packages. We extend the same technique to integro-differential equations, where a discretized integral can be written as a sum in matrix form. This method is compared against experiment and standard atomic structure calculations. We also can use this method for free-electron wavefunctions. This technique is important for spectral line broadening in two ways: improving the atomic structure calculations, and improving the motion of the plasma electrons that collide with the atom.

In this work we outline a general method for calculating Auger spectra in molecules, which accounts for the underlying symmetry of the system. This theory starts from Fano’s formulation of the interaction between discrete and continuum states and generalizes this formalism to deal with the simultaneous presence of several intermediate quasi-bound states and several non-interacting decay channels. Our theoretical description is specifically tailored to resonant autoionization and Auger processes, and includes explicitly the incoming wave boundary conditions for the continuum states and an accurate treatment of the Coulomb repulsion. This approach is implemented and applied to the calculation of the K−LL Auger and autoionization spectra of ozone, which is a C2v symmetric molecule, whose importance in our atmosphere to filter out radiation has been widely confirmed. We also show the effect that the molecular point group, and in particular the localization of the core-hole in the oxygen atoms related by symmetry operations, has on the electronic structure of the Auger states and on the spectral lineshape by comparing our results with experimental data.

A new framework in quantum chemistry has been proposed recently ("An approach to first principles electronic structure calculation by symbolic-numeric computation'' by A. Kikuchi). It is based on the modern technique of computational algebraic geometry, viz. the symbolic computation of polynomial systems. Although this framework belongs to molecular orbital theory, it fully adopts the symbolic method. The analytic integrals in the secular equations are approximated by the polynomials. The indeterminate variables of polynomials represent the wave-functions and other parameters for the optimization, such as atomic positions and contraction coefficients of atomic orbitals. Then the symbolic computation digests and decomposes the polynomials into a tame form of the set of equations, to which numerical computations are easily applied. The key technique is Grröbner basis theory, by which one can investigate the electronic structure by unraveling the entangled relations of the involved variables. In this article, at first, we demonstrate the featured result of this new theory. Next, we expound the mathematical basics concerning computational algebraic geometry, which are necessitated in our study. We will see how highly abstract ideas of polynomial algebra would be applied to the solution of the definite problems in quantum mechanics. We solve simple problems in "quantum chemistry in algebraic variety'' by means of algebraic approach. Finally, we review several topics related to polynomial computation, whereby we shall have an outlook for the future direction of the research.

The theory of one-photon ionization and two-photon above-threshold ionization is formulated for applications to heavy atoms in attosecond science using the Dirac-Fock formalism. A direct comparison of the Wigner-Smith-Eisenbud delays for photoionization is made with delays from the Reconstruction of Attosecond Beating By Interference of Two-photon Transitions (RABBIT) method. Photoionization by an attosecond pulse train, consisting of monochromatic fields in the extreme ultraviolet range, is computed with many-body effects at the level of the Relativistic Random Phase Approximation (RRPA). Subsequent absorption and emission processes of infrared laser photons in RABBIT are evaluated using static ionic potentials as well as asymptotic properties of relativistic Coulomb functions. As expected, light elements, such as Argon, show negligible relativistic effects, while heavier elements, such a Krypton and Xenon, exhibit delays that depend on the fine-structure of the ionic target. The relativistic effects are notable close to ionization thresholds and Cooper minima with differences in fine-structure delays predicted to be as large as tens of attoseconds. The separability of relativistic RABBIT delays into a Wigner-Smith-Eisenbud delay and a universal continuum--continuum delay is studied with reasonable separability found for photoelectrons emitted along the laser polarization axis in agreement with prior non-relativistic results.