24小时热门版块排行榜    

查看: 11923  |  回复: 274
当前主题已经存档。

wuli8

荣誉版主 (知名作家)

…………

优秀版主优秀版主

★ ★ ★ ★ ★
wuchenwf(金币+5,VIP+0):xiexie
Ask
引用回帖:
I am running wien2k version 03 on a machine with intel Pentium 4 processor, redhat linux 7 operating system, fortran compiler (intel) ifort version 9.1.036 and math library mkl 9.1.021. I am working on standard test case TiC. After successful compilation I did the initialization procedure.The master input files e.g TiC.struct & TiC.inst files were generated. Problem started when I run the SCF cycle. The case TiC is ano magnetic no spin one. During the run_lapw the following messages appeared:
cd/home/aruna/TiC;run_lapw-ec0.0001
LAPWO END
/home/aruna/wien2k/lapw1:error while loading shared libraries:libmkl_lapack.so: cannot open shared object file:no such file or directory
I used compiler options as follows
o:-FR -w mp1 prec_div -pad -ipo -i -static
l:-L/opt/intel/fc/9.1.021/lib -i-static
P:'-DParallel'R:-L/opt/intel/mkl/9.1.021/lib/32 -lmkl_lapack -lmkl_ia32 -lguide -lpthread
Would anyone please help to resolve this problem ?

Amswer
引用回帖:
Did you check that the mkl libraries are realy in opt/intel/mkl/9.1.021/lib/32
Most probably you need to use ldconfig to tell youre system where the shared libraries are

…………
141楼2009-01-04 10:12:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wax54

铜虫 (小有名气)

斑竹,用ifort10.1和mkl10.0.2.018编译08.3好像是可以用程序默认设置的,不过mkl一定要用10.0.2.018,我用10.0.5.025编译也会出错。
142楼2009-01-04 19:02:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wax54

铜虫 (小有名气)

麻烦大家帮我看看我试验的k点数和总能量的关系吧,我感觉不怎么对劲
500k(42)-53187.869514
800k(56)-53187.869545
1000k(56)-53187.869236
1500k(108)-53187.882254
2000k(135)-53187.865566
2500k(165)-53187.869175
括号中是ibz中的k点数,外边的是bz中的k点数
按照我的理解感觉大到一定程度就应该没什么区别了,可是为什么会出现波动的情况呢?
我的k点设置不是在较低精度基础上算的,是另算的,我看收敛也很快,这样会有什么问题吗?是不是检验精度高低,一定要在低精度基础上提升精度,那么提高到什么程度是可以信赖的呢?是不是也会出现我这样波动的情况呢?
143楼2009-01-05 09:30:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wuli8

荣誉版主 (知名作家)

…………

优秀版主优秀版主

★ ★ ★ ★ ★
wuchenwf(金币+5,VIP+0):xiexie
引用回帖:
麻烦大家帮我看看我试验的k点数和总能量的关系吧,我感觉不怎么对劲
500k(42)-53187.869514
800k(56)-53187.869545
1000k(56)-53187.869236
1500k(108)-53187.882254
2000k(135)-53187.865566
2500k(165)-53187.869175
括号中是ibz中的k点数,外边的是bz中的k点数
按照我的理解感觉大到一定程度就应该没什么区别了,可是为什么会出现波动的情况呢?
我的k点设置不是在较低精度基础上算的,是另算的,我看收敛也很快,这样会有什么问题吗?是不是检验精度高低,一定要在低精度基础上提升精度,那么提高到什么程度是可以信赖的呢?是不是也会出现我这样波动的情况呢?

总能量和k点有关还和你的收敛精度有关。
你的收敛精度低k点变化后总能量波动就很容易理解了。先提高精度再变化k点。
…………
144楼2009-01-05 13:47:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ylyan

木虫 (小有名气)

请教

我在使用wien2k的过程中,需要考虑自旋轨道耦合,可在计算过程中出现了下面的问题:
lapwso -p   (12:12:02) running LAPWSO in parallel mode
**  LAPWSO crashed!
1.260u 1.712s 0:43.32 6.8%      0+0k 0+0io 0pf+0w
error: command   /home/wang/wien2k/lapwsopara lapwso.def   failed
stop error
我昨天问了一下,有人回答说要在计算过程中修改inso,我不知怎么修改,也没有修改,请问是不是由于这个原因呀?请高手指点!先谢啦!!!
145楼2009-01-08 08:53:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wax54

铜虫 (小有名气)

求助
如何解决GGA+U叠代过程中总能量波动过大的问题?UG中提到过要加SO,可是我想看看加不加SO对体系电子结构有何影响。
146楼2009-01-08 09:48:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wysong

金虫 (小有名气)

★ ★
wuchenwf(金币+2,VIP+0):xiexie
引用回帖:
Originally posted by wax54 at 2009-1-8 09:48:
求助
如何解决GGA+U叠代过程中总能量波动过大的问题?UG中提到过要加SO,可是我想看看加不加SO对体系电子结构有何影响。

你怎老是说能量波动过大呢 那样你继续提高精度 不可以解决吗?
147楼2009-01-08 14:58:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wax54

铜虫 (小有名气)


wuchenwf(金币+1):xiexie讨论 1-17 19:57
谢谢各位了,我重新调整了RKmax和K点,现在已经解决了。
148楼2009-01-10 14:57:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wax54

铜虫 (小有名气)

[求助]关于UG中case.int中部分参数的设置
emin,
delta,
emax
specifies the energy mesh (in Ry) where the DOS is calculated. (emin
should be set slightly below the lowest valence band; emax will be
checked against the lowest energy of the highest band in case.qtl,
and set to the minimum of these two values; delta is the energy increment.
其中emin和emax到底应该怎样设置?
149楼2009-01-14 21:22:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lineng0724

木虫 (小有名气)

请教:wien2k计算如何编辑case.inop

case.inop

1200 1        number of k-points, first k-point
-5.0 2.2      Emin, Emax for matrix elements
2             number of choices (columns in *outmat) - 0: MME into case.mme
1             Re xx
3             Re zz
OFF           ON/OFF   writes MME to unit 4

Choices:
1......Re <;x>;<;x>;
2......Re <;y>;<;y>;
3......Re <;z>;<;z>;
4......Re <;x>;<;y>;
5......Re <;x>;<;z>;
6......Re <;y>;<;z>;
7......Im <;x>;<;y>;
8......Im <;x>;<;z>;
9......Im <;y>;<;z>;
150楼2009-01-16 10:27:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wuli8 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +7 ASdOkHsho7FD 2026-08-28 10/500 2026-08-31 00:42 by hZiFeudZyoGR
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +6 ASdOkHsho7FD 2026-08-28 8/400 2026-08-31 00:21 by hZiFeudZyoGR
[基金申请] 投票:  有多少人是今天查系统知道结果的? +17 爱看书的可乐 2026-08-26 19/950 2026-08-30 22:24 by xiangy672
[基金申请] 能否申诉? +4 echo8914667 2026-08-30 5/250 2026-08-30 21:11 by echo8914667
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 4FFAWE8HcgUD 2026-08-29 5/250 2026-08-30 20:41 by hZiFeudZyoGR
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 jCd0dEvKHShX 2026-08-29 5/250 2026-08-30 19:46 by hZiFeudZyoGR
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 jCd0dEvKHShX 2026-08-29 4/200 2026-08-30 19:34 by hZiFeudZyoGR
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 G6APbkg8SA6w 2026-08-29 4/200 2026-08-30 18:51 by hZiFeudZyoGR
[考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 G6APbkg8SA6w 2026-08-29 5/250 2026-08-30 18:40 by hZiFeudZyoGR
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 gy1nBQXYQJqL 2026-08-29 5/250 2026-08-30 17:56 by hZiFeudZyoGR
[论文投稿] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 ASdOkHsho7FD 2026-08-28 6/300 2026-08-30 17:25 by hZiFeudZyoGR
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 gy1nBQXYQJqL 2026-08-29 8/400 2026-08-30 12:14 by l0VvVHGBGRLv
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +6 gdd2018 2026-08-28 11/550 2026-08-30 08:58 by 超级无敌华子
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 06:57 by ZPa0EcMwuECS
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 06:10 by ZPa0EcMwuECS
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
信息提示
请填处理意见