ESP (Electrostatic Potential in NWChem)

From NWChem

Viewed 192 times, With a total of 1 Posts
Jump to: navigation, search

Just Got Here
Threads 1
Posts 2
I am concerned about the units of the ESP in NWChem. Consider the data on HF in the output file:

1 a.u. = 9.07618 esu/cm ( or statvolts ) 
Point X Y Z Total Potential(a.u.) Diamagnetic shielding(a.u.)
   1 F    0.00000   0.00000  -0.08806     -26.551580            27.123542
2 H 0.00000 0.00000 1.66031 -0.905417 6.053075



Now consider the same data from Gamess on the same molecule:


electric potential 1 a.u. = 9.076814E-02 esu/cm

                                 Total Potential (a.u.)             Diamagnetic shielding (a.u.)
F -26.554587 27.126252
H -0.902754 6.047740


Is it possible that NWChem has an error in the conversion factor from a.u. to esu/cm?
Edited On 1:20:11 PM PDT - Tue, Jul 18th 2017 by Matrix2017

Just Got Here
Threads 1
Posts 2
Proof of an error using an H atom
The potential energy of H in the ground state is exactly 27.2116 Volts
The calculated potential energy = constant * 299.7925 * the calculated value in A.U.
The value, 299.7925 is the conversion from esu/cm to volts where esu / cm are the units of the electrostatic potential energy in NWChem or Gamess.
constant(NWChem) = 9.07618, constant(Gamess)=0.09076814

V(Gamess) = 0.09076814 * 299.7925 * (-0.990160) = -26.943845 Volts

V(NWChem) = 9.07618 * 299.7925 * (-0.990745) = -2695.788108 Volts

V(Exact) = 27.2116 Volts


The results suggest that NWChem may have an error in the conversion constant.


Forum >> NWChem's corner >> Feedback



Who's here now Members 0 Guests 0 Bots/Crawler 1


AWC's: 2.5.10 MediaWiki - Stand Alone Forum Extension
Forum theme style by: AWC