24小时热门版块排行榜    

查看: 3405  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 明天应该可查了!? +5 chengyan1220 2026-08-23 5/250 2026-08-24 23:28 by 我4大白菜
[基金申请] 没有任何消息-是不是就凉了 +7 图啦图啦 2026-08-24 8/400 2026-08-24 22:00 by maomao_da
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 0XLacIJUOj8D 2026-08-24 3/150 2026-08-24 21:59 by BZKMTicpDhFj
[基金申请] 人气不行了 +10 fansofjerry 2026-08-21 10/500 2026-08-24 21:03 by zhanghaozhu
[基金申请] 能否退出参与的面上项目解除限项 +21 koalala 2026-08-24 24/1200 2026-08-24 19:25 by 家与远方
[基金申请] 2026国自然函评费到账 +17 羊腰板 2026-08-21 19/950 2026-08-24 16:52 by iaeyuan
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 21/1050 2026-08-24 14:21 by refreshing11
[基金申请] 朋友圈看到的 +7 wangzilk 2026-08-18 9/450 2026-08-24 10:50 by cmrandy
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 什么时候开奖? +10 CrisMessi 2026-08-18 11/550 2026-08-24 06:50 by 开心的小狮子
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +4 process2012 2026-08-23 5/250 2026-08-23 19:58 by jurkat.1640
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[基金申请] 今天基金会出结果吗?20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
请填处理意见