Just a simple question:
In my calculation the total-drift in z-direction isnt the same as the sum over all z-components of the forces (x and y likewise).
Why is that?
If I understand the source code right, this should however be the case.
What am I missing?
Thx for answering.
drfit and forces dont seem to match
Moderators: Global Moderator, Moderator
-
Franky
- Newbie
- Posts: 33
- Joined: Mon Apr 10, 2006 1:05 pm
- License Nr.: Research Group E. Pehlke
drfit and forces dont seem to match
Last edited by Franky on Thu Oct 05, 2006 9:27 am, edited 1 time in total.
-
admin
- Administrator
- Posts: 2922
- Joined: Tue Aug 03, 2004 8:18 am
- License Nr.: 458
drfit and forces dont seem to match
for the forces on the atoms, the drift is explicitely removed before they are written (in SUBROUTINE SYMVEC), whereas it is not done for the total drift. If the total drift is significantly larger than the foces on the atoms, your basis set most probably is too small.
Last edited by admin on Wed Oct 18, 2006 3:12 pm, edited 1 time in total.