<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 3.2//EN">
<HTML>
<HEAD>
<META HTTP-EQUIV="Content-Type" CONTENT="text/html; charset=iso-8859-1">
<META NAME="Generator" CONTENT="MS Exchange Server version 6.5.7638.1">
<TITLE>PBE Hybrid funcitonals and LDA pseudopotentials</TITLE>
</HEAD>
<BODY>
<!-- Converted from text/plain format -->

<P><FONT SIZE=2>Hi,<BR>
<BR>
I have been attempting to use the hybrid functional PBE0 to tune the bandgap of InGaAs, while using Virtual Crystal Approximation to obtain In(0.57)Ga(0.43)As. I initially tried mixing the PBE pseudos of Ga and In, but ran into a number of problems (unequal nqf mesh when mixing, and using ultrasoft pseudos with hybrid functionals).<BR>
<BR>
I then tried the Norm-conserving LDA pseudos (which also have no nqf information to worry about) for mixing In and Ga (also using LDA for As pseudo, obviously). I then found I could tune the bandgap by specifying "PBE0" in the &SYSTEM card of the input file by changing the proportion of Fock exchange (alpha, or "exx_fraction" in funct.f90), while specifying LDA pseudos for each atomic species.<BR>
<BR>
I am surprised this would work because my understanding is that LDA contains no exchange term, yet modifying the value of alpha has an effect on the bandgap just like if I had specified PBE pseudopotentials in the ATOMIC_SPECIES card. I initially set up this calculation as a test to just prove to myself that the results obtained using LDA pseudopotentials are independant of the amount of exchange, yet this is not so.<BR>
<BR>
My understanding of hybrid functional implementation is obviously incomplete, so can someone help me understand these results?<BR>
<BR>
Thank you very much for the help,<BR>
<BR>
Gabriel Greene,<BR>
PhD student, Electronics Theory Group,<BR>
Tyndall National Institute,<BR>
Cork Ireland<BR>
</FONT>
</P>

</BODY>
</HTML>