| 查看: 1975 | 回复: 8 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
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 |
回复此楼» 猜你喜欢
博士读完未来一定会好吗
已经有6人回复
-
小论文投稿
已经有3人回复
-
Bioresource Technology期刊，第一次返修的时候被退回好几次了
已经有9人回复
-
心脉受损
已经有3人回复
-
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有8人回复
-
申请2026年博士
已经有6人回复
-
请问哪里可以有青B申请的本子可以借鉴一下。
已经有5人回复
-
2025冷门绝学什么时候出结果
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
我的输入文件 一开始就出现错误 怎么办呢、。。。。。
已经有6人回复- 请教高斯g03输入文件错误
已经有5人回复
- 【求助】程序运行已经无错误，但是输入文件有问题，求教。
已经有23人回复
- 【求助】新手学习siesta，出现错误，现求助
已经有4人回复
- 【求助】siesta 计算DOS应该注意什么问题
已经有3人回复
- 【求助】siesta输入文件中赝势文件的作用？
已经有8人回复
- 【求助】帮忙解决晶体CIF文件检查中的错误，非常感谢！
已经有5人回复
- 【求助】siesta3.0 安装，错误一箩筐
已经有11人回复
- 【求助】求助解决用mpich2.1 编译siesta，计算出现错误，
已经有12人回复
- 【求助】siesta使用B3lyp但是没有输出文件，所以找不出错误
已经有6人回复
- 【求助】siesta输入文件错误，昨晚到现在没有弄好，请求帮助
已经有7人回复
zhangguangping
木虫 (著名写手)
- 1ST强帖: 23
- 应助: 71 (初中生)
- 贵宾: 0.031
- 金币: 1973.1
- 散金: 9059
- 红花: 79
- 帖子: 2681
- 在线: 2128.8小时
- 虫号: 529624
- 注册: 2008-03-20
- 性别: GG
- 专业: 原子和分子物理
4楼2012-09-24 21:04:14
zhangguangping
木虫 (著名写手)
- 1ST强帖: 23
- 应助: 71 (初中生)
- 贵宾: 0.031
- 金币: 1973.1
- 散金: 9059
- 红花: 79
- 帖子: 2681
- 在线: 2128.8小时
- 虫号: 529624
- 注册: 2008-03-20
- 性别: GG
- 专业: 原子和分子物理
【答案】应助回帖
★ ★
yangxuezhang: 金币+2, ★有帮助 2012-09-30 15:17:02
yangxuezhang: 金币+2, ★有帮助 2012-09-30 15:17:02
|
从输出文件中看出来，你计算中的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