# [molpro-user] DFT calculations produced different energy with a different guesses

Kirk Peterson kipeters at wsu.edu
Wed Nov 3 15:01:21 GMT 2010

```David,

it looks to me that you're getting stuck into a few local minima.  Do these different calculations converge to different structures?  How different are the energies at the starting geometries?  I would strongly recommend printing the orbitals and taking a careful look.   You could also try doing a CASSCF calculation (print the orbitals and CI vector)  to see how things should really be occupied.

regards,

Kirk

On Nov 3, 2010, at 6:19 AM, David Danovich wrote:

> Hello,
>
> I am calculating CuAu2 molecule in high spin state (4B2 state) using B3LYP method. Below you can find three different possibilities I have used in the calculations. In the first one 1) I have calculated first HF wave function and used it as a guess for DFT calculation. Energy I got was -468.86467917 au. In the second calculation 2) I directly have calculated B3LYP without HF wave function. The energy I got was -468.83995914 au. In the third calculation 3) I first have done B3LYP calculation for state (4A2) and then used it as a guess for the calculations of 4B2. The energy I got was
> -468.85824088 au. As you can see there is quite large difference in energy for the same state. I was trying to use different grids but it does not solved the problem. In my opinion result should not depend on the guess so drastically. What can be a solution for this problem?
>
> Thank you in advance  David
>
> 1)
>
> ***, Peterson PP tz basis set
> memory,250,m
>  r =   2.72684681 ang;
>  a = 60.65125012 degree;
> geometry={Cu1;              !z-matrix geometry input
>           Au2,Cu1,r;
>           Au3,Cu1,r,Au2,a;
>           }
> basis=cc-pVTZ-PP
>
> {hf,maxit=500;
> wf,57,3,3}
>
> {rks,b3lyp3,maxit=500;
> wf,57,3,3}
>
> optg
>
> Total energy  -468.86467917 au
> ____________________________________________________
>
> 2)
> ***, Peterson PP tz basis set
> memory,250,m
>  r =   2.72684681 ang;
>  a = 60.65125012 degree;
> geometry={Cu1;              !z-matrix geometry input
>           Au2,Cu1,r;
>           Au3,Cu1,r,Au2,a;
>           }
> basis=cc-pVTZ-PP
>
> {ks,b3lyp3,maxit=500;
> shift,-0.1,0.0;
> wf,57,3,3}
>
> optg
>
> Total energy -468.83995914 au
>
> ______________________________________________________
>
> 3)
> ***, Peterson PP tz basis set
> memory,250,m
>  r =   2.72684681 ang;
>  a = 60.65125012 degree;
> geometry={Cu1;              !z-matrix geometry input
>           Au2,Cu1,r;
>           Au3,Cu1,r,Au2,a;
>           }
>
> ! cc-pVTZ-PP
>
> basis=cc-pVTZ-PP
>
> {ks,b3lyp3,maxit=500;
> wf,57,4,3}
>
> {ks,b3lyp3,maxit=500;
> shift,-1.0,0.0;
> wf,57,3,3}
>
> optg
>
> Total energy -468.85824088 au
>
> _______________________________________________
> Molpro-user mailing list
> Molpro-user at molpro.net
> http://www.molpro.net/mailman/listinfo/molpro-user

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.molpro.net/pipermail/molpro-user/attachments/20101103/4bf4c02b/attachment.html>
```