版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

返回列表

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

minmin_0082003

金虫 (正式写手)

[交流] 【求助】优化的结构非常奇怪，Au10 簇被优化成了4组小簇，为什么，输入及结构已贴出

用的网上贴出来的pesudopotential
输入文件如下
&CONTROL
  calculation = "relax",
  prefix       = "iso0-Au10-t"
  nstep       = 150,
  pseudo_dir = "/home/***/PWscf/Pseudo",
  outdir    = "/home/***/PWscf/outdir",
/
&SYSTEM
  nosym    = .TRUE.,
  ibrav    = 0,
  celldm(1) = 1.0D0,
  nat       = 10,
  ntyp       = 1,
  ecutwfc    = 25.0D0,
  ecutrho    = 160.0D0,
  occupations = "smearing",
  degauss    = 0.01D0,
  smearing = "gauss",
/
&ELECTRONS
  conv_thr = 1.D-6,
  mixing_beta = 0.5D0,
/
&IONS
  bfgs_ndim       = 3,
  upscale          = 100.D0,
  pot_extrapolation = "second_order",
  wfc_extrapolation = "second_order",
/
ATOMIC_SPECIES
Au  196.96655  Au.blyp-mt.UPF
ATOMIC_POSITIONS { angstrom }
Au  0.351644000    0.323076000    0.020090000
Au  0.353167000    -2.459480000    0.001629000
Au  0.325466000    -5.209497000    0.021409000
Au -1.052830000    -1.054437000    1.988954000
Au -1.026740000    1.696543000    1.965435000
Au -1.051203000    -3.836601000    1.970273000
Au  1.727775000    -1.054277000    -1.924117000
Au  1.723939000    -3.829706000    -1.938591000
Au -2.418980000    0.317707000    3.931133000
Au -2.429430000    -2.457304000    3.913504000
K_POINTS {Gamma}
CELL_PARAMETERS { cubic }
  15.00000  0.00000  0.00000
0.00000  15.00000  0.00000
0.00000  0.00000  15.00000

输入文件是从高斯优化得到的

输出文件时PWscf得到的

回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

优化聚合物结构时，是摆成环形还是应该摆成直线形已经有5人回复
ubuntu下载下来提示安装后什么都没安装就可以启动，感觉奇怪啊。。。已经有21人回复
用w2web界面计算性质之前的结构优化已经有7人回复
求助解决活性炭浓硫酸改性发生的奇怪现象已经有12人回复
正交优化实验是怎么设计的啊？已经有16人回复
B2被要求补交research plan，好奇怪，求交流！已经有13人回复
castep固定体积的结构优化和计算声子谱的问题已经有25人回复
求助链霉菌培养及优化已经有5人回复
CASTEP对优化后的晶体进行分态密度分析，应该怎么操作已经有14人回复
【关注SCI期刊点评】系统优化意见征集中~~~~~~~~~~ 已经有220人回复
【求助】MS4.4 优化失败，但是此结构放在MS2.2却能优化成功已经有14人回复
【求助】求解：geom=connectivity对过渡金属簇合物优化结果有什么影响已经有8人回复
【求助】dmol3优化计算团簇的温度应该是多少已经有4人回复

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2010-12-17 11:39:54

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

minmin_0082003

金虫 (正式写手)

★
ellsaking(金币+1):写一写出了什么问题也让大家借鉴一下吧，哈哈~ 2010-12-22 09:57:08

原因已经找到，问题初步解决

赞一下(1人)

回复此楼

3楼2010-12-19 14:04:25

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

查看全部 13 个回答

hakuna

木虫 (知名作家)

1ST强帖: 3
应助: 82 (初中生)
金币: 8735.9
帖子: 5599
在线: 1367.2小时
虫号: 167622

★ ★
zzy870720z(金币+2):鼓励交流 2010-12-19 11:41:23
minmin_0082003(金币+3):谢谢，原因已经解决，是因为pseudopotential的问题，但是为什么configure 不同得到结果不同还不清楚 2010-12-19 14:04:01

个人感觉初始结构不合理，你把六边形外面的两个Au放在六边形里两个相对空隙上面试试看，如果一定要优化出一个平面结构，建议把左侧Au原子位置下移，使之与左下侧连个Au原子形成三角形结构，这样对称性高一些.....

赞一下(1人)

回复此楼

2楼2010-12-19 10:26:47

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

yoghurt117

木虫 (正式写手)

应助: 3 (幼儿园)
金币: 8018.3
帖子: 576
在线: 134.9小时
虫号: 761174

★ ★
zzy870720z(金币+2):谢谢分享 2010-12-21 08:16:55
minmin_0082003(金币+1):谢谢，Au10 有很多的异构体，我算的是最稳定的一个 2010-12-21 09:31:22

曾经从网上找了个Au10的结构，跟你的好像不一样
Au -.79776090 -.69523418 2.61154030
Au    .21885458 1.44933212 1.53546771
Au 1.72079748 -.80162079 1.64308518
Au -2.27547260 .67789930 .81619191
Au -.38520612  -1.04247831 -.00370155
Au 1.88256086 .84182724 -.50326931
Au -.60370245 1.76230683  -1.01776187
Au 2.04405802  -1.73293814 -.87099203
Au -2.20281778 -.16059519  -1.74977567
Au    .39868891 -.29849888  -2.46078469
可以算算看看

赞一下(1人)

回复此楼

4楼2010-12-21 00:07:57

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