24小时热门版块排行榜

返回列表

当前主题已经存档。

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

minmin_0082003

金虫 (正式写手)

莫愁莫愁

1ST强帖: 1
应助: 2 (幼儿园)
贵宾: 0.08
金币: 921.3
散金: 149
帖子: 910
在线: 218.4小时
虫号: 366989
注册: 2007-05-10
性别: MM
专业: 凝聚态物性I:结构、力学和

[交流] 【求助】吸附的分子又分解了，而且还不收敛，郁闷了，到底哪里错了

输入文件
NumberOfAtoms 33

NumberOfSpecies 5

%block ChemicalSpeciesLabel
      1    79    Au
      2    1    H
      3    6    C
      4    7    N
      5    8    O
%endblock ChemicalSpeciesLabel

PAO.BasisType          split
PAO.BasisSize SZ
PAO.EnergyShift 200 meV

LatticeConstant 8.6514 Ang

%block Latticeparameters
   1  1  3.128 90. 90. 120.
%endblock

AtomicCoordinatesFormat Ang

%block AtomicCoordinatesAndAtomicSpecies
  0.00000000 0.00000000 2.35460760 1
  1.44189679 0.83247950 0.00000000 1
  2.88379359 0.00000000 2.35460760 1
  4.32569038 0.83247950 0.00000000 1
  5.76758717 0.00000000 2.35460760 1
  7.20948396 0.83247950 0.00000000 1
-1.44189679 2.49743850 2.35460760 1
  0.00000000 3.32991801 0.00000000 1
  1.44189679 2.49743850 2.35460760 1
  2.88379359 3.32991801 0.00000000 1
  4.32569038 2.49743850 2.35460760 1
  5.76758717 3.32991801 0.00000000 1
-2.88379359 4.99487701 2.35460760 1
-1.44189679 5.82735651 0.00000000 1
  0.00000000 4.99487701 2.35460760 1
  1.44189679 5.82735651 0.00000000 1
  2.88379359 4.99487701 2.35460760 1
  4.32569038 5.82735651 0.00000000 1
  1.79669955 2.39693886 4.45247644 3
  1.97925927 2.34181771 5.59667128 4
  3.46881163 2.27267535 7.62901694 3
  4.32556961 2.32798369 6.97379283 2
  3.65112312 2.20500939 9.01545892 3
  4.64846405 2.20797205 9.42973770 2
  2.54058155 2.13331871 9.86487387 3
  1.24763727 2.12943163 9.32779258 3
  0.39087933 2.07412329 9.98301668 2
  1.06534096 2.19719237 7.94137912 3
  0.06798500 2.19413491 7.52707179 2
  2.17597398 2.26874631 7.09201869 3
  2.73715882 2.06026474 11.36024020 4
  1.77718841 1.99830441 12.09429137 5
  3.85452871 2.06370875 11.82451839 5
%endblock AtomicCoordinatesAndAtomicSpecies

%block GeometryConstraints
position 1
position  4
position  6
position  8
position  10
position  12
position  14
position  16
position  18
%endblock GeometryConstraints
kgrid_cutoff          30 Ang
%block kgrid_Monkhorst_pack
5 0 0 0.0
0 5 0 0.0
0 0 1 0.0
%endblock

XC.functional GGA
XC.authors    PBE

MeshCutoff             450 Ry
MaxSCFIterations       150
DM.NumberPulay 8
DM.NumberKick          15
DM.MixingWeight    0.001
DM.tolerance    5.d-5
SolutionMethod  diagon
ElectronicTemperature 100 K
MD.TypeOfRun          CG
MD.NumCGsteps    800
MD.MaxCGDispl          0.1 Ang
MD.MaxForceTol  0.003 eV/Ang
PAO.energy shift 50 meV
DM.AllowReuse             .true. # .true. is default
DM.AllowExtrapolation    .false. # .true. is default
DM.UseSaveDM             .false. # .false. is default
MD.UseSaveXV             .false.    # .false. is default

WriteCoorStep          .true.
WriteForces          .true.
WriteCoorXmol          .true.
WriteMDXmol          .true.

输出结构如下
  79    0.000000    0.000000    2.354608
  79    1.477225    0.843030 -0.278061
  79    2.896309 -0.025734    2.733369
  79    4.325690    0.832480    0.000000
  79    5.766163 -0.004577    2.765347
  79    7.209484    0.832480    0.000000
  79 -1.494599    2.500183    2.745271
  79    0.000000    3.329918    0.000000
  79    1.311405    2.540214    2.241525
  79    2.883794    3.329918    0.000000
  79    4.368079    2.492884    2.746128
  79    5.767587    3.329918    0.000000
  79 -2.885823    4.999363    2.762316
  79 -1.441897    5.827357    0.000000
  79 -0.018963    5.041639    2.738072
  79    1.441897    5.827357    0.000000
  79    2.902356    5.043507    2.739694
  79    4.325690    5.827357    0.000000
6    1.883956    2.365830    4.387311
7    1.945563    2.347358    5.485790
6    3.612296    2.278235    7.604465
1    4.584955    2.348022    6.731806
6    3.828853    2.211916    8.967603
1    4.907387    2.203624    9.650012
6    2.546443    2.134439    9.859743
6    1.100608    2.124144    9.353928
1    0.182998    2.057351 10.210496
6    0.908787    2.191538    7.968062
1 -0.207907    2.193948    7.409203
6    2.172007    2.269452    7.078789
7    2.724494    2.049535 11.577601
8    1.456706    1.981410 12.252734
8    4.204885    2.069137 11.891022

回复此楼

» 猜你喜欢

请问哪里可以有青B申请的本子可以借鉴一下。已经有3人回复
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗已经有3人回复
孩子确诊有中度注意力缺陷已经有14人回复
三甲基碘化亚砜的氧化反应已经有4人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
请问有评职称，把科研教学业绩算分排序的高校吗已经有5人回复
2025冷门绝学什么时候出结果已经有3人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复
AI论文写作工具：是科研加速器还是学术作弊器？已经有3人回复

» 本主题相关商家推荐: (我也要在这里推广)

1楼2010-01-22 12:26:35

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

hyb99

金虫 (初入文坛)

应助: 0 (幼儿园)
金币: 800.5
帖子: 11
在线: 9.4小时
虫号: 940360
注册: 2010-01-10
专业: 凝聚态物性 II ：电子结构

★ ★ ★ ★ ★
aylayl08(金币+5):谢谢指点 2010-01-28 14:29

引用回帖:

Originally posted by minmin_0082003 at 2010-01-28 09:43:02:

Hi，请问，建立MS的时候会出现类似的一行
PBC 5.7676 2.8838 17.3576 90 90 120（P1）
该如何转换成lattice vector的形式呢？
我转换成
Lattice constant 2.8838 ang
Latticeparameters 2 1 6 ...

有两个方法

1）你在ms中导出坐标是用res格式，这是一种内坐标。然后在siesta里面
把AtomicCoordinatesFormat设成ScaledByLatticeVectors就行了。

2）你可以在ms中选中晶格，然后在属性栏里面的最底下可以找到三个晶格矢量
然后在siesta里面写成如下形式
LatticeConstant    4.0655 Ang
%block LatticeVectors
   0.5000 0.50000 0.00000
   0.5000 0.00000 0.50000
   0.0000 0.50000 0.50000
%endblock LatticeVectors