Hybrid functionals: Difference between revisions

From Wiki
Jump to navigation Jump to search
Fraenzine (talk | contribs)
New page: go back to Main Page, Group Pages, Núria López and Group This page is under edition, it may not make sense yet :) (feel free to ask, if you have questions) =Hybrid function...
 
Fraenzine (talk | contribs)
No edit summary
Line 11: Line 11:
==VASP==  
==VASP==  


* For [all] hybrid functional calculations, set:  
* For all hybrid functional calculations, set:  
   LHFCALC = .TRUE.
   LHFCALC = .TRUE.


(...)
* For HSE calculations, set the HFSCREEN parameter:
  HFSCREEN = 0.2 # HSE06
  HFSCREEN = 0.3 # HSE03


The HFSCREEN parameter specifies the distance (in 1/Å) at which the exact (Hartree-Fock) exchange is screened.
Therefore HF exact exchange is screened at a distance larger than 5 Å (HSE06), or 3.333 Å (HSE03), respectively, the way it is implemented in VASP. The difference between HSE03 and HSE06 is discussed in detail by Krukau et al. J. Chem. Phys. 125, 224106 (2006): https://aip.scitation.org/doi/abs/10.1063/1.2404663


* Defining the amount of introduced exact exchange:
  AEXX  = 0.25  # Amount of exact HF exchange, ie. 25%
  AGGAX = 0.75  # Amount of GGA exchange, i.e. 75 %


(...)
This is the default setting with 25 % exact exchange, you may reduce that number, in particular if you have compounds with transition metals and/or degenerate electronic states (see theory above).


(see also: https://cms.mpi.univie.ac.at/wiki/index.php/Specific_hybrid_functionals)
* HSE03-13: This functional was optimised to fit the band gap of CoFe Prussian blue compounds; it generally presents a good compromise for calculating ionic solid compounds, which contain coordinated transition metal centres
  LHFCALC = .TRUE.
  HFSCREEN = 0.3 # HSE03
  AEXX  = 0.13  # Amount of exact HF exchange, ie. 13%
  AGGAX = 0.87  # Amount of GGA exchange, i.e. 87 %


There are of course many other hybrid functionals, which require different settings, but generally they can be described by defining the GGA functional and the above-listed parameters.


(see also: https://cms.mpi.univie.ac.at/wiki/index.php/Specific_hybrid_functionals)


Install the full 'Demeter' package, that includes the 'Athena' program. It only works well in Windows 64-bits. https://bruceravel.github.io/demeter/


==Quantum Espresso==
==Quantum Espresso==


For each structure you want to analyze:
* Create a folder.
* Put the VASP output in that folder.
* Convert the structure to a .cif file.
* Convert the .cif to feff.inp using this page: http://millenia.cars.aps.anl.gov/webatoms/. If that page is offline, open Artemis (part of Demeter package), file, open project or data, pick your cif file.
* In the webatoms page or in Artemis, select the edge, select feff6_element style, and select the atom that will be the center of the analysis. In Artemis, clic in "Run Atoms". It will generate a file that contains something like this:
  (...)
  (...)

Revision as of 13:16, 2 July 2019

go back to Main Page, Group Pages, Núria López and Group

This page is under edition, it may not make sense yet :) (feel free to ask, if you have questions)

Hybrid functional theory

tbd

Implementation

VASP

  • For all hybrid functional calculations, set:
 LHFCALC = .TRUE.
  • For HSE calculations, set the HFSCREEN parameter:
 HFSCREEN = 0.2 # HSE06
 HFSCREEN = 0.3 # HSE03

The HFSCREEN parameter specifies the distance (in 1/Å) at which the exact (Hartree-Fock) exchange is screened. Therefore HF exact exchange is screened at a distance larger than 5 Å (HSE06), or 3.333 Å (HSE03), respectively, the way it is implemented in VASP. The difference between HSE03 and HSE06 is discussed in detail by Krukau et al. J. Chem. Phys. 125, 224106 (2006): https://aip.scitation.org/doi/abs/10.1063/1.2404663

  • Defining the amount of introduced exact exchange:
 AEXX  = 0.25   # Amount of exact HF exchange, ie. 25%
 AGGAX = 0.75   # Amount of GGA exchange, i.e. 75 %

This is the default setting with 25 % exact exchange, you may reduce that number, in particular if you have compounds with transition metals and/or degenerate electronic states (see theory above).

  • HSE03-13: This functional was optimised to fit the band gap of CoFe Prussian blue compounds; it generally presents a good compromise for calculating ionic solid compounds, which contain coordinated transition metal centres
 LHFCALC = .TRUE.
 HFSCREEN = 0.3 # HSE03
 AEXX  = 0.13   # Amount of exact HF exchange, ie. 13%
 AGGAX = 0.87   # Amount of GGA exchange, i.e. 87 %

There are of course many other hybrid functionals, which require different settings, but generally they can be described by defining the GGA functional and the above-listed parameters.

(see also: https://cms.mpi.univie.ac.at/wiki/index.php/Specific_hybrid_functionals)


Quantum Espresso

(...)