24小时热门版块排行榜    

查看: 1659  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 01:49 by ZPa0EcMwuECS
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 01:45 by ZPa0EcMwuECS
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 gy1nBQXYQJqL 2026-08-29 7/350 2026-08-30 01:36 by ZPa0EcMwuECS
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:36 by ZPa0EcMwuECS
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:25 by ZPa0EcMwuECS
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 01:14 by ZPa0EcMwuECS
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 ASdOkHsho7FD 2026-08-28 10/500 2026-08-30 00:53 by ZPa0EcMwuECS
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 ASdOkHsho7FD 2026-08-28 5/250 2026-08-30 00:52 by ZPa0EcMwuECS
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 ASdOkHsho7FD 2026-08-28 12/600 2026-08-30 00:31 by ZPa0EcMwuECS
[基金申请] 国社科申报系统有变化 +3 kynobel 2026-08-23 3/150 2026-08-29 21:57 by 余韵清
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 gy1nBQXYQJqL 2026-08-29 5/250 2026-08-29 21:44 by 4qO8zhfjr2AX
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[论文投稿] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 gy1nBQXYQJqL 2026-08-29 3/150 2026-08-29 18:26 by 4FFAWE8HcgUD
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
信息提示
请填处理意见