<div dir="ltr"><div><div>I hope this isn't a terribly pedestrian question, but is there a still a low limit to the number of points at which one can compute the electrostatic potential around a molecule with MolPro?  I found this as about the most helpful response, but it is for old code:<br><br><a href="http://www.molpro.net/pipermail/molpro-user/2008-September/002606.html" target="_blank">http://www.molpro.net/pipermail/molpro-user/2008-September/002606.html</a><br><br></div>I wasn't able to get very far with the example in that thread, but I was able to run the following input:<br><br>>>><br>memory,500,m<br><br>! AAA2<br>geometry={<br>  H1    0.9572000000000   0.0000000000000   0.0000000000000<br>  O2    0.0000000000000   0.0000000000000   0.0000000000000<br>  H3   -0.2399900000000   0.9266300000000   0.0000000000000<br>}<br>gdirect<br>basis=cc-pvtz<br>hf<br>df-mp2<br><br>!{matprop<br>!read, savmat, DEN, charge, file=density.mat<br>!save, savmat, 2140.2, DENSITY}<br><br>{property<br>density<br>pot,,   4.00000,   0.00000,   0.00000<br>pot,,  -1.00288,   3.86820,   0.00000<br>pot,,   0.00000,   0.00000,   4.00000<br>}<br><<<<br><br></div><div>Once I figured out that geometry={ } gets specified in Angstroms and {property } coordinates get specified in Bohr I was in good shape.  However, it seems that the limit of thirty points using this method is still in effect.  Is there a better way that I can compute the electrostatic potential for perhaps 5000 points following one MP2 calculation?<br><br></div><div>Thanks,<br><br>Dave<br><br></div></div>