24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1422  |  回复: 9

zhengbohan

木虫 (小有名气)

[求助] 计算弹性常数的问题

在用vasp5.2算单斜YVO4的弹性常数,得出结果都是负值,是不是哪里设置的问题?求助
TOTAL ELASTIC MODULI (kBar)
Direction    XX          YY          ZZ          XY          YZ          ZX
--------------------------------------------------------------------------------
XX       -4478.7063  -4089.5812  -7210.0791      0.0000  -1034.4726    900.4692
YY       -4089.5812  -1886.8150  -5325.5829      0.0000  -1045.9579    217.1807
ZZ       -7210.0791  -5325.5829  -5696.3147      0.0000  -1221.4168    989.1406
XY           0.0000      0.0000      0.0000    478.8831      0.0000      0.0000
YZ       -1034.4726  -1045.9579  -1221.4168      0.0000   -782.3298   -623.3351
ZX         900.4692    217.1807    989.1406      0.0000   -623.3351  -2280.0065
--------------------------------------------------------------------------------
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

成功不必在我,而功力必不唐捐
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (职业作家)

本帖仅楼主可见
2楼2015-01-06 20:44:01
已阅   申请1ST强帖   回复此楼   编辑   查看我的主页

infrared123

铁杆木虫 (著名写手)

To be excellent

不要怕,不要悔
3楼2015-01-07 11:05:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhengbohan

木虫 (小有名气)

Start parameter for this run:
LWAVE   = False
LCHARG  = False
Electronic relaxation:
ENCUT   = 400.0     ! Cutoff energy for plane waves basis set in eV
#ENAUG   = 750.0     ! Cutoff energy for the augmentation charge
ALGO    = Fast    ! A Davidson blocked iteration scheme
NELMIN  = 8         ! Minimum number of eletronic selfconsistency (SC) steps
NELM    = 1000        ! Maximum number of electronic SC steps
NELMDL  = -4        ! Number of NON-selfconsistency steps
EDIFF   = 1.0E-4    ! Global-break condition for the electronic SC-loop (ELM)
AMIX    = 0.4
BMIX    = 1.000
Calculation mode:
PREC = Accurate  ! Calcululation level (Changes FFT-grids)
#ISPIN   = 2         ! Non-spin-polarized calculations
ADDGRID = .TRUE.    ! PAW flag
#MAGMOM =  2*0 3*1 3*-1 8*0 24*0
#PREC    = Normal    ! Should not be used for final total energy calculations
#LREAL   = AUTO      ! Should not be used for final total energy calculations
#####LMAXMIX = 6         ! PAW flag

Ionic relaxation:
NSW     = 100        ! Number of ionic steps
EDIFFG  = -0.0200   ! stop if all forces are smaller than |EDIFFG|
IBRION  = 6          !
ISIF    = 3          ! Force, stress, relax ions, change cell shape, change cell volume
NFREE   = 4
POTIM   = 0.30

Integration over the Brillouin zone (BZ):
ISMEAR  = -5
SIGMA   = 0.05
#LOPTICS = .TRUE.
#NBANDS = 400
DOS calculation:
#LORBIT  = 10        ! Calculate the DOS without providing the Wigner Seitz radius
#NEDOS   = 4001      ! Number of points to calculate the DOS
#EMin=-25
#Emax=15

OUTCAR size:
NWRITE  = 2         ! Determines how much information will be written in OUTCAR

#LVTOT = .TRUE.
#IDIPOL =      3 ! dipole correction
#DIPOL = 0.5 0.5 0.5
#LDIPOL = .TRUE.

Key for parallel mode calculation:
#NPAR = 1
LPLANE = .TRUE.
#NSIM = 4
EOF
成功不必在我,而功力必不唐捐
4楼2015-01-07 11:43:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (职业作家)

本帖仅楼主可见
5楼2015-01-09 11:36:42
已阅   申请1ST强帖   回复此楼   编辑   查看我的主页

zhengbohan

木虫 (小有名气)

引用回帖:
5楼: Originally posted by zyl雨田 at 2015-01-09 11:36:42
我算的体系是六角结构 没有对称性的时候  弹性常数里面有 负值  提高对称性后 没有了

怎么提高的对称性?小白本科僧一枚

[ 发自小木虫客户端 ]
成功不必在我,而功力必不唐捐
6楼2015-01-09 23:15:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (职业作家)

★ ★
fzx2008: 金币+2, 谢谢指导 2015-01-14 08:37:41
本帖仅楼主可见
7楼2015-01-10 10:22:01
已阅   申请1ST强帖   回复此楼   编辑   查看我的主页

zhengbohan

木虫 (小有名气)

引用回帖:
7楼: Originally posted by zyl雨田 at 2015-01-10 10:22:01
1、你先在OUTCAR里查看,VASP计算时找到的体系的对称性是否与体系应该具有的对称性一致
2、如果不一致,将参数SYMPREC设成 1e-3 (默认值1e-5) (这个参数在OUTCAR里没有显示)...

我试试哈

[ 发自小木虫客户端 ]
成功不必在我,而功力必不唐捐
8楼2015-01-10 12:22:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhengbohan

木虫 (小有名气)

引用回帖:
7楼: Originally posted by zyl雨田 at 2015-01-10 10:22:01
1、你先在OUTCAR里查看,VASP计算时找到的体系的对称性是否与体系应该具有的对称性一致
2、如果不一致,将参数SYMPREC设成 1e-3 (默认值1e-5) (这个参数在OUTCAR里没有显示)...

我怎么没有找到这个参数呢?它在哪里?

[ 发自小木虫客户端 ]
成功不必在我,而功力必不唐捐
9楼2015-01-12 19:15:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (职业作家)

本帖仅楼主可见
10楼2015-01-12 21:03:57
已阅   申请1ST强帖   回复此楼   编辑   查看我的主页
相关版块跳转 我要订阅楼主 zhengbohan 的主题更新
信息提示
请填处理意见