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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Siesta&Smeagol&Atk » 请大家帮忙看看siesta输入文件有哪些错误
51/1返回列表
查看: 1973  |  回复: 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

金虫 (正式写手)
引用回帖:
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的回帖
查看全部 9 个回答

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的回帖
查看全部 9 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭