24小时热门版块排行榜    

查看: 1502  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 南昌大学材料专硕311分求调剂 +5 77chaselx 2026-03-20 5/250 2026-03-20 23:42 by lovewei0727
[考研] 317求调剂 +5 申子申申 2026-03-19 9/450 2026-03-20 22:26 by JourneyLucky
[考研] 290求调剂 +7 ^O^乜 2026-03-19 7/350 2026-03-20 21:43 by JourneyLucky
[考研] 北科281学硕材料求调剂 +5 tcxiaoxx 2026-03-20 5/250 2026-03-20 21:35 by laoshidan
[考研] 319求调剂 +3 小力气珂珂 2026-03-20 3/150 2026-03-20 19:47 by JourneyLucky
[考研] 261求B区调剂,科研经历丰富 +3 牛奶很忙 2026-03-20 4/200 2026-03-20 19:34 by JourneyLucky
[考研] 086500 325 求调剂 +3 领带小熊 2026-03-19 3/150 2026-03-20 18:38 by 尽舜尧1
[考研] 265求调剂 +8 梁梁校校 2026-03-17 8/400 2026-03-20 14:40 by 27道科特
[考研] 材料学硕318求调剂 +5 February_Feb 2026-03-19 5/250 2026-03-19 23:51 by 23Postgrad
[考研] 266求调剂 +5 阳阳哇塞 2026-03-14 10/500 2026-03-19 15:08 by 阳阳哇塞
[考研] 一志愿福大288有机化学,求调剂 +3 小木虫200408204 2026-03-18 3/150 2026-03-19 13:31 by houyaoxu
[考研] 0703化学 305求调剂 +4 FY_yy 2026-03-14 4/200 2026-03-19 05:54 by anny19840123
[考研] 085601专硕,总分342求调剂,地区不限 +5 share_joy 2026-03-16 5/250 2026-03-18 14:48 by haxia
[考研] 297求调剂 +8 戏精丹丹丹 2026-03-17 8/400 2026-03-18 14:30 by laoshidan
[考研] 293求调剂 +11 zjl的号 2026-03-16 16/800 2026-03-18 08:10 by zhukairuo
[考研] 290求调剂 +3 p asserby. 2026-03-15 4/200 2026-03-17 16:35 by wangkm
[考研] 一志愿,福州大学材料专硕339分求调剂 +3 木子momo青争 2026-03-15 3/150 2026-03-17 07:52 by laoshidan
[考研] 求老师收留调剂 +4 jiang姜66 2026-03-14 5/250 2026-03-15 20:11 by Winj1e
[考研] 085601材料工程315分求调剂 +3 yang_0104 2026-03-15 3/150 2026-03-15 10:58 by peike
[考研] 本科南京大学一志愿川大药学327 +3 麦田耕者 2026-03-14 3/150 2026-03-14 20:04 by 外星文明
信息提示
请填处理意见