otter.bound_wavefunctions

otter.bound_wavefunctions(aa: Mapping[str, Any], *, multiply_fd: bool = False) → ndarray[source]

Return radial \(R_{nl}(r)\) on aa['r_bound'].

The (l, n_index, r_bound) result is unweighted by default. If requested, multiply the amplitude by f_FD(E), not sqrt(f_FD), as a display/analysis convention. Neither variant applies M(E), f_cut, or degeneracy. The FD variant is not a new orbital normalization and must not be used to reconstruct the ionic density (which already includes f_FD once).