应《网络安全法》要求，自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下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 生物学308分求调剂（一志愿华东师大） +9 相信必会光芒万� 2026-04-06 12/600 2026-04-06 21:49 by 52305043001
[考研] 调剂求助（生物与医药） +3 @6952 2026-04-06 3/150 2026-04-06 21:48 by Ecowxq666！
[考研] 270分求调剂 +5 maxjxbsk 2026-04-01 5/250 2026-04-06 15:32 by lin-da
[考研] 求调剂到材料 +5 程9915 2026-04-06 5/250 2026-04-06 15:21 by yulian1987
[考研] 298求调剂 +3 manman511 2026-04-05 3/150 2026-04-05 18:09 by kk112233
[考研] 085600调剂 +9 东照照照 2026-04-04 9/450 2026-04-05 13:44 by ujn_zhuj
[考研] 一志愿华中农业大学0710（A）初试329分 求调剂 +4 一名26考研生 2026-04-04 4/200 2026-04-05 10:01 by barlinike
[考研] 求生物学调剂 +14 15172915737 2026-04-01 14/700 2026-04-04 20:13 by babysonlkd
[考研] 求调剂 +3 农业工程与信息� 2026-04-04 3/150 2026-04-04 12:19 by 舍而后得
[考研] 求调剂，一志愿北京中医药大学 +3 小小达不溜 2026-04-02 3/150 2026-04-03 22:55 by 冲矢昴星团
[考研] 266求调剂 +3 08电气工程 2026-04-03 3/150 2026-04-03 14:05 by 1753564080
[考研] 材料考研调剂 +10 Gs大王 2026-04-02 10/500 2026-04-03 09:47 by 遗忘消失的灆
[考研] 312求调剂 +4 赊月色 2026-04-02 5/250 2026-04-03 08:21 by fangshan711
[考研] 362求调剂 +14 西南交材料专硕3 2026-03-31 14/700 2026-04-02 17:50 by yunlongyang
[考研] 一志愿北京科技大学材料学硕328分求调剂 +6 1段时间 2026-03-31 7/350 2026-04-02 13:57 by 3041
[考研] 354求调剂 +4 lxb598 2026-03-31 5/250 2026-04-02 09:55 by Jaylen.
[考研] 材料调剂 +12 一样YWY 2026-04-01 12/600 2026-04-02 09:15 by olim
[考研] 292求调剂 +17 木虫er12138 2026-04-01 17/850 2026-04-01 21:37 by 七度不信任
[考研] 08工科，295，接受跨专业调剂 +6 lmnlzy 2026-03-31 6/300 2026-04-01 11:02 by 逆水乘风
[考研] 339求调剂 +5 zjjkt 2026-03-31 5/250 2026-04-01 09:18 by JourneyLucky
信息提示
关闭
请填处理意见
关闭