[molpro-user] Expectation of second moment and kinetic energy

Yannis Petousis petousis at stanford.edu
Thu Feb 27 19:04:53 GMT 2014

Dear molpro users,

In a simple rhf calculation of around 50 atoms, I am trying to find the expectation values of kinetic energy and second moment of position.
My understanding is that this is done using the respective operators, EKIN and SM.
However, the expectation of the kinetic energy I am getting, violates the virial theorem and the expectation of the second moment components are negative!

Here is my input:

And here is the output:
!RHF expec            <1.1|XX|1.1>  -532.748689831615
!RHF expec            <1.1|XX|1.1>  -532.748689831615
!RHF expec            <1.1|ZZ|1.1>  -522.658296533770
!RHF expec          <1.1|EKIN|1.1>   122.891738342572

Note: I am using no symmetry so my only state is 1.1 

Can anyone help?

Thank you,

Ioannis Petousis
PhD Student
Nanoscale Prototyping Laboratory
Mechanical Engineering dpt.
Stanford University
email: petousis at stanford.edu

More information about the Molpro-user mailing list