Hartree-Fock high-precision analytical functions of open-shell atoms

The algorithm of high-precision optimization of basis functions suggested previously for calculating the analytical Hartree-Fock orbitals of closed-shell atoms is generalized to open-shell systems described by the Roothaan method (1960). Expressions for the first (free gradient) and second (Hesse matrix) derivatives of the system's energy with respect to the nonlinear parameters (orbital exponents) of the basis functions are derived in terms of density matrices for the filled and open shells. An algorithm is proposed for high-precision optimization of the nonlinear parameters using these equations based on Murtagh-Sargent and Newton minimization procedures. To illustrate the application of this algorithm, we give optimization of the basis sets of Slater type functions for atoms from the second row, as well as for Al, Si, P, K, Sc, and Fe atoms. The analytical Hartree-Fock orbitals giving nearly Hartree-Fock energies are calculated with a high degree of accuracy.