24小时热门版块排行榜    

Znn3bq.jpeg
查看: 816  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料化工总分334求调剂 +15 Riot2025 2026-04-08 15/750 2026-04-09 19:15 by 诗与自由
[考研] 269求调剂 +7 跪求收留。 2026-04-04 7/350 2026-04-09 19:06 by 探123
[考研] 初试分332,一志愿报考西北工业大学, +8 故人?? 2026-04-09 8/400 2026-04-09 18:28 by l_paradox
[考研] 0703化学调剂325分 +13 15771691647 2026-04-04 15/750 2026-04-09 16:55 by 15771691647
[考研] 288求调剂 +15 没有答案_ 2026-04-05 15/750 2026-04-09 10:22 by 5268321
[考研] 302分求调剂 +4 凡语祈愿 2026-04-08 5/250 2026-04-08 22:03 by 土木硕士招生
[考研] 263分B区求调剂 +6 李nihao 2026-04-08 6/300 2026-04-08 09:38 by 南开小綦
[考研] 求考研材料调剂 +3 材化李可 2026-04-07 3/150 2026-04-08 00:21 by JourneyLucky
[考研] 316求调剂 +4 15318418673 2026-04-07 4/200 2026-04-07 22:12 by hemengdong
[考研] 22408 一志愿双一流人工智能300分 四六级,数据分析国奖 +4 zzfeng123 2026-04-06 6/300 2026-04-07 21:02 by zzfeng123
[考研] 315求调剂 +3 TUZEIQAQ 2026-04-02 3/150 2026-04-07 17:32 by chenp123
[考研] 一志愿华中农业大学0710(A)初试329分 求调剂 +5 一名26考研生 2026-04-04 5/250 2026-04-07 08:54 by 18828373951
[考研] 生物与医药求调剂 +7 heguanhua 2026-04-05 8/400 2026-04-06 18:41 by macy2011
[考研] 08专硕275调剂 +5 AaAa7420 2026-04-05 5/250 2026-04-05 18:01 by jkddd
[考研] 288求调剂,一志愿华南理工大学071005 +6 ioodiiij 2026-04-04 6/300 2026-04-05 10:09 by guoweigw
[考研] 342求调剂 +3 Liang7111 2026-04-04 5/250 2026-04-04 19:47 by dongzh2009
[考研] 368求调剂 +5 今华习 2026-04-03 7/350 2026-04-04 18:47 by imissbao
[考研] 调剂0855-288 +5 x熊二a 2026-04-03 5/250 2026-04-04 00:19 by 猪会飞
[考研] 本科985,专业0812分336求调剂 +4 莫莫很行 2026-04-03 4/200 2026-04-03 21:31 by zhq0425
[考研] 求调剂 +3 usbdndj 2026-04-03 3/150 2026-04-03 14:10 by dxiaoxin
信息提示
请填处理意见