Author Topic: Bader charge analysis  (Read 8828 times)

0 Members and 1 Guest are viewing this topic.

Offline bensg

  • New QuantumATK user
  • *
  • Posts: 2
  • Country: fr
  • Reputation: 0
    • View Profile
Bader charge analysis
« on: April 25, 2017, 10:51 »
Dear All,

I use the bader charge analysis tool in Virtual Nanolab to analyze VASP calculations. Would you know what version of the bader software is implemented in VNL ?
When using the bader software directly (version 0.95 11/17/15, as technical issues prevented me from using the version 1.00 01/11/17), I obtain different charges than the ones given by VNL.
Also, I face an issue concerning periodic boundary conditions: for some calculations, the bader charges given by VNL are completely incoherent for some of the atoms on the edge of the simulation cell. For instance, with a hexagonal boron nitride single layer, although most of the atomic charges makes sense, on the edge of the cell I could have a nitrogen with a +5e charge and a boron with a -5e charge. Would you know a possible origin for this issue and a possible fix ?

Thank you very much,

Benoit

Offline Daniele Stradi

  • Supreme QuantumATK Wizard
  • *****
  • Posts: 286
  • Country: dk
  • Reputation: 3
    • View Profile
Re: Bader charge analysis
« Reply #1 on: April 25, 2017, 12:46 »
Dear Benoit,

We do not use the Bader software by the Henkelman group, we have an original implementation. Could it be possible for you to send us or attach an small example of the calculations which gives you problems for debugging, including the output files? It should also include the results that you get with the original Bader code.

Regards,
Daniele.

Offline bensg

  • New QuantumATK user
  • *
  • Posts: 2
  • Country: fr
  • Reputation: 0
    • View Profile
Re: Bader charge analysis
« Reply #2 on: April 25, 2017, 16:39 »
Dear Daniele,

Thank you very much for your fast answer. By which mean would you like me to send you an example of calculations results (I have a 372 Mo zip file to provide you)?

Best wishes,

Benoit

Offline Daniele Stradi

  • Supreme QuantumATK Wizard
  • *****
  • Posts: 286
  • Country: dk
  • Reputation: 3
    • View Profile
Re: Bader charge analysis
« Reply #3 on: April 25, 2017, 22:58 »
Dear Benoit,

I have replied you in a private message. We will investigate if this is an issue with the bader analysis tool for VASP implemented in ATK.

Regards,
Daniele.

Offline Kim_W

  • Heavy QuantumATK user
  • ***
  • Posts: 59
  • Country: cn
  • Reputation: 0
    • View Profile
Re: Bader charge analysis
« Reply #4 on: October 26, 2017, 10:00 »
How to know the result is resonable? I done it with FHI-aims.  Besides, how to get the Bader charges transferred from the Al atom to N atom with below txt scprit.

Offline Ulrik G. Vej-Hansen

  • QuantumATK Staff
  • Supreme QuantumATK Wizard
  • *****
  • Posts: 426
  • Country: dk
  • Reputation: 9
    • View Profile
Re: Bader charge analysis
« Reply #5 on: October 26, 2017, 13:35 »
Please provide the script you have used to get the output you have posted.

Offline Daniele Stradi

  • Supreme QuantumATK Wizard
  • *****
  • Posts: 286
  • Country: dk
  • Reputation: 3
    • View Profile
Re: Bader charge analysis
« Reply #6 on: October 26, 2017, 16:22 »
Dear Wang,

How to know the result is resonable?
The reliability of the result is based on chemical intuition. If there is an error, you would normally see an unphysically large charge on some atoms, which is not the case in your calculation.

Besides, how to get the Bader charges transferred from the Al atom to N atom with below txt scprit
To get the charge transferred to a certain atom, you should compare the charge calculated using Bader analysis with the atomic number of the corresponding isolated atom. In your case, for example, you can see that 9.4122 - 7 = 2.4122 electrons have been donated to N. However notice the two WARNING messages at the beginning of the log file. I would suggest you to increase the fineness of the grid to obtain more converged results.

Daniele

Offline Kim_W

  • Heavy QuantumATK user
  • ***
  • Posts: 59
  • Country: cn
  • Reputation: 0
    • View Profile
Re: Bader charge analysis
« Reply #7 on: October 27, 2017, 05:03 »
How to increase the fineness of the grid? Is that k-point or Accuracy or else? I calculated with FHI-aims calculator.  Thank you very much.

Offline Ulrik G. Vej-Hansen

  • QuantumATK Staff
  • Supreme QuantumATK Wizard
  • *****
  • Posts: 426
  • Country: dk
  • Reputation: 9
    • View Profile
Re: Bader charge analysis
« Reply #8 on: October 27, 2017, 10:11 »
Again, please provide the script, it makes it much easier to help you. The text-file you have posted indicates some obvious mistakes, but we cannot help you with correcting them before you have provided the scripts.

Offline Kim_W

  • Heavy QuantumATK user
  • ***
  • Posts: 59
  • Country: cn
  • Reputation: 0
    • View Profile
Re: Bader charge analysis
« Reply #9 on: October 30, 2017, 04:35 »
 :)Thank you!

Offline Ulrik G. Vej-Hansen

  • QuantumATK Staff
  • Supreme QuantumATK Wizard
  • *****
  • Posts: 426
  • Country: dk
  • Reputation: 9
    • View Profile
Re: Bader charge analysis
« Reply #10 on: October 30, 2017, 10:43 »
Please provide also the script which actually gave the output you posted previously. This script will not give the output you posted before.

Offline Kim_W

  • Heavy QuantumATK user
  • ***
  • Posts: 59
  • Country: cn
  • Reputation: 0
    • View Profile
Re: Bader charge analysis
« Reply #11 on: November 1, 2017, 03:04 »
Firstly, I calculated the ED with this script. Then, I dealed with "Bader Charge Analysis" tool.

Offline Ulrik G. Vej-Hansen

  • QuantumATK Staff
  • Supreme QuantumATK Wizard
  • *****
  • Posts: 426
  • Country: dk
  • Reputation: 9
    • View Profile
Re: Bader charge analysis
« Reply #12 on: November 1, 2017, 15:16 »
Thanks for the update. I cannot reproduce the error with my version of VNL, which is 2017.1. Which version did you use?

I think the problem is due to a previous bug in our FHI-AIMS interface, which gave incorrect electron densities, and was fixed in 2017.1. However, it appears that the grid is much too coarse, so I suggest that you try increasing it quite a bit, if you want accurate results.
« Last Edit: November 1, 2017, 15:19 by Ulrik G. Vej-Hansen »

Offline Kim_W

  • Heavy QuantumATK user
  • ***
  • Posts: 59
  • Country: cn
  • Reputation: 0
    • View Profile
Re: Bader charge analysis
« Reply #13 on: November 3, 2017, 10:28 »
The question is: How to increase the fineness of the grid? Is that k-point or Accuracy or else? Thank you a lot.

Offline Ulrik G. Vej-Hansen

  • QuantumATK Staff
  • Supreme QuantumATK Wizard
  • *****
  • Posts: 426
  • Country: dk
  • Reputation: 9
    • View Profile
Re: Bader charge analysis
« Reply #14 on: November 3, 2017, 11:23 »
When we talk about the grid, we usually mean the real-space grid which is used to represent the density. The relevant parameter is called density_mesh_cutoff. See more here: https://docs.quantumwise.com/manuals/Types/NumericalAccuracyParameters/NumericalAccuracyParameters.html
« Last Edit: November 3, 2017, 12:06 by Petr Khomyakov »