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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Siesta&Smeagol&Atk » 【求助】吸附的分子又分解了，而且还不收敛，郁闷了，到底哪里错了
北京石油化工学院2026年研究生招生接收调剂公告
212/3返回列表上一页123下一页
查看: 1003  |  回复: 20
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

fanchen021

铁杆木虫 (正式写手)
★
ice_rain(金币+1):谢谢分享 1-27 18:28
给你简单生成了个VASP的POSCAR格式的坐标，你可以稍微改下。另外你这个分子比较大，3×3的表面可能有点小，lateral interaction 影响会比较大。
Au C N C H C H C H C H C N O
   1.000000000000000
   7.492315513494470  -4.325690378569650   0.000000000000001
   0.000000000000000   8.651380757139290   0.000000000000000
  -0.000000000000000   0.000000000000000  14.354607602836101
   18   1   1   1   1   1   1   2   1   1   1   1   1   2
Selective dynamics
Direct
   0.222074071569498   0.111037035784749   0.000000000000000   T   T   T
   0.333185182680609   0.333259258006971   0.164031485080153   T   T   T
   0.555407404902831   0.111037035784749   0.000000000000000   T   T   T
   0.666518516013942   0.333259258006971   0.164031485080153   T   T   T
   0.888740738236165   0.111037035784749   0.000000000000000   T   T   T
   0.999851849347275   0.333259258006971   0.164031485080153   T   T   T
   0.222074071569498   0.444370369118082   0.000000000000000   T   T   T
   0.333185182680609   0.666592591340304   0.164031485080153   T   T   T
   0.555407404902831   0.444370369118082   0.000000000000000   T   T   T
   0.666518516013942   0.666592591340304   0.164031485080153   T   T   T
   0.888740738236165   0.444370369118082   0.000000000000000   T   T   T
   0.999851849347275   0.666592591340304   0.164031485080153   T   T   T
   0.222074071569498   0.777703702451416   0.000000000000000   T   T   T
   0.333185182680609   0.999925924673638   0.164031485080153   T   T   T
   0.555407404902831   0.777703702451416  -0.000000000000000   T   T   T
   0.666518516013942   0.999925924673638   0.164031485080153   T   T   T
   0.888740738236165   0.777703702451416   0.000000000000000   T   T   T
   0.999851849347275   0.999925924673638   0.164031485080153   T   T   T
   0.230388362214848   0.846886688666260   0.307663266664798   T   T   T
   0.307045875199473   0.768725150766002   0.307928766992413   T   T   T
   0.588741989191645   0.750274734436410   0.308976552854301   T   T   T
   0.664903131083540   0.894343030402537   0.309300780959424   T   T   T
   0.676771031510386   0.651546923623279   0.309284982668288   T   T   T
   0.820831391286255   0.719450129928511   0.309836898562867   T   T   T
   0.578054814387281   0.464794790500489   0.308865160650123   T   T   T
   0.391295558854277   0.376766522746398   0.308146388321157   T   T   T
   0.315134421487918   0.232698231750831   0.307822160129361   T   T   T
   0.303269720279069   0.475494097167495   0.307844639715022   T   T   T
   0.159206174090524   0.407591145132435   0.307286040924886   T   T   T
   0.401999974129583   0.662250183259506   0.308255045526763   T   T   T
   0.672992630817064   0.358304159781186   0.309192890934224   T   T   T
   0.587651430980534   0.196889254717610   0.308829959865763   T   T   T
   0.834397114682275   0.434368304899188   0.309820443162921   T   T   T
赞一下(10人) 回复此楼
11楼2010-01-27 16:08:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

minmin_0082003

金虫 (正式写手)

莫愁莫愁
引用回帖:
Originally posted by fanchen021 at 2010-01-27 16:02:02:
分子构型没错的话，根据化学直觉，你这个分子吸附的时候应该平躺在表面，就是六元环和表面平行，而且O原子和表面形成化学键的可能性比较大

这个我没试过，这个已经有人做过了，的确是竖着的。我只是修改了一下对位的取代基...
赞一下 回复此楼
12楼2010-01-28 09:30:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

minmin_0082003

金虫 (正式写手)

莫愁莫愁
引用回帖:
Originally posted by fanchen021 at 2010-01-27 16:08:03:
给你简单生成了个VASP的POSCAR格式的坐标，你可以稍微改下。另外你这个分子比较大，3×3的表面可能有点小，lateral interaction 影响会比较大。
Au C N C H C H C H C H C N O
   1.000000000000000
   7.492 ...

非常感谢你，但是有点看不懂，我试试用4X4的做做看看。
赞一下 回复此楼
13楼2010-01-28 09:31:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

minmin_0082003

金虫 (正式写手)

莫愁莫愁
引用回帖:
Originally posted by fanchen021 at 2010-01-27 16:08:03:
给你简单生成了个VASP的POSCAR格式的坐标，你可以稍微改下。另外你这个分子比较大，3×3的表面可能有点小，lateral interaction 影响会比较大。
Au C N C H C H C H C H C N O
   1.000000000000000
   7.492 ...

Hi，请问，建立MS的时候会出现类似的一行
PBC 5.7676  2.8838   17.3576   90    90   120（P1）
该如何转换成lattice vector的形式呢？
我转换成
Lattice constant 2.8838 ang
Latticeparameters 2 1 6 90 90 120
前面有人回帖说可能是这个的错误。转换成vector会比较好？
请问该如何转换呢？不胜感激
赞一下 回复此楼
14楼2010-01-28 09:43:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

minmin_0082003

金虫 (正式写手)

莫愁莫愁
引用回帖:
Originally posted by hyb99 at 2010-01-27 10:45:43:
这个问题有可能是晶格矢量的设置问题

我看你在输入文件中是用Latticeparameters，siesta会根据6个参数自动生成晶格，但是这个晶格可能与你想要的晶格不一样（有旋转）。这种情况下最好要用内坐标，这样不会搞错 ...

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

请问该如何转换成vector呢？不胜感激
赞一下 回复此楼
15楼2010-01-28 10:48:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hyb99

金虫 (初入文坛)
★ ★ ★ ★ ★
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
赞一下(10人) 回复此楼
16楼2010-01-28 11:38:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

minmin_0082003

金虫 (正式写手)

莫愁莫愁
引用回帖:
Originally posted by hyb99 at 2010-01-28 11:38:29:


有两个方法

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

2）你可以在ms中选中晶格，然后在属性栏里面的最底下 ...

谢谢你了，这两个法子我都用不了，我没有MS。额。用的话会被罚款18，000，000日元......
赞一下 回复此楼
17楼2010-01-28 12:45:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fanchen021

铁杆木虫 (正式写手)
minmin_0082003(金币+9):x谢谢，我距离解决问题又近了一步 2010-01-28 18:13
引用回帖:
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  ...

只根据a b c 和夹角，并不能确定lattice vector.
还跟你的原子位置有关系。
赞一下 回复此楼
18楼2010-01-28 14:52:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fanchen021

铁杆木虫 (正式写手)
引用回帖:
Originally posted by minmin_0082003 at 2010-01-28 09:31:22:

非常感谢你，但是有点看不懂，我试试用4X4的做做看看。

1.000000000000000
   7.492315513494470  -4.325690378569650   0.000000000000001
   0.000000000000000   8.651380757139290   0.000000000000000
  -0.000000000000000   0.000000000000000  14.354607602836101
这部分就是lattice vector, 下面的是相对坐标 （fractional coordinates）

如果你确定是竖着的话，3×3应该也够了。
我做过苯的吸附，是平躺的。
赞一下 回复此楼
19楼2010-01-28 14:55:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

minmin_0082003

金虫 (正式写手)

莫愁莫愁
引用回帖:
Originally posted by fanchen021 at 2010-01-28 14:55:21:

   1.000000000000000
   7.492315513494470  -4.325690378569650   0.000000000000001
   0.000000000000000   8.651380757139290   0.000000000000000
  -0.000000000000000   0.000000000000000  14.35 ...

额，我查看了我的输出文件，发现 晶矢竟然是
outcell: Unit cell vectors (Ang):
        8.651400    0.000000    0.000000
       -4.325700    7.492332    0.000000
        0.000000    0.000000   27.061579
赞一下 回复此楼
20楼2010-01-28 18:05:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 minmin_0082003 的主题更新
212/3返回列表上一页123下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿河北工业大学材料工程，初试344求专硕调剂 +6 15933906766 2026-04-05 6/300 2026-04-06 13:21 by 无际的草原
[考研] 336材料与化工085600求调剂 +9 水星记infp 2026-04-05 12/600 2026-04-06 11:46 by 水星记infp
[考研] 288求调剂 一志愿哈工大 材料与化工 +13 洛神哥哥 2026-04-03 13/650 2026-04-05 17:27 by zzx2138
[考研] 工科277分求调剂材料 +8 上了上了上哦 2026-04-05 9/450 2026-04-05 13:05 by wwytracy
[考研] 083200 333求调剂 +3 十二！！ 2026-04-04 3/150 2026-04-05 08:28 by barlinike
[考研] 材料调剂 +12 一样YWY 2026-04-04 12/600 2026-04-05 08:24 by 544594351
[考研] 085601，一志愿厦大334复试被刷求调剂 +13 曾仰之 2026-04-03 15/750 2026-04-04 20:13 by dongzh2009
[考研] 一志愿东北大学085901土木专硕345求调剂 +3 zxt11111 2026-04-04 3/150 2026-04-04 14:21 by 土木硕士招生
[考研] 求调剂 +4 压力??大 2026-04-03 4/200 2026-04-03 21:36 by 啵啵啵0119
[考研] 326分求调剂 +3 于是乎呢 2026-04-01 5/250 2026-04-03 14:23 by 于是乎呢
[考研] 085600专硕材料与化工348分求调剂 +10 上学啦！ 2026-04-01 11/550 2026-04-03 14:13 by 百灵童888
[硕博家园] 求老师收留 +9 lllq123 2026-04-03 9/450 2026-04-03 13:48 by 呼吸都是减肥
[考研] 一志愿山东大学化学与化工学院材料与化工专硕，360分求调剂 +4 不愿透露姓名的� 2026-04-02 4/200 2026-04-03 09:29 by 遗忘消失的灆
[考研] 318求调剂，计算材料方向 +10 吸喵有害笙命 2026-04-01 11/550 2026-04-02 16:29 by oooqiao
[考研] 266求调剂 +4 学员97LZgn 2026-04-02 4/200 2026-04-02 13:03 by yulian1987
[考研] 一志愿北京科技，085601总分305求调剂 +9 半生瓜！ 2026-04-01 11/550 2026-04-02 08:28 by Wang200018
[考研] 生物与医药考研调剂 +5 铁憨憨123425 2026-03-31 5/250 2026-04-01 18:01 by syh9288
[考研] 求调剂，一志愿北林食品与营养095500，301分，已过六级，有科研经历 +4 快乐储蓄罐 2026-03-31 4/200 2026-04-01 09:26 by JourneyLucky
[考研] 英一数一总分334求调剂 +4 陈阳坤 2026-03-31 4/200 2026-03-31 14:22 by 记事本2026
[考研] 323分 食品与营养调剂 +3 嘿ooo 2026-03-31 3/150 2026-03-31 09:38 by longlotian
信息提示
关闭
请填处理意见
关闭