quemb.molbe.sparse_2el_integral.get_sparse_DF_integrals

quemb.molbe.sparse_2el_integral.get_sparse_DF_integrals(mol, auxmol, screening_cutoff=None)

Return the 3-center 2-electron integrals in a sparse format with density fitting

One can obtain g[p, q, r, s] == ints_3c2e[p, q, :] @ df_coef[r, s, :].

Parameters:
  • mol (Mole) – The molecule.

  • auxmol (Mole) – The molecule with auxiliary basis functions.

  • screening_cutoff (int | float | floating | Callable[[int | float | floating], int | float | floating] | Mapping[str, int | float | floating] | None) – The screening cutoff is given by the overlap of van der Waals radii. By default, the radii are determined by find_screening_radius(). Alternatively, a fixed radius, callable or a dictionary can be passed. The callable is called with the tabulated van der Waals radius of the atom as argument and can be used to scale it up. The dictionary can be used to define different van der Waals radii for different elements. Compare to the modify_element_data argument of get_bonds().

Return type:

tuple[SemiSparseSym3DTensor, SemiSparseSym3DTensor]