Is there a relationship between the momentum operator matrix elements and the following:
<φ|dH/dkx|ψ>
where kx is the Bloch wave number
such that if I have the latter calculated for the x direction as a matrix, I can get the momentum operator matrix elements from it?