fn open_shell_pt2(scf_data: &SCF) -> Result<[f64; 3]>