应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】如何在以下高斯03输出文件中得到二阶超极化率
161/2返回列表12下一页
查看: 3406  |  回复: 15
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2010-04-12 11:00:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

qzhaosdu

金虫 (著名写手)
gaussian输出不了二阶超极化率
赞一下(1人) 回复此楼
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）
具体算法可以去查资料
赞一下(1人) 回复此楼
谋事在人成事在天
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 的主题更新
161/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 0XLacIJUOj8D 2026-08-24 3/150 2026-08-25 01:47 by BZKMTicpDhFj
[硕博家园] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 0XLacIJUOj8D 2026-08-24 6/300 2026-08-25 00:11 by BZKMTicpDhFj
[基金申请] 明天应该可查了！？ +5 chengyan1220 2026-08-23 5/250 2026-08-24 23:28 by 我4大白菜
[基金申请] 人气不行了 +10 fansofjerry 2026-08-21 10/500 2026-08-24 21:03 by zhanghaozhu
[基金申请] 我面上完蛋了 +9 且听虎啸 2026-08-20 10/500 2026-08-24 20:08 by 苦寒来香
[基金申请] 能否退出参与的面上项目解除限项 +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
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 今日不放榜？网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 21/1050 2026-08-24 14:21 by refreshing11
[基金申请] 估计是周四 +4 archvillain 2026-08-18 4/200 2026-08-24 13:53 by zzuzxg
[基金申请] 科研孤儿太难了 +18 我4大白菜 2026-08-20 19/950 2026-08-24 09:47 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
[基金申请] 今天基金会出结果吗？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
[基金申请] 看来今天不会放榜了？ +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 应该是下周三26日公布了吧？ +4 哈哈蛤？ 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 重要消息，中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 明天放榜？ +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 家与远方
信息提示
关闭
请填处理意见
关闭