24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1562  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 今年审到国自然15份,谈谈感受 +15 国自然国社科中 2026-05-17 15/750 2026-05-18 14:17 by 能丫
[硕博家园] 考博自荐 +5 科研狗111 2026-05-13 6/300 2026-05-18 11:22 by 糊糊涂涂好
[基金申请] 青C资助名额大幅增加! +12 西葫芦炒鸡蛋 2026-05-13 16/800 2026-05-18 10:02 by Equinoxhua
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 xx7gd5zq4e 2026-05-15 6/300 2026-05-17 19:36 by Equinoxhua
[公派出国] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 v9tggjlwd0 2026-05-15 5/250 2026-05-17 18:39 by Equinoxhua
[考博] 光量子物理方向 博士招生 1人(2026.09) +3 sandyworld 2026-05-15 4/200 2026-05-17 14:38 by sandyworld
[基金申请] 这年头没有找到涵评专家,还有中面上的可能吗 +13 dd921ww 2026-05-12 16/800 2026-05-17 12:38 by 于轩
[考博] 26/27博士推荐 +3 1木头人13949 2026-05-13 3/150 2026-05-17 09:41 by YuY66
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 v9tggjlwd0 2026-05-15 4/200 2026-05-17 08:06 by 11n4dfd8yn
[考博] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 7hlccowb3h 2026-05-15 4/200 2026-05-17 07:46 by 11n4dfd8yn
[考博] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +6 l7k6xnh0yc 2026-05-14 6/300 2026-05-17 07:11 by 11n4dfd8yn
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +5 cjf4bx70cj 2026-05-14 7/350 2026-05-17 06:55 by 11n4dfd8yn
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 k37jurhrau 2026-05-16 3/150 2026-05-17 01:25 by ue3ir18jc3
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 ky2p12rrjj 2026-05-15 4/200 2026-05-17 00:50 by ue3ir18jc3
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 x0mp7owy2b 2026-05-15 4/200 2026-05-17 00:25 by ue3ir18jc3
[有机交流] 求有机合成大神指点三硫酸乙烯酯(CAS:2793408-99-6)的合成路线 30+3 Leekmid 2026-05-13 10/500 2026-05-16 16:37 by czyzsu
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 k37jurhrau 2026-05-16 3/150 2026-05-16 13:57 by vcdazktkjx
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 l7k6xnh0yc 2026-05-14 6/300 2026-05-16 11:29 by h3oerqvkv9
[硕博家园] 申请博士 +3 呃?呃 2026-05-15 3/150 2026-05-16 11:01 by a4742549
[文学芳草园] 风把牡丹吹跑了 +5 myrtle 2026-05-12 9/450 2026-05-15 15:27 by myrtle
信息提示
请填处理意见