24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 3084  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0703总分331求调剂 +12 ZY-05 2026-04-04 16/800 2026-04-06 17:54 by lijunpoly
[考研] 求调剂到材料 +5 程9915 2026-04-06 5/250 2026-04-06 15:21 by yulian1987
[考研] 362求调剂一志愿中国石油大学 +4 我要考大 2026-04-06 6/300 2026-04-06 14:11 by 无际的草原
[考研] 一志愿北交大材料工程总分358求调剂 +6 cs0106 2026-04-05 6/300 2026-04-05 16:34 by imissbao
[考研] 325求调剂 +4 春风不借意 2026-04-04 4/200 2026-04-04 22:08 by 啵啵啵0119
[考研] 求调剂 +3 ffyyu 2026-04-02 3/150 2026-04-04 19:03 by 蓝云思雨
[考研] 一志愿南农090401,268,求调剂 +5 一木鸟然 2026-04-04 5/250 2026-04-04 17:07 by babysonlkd
[考研] 一志愿东北大学085901土木专硕345求调剂 +3 zxt11111 2026-04-04 3/150 2026-04-04 14:21 by 土木硕士招生
[考研] 280求调剂 +21 咕噜晓晓 2026-04-02 22/1100 2026-04-04 11:12 by 猪会飞
[考研] 求调剂 +8 akdhjs 2026-04-03 8/400 2026-04-03 18:17 by 戴维ING
[考研] 专硕 351 086100 也是考的材科基 本科也是材料 +8 202451007219 2026-04-02 8/400 2026-04-03 09:50 by 蓝云思雨
[考研] 一志愿安徽大学0817化学工程与技术,求调剂 +14 我不是只因 2026-04-02 15/750 2026-04-03 09:49 by 蓝云思雨
[考研] 一志愿a区211,085601-307分求调剂 +13 党嘉豪 2026-03-31 26/1300 2026-04-03 08:33 by 495374996
[考研] 321求调剂 一志愿 浙江工业大学生物医药 +5 嘿嘿HC 2026-04-01 6/300 2026-04-02 15:23 by sophie2180
[考研] 081200-11408-276学硕求调剂 +3 崔wj 2026-04-02 3/150 2026-04-02 15:06 by cal0306
[考研] 土木304求调剂 +6 兔突突突, 2026-03-31 7/350 2026-04-02 09:06 by coolminer
[考研] 303分 0807学硕求调剂 +3 TYC3632 2026-04-01 3/150 2026-04-01 19:24 by lwk2004
[考研] 353求调剂 +4 拉钩不许变 2026-04-01 4/200 2026-04-01 18:10 by 记事本2026
[考研] 材料调剂 +10 Eujd1 2026-03-31 11/550 2026-04-01 11:23 by ivanqyq
[考研] 一志愿西电085401数一英一299求调剂 六级521 +4 爱吃大鸭梨 2026-03-31 4/200 2026-03-31 11:51 by 搏击518
信息提示
请填处理意见