24小时热门版块排行榜    

查看: 2305  |  回复: 15
【有奖交流】积极回复本帖子,参与交流,就有机会分得作者 0011608 的 96 个金币

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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 2024国社科通讯评审 +14 qsd10086 2024-06-13 27/1350 2024-06-18 09:48 by robin_work
[找工作] 药学硕士找不到工作,打算去做科研助理了 +6 pom戴墨镜 2024-06-14 16/800 2024-06-18 09:27 by pom戴墨镜
[基金申请] 最新消息:2024国自然下载文件名变了 +16 dream200606 2024-06-16 29/1450 2024-06-18 09:26 by 可7可D
[找工作] 初始合伙人来啦!(生物试剂耗材标准品) +12 欢快的小科研人 2024-06-15 22/1100 2024-06-18 06:50 by swq120
[文学芳草园] 累并快乐着 +15 MYHLD521 2024-06-14 15/750 2024-06-18 01:28 by liyeqik
[硕博家园] 硕博巨婴,也许才刚刚开始 +33 SNaiL1995 2024-06-12 85/4250 2024-06-18 00:06 by 博士求上上岸
[基金申请] 博后基金,博管会会提前知道消息吗? +4 yuyiang 2024-06-13 4/200 2024-06-16 11:40 by yangyuzhong4
[论文投稿] 投稿被一个审稿人恶意评审了怎么样? +5 1chen 2024-06-14 7/350 2024-06-15 23:15 by xy66xy
[基金申请] BO4的YQ答辩通知发布了吗? +6 博学笃行 2024-06-11 6/300 2024-06-15 16:04 by 悲催科研狗
[基金申请] 为什么我的博后基金还在流动站审核中?不会是学院给我卡了吧? +14 王凯12 2024-06-13 26/1300 2024-06-15 15:22 by 好人与坏人
[论文投稿] 投了一篇4区的SCI,审稿人一个拒稿,一个小修,编辑给了大修。 +9 安稳22123 2024-06-13 10/500 2024-06-14 23:45 by jurkat.1640
[考研] 物理化学一对一辅导 +3 林大diao 2024-06-12 5/250 2024-06-14 20:57 by 林大diao
[基金申请] 工材E口JQ有消息了吗 +4 babyduck 2024-06-11 4/200 2024-06-14 17:23 by firepick
[基金申请] E12面上申请 +4 汉风之遗 2024-06-13 4/200 2024-06-14 15:28 by 天外飞去来
[考博] 申博找导师 +4 疏影横斜水清浅3 2024-06-13 6/300 2024-06-14 14:31 by zxl_1105
[基金申请] 75批博后基金 +10 kyukitu 2024-06-13 13/650 2024-06-14 10:31 by kyukitu
[论文投稿] 文章proof要求使用机构的邮箱 5+3 不可不信缘 2024-06-11 11/550 2024-06-14 07:00 by 3001160025
[硕博家园] 机械研究生如何拿到年薪40+w +13 阿巴阿巴哦哦 2024-06-11 15/750 2024-06-13 15:40 by 113745685
[有机交流] 原料反应完了,怎么知道是产物还是中间体 +6 小胡在努力 2024-06-11 8/400 2024-06-13 13:33 by 091602
[硕博家园] 考博应该从哪里开始准备 +6 亮姐2019 2024-06-11 8/400 2024-06-13 12:53 by 帕秋莉GO
信息提示
请填处理意见