应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Siesta&Smeagol&Atk » 请大家帮忙看看siesta输入文件有哪些错误
91/1返回列表
查看: 2323  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

yangxuezhang

金虫 (正式写手)

[求助] 请大家帮忙看看siesta输入文件有哪些错误

我这是第一次写输入文件，以前也没师兄做这个，也每个人问，希望高手不要嫌低级，帮我看一下，哪里有错误，不甚感激
SystemName Test for Ti3SiC2
SystemLabel Ti3SiC2

NumberOfAtoms           11            # Number of atoms
NumberOfSpecies         3           # Number of species

%block Chemical_Species_Label
  1    22  Ti
  2    6   C
  3    14  Si
%endblock Chemical_Species_Label

PAO.BasisSize      SZ   #???

# Lattice, coordinates, k-sampling

AtomicCoordinatesFormat ScaledCartesian # Format for coordinates
AtomicCoorFormatOut     Ang

%block AtomicCoordinatesAndAtomicSpecies
      0.000000000    0.000000000    0.000000000     1
      0.000000000    0.000000000    4.401250000     3
      0.000000000    0.000000000    13.20375        3
      1.541000000    0.889696765    1.271081005     2
      1.541000000    0.889696765    7.531419        2
      0.000000000    1.779393530    10.073581       2
      0.000000000    1.779393530    16.333919       2
      0.000000000    1.779393530    2.385477481     1
      0.000000000    1.779393530    6.4170225       1
      1.541000000    0.889696765    11.1879775      1
      1.541000000    0.889696765    15.2195225      1
%endblock AtomicCoordinatesAndAtomicSpecies

LatticeConstant     3.082   Ang

%block LatticeVectors
  1.000  0.500  0.000
  0.000  0.866  0.000
  0.000  0.000  1.000
%endblock LatticeVectors

kgrid_cutoff        7. Ang   #?????

# DFT, Grid, SCF

XC.functional           LDA         # Exchange-correlation functional type
XC.authors              CA          # Particular parametrization of xc func
SpinPolarized           .false.     # Spin unpolarized calculation
MeshCutoff              200. Ry     # Equivalent planewave cutoff for the grid
MaxSCFIterations        100         # Maximum number of SCF iterations per step
DM.MixingWeight         0.3         # New DM amount for next SCF cycle
DM.Tolerance            1.d-4       # Tolerance in maximum difference
                                    # between input and output DM
DM.NumberPulay          3           # Number of SCF steps between pulay mixing

# Eigenvalue problem: order-N or diagonalization

SolutionMethod          diagon      # OrderN or Diagon
ElectronicTemperature   5 K        # Temp. for Fermi smearing

# Molecular dynamics and relaxations

MD.TypeOfRun            cg          # Type of dynamics:

# Output options

WriteCoorInitial
WriteCoorStep
WriteForces
WriteKpoints            .false.
WriteEigenvalues        .false.
WriteKbands             .false.
WriteBands              .false.
WriteMullikenPop        1            # Write Mulliken Population Analysis
WriteCoorXmol           .false.
WriteMDCoorXmol         .false.
WriteMDhistory          .false.
WriteCoorXmol           .false.

# Options for saving/reading information

DM.UseSaveDM                         # Use DM Continuation files
MD.UseSaveXV            .false.      # Use stored positions and velocities
MD.UseSaveCG            .false.      # Use stored positions and velocities
SaveRho                              # Write valence pseudocharge at the mesh
SaveDeltaRho                         # Write RHOscf-RHOatm at the mesh
SaveElectrostaticPotential .false.   # Write the total elect. pot. at the mesh
SaveTotalPotential      .false.      # Write the total pot. at the mesh
WriteSiestaDim          .false.      # Write minimum dim to siesta.h and stop
WriteDenchar                         # Write information for DENCHAR
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-09-22 21:31:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangxuezhang

金虫 (正式写手)
为什么我的问题都木有人跟帖呢？
赞一下 回复此楼
2楼2012-09-23 15:30:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangxuezhang

金虫 (正式写手)
这是提交后的scrip.err文件里面的
ERROR STOP from Node:    0
--------------------------------------------------------------------------
MPI_ABORT was invoked on rank 0 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.
--------------------------------------------------------------------------
--------------------------------------------------------------------------
mpiexec has exited due to process rank 0 with PID 14863 on
node dell745-2-1.local exiting without calling "finalize". This may
have caused other processes in the application to be
terminated by signals sent by mpiexec (as reported here).
--------------------------------------------------------------------------

real    0m0.866s
user    0m0.041s
sys     0m0.054s
赞一下 回复此楼
3楼2012-09-23 15:46:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangguangping

木虫 (著名写手)
★ ★
franch: 金币+2, 谢谢回帖交流 2012-09-25 09:59:38
给出你的标准输出的最后的出错信息。这些是MPI的，看不出来。
赞一下 回复此楼
弘德明志博学笃行
4楼2012-09-24 21:04:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangxuezhang

金虫 (正式写手)
引用回帖:
4楼: Originally posted by zhangguangping at 2012-09-24 21:04:14
给出你的标准输出的最后的出错信息。这些是MPI的，看不出来。

下面是*.out文件里面的，为了找错误，和前面的输入文件比，我注释掉了很多项，不过还是不行，不知道是什么问题
Siesta Version:                                        siesta-3.1
Architecture  : x86_64-unknown-linux-gnu--Intel
Compiler flags: /opt/openmpi-intel11/bin/mpif90 -g -O2
PARALLEL version

* Running on    2 nodes in parallel
>> Start of run:  24-SEP-2012  15:08:55

                           ***********************
                           *  WELCOME TO SIESTA  *
                           ***********************

reinit: Reading from standard input
************************** Dump of input data file ****************************
SystemName Test for Ti3SiC2
SystemLabel Ti3SiC2
NumberOfAtoms           11            # Number of atoms
NumberOfSpecies         3           # Number of species
%block Chemical_Species_Label
  1    22  Ti
  2    6   C
  3    14  Si
%endblock Chemical_Species_Label
#PAO.BasisSize      SZ   #???
# Lattice, coordinates, k-sampling
AtomicCoordinatesFormat Ang  # Format for coordinates
AtomicCoorFormatOut     Ang
%block AtomicCoordinatesAndAtomicSpecies
      0.000000000    0.000000000    0.000000000     1
      0.000000000    0.000000000    4.401250000     3
      0.000000000    0.000000000    13.20375        3
      1.541000000    0.889696765    1.271081005     2
      1.541000000    0.889696765    7.531419        2
      0.000000000    1.779393530    10.073581       2
      0.000000000    1.779393530    16.333919       2
      0.000000000    1.779393530    2.385477481     1
      0.000000000    1.779393530    6.4170225       1
      1.541000000    0.889696765    11.1879775      1
      1.541000000    0.889696765    15.2195225      1
%endblock AtomicCoordinatesAndAtomicSpecies
LatticeConstant     3.082   Ang
%block LatticeVectors
  1.000  0.500  0.000
  0.000  0.866  0.000
  0.000  0.000  1.000
%endblock LatticeVectors
#kgrid_cutoff        7. Ang   #?????
# DFT, Grid, SCF
#XC.functional           LDA         # Exchange-correlation functional type
#XC.authors              CA          # Particular parametrization of xc func
#SpinPolarized           .false.     # Spin unpolarized calculation
#MeshCutoff              200. Ry     # Equivalent planewave cutoff for the grid
#MaxSCFIterations        100         # Maximum number of SCF iterations per step
#DM.MixingWeight         0.3         # New DM amount for next SCF cycle
#DM.Tolerance            1.d-4       # Tolerance in maximum difference
                                    # between input and output DM
#DM.NumberPulay          3           # Number of SCF steps between pulay mixing
# Eigenvalue problem: order-N or diagonalization
SolutionMethod          diagon      # OrderN or Diagon
#ElectronicTemperature   5 K        # Temp. for Fermi smearing
# Molecular dynamics and relaxations
#MD.TypeOfRun            cg          # Type of dynamics:
# Output options
#WriteCoorInitial
#WriteCoorStep
#WriteForces
#WriteKpoints            .false.
#WriteEigenvalues        .false.
#WriteKbands             .false.
#WriteBands              .false.
#WriteMullikenPop        1            # Write Mulliken Population Analysis
#WriteCoorXmol           .false.
#WriteMDCoorXmol         .false.
#WriteMDhistory          .false.
#WriteCoorXmol           .false.
# Options for saving/reading information
#DM.UseSaveDM                         # Use DM Continuation files
#MD.UseSaveXV            .false.      # Use stored positions and velocities
#MD.UseSaveCG            .false.      # Use stored positions and velocities
#SaveRho                              # Write valence pseudocharge at the mesh
#SaveDeltaRho                         # Write RHOscf-RHOatm at the mesh
#SaveElectrostaticPotential .false.   # Write the total elect. pot. at the mesh
#SaveTotalPotential      .false.      # Write the total pot. at the mesh
#WriteSiestaDim          .false.      # Write minimum dim to siesta.h and stop
#WriteDenchar                         # Write information for DENCHAR
************************** End of input data file *****************************

reinit: -----------------------------------------------------------------------
reinit: System Name: Test for Ti3SiC2
reinit: -----------------------------------------------------------------------
reinit: System Label: Ti3SiC2
reinit: -----------------------------------------------------------------------

initatom: Reading input for the pseudopotentials and atomic orbitals ----------
Species number:            1  Label: Ti Atomic number:          22
Species number:            2  Label: C Atomic number:           6
Species number:            3  Label: Si Atomic number:          14
Ground state valence configuration:   4s02  3d02
Reading pseudopotential information in formatted form from Ti.psf

Pseudopotential generated from a relativistic atomic calculation
There are spin-orbit pseudopotentials available
Spin-orbit interaction is not included in this calculation

Valence configuration for pseudopotential generation:
3s( 2.00) rc: 1.30
3p( 6.00) rc: 1.30
3d( 2.00) rc: 1.30
4f( 0.00) rc: 1.98
Ground state valence configuration:   2s02  2p02
Reading pseudopotential information in formatted form from C.psf

Valence configuration for pseudopotential generation:
2s( 2.00) rc: 1.49
2p( 2.00) rc: 1.52
3d( 0.00) rc: 1.58
Ground state valence configuration:   3s02  3p02
Reading pseudopotential information in formatted form from Si.psf

Valence configuration for pseudopotential generation:
3s( 2.00) rc: 1.89
3p( 2.00) rc: 1.89
3d( 0.00) rc: 1.89
4f( 0.00) rc: 1.89
Semicore shell(s) with  8 electrons included in the valence for
Ti
Ti                   must be in PAO.Basis (it has semicore states)
ERROR STOP from Node:    0
希望大侠多多帮忙
赞一下 回复此楼
5楼2012-09-26 14:53:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangguangping

木虫 (著名写手)

【答案】应助回帖

★ ★
yangxuezhang: 金币+2, ★有帮助 2012-09-30 15:17:02
引用回帖:
5楼: Originally posted by yangxuezhang at 2012-09-26 14:53:20
下面是*.out文件里面的，为了找错误，和前面的输入文件比，我注释掉了很多项，不过还是不行，不知道是什么问题
Siesta Version:                                        siesta-3.1
Architecture  : x86_64-unk ...

从输出文件中看出来，你计算中的Ti元素的赝势将8个本应该属于core的电子放到了价电子的考虑范围，这样可以更好的描述不同环境中的Ti的电子结构。但是你需要对Ti的基函数进行手动设置。而不能用简单的PAO.BasisSize      SZ关键字，使用这个关键字默认只给价电子产生基函数。所以你需要使用
%block PAO.Basis
%endblock PAO.Basis
来设置，比如
PAO.BasisSize      DZP
PAO.EnergyShift       100 meV
对应的C元素的基函数为：
%block PAO.Basis
C                     2                    # Species label, number of l-shells
n=2   0   2                         # n, l, Nzeta
   4.633      3.432   
   1.000      1.000   
n=2   1   2 P   1                   # n, l, Nzeta, Polarization, NzetaPol
   5.659      3.653   
   1.000      1.000
%endblock PAO.Basis
当然，这个与具体的赝势还有点关系。
你可以比照着这个对Ti进行一下设置，将那8个电子的基函数明确的写出来。我没有算过Ti,所以不会设置。  

还有一个办法就是你可以使用一个没有semicore态电子的Ti赝势。你从官网上下载一个就是了。我看了一下官网上的LDA和GGA的Ti赝势就都没有semicore态。
http://icmab.cat/leem/siesta/Dat ... init/Tl_html/Tl.psf
http://icmab.cat/leem/siesta/Dat ... init/Tl_html/Tl.psf
赞一下 回复此楼
弘德明志博学笃行
6楼2012-09-26 15:22:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangxuezhang

金虫 (正式写手)
引用回帖:
6楼: Originally posted by zhangguangping at 2012-09-26 15:22:26
从输出文件中看出来，你计算中的Ti元素的赝势将8个本应该属于core的电子放到了价电子的考虑范围，这样可以更好的描述不同环境中的Ti的电子结构。但是你需要对Ti的基函数进行手动设置。而不能用简单的PAO.BasisSize ...

我是刚入门，我也不会设置Ti，我试了你给的两种赝势，另外还改了
PAO.BasisSize      DZP   #???
但是还是不行，错误文件
Semicore shell(s) with -1 electrons included in the valence for
Ti
Ti                   must be in PAO.Basis (it has semicore states)
ERROR STOP from Node:    0
赞一下 回复此楼
7楼2012-09-27 15:59:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tic2009

金虫 (初入文坛)
★ ★
yangxuezhang: 回帖置顶 2012-09-30 15:18:54
yangxuezhang: 取消置顶 2012-09-30 15:20:06
liliangfang: 金币+2, 谢谢交流 2012-09-30 18:34:55
您使用的赝势含有semi-core states，基组必须显式表达，利用PAO.Basis模块来定义
赞一下 回复此楼
西部二本高校青年教师
8楼2012-09-30 12:00:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangxuezhang

金虫 (正式写手)
引用回帖:
8楼: Originally posted by tic2009 at 2012-09-30 12:00:09
您使用的赝势含有semi-core states，基组必须显式表达，利用PAO.Basis模块来定义

谢谢你的回答，我是刚入门的，PAO.Basis这块还不会定义，所以不用显式，不知道为什么不可以给金币了？
赞一下 回复此楼
9楼2012-09-30 15:16:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yangxuezhang 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[精细化工] 招聘 金属平磨液，抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 有时候，自然基金真的不能太认真 （我的申报经验） +9 majunge000 2026-08-11 11/550 2026-08-15 23:20 by botar
[基金申请] 时间戳又变了8-15 +12 archvillain 2026-08-15 23/1150 2026-08-15 20:47 by foyo404
[基金申请] 小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊 +12 Tide man 2026-08-10 13/650 2026-08-15 16:34 by 氺木
[基金申请] 关于Filecode分析方法 +10 majunge000 2026-08-10 13/650 2026-08-15 12:26 by hanpeng972
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 22/1100 2026-08-14 23:45 by Noways
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +23 Lanmanbaby 2026-08-09 37/1850 2026-08-14 13:40 by Equinoxhua
[硕博家园] 请教兼职经验 +3 是阿文鸭 2026-08-09 3/150 2026-08-14 12:07 by HER12025
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] FileCode能看出啥？ +10 要乐观耀哥 2026-08-10 32/1600 2026-08-14 09:37 by 要乐观耀哥
[基金申请] 我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人 +10 家与远方 2026-08-10 15/750 2026-08-14 02:08 by 绵羊哥哥
[基金申请] 重要来源：本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] Filecode 又变了，巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 结合人工智能，周易传统文化，filecode打分制来了，3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见，大家看看几个A，几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分？ +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 什么时候出结果，有咨询渠道？？？ +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
信息提示
关闭
请填处理意见
关闭