<div dir="ltr"><div class="gmail_default" style="font-family:times new roman,serif;font-size:small">Hi all,<br><br></div><div class="gmail_default" style="font-family:times new roman,serif;font-size:small">I generated and saved a 3D grid in a molpro output file. The x,y,z coordinates and weights (w_k) are printed for each atom.For evaluating the atomic contribution to any molecular integral (say one-electron for now), I use the formula:<br>

<br></div><div class="gmail_default" style="font-family:times new roman,serif;font-size:small">I_k = sum(x,y,z) F(x,y,z) w_k(x,y,z) where k runs over all atomic indices.<br><br></div><div class="gmail_default" style="font-family:times new roman,serif;font-size:small">

If I then sum the atomic contributions, I should get the total number of electrons. My question is:</div><div class="gmail_default" style="font-family:times new roman,serif;font-size:small">Is there a missing factor involved in this formula? <br>

I tried a simple test of evaluating the atomic contribution to the electron ground state density where F stands for<br>F(x,y,z) = sum(AO1,AO2) AO1(x,y,z)* AO2(x,y,z) Denmat(AO1,AO2)<br>, and did not obtain the correct number of electrons. <br>

<br>Using the following molpro input for example, in a propane calculation gives me<br>Density functional                    25.99999792     STEST=25.99999792.<br><br></div><div class="gmail_default" style="font-family:times new roman,serif;font-size:small">

Input file:<br></div><div class="gmail_default" style="font-family:times new roman,serif;font-size:small"></div><div class="gmail_default" style="font-family:times new roman,serif;font-size:small">***,propane print grid<br>

</div><div class="gmail_default" style="font-family:times new roman,serif;font-size:small">memory,128,m<br></div><div class="gmail_default" style="font-family:times new roman,serif;font-size:small">symmetry,nosym<br> angstrom<br>

 gthresh,grid=1.0e-8<br> geometry={<br> 11<br><br> C        ,         1.82738088 ,  0.82142856 ,  0.00000000 ,<br> H        ,         2.18403531 , -0.18738144 ,  0.00000000 ,<br> H        ,         2.18405372 ,  1.32582675 , -0.87365150 ,<br>

 H        ,         0.75738088 ,  0.82144174 ,  0.00000000 ,<br> C        ,         2.34072310 ,  1.54738483 ,  1.25740497 ,<br> H        ,         3.41072308 ,  1.54720226 ,  1.25750243 ,<br> H        ,         1.98389084 ,  1.04309903 ,  2.13105625 ,<br>

 C        ,         1.82763676 ,  2.99939805 ,  1.25726452 ,<br> H        ,         2.18547770 ,  3.50398400 ,  0.38419930 ,<br> H        ,         0.75763749 ,  2.99993878 ,  1.25613883 ,<br> H        ,         2.18363852 ,  3.50325677 ,  2.13150084<br>

}<br>cartesian<br>basis=6-31g<br>{grid,1800.2,new;<br>radial,log,3,1.0;<br>angular,legendre;<br>gridprint,grid=2}<br>dfunc,stest<br>ks-scf<br><br></div><div class="gmail_default" style="font-family:times new roman,serif;font-size:small">

I can see that molpro correctly obtains the total number of electrons. <br></div><div class="gmail_default" style="font-family:times new roman,serif;font-size:small">Thanks for any help and suggestions,<br>Jayashree <br>
</div>
</div>