frequency calculation
Moderators: Global Moderator, Moderator
-
- Newbie
- Posts: 6
- Joined: Fri Jul 09, 2010 7:13 am
frequency calculation
When doing frequency calculation (IBROIN=5),it is suggested to build a supercell large enough.However,if my cell is very small,building a large supercell(eg 3*3*3) would induce too many atoms in the system,and slow down the calculation significantly.What should I do ? Thanks!
Last edited by tssunrise on Wed Jul 21, 2010 10:08 am, edited 1 time in total.
-
- Hero Member
- Posts: 586
- Joined: Tue Nov 16, 2004 2:21 pm
- License Nr.: 5-67
- Location: Germany
frequency calculation
Why do want to calculate a supercell if your system does _not_ require one?
Last edited by alex on Wed Jul 21, 2010 12:04 pm, edited 1 time in total.
-
- Sr. Member
- Posts: 339
- Joined: Mon Apr 24, 2006 9:07 am
- License Nr.: 173
- Location: Gothenburg, Sweden
frequency calculation
@alex: I do not understand your answer. If you want to evaluate the dynamical matrix by using finite difference, you will have to use a larger supercell. Otherwise your force constants will be effected significantly by the displacements in the periodic cells.
@ tssunrise: I would suggest that you first do the calculation using a 1*1*1 cell and a 2*2*2 cell. Then if there is a significant difference and/or to check for convergence, you will have to accept that it is necessary and perform the 3*3*3 supercell calculation as well. (Alternatively you can gradually extent the supercell in one direction at the time to estimate the effect, 2*2*2 -> 2*2*3 -> 2*3*3 -> 3*3*3. )
EDIT: If you have access to Vasp 5.2, I recommend to use IBRION = 6 instead to reduce the number of needed displacements. Otherwise if you are using Vasp 4.6, you can do a corresponding symmetry analysis, and then assemble the dynamical matrix by hand.
Cheers,
/Dan
<span class='smallblacktext'>[ Edited Wed Jul 21 2010, 04:09PM ]</span>
@ tssunrise: I would suggest that you first do the calculation using a 1*1*1 cell and a 2*2*2 cell. Then if there is a significant difference and/or to check for convergence, you will have to accept that it is necessary and perform the 3*3*3 supercell calculation as well. (Alternatively you can gradually extent the supercell in one direction at the time to estimate the effect, 2*2*2 -> 2*2*3 -> 2*3*3 -> 3*3*3. )
EDIT: If you have access to Vasp 5.2, I recommend to use IBRION = 6 instead to reduce the number of needed displacements. Otherwise if you are using Vasp 4.6, you can do a corresponding symmetry analysis, and then assemble the dynamical matrix by hand.
Cheers,
/Dan
<span class='smallblacktext'>[ Edited Wed Jul 21 2010, 04:09PM ]</span>
Last edited by forsdan on Wed Jul 21, 2010 12:41 pm, edited 1 time in total.
-
- Hero Member
- Posts: 586
- Joined: Tue Nov 16, 2004 2:21 pm
- License Nr.: 5-67
- Location: Germany
frequency calculation
@forsdan: right. I was not considering tiny cells
Last edited by alex on Wed Jul 21, 2010 3:21 pm, edited 1 time in total.