help,how can I get the potential energy of each atom in a structure using DFT package


by wkxez
Tags: classical md, dft(vasp), potential energy
wkxez
wkxez is offline
#1
Mar14-12, 04:46 AM
P: 6
I am a freshman in DFT calculation. I don't understand the DFT method clearly,but I want to use DFT to do some calculations.
My question is that how how can I get the potential energy of each atom in a structure using DFT package ,which I have done by using classical MD method in Lammps(the command: compute 1 all pe/atom ,Lammps manual :The per-atom energy is calculated by the various pair, bond, etc potentials defined for the simulation. If no extra keywords are listed, then the potential energy is the sum of pair, bond, angle, dihedral,improper, and kspace energy.).
I want to know how can I do such a calculation.Please give some suggestions.
Phys.Org News Partner Physics news on Phys.org
Global scientific team 'visualizes' a new crystallization process (w/ video)
Modification of structural composite materials to create tailored lenses
Novel technique opens door to better solar cells
wkxez
wkxez is offline
#2
Mar15-12, 12:20 AM
P: 6
I now know that i can get the force F of each atom from the results of vasp relaxation calculation .In the classical theory,F= dU/dr(U is the potential energy),,so if it's possible to get the potential energy from the VAsp or other package calculation?
Useful nucleus
Useful nucleus is offline
#3
Mar15-12, 06:15 AM
Useful nucleus's Avatar
P: 240
The notion of energy per atom is ill defined for a non-crystalline material. IF you are modelling a crystal , you simply divide the total energy by the number of unit formulas of the crysta.

wkxez
wkxez is offline
#4
Mar18-12, 08:22 PM
P: 6

help,how can I get the potential energy of each atom in a structure using DFT package


Quote Quote by Useful nucleus View Post
The notion of energy per atom is ill defined for a non-crystalline material. IF you are modelling a crystal , you simply divide the total energy by the number of unit formulas of the crysta.

Thank you for your reply,Why is it ill defined for non-crystalline material? And you may misunderstand my question. I want to get the relative potential energy of each atom in a structure,not energy per atom.
DrDu
DrDu is offline
#5
Mar19-12, 03:12 AM
Sci Advisor
P: 3,366
Quote Quote by wkxez View Post
Thank you for your reply,Why is it ill defined for non-crystalline material? And you may misunderstand my question. I want to get the relative potential energy of each atom in a structure,not energy per atom.
In an infinite solid this is difficult, in a molecule you can in principle remove that atom from the calculation and compare the energy of the molecule with and without the atom. In practice there is a problem called basis set superposition error, and you may have to replace the atom by a ghost atom with the same basis set but zero nuclear charge.
For a solid you would have to consider probably some large super cell.


Register to reply

Related Discussions
fine structure, energy levels of the hydrogen atom Introductory Physics Homework 0
potential energy in hydrogen atom Quantum Physics 15
Potential energy of an Atom Introductory Physics Homework 3
Potential energy of a helium atom Introductory Physics Homework 8
potential energy of a loaded structure General Physics 2