当前位置:文档之家› 单晶结构解析步骤

单晶结构解析步骤

运行xp
自动调入test.res文件
fmol
proj
选择好的方向,可以观察到很多原子,包含您的分子模型
pick
删除不是你分子模型的其他Q原子(保留:空格键;定原子:直接输入原子编号并回车
去掉:回车;错误更改:回车后再back键)
在闪烁的Q原子上直接输入原子类型和编号,一般新产生的环编号较小
回车
HADD
2、修改ins文件,将X和X‘的原子的占有率先定为10.500000
3、相应的氢的占有率先定为10.500000
4、运行XL和精修到收敛,打开CIF,观察X和X‘的热参数值,如X大,则降低其占有率,但保持X和X‘的占有率为1,例如X为10.450000, X‘为10.550000.氢做相应的修改,
5、运行运行XL和精修到收敛,打开CIF,观察X和X‘的热参数值,重复以上直到X和X‘的热参数值相同或接近。
plane
xp
read file name
fmol
mpln atom1 atom 2.....
enter
angle
xp
read file name
fmol
mpla n(atom number) atom1 atom 2.....
mpla n(atom number) atom1 atom 2.....
edit ins
将cell的三个数据参数分别成以0.99(或0.98, 0.97....)
保存ins
再画图即可。
补充:
产生report
运行xcif
运行默认的命令
File......输入report
File......[ang]输入def
运行[F]时改运行[Q]
是否需要运行消光吸收校正
在运行XL中有一行提示
在ins中加入EXTI命令
mol图的大小
TELP 0 -30
-30表示30%
如-50表示50%
解析晶体的一般顺序
1、先定出结构模型XP(file)
2、运行XL
2、各项异性修正XP(file)(在ins中加入anis命令,运行XL一般r1值大大变小,XP中envi中Q<1)
3、加氢XP(file)(HADD $C HADD $N HADD $O)
SGEN 2676则会产生另外一个分子
CENT/X C3 to C8则会产生X1A
CENT/X C3A to C8A则会产生X1B
LINK X1A X1B
MATR 1 OR 2 OR 3
PACK
PROJ
使用UNDO命令删除多余的键
Telp cell
……
氢键的对称操作
例如C2-H2A…N1
则ENVI C2 3
shelxtl
open new
name
xp
fmol
kill $q
proj
select the good direction
exit
telp 0 -30
plotfile enter file name
draw file name
select file(ps file)
black and white
cell
寻找N1的对称操作,如为7456
SYMM
7555=XYZ
对映值相加。
然后
link X1A C? (C? 就是你需要的那个CH...pi 中的C原子)
然后
envi X1A 3
就能显示离X1A原子距离为4埃以内的所有原子的对称操作键长键角
计算π•••π stacking也是如此
如C3-C8苯环存在堆积
Envi c3 3
观察除1555 C3-C8以外的对称操作,例如C3 2676
解决无序问题最好方法:
1、根据Q峰(一般>1)和键长,命名为X’,除去所有的Q峰,保存文件。
2、修改ins文件,将X和X‘的原子定为part 1
Part 2
Part 0三个部分
例如C26和C26’无序
PART 1
C26 1 0.390625 0.672967 1.176015 21.000000 0.171080 0.118300 =
选中画笔
点出两个点
按ESC
点选择键
选中画笔
鼠标移动至出现小手
拖动到其他角度。
氢键数据的列出
eqiv $1 x, y, z(对称操作)
htab c12 o2_$1(形成氢键的两个重原子)
运行XL命令
在cif和1st文件中即有相关的数据
对称操作的寻找
XP中
cell中扫描NON-B
找到形成氢键的两个重原子
在钻石中找到相应的角
例如C12-H12...O2的角度中的对称操作就是
如有多个氢键
请输入$1, $2....
解析晶体结构的方法
打开shelxtl
选择new
选择路径
选择raw文件
输入你的文件名(随便)
运行xprep
多个回车键
输入分子式(不知道的输入可能的原子类型)
输入另外一个全新的文件名(例如test,将产生test.ins文件)
(一般加入所有H原子)
如有的氢加错误可以去掉所有氢(kill $h
也可以选择性例如去掉h18(kill h18)
单独加氢(hadd氢所连接的非氢原子类型)
类型(1:叔氢
2:亚甲基
3:甲基
4:芳氢
proj
观察分子包括氢是否对
file test(如您保存的ins文件名为test)
exit
上述步骤产生test.ins(如您保存的ins文件名为test)
0.061140 0.023050 -0.000040 -0.062810
AFIX 33
H26A 2 0.340249 0.671478 1.257577 21.000000 -1.500000
H26B 2 0.478542 0.699786 1.170814 21.000000 -1.500000
询问是否保存新的hkl文件
输入y
回车
上述
选择刚产生的test.hkl文件
运行xs
程序自动调入test.ins文件
如在黑色屏幕的最下行显示re<0.3(例如0.289)
则可以解出结构
上述步骤产生test.res文件(如您保存的res文件名为test)
********
运行xl
自动调入test.ins文件
产生test.res, test.cif, test.lst文件
运行xshell
自动新产生的test.res文件
选择edit中的去掉所有的Q
运行refine
选择aniso, bond, acta
cycle任意选择非0数值
运行
产生test.res, test.cif, test.lst文件
AFIX 0
PART 0
3、FVAR 0.699860(这个值原先就有) 加上 0.5200
运行XL和精修到收敛即可。
Xp下用CENT命令做出苯环中心哑原子,然后在与相应的C(或H)链接即可
如:
CENT/x C1 to C6 (苯环上连续6个C,这样写也行C1 C2 C3 C4 C5 C6)
会生成X1A哑原子
AFIX 33
H26D 2 0.062772 0.949197 0.982469 -21.000000 -1.500000
H26E 2 -0.003777 0.938874 1.119104 -21.000000 -1.500000
H26F 2 0.015775 0.824384 1.040314 -21.000000 -1.500000
H26C 2 0.430851 0.581169 1.150402 21.000000 -1.500000
PART 2
AFIX 0
C26' 1 0.060971 0.885965 1.053770 -21.000000 0.077010 0.161380 =
0.087380 -0.055200 0.017380 -0.058090
例如c1 c2 c3 c4 c5 n1中,c1 c2 c4 c5共平面
mpla 4 c1 c2 c4 c5 c3 n1
txt
运行xcif
选择t
两次回车
输入文件名.txt
选择def
回车直到选择q
理论加氢
在ins中输入
HFIX要加氢的原子
保存ins
运行XL
打开RES
拷贝相应的数据到ins中即可。
CHEMICAL DRAW
fmol
kill $q
matr 1=a 2=b 3=c
pbox 5 15
pack
select (space=keep, enter=del)
fmol
telp cell
enter file name
draw file name
select file type(a=psfile)
black and white(enter)
有效体积大小一般接近17最好
如解不出结构
1、晶系的选择M---A
2、空间群P-1---P1;C2/C---CC
定原子模型
1、KILL相差较大的Q峰(例如0.99----0.33,0.33的可以去掉)
2、name q? c?
name q?? c??
再name cxx---o?....
修正
XL运行后R1变小证明刚才的原子定的是正确的
打开test.cif,
观察r1(<0.07)
good of factor(接近1)
max(0.001,或0.000)
相关主题