24小时热门版块排行榜    

查看: 15950  |  回复: 82
【奖励】 本帖被评价43次,作者xi2004增加金币 37.6
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

xi2004

至尊木虫 (职业作家)


[资源] 【分享】教你如何在CIF文件中表达氢键

本来是回答4028...同学的问题的. 但那个主题不知道为何被删除了.
看其他人也老在问, 就在这里说一说.


问题一:
向ACta Cryst. C或E投稿, 需要产生氢键表, 该怎么办呢?

正确次序如下:

1. 先在.ins中加

HTAB

(加在UNIT后面和FVAR之前的任何地方)

2. 从.lst中 查看可能的氢键,

一般会看到如下

Hydrogen bonds with  H..A < r(A) + 2.000 Angstroms  and   110 deg.

D-H           d(D-H)   d(H..A)    O25-H13        0.902    2.343   149.03    3.151    O21
O19-H3         0.864    2.168   154.00    2.969    O4 [ -x, -y, -z ]
O21-H6         0.857    1.949   144.74    2.695    O27 [ -x, -y, -z+1 ]

3. 在 .ins中, 加以下行

ACTA
BOND $H
HTAB
HTAB O25 O21
EQIV $1 -x, -y, -z
EQIV $2 -x, -y, 1-z
HTAB O19 O4_$1
HTAB O21 O27_$2

(加在UNIT后面和FVAR之前的任何地方)

4. 精修n轮到收敛

5. 那么, 在.lst中

你会看到

Specified hydrogen bonds (with esds except fixed and riding H)

  D-H          H...A        D...A        <(DHA)

  0.90         2.34         3.151(8)     149.0        O25-H13...O21
  0.86(2)      2.17(4)      2.969(8)     154(8)       O19-H3...O4_$1
  0.86(2)      1.95(6)      2.695(7)     145(8)       O21-H6...O27_$2

5. 而在 .cif中,

你会看到以下信息. 注意一下 symmtry code的对应关系.
比如3556 表示 -x, -y, 1-z


_symmetry_cell_setting            monoclinic
_symmetry_space_group_name_H-M    P2(1)/n

loop_
_symmetry_equiv_pos_as_xyz
'x, y, z'
'-x+1/2, y+1/2, -z+1/2'
'-x, -y, -z'
'x-1/2, -y-1/2, z-1/2'

(这是P21/n空间群相对应的四个symmetry operations)

loop_
_geom_hbond_atom_site_label_D
_geom_hbond_atom_site_label_H
_geom_hbond_atom_site_label_A
_geom_hbond_distance_DH
_geom_hbond_distance_HA
_geom_hbond_distance_DA
_geom_hbond_angle_DHA
_geom_hbond_site_symmetry_A
O25 H13 O21  0.90 2.34 3.151(8) 149.0 .
O19 H3 O4  0.86(2) 2.17(4) 2.969(8) 154(8) 3
O21 H6 O27  0.86(2) 1.95(6) 2.695(7) 145(8) 3_556


问题二: 如何查看不常见的氢键?,
用别的程序算出来的氢键, 其acceptor的symmetry code对称码与shelx算的不一致怎么办?

象C-H...O这类氢键, Shelx好象不给计算, 至少是我不知道如何算. 那么怎么算呢?

方法一: 用Platon, 或者PWT中的Platon 其中有Calc H-bond或Calc all的命令
方法二: 用Material Studio可以检查, 但要把C设为donor

但是, 根据一些虫子的报道, Platon算出来的氢键, 其对称操作码会和XP和shelx 算出来的不一致.
这个, 据我查看, 好象是Platon计算时, 有时会把你的文件中的有些原子坐标做个平移.(是不是这个原因, 要请熟悉Platon的人说说).
而要注意的是, 有时候, XP中的symmetry code (用SYMM指令)和最终CIF中的 也会有所不同, 比如 P2(1)/n, 在XP中, 就把-x, -y, -z排在2, 而在cif中是在3.

我的看法, 如果你的结构是用Shelx或Shextl解的, 可以用Platon计算非常规氢键, 但不要相信platon所给的对称操作码, 因为platon做计算时可能把有些原子的坐标平移了.
这时候, 要回过头来用XP来查对正确的对称操作码.
比如我举的上述例子, 用Platon算的, 有个氢键是
O(21) --H(6)   ..O(27)  [  2545.07]       0.86       1.95     2.6946        145

那你用XP打开.ins,
fmol后,
arad 0.2 2 2 O21
envi o21

你就会发现与O21距离是2.6946对应的O27的symmetry code是2556,
symm 一下, 发现
2 对应的是-x, -y, -z
所以, O27的对称操作是-x, -y, 1-z
而要注意的是, 在CIF中, -x, -y, -z是排在第三个的, 所以O27的symmetry code 最后在CIF中应为3556.


所以, 关键是要明白symmetry code的含义.
然后要以CIF中的为准.



问题三: 对称操作的含义



关于对称操作码的含义, 举个例子:

比如 P21/n 这个空间群,
在CIF中有
_symmetry_cell_setting            monoclinic
_symmetry_space_group_name_H-M    P2(1)/n

loop_
_symmetry_equiv_pos_as_xyz
'x, y, z'
'-x+1/2, y+1/2, -z+1/2'
'-x, -y, -z'
'x-1/2, -y-1/2, z-1/2'

上面四个就是P21/n这个空间群对应的对称操作.


如果你有个氢键是:

O21 H6 O27  0.86(2) 1.95(6) 2.695(7) 145(8) 3_556

就表明: O27的坐标是在 -x, -y, 1-z的位置.


其中3表示空间群P21/n的第三个对称操作, 即'-x, -y, -z' ; 也就是在_symmetry_equiv_pos_as_xyz 中排第三的那个操作.
而556表示在a, b ,c方向上的平移. 5不平移. 6是往正方向平移一个单胞, 即+1, 4是-1.


注意一下, 在.ins中, 对称操作数目对有心的空间群来说是减半的, 而且x, y, z这个操作也不会出现. 比如
上述P21/n,  在.ins中, 就只出现 0.5-X, 0.5+Y, 0.5-Z

具体举例如下:

TITL NONE in P2(1)/n  
CELL 0.71073   1.0000   1.0000   1.0000  90.000  90.000  90.000
ZERR    1.00   0.0002   0.0002   0.0002   0.000   0.030   0.000
LATT  1
SYMM 0.5-X, 0.5+Y, 0.5-Z

[ Last edited by xi2004 on 2009-7-31 at 09:15 ]
回复此楼

» 收录本帖的淘帖专辑推荐

晶体生长合成专辑 初学晶体 优秀资源 单晶解析
晶体解析相关 解晶体 软件专辑以及解析晶体相关问题 晶体好贴
结构分析软件学习 与晶体相关的好帖子 晶体 crystal
xi2004 cif check cif 先辈精华 我的
测试 晶体 专业知识 单晶解析问题集锦
有机问题解答 专用帖子集 氢键-超分子 杂多酸
晶体解析 晶体方面 晶体相关

» 猜你喜欢

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

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

垂直极限

木虫 (正式写手)


谢谢楼主!学习一下……
7楼2008-09-15 20:27:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 83 个回答
谢谢 非常有帮助 收藏 为什么不能评分·~
2楼2008-09-15 12:29:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

4220824

木虫 (正式写手)



我那个帖子斑竹说求助完了,就要删的,哈哈,正愁没保存下来呢,谢谢LZ
3楼2008-09-15 13:09:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gujm516

金虫 (小有名气)



纳米镍粉(金币+1,VIP+0):
说的好。
不过第5点,3_556后面的对称操作笔误了。
4楼2008-09-15 15:11:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 279分求调剂 一志愿211 +9 chaojifeixia 2026-03-19 10/500 2026-03-20 12:29 by lature00
[考研] 085700资源与环境308求调剂 +10 墨墨漠 2026-03-18 11/550 2026-03-20 11:57 by 学员8dgXkO
[考研] 271材料工程求调剂 +7 .6lL 2026-03-18 7/350 2026-03-20 09:10 by xingguangj
[考研] 一志愿苏州大学材料求调剂,总分315(英一) +3 sbdksD 2026-03-19 3/150 2026-03-19 23:21 by fmesaito
[考研] 0703化学调剂 +4 18889395102 2026-03-18 4/200 2026-03-19 16:13 by 30660438
[考研] 本人考085602 化学工程 专硕 +17 不知道叫什么! 2026-03-15 19/950 2026-03-19 15:06 by 尽舜尧1
[考研] 286求调剂 +6 lemonzzn 2026-03-16 10/500 2026-03-19 14:31 by lemonzzn
[考研] 324分 085600材料化工求调剂 +3 llllkkkhh 2026-03-18 3/150 2026-03-19 14:22 by houyaoxu
[考研] 一志愿 西北大学 ,070300化学学硕,总分287,双非一本,求调剂。 +3 晨昏线与星海 2026-03-19 3/150 2026-03-19 13:36 by houyaoxu
[考研] 一志愿福大288有机化学,求调剂 +3 小木虫200408204 2026-03-18 3/150 2026-03-19 13:31 by houyaoxu
[考研] 0703化学 305求调剂 +4 FY_yy 2026-03-14 4/200 2026-03-19 05:54 by anny19840123
[考研] 354求调剂 +4 Tyoumou 2026-03-18 7/350 2026-03-18 21:45 by Tyoumou
[考研] 材料专业求调剂 +5 hanamiko 2026-03-18 5/250 2026-03-18 20:19 by 星空星月
[考研] 化学工程321分求调剂 +15 大米饭! 2026-03-15 18/900 2026-03-18 14:52 by haxia
[考研] 08工科 320总分 求调剂 +5 梨花珞晚风 2026-03-17 5/250 2026-03-18 14:49 by haxia
[考研] 331求调剂(0703有机化学 +7 ZY-05 2026-03-13 8/400 2026-03-18 14:13 by 007_lilei
[考研] 304求调剂 +12 小熊joy 2026-03-14 13/650 2026-03-18 12:34 by Linda Hu
[考研] 有没有道铁/土木的想调剂南林,给自己招师弟中~ +3 TqlXswl 2026-03-16 7/350 2026-03-17 15:23 by TqlXswl
[考研] 药学383 求调剂 +3 药学chy 2026-03-15 4/200 2026-03-16 20:51 by 元子^0^
[考研] 304求调剂 +5 素年祭语 2026-03-15 5/250 2026-03-16 17:00 by 我的船我的海
信息提示
请填处理意见