ini_wfn = psi4.core.Wavefunction.build(h2o, psi4.core.get_global_option(‘basis’))
print “Alpha Orbital coefficient”
print ini_wfn.basisset()

mints = psi4.core.MintsHelper(ini_wfn.basisset())
S = np.asarray(mints.ao_overlap())
T = np.asarray(mints.ao_kinetic())
V = np.asarray(mints.ao_potential())

H = T + V
print “Overlap matrix”
print S
print “Kinetic matrix”
print T
print “Potential matrix”
print V
print “Core-Hamiltonian matrix”
print H

Thank you for the reply. I’ve already seen that tutorial but couldn’t find how to get kinetic energy matrix without mintshelper. I mean, wonder how to make AO basis integral explicitly with psi4numpy without using mintshelper and only with AO basis information. Sorry for the confusion.

MintsHelper is how you obtain this information. Is there a reason this is not sufficient?

MintsHelper is effectively a shell looper for one and two electron integrals. Is that you want other one electron integrals or is the object too heavy?

Yes, I want to do treat some non-trivial (may be not physically meaningful…) 1e integral term and to do so, I guess it will be a great start if I can make the same energy matrices as mintshelper did.

If its a custom integral you will have to loop over the shells and build these yourself. The following line compute the kinetic integrals for a shell pair:

You can pull Gaussian Shells out of BasisSets by using the basis.shell(n) function where n is a number. The number of shells in a basis can be obtained as basis.nshell(). As you can see its a fairly non-trivial task.