I am trying to develop an implementation of a certain electronic structure method within Julia aimed at modelling strong light-matter interactions. The method is closely related to the widely popular Hartree-Fock method. For this, I would be needing a lot of basis-set related integrals to be computed, namely the one and two electron electron repulsion integrals (which I found can be easily gotten by specifying molecular geometry and basis set using Fermi.jl or GaussianBasis.jl) and dipole moment integrals. What I cannot understand is that how to retrieve are the dipole moment integrals in the atomic orbital (AO) basis, this is essential for the method in question to work.
Given the time constraints, obviously it is not a good idea to write the handling of basis sets and integral libraries from scratch. Is there any way I can get these integrals using some package or from somewhere else?