Author Topic: effective potential  (Read 7130 times)

0 Members and 1 Guest are viewing this topic.

Offline xiaolight

  • Regular QuantumATK user
  • **
  • Posts: 14
  • Country: us
  • Reputation: 0
    • View Profile
effective potential
« on: August 22, 2017, 16:44 »
in http://quantumwise.com/forum/index.php?topic=5114.msg22147#msg22147
It tells me I can use effective potential to calculate the electron affinity but the data is really large than I thought. I use the MoS2 supercell and make the vacuum level really large.

Offline Anders Blom

  • QuantumATK Staff
  • Supreme QuantumATK Wizard
  • *****
  • Posts: 5576
  • Country: dk
  • Reputation: 96
    • View Profile
    • QuantumATK at Synopsys
Re: effective potential
« Reply #1 on: August 22, 2017, 20:05 »
I'm not sure you can calculate the workfunction/affinity of a 2D sheet this way. The method essentially assumes that you can define a "bulk-like" region and compare the potential there to vacuum. In this case, it would mean cutting a surface in the 2D plane. So, in your case, rather, you would repeat the structure many times in B or C (not both), then add vacuum in that direction, and then study the potential in this direction, using the averaged value across the perpendicular plane (and ideally also average in the projection direction to get rid of the atomic variation).

Alternatively, and maybe better, you might try to use Dirichlet/Neumann boundary conditions and use the fact that then the Fermi level represents the work function (for metal); for semiconductors you need to also factor in the band gap.

Offline Kim_W

  • Heavy QuantumATK user
  • ***
  • Posts: 59
  • Country: cn
  • Reputation: 0
    • View Profile
Re: effective potential
« Reply #2 on: October 22, 2017, 15:10 »
By comparing the values reported in the tutorials for the (100), (110) and (111) crystalline faces of Ag, it appears clear that the work function depends on the crystalline face considered. Thus, the method you posted may be not reasonable. Can we use two layers MoS2 for this calculation? And turning the first layer MoS2 into gost atoms.

Offline Petr Khomyakov

  • QuantumATK Staff
  • Supreme QuantumATK Wizard
  • *****
  • Posts: 1290
  • Country: dk
  • Reputation: 25
    • View Profile
Re: effective potential
« Reply #3 on: October 23, 2017, 09:30 »
By comparing the values reported in the tutorials for the (100), (110) and (111) crystalline faces of Ag, it appears clear that the work function depends on the crystalline face considered.
This is a well-known experimentally-confirmed fact, and not for Ag only.

Thus, the method you posted may be not reasonable.
Which method? And why not reasonable?

Can we use two layers MoS2 for this calculation? And turning the first layer MoS2 into gost atoms.
You can certainly do it. Whether this will give an improvement for accuracy of work function or affinity calculation, it is another question.

Offline Kim_W

  • Heavy QuantumATK user
  • ***
  • Posts: 59
  • Country: cn
  • Reputation: 0
    • View Profile
Re: effective potential
« Reply #4 on: October 23, 2017, 10:28 »
Thus, the method you posted may be not reasonable.
Which method? And why not reasonable?

So, in your case, rather, you would repeat the structure many times in B or C (not both), then add vacuum in that direction, and then study the potential in this direction, using the averaged value across the perpendicular plane (and ideally also average in the projection direction to get rid of the atomic variation).

This one.


Offline Petr Khomyakov

  • QuantumATK Staff
  • Supreme QuantumATK Wizard
  • *****
  • Posts: 1290
  • Country: dk
  • Reputation: 25
    • View Profile
Re: effective potential
« Reply #5 on: October 23, 2017, 12:43 »
Why is it not a "reasonable" method for work function calculation?

Offline Jess Wellendorff

  • QuantumATK Staff
  • Supreme QuantumATK Wizard
  • *****
  • Posts: 933
  • Country: dk
  • Reputation: 29
    • View Profile
Re: effective potential
« Reply #6 on: October 25, 2017, 14:45 »
Dear Wang, the work function does in general depend on surface facet.