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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 【求助】优化的结构非常奇怪，Au10 簇 被优化成了4组小簇，为什么，输入及结构已贴出
131/1返回列表
查看: 2988  |  回复: 12
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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得到的
回复此楼

» 猜你喜欢

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

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

查看全部散金贴
1楼 2010-12-17 11:39:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hakuna

木虫 (知名作家)
★ ★
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的回帖

minmin_0082003

金虫 (正式写手)
★
ellsaking(金币+1):写一写出了什么问题也让大家借鉴一下吧，哈哈~ 2010-12-22 09:57:08
原因已经找到， 问题初步解决
赞一下(1人) 回复此楼
3楼2010-12-19 14:04:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yoghurt117

木虫 (正式写手)
★ ★
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的回帖

yoghurt117

木虫 (正式写手)
★
zzy870720z(金币+1):鼓励交流 2010-12-22 07:57:44
10至56个金Au的团簇结构请到下面的网址下载
http://u.115.com/file/f4e82e1827

这些结构也是我当初从网上找的，rar里面有网页介绍
赞一下(1人) 回复此楼
5楼2010-12-22 00:36:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sg18408926

至尊木虫 (著名写手)
好像没有了啊
回复此楼
6楼2010-12-22 08:54:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yoghurt117

木虫 (正式写手)
★
zzy870720z(金币+1):谢谢上传 2010-12-22 14:33:09
http://ishare.iask.sina.com.cn/f/12449694.html
还是ishare好用些，这里吧
赞一下(1人) 回复此楼
7楼2010-12-22 14:07:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hchfox

木虫 (正式写手)
★
zzy870720z(金币+1):谢谢讨论 2010-12-22 14:33:21
觉得参考价值不大，20左右金原子数的Aun团簇应该是中空结构才对，较新的DFT计算和一些实验也证实了这点的。。。
引用回帖:
Originally posted by yoghurt117 at 2010-12-22 14:07:28:
http://ishare.iask.sina.com.cn/f/12449694.html
还是ishare好用些，这里吧

赞一下(1人) 回复此楼
8楼2010-12-22 14:22:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

minmin_0082003

金虫 (正式写手)
个人觉得应该跟金属的结构没有关系
赞一下 回复此楼
9楼2010-12-22 17:23:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

minmin_0082003

金虫 (正式写手)
★
ellsaking(金币+1):谢谢回帖交流~ 2010-12-22 18:39:40
我开始用的
5d45d66s1的超软赝式
后面重新生成的模守恒赝式5d106s1，虽然好多了，但是还是达不到我的要求，原因还不知道，弄清楚再告诉大家，不要在纠结结构了，我觉得结构不是主要因素。Au10有好多个稳定的异构体。
赞一下(1人) 回复此楼
10楼2010-12-22 17:26:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hakuna

木虫 (知名作家)
★
zzy870720z(金币+1):谢谢讨论交流 2010-12-23 16:22:37
引用回帖:
Originally posted by hchfox at 2010-12-22 14:22:30:
觉得参考价值不大，20左右金原子数的Aun团簇应该是中空结构才对，较新的DFT计算和一些实验也证实了这点的。。。


金球是可能存在的，但不应该是最稳定结构，最稳定结构应该是成键最多的结构，显然没心没肺的金球不是成键最多的....
C球能稳定赖于共价键，但它稳定性上恐怕还是赶不上金刚钻.....
声明一下：没有计算过，所以只是一个判断
赞一下(1人) 回复此楼
11楼2010-12-23 09:41:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hchfox

木虫 (正式写手)
★
zzy870720z(金币+1):鼓励讨论交流 2010-12-23 16:22:56
呵呵，对与Au团簇，如果没算过，是不能轻易下判断的，我看到一片文献好像是Au20，是金字塔形状，但中间是中空的。根据我自己的经验，以前对Au13作了检验，用模拟退火方法找出来的结构已经有呈中空状态的趋势了，原先的高对称性结构根本是不稳定的，能量差异巨大。。。
引用回帖:
Originally posted by hakuna at 2010-12-23 09:41:09:



金球是可能存在的，但不应该是最稳定结构，最稳定结构应该是成键最多的结构，显然没心没肺的金球不是成键最多的....
C球能稳定赖于共价键，但它稳定性上恐怕还是赶不上金刚钻.....
声明一下：没有计算过， ...

赞一下(1人) 回复此楼
12楼2010-12-23 09:50:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hakuna

木虫 (知名作家)
引用回帖:
Originally posted by hchfox at 2010-12-23 09:50:38:
呵呵，对与Au团簇，如果没算过，是不能轻易下判断的，我看到一片文献好像是Au20，是金字塔形状，但中间是中空的。根据我自己的经验，以前对Au13作了检验，用模拟退火方法找出来的结构已经有呈中空状态的趋势了，原 ...

总感觉那样与常识相悖，能给出个令人信服的理由吗？
虽然空的金字塔里流荡着法老的幽灵，但那真的是最稳定结构吗？亦或是作者弄出来的玄虚？
不知道你用什么势函数退的火？SA理论上可以找到能量最低点，只是理论而已，相比之下，GA可能更好一些

[ Last edited by hakuna on 2010-12-23 at 19:41 ]
赞一下 回复此楼
13楼2010-12-23 19:32:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 minmin_0082003 的主题更新
131/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 zSbvVX9JFl9J 2026-11-09 4/200 2026-11-10 19:17 by czAiolzR0IJC
[教师之家] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 cqv1zKnzd7hr 2026-11-09 4/200 2026-11-10 12:52 by C6CRn66RYYXS
[考博] 申请环境，微生物方向博士 +3 浅笑09 2026-11-08 3/150 2026-11-10 12:44 by 可淡不可忘
[考研] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 cqv1zKnzd7hr 2026-11-09 4/200 2026-11-10 11:46 by YHiw1gou1Eda
[论文投稿] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 cqv1zKnzd7hr 2026-11-09 4/200 2026-11-10 11:35 by YHiw1gou1Eda
[教师之家] 延安大学回应“副教授留绝笔信”：确已离世，正核实举报内容 +4 yexuqing 2026-11-04 6/300 2026-11-10 11:13 by susizheng
[考研] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 0TEMgxdT1eIw 2026-11-09 4/200 2026-11-10 10:51 by YHiw1gou1Eda
[硕博家园] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 0TEMgxdT1eIw 2026-11-09 4/200 2026-11-10 10:18 by YHiw1gou1Eda
[教师之家] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 zSbvVX9JFl9J 2026-11-09 5/250 2026-11-10 09:23 by 1yv6QACPoY1E
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 zSbvVX9JFl9J 2026-11-09 3/150 2026-11-10 01:40 by FChKX8meu1ds
[考博] 27博士申请 +3 白白白ly 2026-11-09 3/150 2026-11-09 21:39 by szq0819
[文学芳草园] 厉娜 +3 张丹枫1991 2026-11-09 5/250 2026-11-09 16:53 by 张丹枫1991
[考博] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +5 kvewwGtKktqB 2026-11-04 12/600 2026-11-09 13:26 by Equinoxhua
[考研] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +5 ttH3f1oA6ace 2026-11-04 15/750 2026-11-09 13:19 by Equinoxhua
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +5 5sJFDKp2kwuI 2026-11-06 9/450 2026-11-08 22:46 by l48hFd8Nwae2
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +5 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 16:12 by 5sJFDKp2kwuI
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 ttH3f1oA6ace 2026-11-04 11/550 2026-11-06 16:10 by 5sJFDKp2kwuI
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +6 kvewwGtKktqB 2026-11-04 16/800 2026-11-06 05:00 by qL1cU46YIdK0
[考研] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 02:04 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 hNpp4Y2CGvSC 2026-11-04 12/600 2026-11-06 01:53 by qL1cU46YIdK0
信息提示
关闭
请填处理意见
关闭