24小时热门版块排行榜    

查看: 3001  |  回复: 15

0011608

铁虫 (小有名气)

[交流] 【求助】如何在以下高斯03输出文件中得到二阶超极化率 已有6人参与

<γ>的各个分量是多少,在哪里?

Dipole moment (field-independent basis, Debye):
    X=     0.6915    Y=     0.2627    Z=     0.0000  Tot=     0.7397
Quadrupole moment (field-independent basis, Debye-Ang):
   XX=   -72.1916   YY=   -73.6973   ZZ=   -89.9511
   XY=    -2.3909   XZ=     0.0000   YZ=     0.0000
Traceless Quadrupole moment (field-independent basis, Debye-Ang):
   XX=     6.4218   YY=     4.9160   ZZ=   -11.3377
   XY=    -2.3909   XZ=     0.0000   YZ=     0.0000
Octapole moment (field-independent basis, Debye-Ang**2):
  XXX=    40.6015  YYY=    -1.8219  ZZZ=     0.0000  XYY=    -3.9133
  XXY=    15.6956  XXZ=     0.0000  XZZ=    -6.2461  YZZ=    -0.7539
  YYZ=     0.0000  XYZ=    -0.0004
Hexadecapole moment (field-independent basis, Debye-Ang**3):
XXXX= -4562.8469 YYYY=  -492.1221 ZZZZ=  -109.5764 XXXY=   -94.1069
XXXZ=    -0.0017 YYYX=     2.9680 YYYZ=     0.0000 ZZZX=     0.0001
ZZZY=     0.0000 XXYY=  -797.5517 XXZZ=  -936.9767 YYZZ=  -123.8276
XXYZ=     0.0002 YYXZ=     0.0000 ZZXY=     3.4891
N-N= 7.628912261271D+02 E-N=-3.620480346817D+03  KE= 8.897100382751D+02
  Exact polarizability: 240.264  12.853 136.473   0.000   0.000  80.021
Approx polarizability: 180.402  12.521 143.400   0.000   0.000  88.535
Full mass-weighted force constant matrix:
Low frequencies ---   -0.3003   -0.2018    0.0013    0.0024    0.0025    0.2635
Low frequencies ---   26.2117   69.8678   88.0743
Diagonal vibrational polarizability:
       10.5242669       2.3566937      17.4199031
Diagonal vibrational hyperpolarizability:
      -34.2316821      -8.3324288      -0.0198722
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

qzhaosdu

金虫 (著名写手)

gaussian输出不了二阶超极化率
Anewday,anewhour,anewminute,anewpeople.
3楼2010-04-13 13:05:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

dukelr

金虫 (小有名气)

0011608(金币+1): 2010-05-23 12:01:10
我也想知道....
2楼2010-04-13 12:48:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dukelr

金虫 (小有名气)

但是 根据输出文件可以有公式算出来伽马值 从而算出来二阶超极化率的吧?
4楼2010-04-13 13:10:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sculhf

荣誉版主 (知名作家)

断翼天使

优秀版主优秀版主

★ ★ ★
aylayl08(金币+3):感谢指点 2010-04-13 21:00
平均静态极化率=(axx+ayy+azz)/3
对于二阶超极化率计算方法有几种近似
比如有限场法(Finite-Field Method)
具体算法可以去查资料
谋事在人成事在天
5楼2010-04-13 16:31:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xtdut

银虫 (小有名气)

可以用完全态(sos)计算!

0011608(金币+2): 2010-05-23 12:01:40
在关键词TD=(full,sos)可以计算完全态求和(sos)方法计算二阶超极化率,可是我不清楚输出的十八个分量在哪里,如何归属,另外一个问题是,如何计算不同频率下的二阶超极化率?
6楼2010-05-07 13:20:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

faqianliu

银虫 (小有名气)

0011608(金币+1): 2010-05-23 12:01:47
帮顶一下,有没有高手指点
7楼2010-05-11 21:17:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunhaitao

木虫 (小有名气)

二阶超极化率=(xxxx+yyyy+zzzz+2xxyy+2xxzz+2yyzz)/5,你去数据库中找那些超极化率的量化计算就好了,有很多的
一点一点进步!
8楼2010-05-12 10:52:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qinyuzhu5233

金虫 (小有名气)

上帝的指纹

请问高手,你的g03.out文件中,关键字用的是什么呀,就能得到您下面的结果
圣经上说:在信的人凡事都能!我信这话
9楼2010-09-13 15:04:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xcyqyz5233

金虫 (小有名气)

记得给红包啊

polar 极化率和超极化率(2)论坛上的一个讨论
2008-10-20 19:21
问题1.
对一个有机染料在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
中这些数值的含义是什么呀?单位呢?
回答
保留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文件里就能很清晰的找到。
问题2.
对于计算得到的超极化率 共有十个数值,他们分别代表那几个张量呀?
是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

回答
以上分量按顺序分别为:
mx
my
mz
axx
axy
ayy
axz
ayz
azz
βxxx
βxxy
βxyy
βyyy
βxxz
βxyz
βyyz
βxzz
βyzz
βzzz
10楼2010-09-14 17:13:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 0011608 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085601专硕,总分342求调剂,地区不限 +4 share_joy 2026-03-16 4/200 2026-03-17 21:32 by hmn_wj
[考研] 085601材料工程专硕求调剂 +5 慕寒mio 2026-03-16 5/250 2026-03-17 21:31 by hmn_wj
[考研] 296求调剂 +5 大口吃饭 身体健 2026-03-13 5/250 2026-03-17 21:05 by 不惑可乐
[考研] 304求调剂 +8 小熊joy 2026-03-14 8/400 2026-03-17 17:29 by ruiyingmiao
[考研] 311求调剂 +8 冬十三 2026-03-15 8/400 2026-03-17 16:59 by ruiyingmiao
[考研] 梁成伟老师课题组欢迎你的加入 +8 一鸭鸭哟 2026-03-14 10/500 2026-03-17 15:07 by 一鸭鸭哟
[考研] 290求调剂 +6 孔志浩 2026-03-12 11/550 2026-03-17 14:41 by 周舟舟77
[考研] 302求调剂 +4 小贾同学123 2026-03-15 8/400 2026-03-17 10:33 by 小贾同学123
[考研] 一志愿985,本科211,0817化学工程与技术319求调剂 +5 Liwangman 2026-03-15 5/250 2026-03-16 17:10 by 我的船我的海
[基金申请] 今年的国基金是打分制吗? 50+3 zhanghaozhu 2026-03-14 3/150 2026-03-16 17:07 by 北京莱茵润色
[考研] 0703化学调剂 +6 妮妮ninicgb 2026-03-15 9/450 2026-03-16 16:40 by houyaoxu
[考研] 285求调剂 +6 ytter 2026-03-12 6/300 2026-03-16 15:05 by njzyff
[考研] 一志愿华中师范071000,325求调剂 +6 RuitingC 2026-03-12 6/300 2026-03-16 14:50 by 可淡不可忘
[考研] 288求调剂 +4 奇点0314 2026-03-14 4/200 2026-03-14 23:04 by JourneyLucky
[考研] 【0703化学调剂】-一志愿华中师范大学-六级475 +5 Becho359 2026-03-11 5/250 2026-03-14 11:35 by 哦哦123
[考研] 0703,333分求调剂 一志愿郑州大学-物理化学 +3 李魔女斗篷 2026-03-11 3/150 2026-03-13 22:24 by JourneyLucky
[考研] 290求调剂 +7 ADT 2026-03-12 7/350 2026-03-13 15:17 by JourneyLucky
[考研] 工科调剂 +4 Jiang191123! 2026-03-11 4/200 2026-03-13 15:15 by Miko19
[考博] 福州大学杨黄浩课题组招收2026年专业学位博士研究生,2026.03.20截止 +3 Xiangyu_ou 2026-03-12 3/150 2026-03-13 09:36 by duanwu655
[考研] 283求调剂,材料、化工皆可 +8 苏打水7777 2026-03-11 10/500 2026-03-13 09:06 by Linda Hu
信息提示
请填处理意见