def analyze_optimized(name, atoms, A_idx, B_idx, basis='cc-pVTZ', xc='B3LYP'): # Get atom coords (in Bohr, convert to Angstrom for find_BCP) A_ang = mol_eq.atom_coord ...