Author Topic: Different results by LDA and GGA  (Read 5953 times)

0 Members and 1 Guest are viewing this topic.

Offline yongjunwinwin

  • Heavy QuantumATK user
  • ***
  • Posts: 36
  • Country: cn
  • Reputation: 1
    • View Profile
Different results by LDA and GGA
« on: July 23, 2011, 13:15 »
I calculated the band structure of monolayer MoTe2 respectively with LDA and GGA, but got the complete different result! The band structure by LDA suggest metal property while GGA suggest semiconductor property with a direct band gap of 1.118 eV. What cause the difference and which result is right? (The results and the script  of the unite cell of monolayer MoTe2 is attached)

Offline zh

  • Supreme QuantumATK Wizard
  • *****
  • Posts: 1141
  • Reputation: 24
    • View Profile
Re: Different results by LDA and GGA
« Reply #1 on: July 23, 2011, 13:31 »
It has been reported in literature that the single atomic sheet of MoS2 is a semiconductor. If the LDA and GGA band structures are calculated for a same lattice of MoTe (i.e., the lattice constants and atomic coordinates in two cases are same), I guess that the possible reason comes from the  LDA pseudoptentials of Mo or Te. 

Now you had better check the LDA and GGA band structure of bulk MoTe2, and compare them with the experimental results and the previous theoretical work:
Th. Böker, et al, Phys. Rev. B 64, 235305 (2001), Band structure of MoS2, MoSe2, and α-MoTe2: Angle-resolved photoelectron spectroscopy and ab initio calculations
http://prb.aps.org/abstract/PRB/v64/i23/e235305 

Offline Anders Blom

  • QuantumATK Staff
  • Supreme QuantumATK Wizard
  • *****
  • Posts: 5576
  • Country: dk
  • Reputation: 96
    • View Profile
    • QuantumATK at Synopsys
Re: Different results by LDA and GGA
« Reply #2 on: August 7, 2011, 20:16 »
The results for LDA look strange, and actually the LDA basis set for Te looks wrong too (you can inspect it in VNL 11.2 by choosing Analyzer>Basis Set Explorer). There are 2 valence states (5s,5p), but if we check the DoubleZeta basis set, for GGA it looks normal, but in LDA we have 3 functions with nodes at r=0. If you zoom in, however, one of them is actually finite at r=0, but much too small. Also, the s state seems completely linear for small r.

Since the GGA results seem reasonable, you should definitely use GGA for this system.
« Last Edit: August 7, 2011, 20:21 by Anders Blom »