[molpro-user] Expectation of second moment and kinetic energy
Peter Knowles
KnowlesPJ at cardiff.ac.uk
Fri Feb 28 08:00:54 GMT 2014
On 27 Feb 2014, at 20:04, Yannis Petousis <petousis at stanford.edu> wrote:
> 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
Without seeing your whole job, it's difficult to comment. Hartree-Fock satisfies the virial theorem, but only in the limit of a complete basis set. You should not expect the kinetic energy to be equal to minus the total energy unless the forces are zero.
> and the expectation of the second moment components are negative!
The second moment of course depends on the origin, and you will in any case get contributions of opposite sign from nuclei and electrons, so it is difficult to comment.
Peter
>
> Here is my input:
> {hf;
> start,2100.2;
> wf,250,1,0;
> expec,pot,ekin,ov,sm,delta,tm;
> }
> show,epot,kin,ov,sm,delta,tm;
>
>
> 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,
> Yannis
>
> Ioannis Petousis
> PhD Student
> Nanoscale Prototyping Laboratory
> Mechanical Engineering dpt.
> Stanford University
> email: petousis at stanford.edu
> _______________________________________________
> Molpro-user mailing list
> Molpro-user at molpro.net
> http://www.molpro.net/mailman/listinfo/molpro-user
--
Prof. Peter J. Knowles
School of Chemistry, Cardiff University, Main Building, Park Place, Cardiff CF10 3AT, UK
Telephone +44 29 208 79182 Email KnowlesPJ at Cardiff.ac.uk
WWW http://www.cardiff.ac.uk/chemy/contactsandpeople/academicstaff/knowles.html
More information about the Molpro-user
mailing list