[CPMD-list] combination of PBE USPP with BLYP DFT
Yi Liu
yliu252 at uwo.ca
Fri Oct 31 20:15:50 CET 2003
----- Original Message -----
From: "Erik Emilio Santiso" <eesantis at unity.ncsu.edu>
To: "Jung Goo Lee" <jglee at nemo.physics.ncsu.edu>
Cc: <cpmd-list at cpmd.org>
Sent: Thursday, October 30, 2003 9:14 PM
Subject: Re: [CPMD-list] combination of PBE USPP with BLYP DFT
>
> Hi! Well, this looks more like a coincidence than anything else. I don't
> think it makes any sense to find a PP with one xc functional and then do
> the calculation with another one. What happens if you use the PBE PP's
> with the PBE functional?
>
> Erik.
>
> On Thu, 30 Oct 2003, Jung Goo Lee wrote:
>
> > Hi everyone,
> >
> > I would like to know whether the combination of PBE uspp in ATOMS option
> > and BLYP in DFT
> > option (the example of partial input is shown below)makes sense or not,
> > especially, for calculations
> > of free energy (constraint). This is because BLYP method underestimate
> > the dissociation free energy
> > compared with experiment. For example, water dissociation in 32 water
> > molecules with MT PP/
> > BLYP method gave 17 kcal/mol in free energy while experimental value is
> > 21 kcal/mol(Trout and
> > Parrinello, Chem. Phys. Lett, 288(1998)343.
> > I tried the same calculation with PBE uspp and BLYP method.
> > I got about 20 kcal/mol althogh I tried with a few hundred fs (instead
> > of 1 ps) MD at each point of
> > O-H bond constraint. So, this combination (PBE uspp/BLYP) seems to be
> > o.k. in free energy
> > calculations.
> >
> > Is this a reasonable and cheap(computationally) way to get out of the
> > problem with free energy
> > calculations with gbl uspp/BLYP?
> > Does any one have any rationalization or contradiction with reasonable
> > explanation?
> >
> > Jung-Goo Lee
> >
> > Example
> > --
> > --
> > &ATOMS
> > *008-O-gpbe--bm.uspp NEWF BINARY
> > LMAX=P
> > 1
> > 4.3728 4.5000 4.5000 1 8
> >
> > *001-H-gpbe--bm.uspp NEWF BINARY
> > LMAX=S
> > 2
> > 5.5172 4.5000 5.9238 1 1
> > 5.5172 4.5000 3.0761 1 1
> > &END
> >
> > &DFT
> > OLDCODE
> > FUNCTIONAL BLYP
> > GRADIENT CORRECTION BECKE88 LYP
> > LDA CORRELATION LYP
> > &END
> > --
> > --
> >
> >
> > _______________________________________________
> > CPMD-list mailing list
> > CPMD-list at cpmd.org
> > http://www.cpmd.org/mailman/listinfo/cpmd-list
> >
> _______________________________________________
> CPMD-list mailing list
> CPMD-list at cpmd.org
> http://www.cpmd.org/mailman/listinfo/cpmd-list
More information about the CPMD-list
mailing list