24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1280  |  回复: 4
当前主题已经存档。

zhugezilong

木虫 (著名写手)

[交流] [求助]Gaussian计算中怎么来求超极化率呀

本人是个菜鸟级的。对一个有机染料在Chem3D中进行超极化率的计算,
3D显示结果为
Gaussian Job: # UHF/6-31G Polar Test
Hyper Polarizability =  87.8939 797.5138 562.1404 230.6036 52.0865 82.6628 39.709 -60.3895 46.4228 -10.2369

而gaussian的输出结果为

Dipole moment (field-independent basis, Debye):
    X=   -11.1536    Y=     4.1110    Z=     1.0181  Tot=    11.9306
Quadrupole moment (field-independent basis, Debye-Ang):
   XX=  -155.1729   YY=  -141.7660   ZZ=  -154.0228
   XY=    -7.5755   XZ=    -9.4532   YZ=     3.2320
Traceless Quadrupole moment (field-independent basis, Debye-Ang):
   XX=    -4.8524   YY=     8.5546   ZZ=    -3.7022
   XY=    -7.5755   XZ=    -9.4532   YZ=     3.2320
Octapole moment (field-independent basis, Debye-Ang**2):
  XXX=  -143.3295  YYY=     3.6014  ZZZ=    38.5871  XYY=   -52.9743
  XXY=    58.8708  XXZ=   104.6381  XZZ=   -94.9501  YZZ=    14.0425
  YYZ=    -1.4441  XYZ=   -25.8536
Hexadecapole moment (field-independent basis, Debye-Ang**3):
XXXX=-23081.9677 YYYY= -1743.3300 ZZZZ=  -595.0245 XXXY=  -599.7630
XXXZ=  -786.5680 YYYX=    25.3704 YYYZ=     4.0096 ZZZX=  -303.9212
ZZZY=    48.6026 XXYY= -3860.9585 XXZZ= -3800.9510 YYZZ=  -439.8183
XXYZ=   146.9746 YYXZ=   -11.4379 ZZXY=   -56.3928
N-N= 2.075309850846D+03 E-N=-6.857772005440D+03  KE= 1.158057797928D+03
  Exact polarizability: 551.479   9.063 212.191  37.252   5.144 104.035
Approx polarizability: 339.761  11.998 205.778  26.557   4.962  98.392


我想问的是,那个才是我要求的一阶超极化率的数值呀。具体要怎么去求呢?
Hyper Polarizability =  87.8939 797.5138 562.1404 230.6036 52.0865 82.6628 39.709 -60.3895 46.4228 -10.2369
中这些数值的含义是什么呀?单位呢?

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

yuhuobuku

木虫 (正式写手)

★ ★ ★ ★ ★ ★ ★ ★ ★
zzgyb(金币+4,VIP+0):谢谢你的参与,欢迎再次光临计算模拟版!
zhugezilong(金币+5,VIP+0):谢谢解答
保留chk文件了吗?如果保留了,转化成fchk文件,打开后看最后那部分,有如下:
Dipole Moment                              R   N=           3
-1.78977816E-01  3.21669792E-01 -4.55191440E-14
Polarizability                             R   N=           6
  1.01401502E+03  1.86464510E+02  7.59784978E+02 -2.64393593E-05  3.99348314E-05
  6.83520941E+02

polarizability的那些值就是各个分量,就是一阶非线性系数,就是极化率,也就是alpha。总的和平均值都需要带入公式,公式就看文献吧。超极化率是指beta和gama。用polar只能算极化率,如果用polar=enonly能计算极化率和第一超极化率,就是alpha和beta。gama得不能用polar算出。得用其他方法,比如态求和(SOS)或者有现场(FF)方法。上面给你的那个例子就是用polar算出来的结果,直接在fchk文件里就能很清晰的找到。

[ Last edited by yuhuobuku on 2008-10-12 at 13:44 ]
薄酒可以忘忧,丑妻可以白头,徐行不必驷马,称身不必狐裘平淡才是生活的主旋律
2楼2008-10-12 13:43:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhugezilong

木虫 (著名写手)

真是非常感谢。我还想问问
对于计算得到的超极化率 共有十个数值,他们分别代表那几个张量呀?
是xxx,xxy,xxz之类的吗?
Dipole Moment                              R   N=           3
-4.38816493E+00  1.61739620E+00  4.00534124E-01
Polarizability                             R   N=           6
  5.51479003E+02  9.06290697E+00  2.12191078E+02  3.72516771E+01  5.14358937E+00
  1.04035192E+02
HyperPolarizability                        R   N=          10
  8.78939450E+01  7.97513786E+02  5.62140446E+02  2.30603633E+02  5.20865354E+01
  8.26627679E+01  3.97090078E+01 -6.03894687E+01  4.64228105E+01 -1.02368810E+01
3楼2008-10-12 20:42:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yuhuobuku

木虫 (正式写手)

★ ★
zzgyb(金币+2,VIP+0):谢谢你的参与,欢迎再次光临计算模拟版!
以上分量按顺序分别为:
mx
my
mz
axx
axy
ayy
axz
ayz
azz
βxxx
βxxy
βxyy
βyyy
βxxz
βxyz
βyyz
βxzz
βyzz
βzzz
薄酒可以忘忧,丑妻可以白头,徐行不必驷马,称身不必狐裘平淡才是生活的主旋律
4楼2008-10-13 13:08:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

greatsea1212

铜虫 (小有名气)

请问,用FF方法计算超极化率应该用什么关键字啊?急切
5楼2008-11-14 10:55:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zhugezilong 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿A区211,22408 321求调剂 +7 随心所欲☆ 2026-04-15 8/400 2026-04-21 08:22 by Equinoxhua
[考研] 295分求调剂 +6 ?要上岸? 2026-04-17 6/300 2026-04-21 08:18 by Equinoxhua
[考研] 一志愿中科大材料与化工,353分还有调剂学校吗 +11 否极泰来2026 2026-04-15 13/650 2026-04-20 22:31 by Equinoxhua
[教师之家] 又一批高校组建人工智能学院 师资行吗 不是骗人吗 +4 yexuqing 2026-04-19 4/200 2026-04-20 14:47 by brantleo
[考研] 297,工科调剂?河南农业大学本科 +15 河南农业大学-能 2026-04-14 15/750 2026-04-20 12:39 by Equinoxhua
[考研] 337求调剂 +3 jyz04 2026-04-18 3/150 2026-04-20 12:24 by 研可安
[考博] 湖南大学刘巧玲课题组2026年第二批次博士研究生招生信息 +3 南风观火 2026-04-18 5/250 2026-04-20 10:13 by 南风观火
[考研] 求计算机方向调剂 +3 Toffee2 2026-04-16 6/300 2026-04-19 22:37 by ll叶
[考研] 通信工程求调剂!!! +7 zlb770521 2026-04-14 7/350 2026-04-19 20:56 by Equinoxhua
[考研] 291求调剂 +11 关忆北. 2026-04-14 11/550 2026-04-19 17:16 by 中豫男
[考研] 304求调剂 +8 castLight 2026-04-16 8/400 2026-04-19 17:14 by 中豫男
[考研] 求调剂推荐 +9 小聂爱学习 2026-04-14 9/450 2026-04-19 17:03 by 中豫男
[考研] 求调剂 +10 小聂爱学习 2026-04-16 12/600 2026-04-19 16:51 by 中豫男
[考研] 085404 22408 309分求调剂 +10 lzmk 2026-04-14 11/550 2026-04-19 16:42 by 中豫男
[考研] 求调剂 +6 苦命人。。。 2026-04-18 7/350 2026-04-19 16:27 by 中豫男
[考研] 307中医考研调剂 +9 于以采蘩 2026-04-14 9/450 2026-04-19 08:41 by 烟雨流涯
[考研] 294求调剂 +15 淡然654321 2026-04-15 15/750 2026-04-19 08:20 by cuisz
[考研] 0854求调剂 +23 门路摸摸 2026-04-15 27/1350 2026-04-19 01:59 by 烟雨流涯
[考研] 300求调剂 +12 橙a777 2026-04-15 12/600 2026-04-18 23:51 by 路病情
[考研] 322求调剂 +6 tekuzu 2026-04-17 6/300 2026-04-17 13:48 by Espannnnnol
信息提示
请填处理意见