Hello everyone! I would like to calculate surface stress. In the OUTCAR stress in x-, y- and z-directions is given. Are these stress components (Total) in eV/Ã… or eV/(Ã…*Ã…*Ã…)? Often the surface stress is presented in the form eV/(Ã…*Ã…). How is surface stress obtained in this unit (eV/(Ã…*Ã…)) from the x- and y-components given in the OUTCAR?
I have been pondering these questions, because I get significantly different results (x,y,z stress components) for the same system, when I use slabs with different surface sizes. Therefore, I have been thinking, whether these components should be multiplied or divided by the lattice vectors to get comparable values.
Any advice would be great and most welcome. Thank you!
surface stress
Moderators: Global Moderator, Moderator
-
- Newbie
- Posts: 29
- Joined: Mon Jun 19, 2006 7:03 am
- License Nr.: 853
- Location: FAMU
surface stress
Last edited by vasp on Fri Nov 09, 2007 2:33 pm, edited 1 time in total.
-
- Newbie
- Posts: 29
- Joined: Mon Jun 19, 2006 7:03 am
- License Nr.: 853
- Location: FAMU
surface stress
I realize that the surface stress is very sensitive to the details of the calculation, and that the results are comparable between systems of different surface areas. Furthermore, the unit surface area and strain should be taken into account in the units.
Last edited by vasp on Tue Nov 13, 2007 8:47 am, edited 1 time in total.
-
- Newbie
- Posts: 7
- Joined: Thu Jul 26, 2007 2:32 pm
surface stress
To my understanding, stresses are in units eV/reduced length, as stated in the OUTCAR file. Therefore, they are dependent on your unit cell size and have to be divided by the lattice vector in the direction in question to obtain results in eV/Ã….
Last edited by oja on Thu Dec 13, 2007 1:09 pm, edited 1 time in total.
-
- Newbie
- Posts: 29
- Joined: Mon Jun 19, 2006 7:03 am
- License Nr.: 853
- Location: FAMU
surface stress
'To my understanding, stresses are in units eV/reduced length, as stated in the OUTCAR file. Therefore, they are dependent on your unit cell size and have to be divided by the lattice vector in the direction in question to obtain results in eV/Ã….'
Yes, I understood this. Thank you for your answer!
Yes, I understood this. Thank you for your answer!
Last edited by vasp on Thu Dec 27, 2007 8:03 am, edited 1 time in total.