short range + long range wPBE


Just Got Here
Dear NWchem center

May I ask some help to build exchange functional?

if I want use wPBE as DFT part and set alpha=0.2, beta=0.7

can you please help me for that?
thanks.
Best wishes,
Zilong

Forum Regular
Hi Zilong,

I am not sure I completely understand your entire question. I'll try to answer it as clearly as possible.

>>if I want use wPBE as DFT part and set omega=0.160 bohr-1, alpha=0.2, beta=0.7, which following definition is correct?

You could try the following:
xc xwpbe 1.0 cpbe96 1.0 hfexch 1.00
cam 0.160 cam_alpha 0.20 cam_beta 0.70

Here, the long-range aymptotic behavior of the HF term would be 0.9 * HF. If you want to use a different coefficient for the xwpbe, then you can change the 1.0 next to it to the appropriate value.


>>I do not fully understand the meaning of the number, such as in (xwpbe 0.8) or (xwpbe 0.7 xpbe96 0.1)
The 0.8, 0.7 etc above correspond to the mixing coefficient of the exchange or correlation part of the functional

Hope this helps.

Best,
-Niri
niri.govind@pnnl.gov

Just Got Here
Hello, Niri, thanks so much for you message.

Just Got Here


Quote:Niri Jul 18th 10:22 am
Hi Zilong,

I am not sure I completely understand your entire question. I'll try to answer it as clearly as possible.

>>if I want use wPBE as DFT part and set omega=0.160 bohr-1, alpha=0.2, beta=0.7, which following definition is correct?

You could try the following:
xc xwpbe 1.0 cpbe96 1.0 hfexch 1.00
cam 0.160 cam_alpha 0.20 cam_beta 0.70

Here, the long-range aymptotic behavior of the HF term would be 0.9 * HF. If you want to use a different coefficient for the xwpbe, then you can change the 1.0 next to it to the appropriate value.


>>I do not fully understand the meaning of the number, such as in (xwpbe 0.8) or (xwpbe 0.7 xpbe96 0.1)
The 0.8, 0.7 etc above correspond to the mixing coefficient of the exchange or correlation part of the functional

Hope this helps.

Best,
-Niri
niri.govind@pnnl.gov





Hi Niri,

How do one incorporate the long-range DFT exchange fraction i.e., $(1- \alpha - \beta)$ in the above input. It is zero when \alpha+\beta=1. Do we have to set this value as well in the input or the code will pick up the alpha, beta values for this term?


Thank you,
Chandra

Forum Regular
Hi Chandra,

The above definition is for the traditional range-separated xc functionals, where you have short-range DFT and long-range HF.

For long-range DFT and short-range HF, as in the HSE functional, you can define it as follows, where we use the srhfexch.

HSE03:
xc xpbe96 1.0 xcampbe96 -0.25 cpbe96 1.0 srhfexch 0.25
cam 0.33 cam_alpha 0.0 cam_beta 1.0

HSE06:
xc xpbe96 1.0 xcampbe96 -0.25 cpbe96 1.0 srhfexch 0.25
cam 0.11 cam_alpha 0.0 cam_beta 1.0

Hope this helps.

Let me know if there are any issues.

Best,
-Niri
niri.govind@pnnl.gov

Just Got Here
Quote:Niri Apr 24th 2:08 pm
Hi Chandra,

The above definition is for the traditional range-separated xc functionals, where you have short-range DFT and long-range HF.

For long-range DFT and short-range HF, as in the HSE functional, you can define it as follows, where we use the srhfexch.

HSE03:
xc xpbe96 1.0 xcampbe96 -0.25 cpbe96 1.0 srhfexch 0.25
cam 0.33 cam_alpha 0.0 cam_beta 1.0

HSE06:
xc xpbe96 1.0 xcampbe96 -0.25 cpbe96 1.0 srhfexch 0.25
cam 0.11 cam_alpha 0.0 cam_beta 1.0

Hope this helps.

Let me know if there are any issues.

Best,
-Niri
niri.govind@pnnl.gov




Hi Niri,

Thank you very much for your reply. I have tried using above inputs to incorporate long-range DFT (xcampbe96), short-range DFT (xwpbe) and short-range HF (srhexch). But, using above method I cannot use long-range HF.

Is there a way that I can use both long range HF exchange and short range HF exchange together.



Best,
Chandra


Forum >> NWChem's corner >> NWChem functionality