I found the elastic constants and bulk modules of diamond is wrong by a VASP job
the input files below
==========================
SYSTEM = Di
ICHARG = 2
ENCUT = 700
ISMEAR = -5
IBRION = 6
NSW = 10111
ISIF = 4
EDIFF = 0.00000001
NFREE = 2
LCHARG = F
LWAVE = F
POTIM = 0.015
=================================
K
0
M
15 15 15
0 0 0
=============================
here are some outputs:
===================================
TOTAL ELASTIC MODULI (kBar)
Direction XX YY ZZ XY YZ ZX
--------------------------------------------------------------------------------
XX 11070.1035 1500.6893 1500.6893 0.0000 0.0000 0.0000
YY 1500.6893 11070.1035 1500.6893 0.0000 0.0000 0.0000
ZZ 1500.6893 1500.6893 11070.1035 0.0000 0.0000 0.0000
XY 0.0000 0.0000 0.0000 5911.9935 0.0000 0.0000
YZ 0.0000 0.0000 0.0000 0.0000 5911.9935 0.0000
ZX 0.0000 0.0000 0.0000 0.0000 0.0000 5911.9935
--------------------------------------------------------------------------------
--------------------------------------------------------------------------------------------------------
LATTYP: Found a simple cubic cell.
ALAT = 3.5329679502
Lattice vectors:
A1 = ( 3.5329679502, 0.0000000000, 0.0000000000)
A2 = ( 0.0000000000, 3.5329679502, 0.0000000000)
A3 = ( 0.0000000000, 0.0000000000, 3.5329679502)
Subroutine PRICEL returns following result:
LATTYP: Found a face centered cubic cell.
ALAT = 3.5329679502
Lattice vectors:
A1 = ( 1.7664839751, 1.7664839751, 0.0000000000)
A2 = ( 1.7664839751, 0.0000000000, -1.7664839751)
A3 = ( 0.0000000000, 1.7664839751, -1.7664839751)
4 primitive cells build up your supercell.
========================================
<span class='smallblacktext'>[ Edited ]</span>
wrong elstic constans of diamone
Moderators: Global Moderator, Moderator
-
- Newbie
- Posts: 13
- Joined: Tue Mar 23, 2010 5:14 am
wrong elstic constans of diamone
Last edited by ntuxi on Sun Mar 28, 2010 11:19 am, edited 1 time in total.
-
- Newbie
- Posts: 13
- Joined: Tue Mar 23, 2010 5:14 am
wrong elstic constans of diamone
btw: the bulk modules of diamond by the vasp results is >700Gpa, it should be under 500
Last edited by ntuxi on Sun Mar 28, 2010 11:20 am, edited 1 time in total.
-
- Administrator
- Posts: 2921
- Joined: Tue Aug 03, 2004 8:18 am
- License Nr.: 458
wrong elstic constans of diamone
which potential did you use,
how large was the remaining external pressure on the unit cell,
which procedure did you use to determine the bulk lattice constant of your cell?
which convergence criterium did you set for the remaining forces/stress tensors
how large was the remaining external pressure on the unit cell,
which procedure did you use to determine the bulk lattice constant of your cell?
which convergence criterium did you set for the remaining forces/stress tensors
Last edited by admin on Mon Mar 29, 2010 2:45 pm, edited 1 time in total.