24小时热门版块排行榜    

查看: 1653  |  回复: 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的回帖

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的回帖
查看全部 10 个回答

匿名

用户注销 (职业作家)

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

匿名

用户注销 (职业作家)

本帖仅楼主可见
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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 ASdOkHsho7FD 2026-08-28 4/200 2026-08-29 09:48 by G6APbkg8SA6w
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 ASdOkHsho7FD 2026-08-28 5/250 2026-08-29 09:39 by G6APbkg8SA6w
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 ASdOkHsho7FD 2026-08-28 7/350 2026-08-29 09:38 by G6APbkg8SA6w
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 ASdOkHsho7FD 2026-08-28 4/200 2026-08-29 05:45 by gy1nBQXYQJqL
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +13 晴天加油 2026-08-26 14/700 2026-08-28 18:36 by huagongfeihu
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 有没有仍没收到信息的 +3 德尚中行 2026-08-27 3/150 2026-08-28 15:24 by rongshuxia
[基金申请] 基金不中,共勉 +11 eulota 2026-08-26 11/550 2026-08-28 14:22 by 火星超人xi
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +4 gdd2018 2026-08-28 9/450 2026-08-28 14:22 by jklily
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
信息提示
请填处理意见