Author Topic: Electron density difference integral?  (Read 333 times)

0 Members and 1 Guest are viewing this topic.

Offline Gabos

  • New ATK user
  • *
  • Posts: 1
  • Country: fr
  • Reputation: 0
    • View Profile
Electron density difference integral?
« on: April 15, 2021, 11:10 »
Hello ATK users, I am currently studying with DFT the interaction between a molecule and a surface, and I am interested to obtain a (possibly) quantitative analysis of the transferred charge between the two fragments.
Mulliken charges would be the obvious method of choice, however I would also like to explore other possibilities, such as the "charge displacement function". I would then need to compute some integrals on the electron density difference, which I already have. I have seen that there is this possibility within ATK, at least for the electron density (https://docs.quantumatk.com/manual/Types/ElectronDensity/ElectronDensity.html), however I am encountering some issues when I feed in the provided script my EDD file.
If I do
density = nlread('my_EDD.hdf5', ElectronDensity)[0]
I get an error message saying that the list index is out of range.

My question is, is it possible to perform integrals as shown in the manual also for an EDD file? If so, am I missing something trivial given the error message that I get?

Offline Anders Blom

  • QuantumATK Staff
  • Supreme ATK Wizard
  • *****
  • Posts: 4999
  • Country: dk
  • Reputation: 85
    • View Profile
    • QuantumATK at Synopsys
Re: Electron density difference integral?
« Reply #1 on: May 5, 2021, 22:57 »
The error message means your file does not contain any ElectronDensity. Maybe you have only computed and saved the ElectronDifferenceDensity?