Dear All, <div> </div><div>I was wondering if I was doing the right calculations using LRC-μBOP (HYB_GGA_XC_LC_BOP), and LRC-ωPBEh (HYB_GGA_XC_LRC_WPBEH) functionals in cp2k-8.2. Here are the &XC_FUNCTIONAL and &HF sections I used.</div><div> </div><div>For LRC-μBOP:</div><div>&XC</div><div>&XC_FUNCTIONAL </div><div>&HYB_GGA_XC_LC_BOP</div><div>_OMEGA 0.47</div><div>&END HYB_GGA_XC_LC_BOP</div><div>&END XC_FUNCTIONAL</div><div> </div><div>&HF</div><div>&SCREENING        </div><div>EPS_SCHWARZ 1.0E-6</div><div>SCREEN_ON_INITIAL_P TRUE</div><div>&END SCREENING</div><div>&INTERACTION_POTENTIAL</div><div>POTENTIAL_TYPE LONGRANGE </div><div>OMEGA 0.47  </div><div>&END INTERACTION_POTENTIAL</div><div>FRACTION 1.0</div><div>&END HF</div><div> </div><div>&VDW_POTENTIAL</div><div>DISPERSION_FUNCTIONAL PAIR_POTENTIAL</div><div>        &PAIR_POTENTIAL</div><div>          TYPE DFTD3</div><div>          CALCULATE_C9_TERM .TRUE.</div><div>          REFERENCE_C9_TERM .TRUE.</div><div>          LONG_RANGE_CORRECTION .TRUE.</div><div>          VERBOSE_OUTPUT .FALSE.</div><div>          REFERENCE_FUNCTIONAL BLYP</div><div>          PARAMETER_FILE_NAME dftd3.dat</div><div>          R_CUTOFF  9</div><div>          EPS_CN 1.0E-6</div><div>        &END PAIR_POTENTIAL</div><div>      &END VDW_POTENTIAL</div><div>    &END XC</div><div> </div><div>Question: </div><div>Are the two OMEGA parameters (i.e. _OMEGA in &HYB_GGA_XC_LC_BOP and OMEGA in &HF) independent with each other, or should be always the same as I set here? I really want to tune these long range parameters according to some references (e.g. <a href="http://doi.org/10.1021/jp5004243">doi.org/10.1021/jp5004243</a>).</div><div> </div><div>For LRC-ωPBEh where 20% HF is contained in short range:</div><div> </div><div>&XC</div><div>&XC_FUNCTIONAL </div><div>&HYB_GGA_XC_LRC_WPBEH</div><div>_OMEGA 0.47</div><div>&END HYB_GGA_XC_LRC_WPBEH</div><div>&END XC_FUNCTIONAL</div><div> </div><div>&HF</div><div>&SCREENING        </div><div>EPS_SCHWARZ 1.0E-6</div><div>SCREEN_ON_INITIAL_P TRUE</div><div>&END SCREENING</div><div>&INTERACTION_POTENTIAL</div><div>        POTENTIAL_TYPE MIX_CL_TRUNC </div><div>        OMEGA 0.47  </div><div>        SCALE_COULOMB 0.0</div><div>        SCALE_LONGRANGE 1.0</div><div>        CUTOFF_RADIUS  9.4</div><div>&END INTERACTION_POTENTIAL</div><div>FRACTION 1.0</div><div>&END HF</div><div> </div><div>Questions:</div><div>(1) Is such hybrid LRC functional available in the code? </div><div>(2) How can I implement the 20% HF exchange for the short range? It looks to me only LONGRANGE HF is defined here.</div><div>(3) Are the two OMEGA parameters independent with each other in &HYB_GGA_XC_LRC_WPBEH and &HF sections?</div><div> </div><div> </div><div>Thanks for advance!</div><div>Haibei</div>

<p></p>

-- <br />
You received this message because you are subscribed to the Google Groups "cp2k" group.<br />
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>.<br />
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/CAEqr0hOaLmnD6a19kcRyNDOWokATZo%2BcR7LAYh7JEukaAVfgtQ%40mail.gmail.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/CAEqr0hOaLmnD6a19kcRyNDOWokATZo%2BcR7LAYh7JEukaAVfgtQ%40mail.gmail.com</a>.<br />