I am interested in representing an arbitrary charge density (say, of atoms in a molecule) rho(r),;rinmathbbR3 by a finite linear combination of basis functions
rho(r)=sumNi=1qiphii(r)
where phii(r) is normalized to intmathbbR3phii(r)dr=1 and has the interpretation of being the shape of some charge distribution (shape) of a unit charge. rho and the phiis are real-valued but may be positive in some regions and negative in others. The basis functions are nonorthogonal and local in space but not strictly compact. Let's say for now that we use spherical Gaussians of the form phii(r)proptoexp(−alphai|r−Ri|2), where Ri is where the basis function is centered around. The number of basis functions chosen scales approximately as the number of atoms, as we expect charge to concentrate around atomic nuclei. (We may add additional basis function per atom of different shapes until we achieve a reasonable approximation to the desired shape of the charge distribution around an atom.)
The energy of the system can then be given by
E=frac12sumNi,j=1qiqjJij
where the matrix J has elements
Jij=intmathbbR3times2fracphii(r1)phij(r2)|r1−r2|dr1dr2
and represents the Coulomb interaction between the unit charge distributions phii and phij.
One way to look at the matrix J is as a finite dimensional (approximate) representation of the Coulomb operator hatJ=1/|r1−r2|. We know that hatJ has certain nice properties such as positivity, so we expect a "good" representation of hatJ should be a symmetric positive definite matrix.
My question is this: are there conditions on the discrete representation (possibly expressible as conditions on the {phii} basis) to detect whether or not a given claimed representation J is "good" in that it preserves such properties? Or asked another way, if I have some matrix J which is claimed to represent hatJ, what are necessary and sufficient conditions on its matrix elements for it to be a "good" representation of hatJ?
I hope the question makes sense, and that I am not misusing too much terminology.
No comments:
Post a Comment