Author Topic: How to calculate the change in the electrostatic energy  (Read 4383 times)

0 Members and 1 Guest are viewing this topic.

Offline guyue03

  • Regular QuantumATK user
  • **
  • Posts: 11
  • Country: cn
  • Reputation: 0
    • View Profile
Hi,
 
  I wonder how to calculate the change in the electrostatic energy. I've already done a script to calculate the electrostatic potential of the device, but an error occurred "The parameter 'configuration' must be an instance of one of the following: BulkConfiguration, MoleculeConfiguration".
  My goal is to draw a picture like the following picture.
  Can you help me solve this problem?

Offline Petr Khomyakov

  • QuantumATK Staff
  • Supreme QuantumATK Wizard
  • *****
  • Posts: 1290
  • Country: dk
  • Reputation: 25
    • View Profile
Re: How to calculate the change in the electrostatic energy
« Reply #1 on: December 7, 2018, 11:43 »
The Electrostatic Potential analysis object can only be computed for Bulk and Molecule Configuration, see https://docs.quantumwise.com/manual/Types/ElectrostaticPotential/ElectrostaticPotential.html.

For Device Configuration, you may compute Electrostatic Difference Potential , see https://docs.quantumwise.com/manual/Types/ElectrostaticDifferencePotential/ElectrostaticDifferencePotential.html, or actually Hartree Difference Potential energy if you see changes in the electrostatic potential energy.

If you have computed the potential energy for the entire system and different sub-systems of your molecular junction, then you may visualize the potential change by computing dE = E(entire system) - E(sub-system 1) - E(sub-system 2) - ... . The 2 subsystems can be, for instance, the electrodes without the molecules (this should be done using the same device configuration, just remove the molecule)  and the molecule without electrodes (you have to use bulk configuration with the unit cell size as that of the central region of the device, and use the multigrid Poisson solver with the Dirichlet and Neumann boundary conditions imposed on the left and right side of the slab, respectively).

To compute the dE, you can use Grid Operation plugin available in the QuantumATK NanoLab, just select all the potential objects on the LabFloor and click on Grid Operation on the left panel of the GUI.

The visualization of the new grid object can be done in the Viewer, and you may also drag-and-drop the configuration from the LabFloor to the Viewer to superimpose the potential and the atomic structure of the molecular junction.

Note that the same operation can be done for electron density or any object on the real-space grid.

Offline guyue03

  • Regular QuantumATK user
  • **
  • Posts: 11
  • Country: cn
  • Reputation: 0
    • View Profile
Re: How to calculate the change in the electrostatic energy
« Reply #2 on: December 8, 2018, 01:42 »
Hi,

  If I have computed the dE, then is that mean I could get the energy jump at the interface, as shown in the picture?
  Thank you!

Offline guyue03

  • Regular QuantumATK user
  • **
  • Posts: 11
  • Country: cn
  • Reputation: 0
    • View Profile
Re: How to calculate the change in the electrostatic energy
« Reply #3 on: December 9, 2018, 04:59 »
Hi,

  I've followed your guidance to compute the dE. I am faced with a new problem now. I have attached all files, including the error, and py files. Could you help me?
  Thank you so much!!!!

Offline Petr Khomyakov

  • QuantumATK Staff
  • Supreme QuantumATK Wizard
  • *****
  • Posts: 1290
  • Country: dk
  • Reputation: 25
    • View Profile
Re: How to calculate the change in the electrostatic energy
« Reply #4 on: December 10, 2018, 14:49 »
Please read my previous post carefully.