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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Siesta&Smeagol&Atk » 请大家帮忙看看siesta输入文件有哪些错误
91/1返回列表
查看: 2120  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 312求调剂 +8 陌宸希 2026-03-16 9/450 2026-03-18 12:39 by Linda Hu
[考研] 一志愿天津大学化学工艺专业（081702）315分求调剂 +9 yangfz 2026-03-17 9/450 2026-03-18 12:38 by 尽舜尧1
[考研] 【0703化学调剂】-一志愿华中师范大学-六级475 +6 Becho359 2026-03-11 6/300 2026-03-18 12:09 by djl2006
[考研] 0703化学调剂 +3 妮妮ninicgb 2026-03-17 3/150 2026-03-18 10:29 by macy2011
[考研] 生物学071000 329分求调剂 +3 我爱生物生物爱� 2026-03-17 3/150 2026-03-18 10:12 by macy2011
[考研] 290求调剂 +3 p asserby. 2026-03-15 4/200 2026-03-17 16:35 by wangkm
[考研] 材料与化工专硕调剂 +5 heming3743 2026-03-16 5/250 2026-03-17 14:03 by 勇敢太监王公公
[考研] 11408 一志愿西电，277分求调剂 +3 zhouzhen654 2026-03-16 3/150 2026-03-17 07:03 by laoshidan
[考研] 333求调剂 +3 文思客 2026-03-16 7/350 2026-03-16 18:21 by 文思客
[考研] 中科院材料273求调剂 +4 yzydy 2026-03-15 4/200 2026-03-16 15:59 by Gaodh_82
[考研] 0703化学调剂，求各位老师收留 +8 秋有木北 2026-03-14 8/400 2026-03-16 15:21 by 哦哦123
[考研] 070305求调剂 +3 mlpqaz03 2026-03-14 4/200 2026-03-15 11:04 by peike
[考研] 085601材料工程315分求调剂 +3 yang_0104 2026-03-15 3/150 2026-03-15 10:58 by peike
[考研] 297求调剂 +4 学海漂泊 2026-03-13 4/200 2026-03-14 11:51 by 热情沙漠
[考研] 招收0805（材料）调剂 +3 18595523086 2026-03-13 3/150 2026-03-14 00:33 by 123%、
[考研] 315求调剂 +9 小羊小羊_ 2026-03-11 10/500 2026-03-13 21:13 by SXNU李老师
[考研] 材料与化工085600调剂求老师收留 +9 jiaanl 2026-03-11 9/450 2026-03-13 20:22 by JourneyLucky
[考研] 0703化学求调剂 +7 绿豆芹菜汤 2026-03-12 7/350 2026-03-13 17:25 by njzyff
[考研] 求b区学校调剂 +3 周56 2026-03-11 3/150 2026-03-13 16:20 by JourneyLucky
[考研] 材料专硕350 求调剂 +4 王金科 2026-03-12 4/200 2026-03-13 16:02 by ruiyingmiao
信息提示
关闭
请填处理意见
关闭