24小时热门版块排行榜    

查看: 1237  |  回复: 4

小明好开心

新虫 (初入文坛)

[求助] 关于使用siesta计算ordern的参数设置和问题解决 已有1人参与

siesta新手想试一下做ordern的计算,修改参数成 SolutionMethod=ordern后,计算了TiN体系的能量。返回了这样的错误。求大神指导!!需要设置其他的参数吗?ordern计算需要注意什么问题呢。手册上有句话是说ordern subsystem is quite fragile and only works for systems with clearly separated occupied and empty states. Note also that the option to compute the potential automatically does not yet work in parallel.
下面是返回的错误提示。
Cannot do O(N) with k-points.
Stopping Program from Node:    3
Cannot do O(N) with k-points.
Stopping Program from Node:    3
Cannot do O(N) with k-points.
Stopping Program from Node:    2
Cannot do O(N) with k-points.
Stopping Program from Node:    2
--------------------------------------------------------------------------
MPI_ABORT was invoked on rank 3 in communicator MPI_COMM_WORLD
with errorcode 1.

NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
You may or may not see output from other processes, depending on
exactly when Open MPI kills them.
--------------------------------------------------------------------------
forrtl: error (78): process killed (SIGTERM)
Image              PC                Routine            Line        Source            
libpthread.so.0    00002B32BCC635AD  Unknown               Unknown  Unknown
siesta             00000000014560DB  Unknown               Unknown  Unknown
siesta             0000000001455903  Unknown               Unknown  Unknown
siesta             000000000056E409  Unknown               Unknown  Unknown
siesta             00000000006A06A8  Unknown               Unknown  Unknown
siesta             000000000054223B  Unknown               Unknown  Unknown
siesta             000000000056D501  Unknown               Unknown  Unknown
siesta             000000000079A07F  Unknown               Unknown  Unknown
siesta             0000000000418ECC  Unknown               Unknown  Unknown
libc.so.6          00002B32BDA97CDD  Unknown               Unknown  Unknown
siesta             0000000000418DC9  Unknown               Unknown  Unknown
forrtl: severe (28): CLOSE error, unit 13, file "Unknown"
Image              PC                Routine            Line        Source            
siesta             00000000014FE11A  Unknown               Unknown  Unknown
siesta             00000000014FCC16  Unknown               Unknown  Unknown
siesta             00000000014A9D70  Unknown               Unknown  Unknown
siesta             000000000145EE0E  Unknown               Unknown  Unknown
siesta             00000000014625C1  Unknown               Unknown  Unknown
siesta             00000000014641F2  Unknown               Unknown  Unknown
libpthread.so.0    00002B32BCC64500  Unknown               Unknown  Unknown
libpthread.so.0    00002B32BCC635AD  Unknown               Unknown  Unknown
siesta             00000000014560DB  Unknown               Unknown  Unknown
siesta             0000000001455903  Unknown               Unknown  Unknown
siesta             000000000056E409  Unknown               Unknown  Unknown
siesta             00000000006A06A8  Unknown               Unknown  Unknown
siesta             000000000054223B  Unknown               Unknown  Unknown
siesta             000000000056D501  Unknown               Unknown  Unknown
siesta             000000000079A07F  Unknown               Unknown  Unknown
siesta             0000000000418ECC  Unknown               Unknown  Unknown
libc.so.6          00002B32BDA97CDD  Unknown               Unknown  Unknown
siesta             0000000000418DC9  Unknown               Unknown  Unknown
Cannot do O(N) with k-points.
Stopping Program from Node:    1
Cannot do O(N) with k-points.
Stopping Program from Node:    1
[node08:16898] 2 more processes have sent help message help-mpi-api.txt / mpi-abort
[node08:16898] Set MCA parameter "orte_base_help_aggregate" to 0 to see all help / error messages
下面试我的输入文件
SystemName             TiN   
SystemLabel            TiN           
NumberOfSpecies        2      
NumberOfAtoms          8      
%block ChemicalSpeciesLabel
  1  22  Ti
  2  7   N
%endblock ChemicalSpeciesLabel
%block     PAO.BasisSize         
  N      DZP
%endblock  PAO.BasisSize
%block PAO.Basis
Ti    5      1.91
n=3    0    1   E     93.95      5.20
   5.69946662616249
   1.00000000000000
n=3    1    1   E     95.47      5.20
   5.69941339465994
   1.00000000000000
n=4    0    2   E     96.47      5.60
   6.09996398975307        5.09944363262274
   1.00000000000000        1.00000000000000
n=3    2    2   E     46.05      4.95
   5.94327035784617        4.70009988294302
   1.00000000000000        1.00000000000000
n=4    1    1   E      0.50      1.77
   3.05365979938936
   1.00000000000000
%endblock PAO.Basis
LatticeConstant     4.30   Ang
%block LatticeVectors
   1.000  0.000  0.000
   0.000  1.000  0.000
   0.000  0.000  1.000
%endblock LatticeVectors
MeshCutoff         150.0 Ry
DM.MixingWeight      0.3
DM.NumberPulay       3  
DM.Tolerance         1.d-4
kgridcutoff         10. Ang
ElectronicTemperature  25 meV
xc.functional         GGA           
xc.authors            PBE   
MD.TypeOfRun          CG        
SolutionMethod        OrderN
WriteEigenvalues      true
AtomicCoordinatesFormat   Fractional   
%block AtomicCoordinatesAndAtomicSpecies
     0.000000000         0.000000000         0.000000000     1
     0.000000000         0.500000000         0.500000000     1
     0.500000000         0.000000000         0.500000000     1
     0.500000000         0.500000000         0.000000000     1
     0.500000000         0.500000000         0.500000000     2
     0.500000000         0.000000000         0.000000000     2
     0.000000000         0.500000000         0.000000000     2
     0.000000000         0.000000000         0.500000000     2
%endblock AtomicCoordinatesAndAtomicSpecies
COOP.write T
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

charleslian

木虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
OrderN只能用一个gamma点,说明书应该挺清楚了。

发自小木虫IOS客户端
2楼2015-09-24 09:54:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小明好开心

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by charleslian at 2015-09-24 09:54:06
OrderN只能用一个gamma点,说明书应该挺清楚了。

那请问是我的输入文件里设置多于一个gamma点才出现错误吗?怎么知道我设置了几个gamma点?
3楼2015-09-24 16:24:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

charleslian

木虫 (小有名气)

【答案】应助回帖

Kgridcutoff 那一项,实际上对应了多个k点

发自小木虫IOS客户端
4楼2015-09-25 12:14:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小明好开心

新虫 (初入文坛)

引用回帖:
4楼: Originally posted by charleslian at 2015-09-25 12:14:46
Kgridcutoff 那一项,实际上对应了多个k点

请问如何修改呢?把Kgridcutoff改成1就行吗?听导师说可能单胞太小了。但我看看官方给出的例子有一个计算H2O ordern的,那个原子数也少。还有kridcutoff和kgridMonkhorstpack 有什么不一样 ?请高人指教。下面附上H20 ORDERN 的例子fdf文件。这个可以算出来。

SystemName          Water molecule, using order-N method
SystemLabel         h2o
NumberOfAtoms       3
NumberOfSpecies     2

%block ChemicalSpeciesLabel
1  8  O      # Species index, atomic number, species label
2  1  H
%endblock ChemicalSpeciesLabel

AtomicCoordinatesFormat  Ang
%block AtomicCoordinatesAndAtomicSpecies
0.000  0.000  0.000  1
0.757  0.586  0.000  2
-0.757  0.586  0.000  2
%endblock AtomicCoordinatesAndAtomicSpecies

SolutionMethod   ordern
5楼2015-09-25 15:12:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 小明好开心 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 9g0rmhtq5w 2026-06-27 4/200 2026-06-28 03:28 by jjz01vuwzk
[博后之家] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 9g0rmhtq5w 2026-06-27 5/250 2026-06-28 03:27 by jjz01vuwzk
[公派出国] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 9g0rmhtq5w 2026-06-27 3/150 2026-06-28 03:27 by jjz01vuwzk
[考博] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 9g0rmhtq5w 2026-06-27 4/200 2026-06-28 03:25 by jjz01vuwzk
[公派出国] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 dheqf6rz6n 2026-06-27 3/150 2026-06-28 03:09 by jjz01vuwzk
[硕博家园] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +4 dheqf6rz6n 2026-06-27 5/250 2026-06-28 03:08 by jjz01vuwzk
[找工作] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 dheqf6rz6n 2026-06-27 5/250 2026-06-28 03:07 by jjz01vuwzk
[公派出国] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 dheqf6rz6n 2026-06-27 3/150 2026-06-28 02:49 by jjz01vuwzk
[考研] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 dheqf6rz6n 2026-06-26 7/350 2026-06-28 02:12 by jjz01vuwzk
[考研] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 7s8du2bt8y 2026-06-26 6/300 2026-06-28 01:46 by jjz01vuwzk
[考研] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +4 7s8du2bt8y 2026-06-26 8/400 2026-06-28 01:29 by jjz01vuwzk
[基金申请] 时间戳又变了 +4 Vivilian 2026-06-27 4/200 2026-06-28 00:37 by jnhyjjm
[基金申请] 无聊看看filecode +7 流流伤 2026-06-23 10/500 2026-06-27 16:19 by 北京中袖冉老师
[论文投稿] 求推荐期刊,重谢 +4 girlbaby 2026-06-23 4/200 2026-06-26 16:52 by 不打工牛马
[有机交流] 求助!! 5+3 我啥都没看见 2026-06-24 4/200 2026-06-26 09:35 by 951037019
[论文投稿] 职称论文投稿 170+4 guoj5292 2026-06-22 12/600 2026-06-25 22:27 by Yanyanoo
[有机交流] 反应求助 10+3 slz_1986 2026-06-24 6/300 2026-06-25 21:38 by nBu锂
[基金申请] 2026年WR青拔进展 +5 chs564851482 2026-06-24 7/350 2026-06-24 18:15 by chs564851482
[基金申请] 咨询 +3 _xyan818 2026-06-24 3/150 2026-06-24 08:12 by Equinoxhua
[基金申请] 评委有多少概率知道其他专家手中有哪些人的本子? +6 huitong441 2026-06-22 6/300 2026-06-23 15:45 by 新城子曾
信息提示
请填处理意见