Hi,
I was wondering if the following is possible in vasp:
A static calculation of a cell at a non-zero temperature.
I want to use this to check if a defect atom becomes unstable/less stable in favor of the clean system with regard to a certain temperature.
For this I would do calculations at different temperatures for both systems, and their difference in E0 or F(at certain T) should change sign at some point.
Can this be done? or do I have to do some sort of MD calculation for this?
thx,
Danny
Temperature and energy
Moderators: Global Moderator, Moderator
-
- Full Member
- Posts: 201
- Joined: Thu Nov 02, 2006 4:35 pm
- License Nr.: 5-532
- Location: Ghent, Belgium
- Contact:
Temperature and energy
Last edited by Danny on Mon Oct 20, 2008 3:19 pm, edited 1 time in total.
-
- Administrator
- Posts: 2921
- Joined: Tue Aug 03, 2004 8:18 am
- License Nr.: 458
Temperature and energy
please note that vasp only accounts for the electronic T, i.e. the excitation of certain phonon modes at elevated T is not taken into account for. Therefore I doubt whether your problem can be solved easily with vasp.
Last edited by admin on Tue Oct 21, 2008 11:36 am, edited 1 time in total.
-
- Full Member
- Posts: 201
- Joined: Thu Nov 02, 2006 4:35 pm
- License Nr.: 5-532
- Location: Ghent, Belgium
- Contact:
Temperature and energy
(a difficult solution with vasp would also be usefull )
So, just to be sure, the E given for a MD calculation(with IBRION=0, SMASS= -3 and POTIM=5.0) at the end of the first ionic step can not be considered to include the thermal temperature of the system?
thx
Danny
So, just to be sure, the E given for a MD calculation(with IBRION=0, SMASS= -3 and POTIM=5.0) at the end of the first ionic step can not be considered to include the thermal temperature of the system?
thx
Danny
Last edited by Danny on Wed Oct 22, 2008 12:47 am, edited 1 time in total.
-
- Administrator
- Posts: 2921
- Joined: Tue Aug 03, 2004 8:18 am
- License Nr.: 458
Temperature and energy
if you do MD, T includes the electronic T (via the semaring width of the electronic levels) and the contribution due to the kinetic energy of the ions ( T=2*Ekinetic/(the degrees of freedom) ).
Last edited by admin on Wed Oct 22, 2008 8:44 am, edited 1 time in total.