24小时热门版块排行榜    

查看: 2456  |  回复: 17

alexwpch

木虫 (正式写手)

引用回帖:
9楼: Originally posted by beefly at 2012-06-13 00:34:56
多半是自旋多重度不合适,计算了一些的高激发态

以Zr和Ti为例,去掉4个电子成四价了以后,最外层是8个电子,sp3刚好填满,这种情况也是高自旋吗?
11楼2012-06-13 08:45:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

meteoric30

木虫 (正式写手)

【答案】应助回帖

引用回帖:
10楼: Originally posted by alexwpch at 2012-06-13 08:43:53
加了maxcyc=300,xqc,还是L9999错误,并且金属离子离得很远很远到了18A。...

帖输入文件 我看下~~报错的时候几个收敛了?
12楼2012-06-13 09:05:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

meteoric30

木虫 (正式写手)

【答案】应助回帖

引用回帖:
10楼: Originally posted by alexwpch at 2012-06-13 08:43:53
加了maxcyc=300,xqc,还是L9999错误,并且金属离子离得很远很远到了18A。...

我建议你,先放的位置近一些,优化,看long的收敛,和结构,如果发现一直变远,就kill掉,调整初始构型,如果发现结构比较合适报999的错,这是scf不收敛,拉出最后一个结构,继续优化,我感觉你的能量是不是跳跃的,跳的话在木虫搜一下,出这个问题的很多回复~~
13楼2012-06-13 09:17:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

alexwpch

木虫 (正式写手)

引用回帖:
12楼: Originally posted by meteoric30 at 2012-06-13 09:05:33
帖输入文件 我看下~~报错的时候几个收敛了?...

Item               Value     Threshold  Converged?
Maximum Force            0.040574     0.000450     NO
RMS     Force            0.013112     0.000300     NO
Maximum Displacement     0.064171     0.001800     NO
RMS     Displacement     0.014782     0.001200     NO
Predicted change in Energy=-1.206551D-03
四个都不收敛。没法上传附件。
最后一段是
Isotropic Fermi Contact Couplings
        Atom                 a.u.       MegaHertz       Gauss      10(-4) cm-1
     1  C(13)              0.00000       0.00000       0.00000       0.00000
     2  C(13)              0.00000       0.00000       0.00000       0.00000
     3  C(13)              0.00000       0.00000       0.00000       0.00000
     4  C(13)              0.00000       0.00000       0.00000       0.00000
     5  C(13)              0.00000       0.00000       0.00000       0.00000
     6  C(13)              0.00000       0.00000       0.00000       0.00000
     7  H(1)               0.00000       0.00000       0.00000       0.00000
     8  H(1)               0.00000       0.00000       0.00000       0.00000
     9  H(1)               0.00000       0.00000       0.00000       0.00000
    10  H(1)               0.00000       0.00000       0.00000       0.00000
    11  H(1)               0.00000       0.00000       0.00000       0.00000
    12  Cl(35)             0.00000       0.00000       0.00000       0.00000
    13  Zr(91)             0.00000       0.00000       0.00000       0.00000
--------------------------------------------------------
       Center         ----  Spin Dipole Couplings  ----
                      3XX-RR        3YY-RR        3ZZ-RR
--------------------------------------------------------
     1   Atom        0.000000      0.000000      0.000000
     2   Atom        0.000000      0.000000      0.000000
     3   Atom        0.000000      0.000000      0.000000
     4   Atom        0.000000      0.000000      0.000000
     5   Atom        0.000000      0.000000      0.000000
     6   Atom        0.000000      0.000000      0.000000
     7   Atom        0.000000      0.000000      0.000000
     8   Atom        0.000000      0.000000      0.000000
     9   Atom        0.000000      0.000000      0.000000
    10   Atom        0.000000      0.000000      0.000000
    11   Atom        0.000000      0.000000      0.000000
    12   Atom        0.000000      0.000000      0.000000
    13   Atom        0.000000      0.000000      0.000000
--------------------------------------------------------
                        XY            XZ            YZ
--------------------------------------------------------
     1   Atom        0.000000      0.000000      0.000000
     2   Atom        0.000000      0.000000      0.000000
     3   Atom        0.000000      0.000000      0.000000
     4   Atom        0.000000      0.000000      0.000000
     5   Atom        0.000000      0.000000      0.000000
     6   Atom        0.000000      0.000000      0.000000
     7   Atom        0.000000      0.000000      0.000000
     8   Atom        0.000000      0.000000      0.000000
     9   Atom        0.000000      0.000000      0.000000
    10   Atom        0.000000      0.000000      0.000000
    11   Atom        0.000000      0.000000      0.000000
    12   Atom        0.000000      0.000000      0.000000
    13   Atom        0.000000      0.000000      0.000000
--------------------------------------------------------


---------------------------------------------------------------------------------
              Anisotropic Spin Dipole Couplings in Principal Axis System
---------------------------------------------------------------------------------

       Atom             a.u.   MegaHertz   Gauss  10(-4) cm-1        Axes

              Baa     0.0000     0.000     0.000     0.000  1.0000  0.0000  0.0000
     1 C(13)  Bbb     0.0000     0.000     0.000     0.000  0.0000  1.0000  0.0000
              Bcc     0.0000     0.000     0.000     0.000  0.0000  0.0000  1.0000

              Baa     0.0000     0.000     0.000     0.000  1.0000  0.0000  0.0000
     2 C(13)  Bbb     0.0000     0.000     0.000     0.000  0.0000  1.0000  0.0000
              Bcc     0.0000     0.000     0.000     0.000  0.0000  0.0000  1.0000

              Baa     0.0000     0.000     0.000     0.000  1.0000  0.0000  0.0000
     3 C(13)  Bbb     0.0000     0.000     0.000     0.000  0.0000  1.0000  0.0000
              Bcc     0.0000     0.000     0.000     0.000  0.0000  0.0000  1.0000

              Baa     0.0000     0.000     0.000     0.000  1.0000  0.0000  0.0000
     4 C(13)  Bbb     0.0000     0.000     0.000     0.000  0.0000  1.0000  0.0000
              Bcc     0.0000     0.000     0.000     0.000  0.0000  0.0000  1.0000

              Baa     0.0000     0.000     0.000     0.000  1.0000  0.0000  0.0000
     5 C(13)  Bbb     0.0000     0.000     0.000     0.000  0.0000  1.0000  0.0000
              Bcc     0.0000     0.000     0.000     0.000  0.0000  0.0000  1.0000

              Baa     0.0000     0.000     0.000     0.000  1.0000  0.0000  0.0000
     6 C(13)  Bbb     0.0000     0.000     0.000     0.000  0.0000  1.0000  0.0000
              Bcc     0.0000     0.000     0.000     0.000  0.0000  0.0000  1.0000

              Baa     0.0000     0.000     0.000     0.000  1.0000  0.0000  0.0000
     7 H(1)   Bbb     0.0000     0.000     0.000     0.000  0.0000  1.0000  0.0000
              Bcc     0.0000     0.000     0.000     0.000  0.0000  0.0000  1.0000

              Baa     0.0000     0.000     0.000     0.000  1.0000  0.0000  0.0000
     8 H(1)   Bbb     0.0000     0.000     0.000     0.000  0.0000  1.0000  0.0000
              Bcc     0.0000     0.000     0.000     0.000  0.0000  0.0000  1.0000

              Baa     0.0000     0.000     0.000     0.000  1.0000  0.0000  0.0000
     9 H(1)   Bbb     0.0000     0.000     0.000     0.000  0.0000  1.0000  0.0000
              Bcc     0.0000     0.000     0.000     0.000  0.0000  0.0000  1.0000

              Baa     0.0000     0.000     0.000     0.000  1.0000  0.0000  0.0000
    10 H(1)   Bbb     0.0000     0.000     0.000     0.000  0.0000  1.0000  0.0000
              Bcc     0.0000     0.000     0.000     0.000  0.0000  0.0000  1.0000

              Baa     0.0000     0.000     0.000     0.000  1.0000  0.0000  0.0000
    11 H(1)   Bbb     0.0000     0.000     0.000     0.000  0.0000  1.0000  0.0000
              Bcc     0.0000     0.000     0.000     0.000  0.0000  0.0000  1.0000

              Baa     0.0000     0.000     0.000     0.000  1.0000  0.0000  0.0000
    12 Cl(35) Bbb     0.0000     0.000     0.000     0.000  0.0000  1.0000  0.0000
              Bcc     0.0000     0.000     0.000     0.000  0.0000  0.0000  1.0000

              Baa     0.0000     0.000     0.000     0.000  1.0000  0.0000  0.0000
    13 Zr(91) Bbb     0.0000     0.000     0.000     0.000  0.0000  1.0000  0.0000
              Bcc     0.0000     0.000     0.000     0.000  0.0000  0.0000  1.0000


---------------------------------------------------------------------------------



YOU ARE IN A MAZE OF TWISTY LITTLE PASSAGES.
Error termination request processed by link 9999.
Error termination via Lnk1e in d:\gaussian09\l9999.exe at Tue Jun 12 23:56:21 2012.
Job cpu time:  0 days  1 hours 44 minutes 18.0 seconds.
File lengths (MBytes):  RWF=     25 Int=      0 D2E=      0 Chk=      2 Scr=      1
14楼2012-06-13 09:56:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家

引用回帖:
11楼: Originally posted by alexwpch at 2012-06-13 08:45:39
以Zr和Ti为例,去掉4个电子成四价了以后,最外层是8个电子,sp3刚好填满,这种情况也是高自旋吗?...

电子排布不能光靠主观猜测。看Ti的例子。从计算输出的电荷上来看,Ti只能失去两个电子。Ti二价离子基态多重度是闭壳层么?自己去查实验数据:

http://physics.nist.gov/cgi-bin/ASD/levels_pt.pl
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
15楼2012-06-13 12:56:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

alexwpch

木虫 (正式写手)

引用回帖:
15楼: Originally posted by beefly at 2012-06-13 12:56:46
电子排布不能光靠主观猜测。看Ti的例子。从计算输出的电荷上来看,Ti只能失去两个电子。Ti二价离子基态多重度是闭壳层么?自己去查实验数据:

http://physics.nist.gov/cgi-bin/ASD/levels_pt.pl...

这个。。。是我想的太简单了,非常感谢,我先学习下。
16楼2012-06-13 15:44:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zdgclf

银虫 (初入文坛)

【答案】应助回帖

四楼的说法是比较准确的,不过你可以多金属离子采用lanl2dz 那样计算的结果会更好!
%chk=F:\fe.chk
# opt rhf/genecp geom=connectivity pop=mk

Title Card Required

0 1
C                  0.84805700    0.97511100    0.33143200
C                  2.14219300    1.49348500    0.18852800
C                  3.20653100    0.64133500   -0.13542700
C                  2.97673300   -0.72918900   -0.31647800
C                  1.68259700   -1.24756300   -0.17357200
C                  0.61825900   -0.39541300    0.15038300
H                  0.03541100    1.62574500    0.57878200
H                  2.31764900    2.53990900    0.32676700
H                  4.19463300    1.03712500   -0.24453500
H                  3.78937900   -1.37982400   -0.56382300
H                  1.50714200   -2.29398900   -0.31180500
Cl                -1.00703000   -1.04643300    0.32985600
Fe                -2.44503600    0.45514600   -0.21777400

1 2 1.5 6 1.5 7 1.0
2 3 1.5 8 1.0
3 4 1.5 9 1.0
4 5 1.5 10 1.0
5 6 1.5 11 1.0
6 12 1.0
7
8
9
10
11
12 13 1.0
13

C H Cl 0
6-31g
****
Fe 0
lanl2dz
****
这是输入文件
17楼2012-06-14 15:16:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

大雪微蓝

至尊木虫 (正式写手)

【答案】应助回帖

你增大循环次数试试
18楼2013-06-26 00:18:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 alexwpch 的主题更新
信息提示
请填处理意见