24小时热门版块排行榜    

查看: 763  |  回复: 9
当前主题已经存档。

chenyidan

[交流] 【求助】如何对坐标进行"冻结操作"谢谢!

请教大家,如何对一晶体结构,在高斯中进行"冻结操作"啊? 
如果有这方的资料总结,或是相关文献希望能多多指教,谢谢!

[ Last edited by zzgyb on 2008-9-10 at 11:48 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiekf1985

木虫 (正式写手)

★ ★
gromacs(金币+2,VIP+0):谢谢解答
GAussian说明中有,把参数分为Constant和Variable两部分,Constant部分参数被冻结
2楼2008-09-09 13:09:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sandycug

金虫 (小有名气)

★ ★
gromacs(金币+2,VIP+0):谢谢解答
直接把相应的院子序号在坐标后列表
比如将10号原子冻结

10 F
3楼2008-09-09 13:13:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kylin78

金虫 (小有名气)

3楼的说的不是gaussian 里的用法吧
4楼2008-09-09 17:16:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sandycug

金虫 (小有名气)

是gaussian的做法
5楼2008-09-09 20:55:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hbnuwangshuling

铜虫 (小有名气)

加F,文件为什么没冻结?

我的输入是# RHF/6-31G  OPT

frozen 1  2

0 1
C        -14.284202000      -0.482006000       0.000000000
C        -13.006781000       0.330542000       0.000000000
C        -11.687822000      -0.273042000       0.000000000
C        -10.528298000       0.395763000       0.000000000
C         -9.223995000      -0.240691000       0.000000000
C         -8.053933000       0.409834000       0.000000000
C         -6.757951000      -0.243110000       0.000000000
C         -5.582374000       0.397316000       0.000000000
C         -4.291199000      -0.265011000       0.000000000
C         -3.112568000       0.369648000       0.000000000
C         -1.823934000      -0.297610000       0.000000000
C         -0.643918000       0.334402000       0.000000000
C          0.643918000      -0.334402000       0.000000000
C          1.823934000       0.297610000       0.000000000
C          3.112568000      -0.369648000       0.000000000
C          4.291199000       0.265011000       0.000000000
C          5.582374000      -0.397316000       0.000000000
C          6.757951000       0.243110000       0.000000000
C          8.053933000      -0.409834000       0.000000000
C          9.223995000       0.240691000       0.000000000
C         10.528298000      -0.395763000       0.000000000
C         11.687822000       0.273042000       0.000000000
C         13.006781000      -0.330542000       0.000000000
C         14.135655000       0.387518000       0.000000000
H        -15.043347000      -0.035115000       0.000000000
H         15.043347000       0.035115000       0.000000000
H        -14.084629000      -1.461196000       0.000000000
H        -13.060329000       1.405870000       0.000000000
H        -11.658248000      -1.351765000       0.000000000
H        -10.541126000       1.472819000       0.000000000
H         -9.218259000      -1.319030000       0.000000000
H         -8.051725000       1.487135000       0.000000000
H         -6.764352000      -1.321038000       0.000000000
H         -5.571629000       1.474636000       0.000000000
H         -4.304395000      -1.342665000       0.000000000
H         -3.097131000       1.446976000       0.000000000
H         -1.840508000      -1.375086000       0.000000000
H         -0.626647000       1.411773000       0.000000000
H          0.626647000      -1.411773000       0.000000000
H          1.840508000       1.375086000       0.000000000
H          3.097131000      -1.446976000       0.000000000
H          4.304395000       1.342665000       0.000000000
H          5.571629000      -1.474636000       0.000000000
H          6.764352000       1.321038000       0.000000000
H          8.051725000      -1.487135000       0.000000000
H          9.218259000       1.319030000       0.000000000
H         10.541126000      -1.472819000       0.000000000
H         11.658248000       1.351765000       0.000000000
H         13.060329000      -1.405870000       0.000000000
H         14.084629000       1.461196000       0.000000000

1 2 F
运行结果没有冻结1 ,2原子它们之间的键长变化了,这是为什么?哪里出错了呢?
6楼2008-09-10 10:32:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mscic

金虫 (正式写手)

★ ★
lei0736(金币+2,VIP+0):谢谢
1,2 后面要写键参数吧?然后加个F表示冻结,还有S可是扫描等等。。。GAUSSIAN自带例子中有这个的
7楼2008-09-10 10:43:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

布丁3835

金虫 (小有名气)


jghe(金币+1,VIP+0):谢谢你参与讨论!
在分子描述部分使用直角坐标,在元素符号与坐标之间有0表示自由优化,-1表示冻结。
比如
C   0   x   y   z
C    -1 x  y   z
N   -1  x  y   z
8楼2008-10-05 15:17:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yytsnake

至尊木虫 (职业作家)


zzgyb(金币+1,VIP+0):谢谢你的参与,欢迎再次光临计算模拟版!
看看下面的文章就知道了:
http://new.quantumchemistry.net/ ... ew/200602/1021.html
9楼2008-10-05 19:55:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mscic

金虫 (正式写手)

★ ★
zzgyb(金币+2,VIP+0):谢谢你的参与,欢迎再次光临计算模拟版!
在ROUTE SECTION里面的关键词加一个MODREDUNDANT
然后在分子结构后面单独列出需要冻结的原子编号,加个F就可以.在GAUSSIAN的OPT关键字中有介绍.
10楼2008-10-07 09:17:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 chenyidan 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见