Author Topic: Single zeta polarized and Double Zeta polarised  (Read 3852 times)

0 Members and 1 Guest are viewing this topic.

Offline abhishek77492652

  • QuantumATK Guru
  • ****
  • Posts: 111
  • Country: in
  • Reputation: 0
    • View Profile
Single zeta polarized and Double Zeta polarised
« on: November 7, 2017, 15:17 »
Hello sir,
I have calculated the chemical potential for silicon crystal (alpha silicon) using both the DZP ansd SZP basis set and get the value as -4.081992 eV and -4.265549 eV  respectively. Please tell me which one is more accurate .
Also mention how to select basis set  (like DZP, SZP,...) in different situation .

THANKS in advance

Offline Petr Khomyakov

  • QuantumATK Staff
  • Supreme QuantumATK Wizard
  • *****
  • Posts: 1290
  • Country: dk
  • Reputation: 25
    • View Profile
Re: Single zeta polarized and Double Zeta polarised
« Reply #1 on: November 7, 2017, 22:10 »
You may look at the total energy, not the chemical potential. Lower energy usually means more accurate basis set, because the total energy is a variational  quantity.

However, other physical quantities may not converge with respect to the LCAO-based basis set monotonously, so in that regard it is hard to say whether your electronic structure is already converged with respect the basis. You may want to use even more accurate basis set such as SG15-based one to have a better idea on this.

Another option is to use our plane-wave calculator in ATK-2017.2. Then you may systematically check for convergence of physical quantities with respect to the plane-wave bases basis set, by varying a single parameter, which is the kinetic energy cut-off.
« Last Edit: November 7, 2017, 22:15 by Petr Khomyakov »

Offline Jess Wellendorff

  • QuantumATK Staff
  • Supreme QuantumATK Wizard
  • *****
  • Posts: 933
  • Country: dk
  • Reputation: 29
    • View Profile
Re: Single zeta polarized and Double Zeta polarised
« Reply #2 on: November 8, 2017, 07:35 »
In addition, the calculated chemical potential is not very meaningful in itself; you need to compare it to some other energy-quantity computed using the same pseudopotentials+basis set in order to use it.