24小时热门版块排行榜

返回列表

当前主题已经存档。

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

zhugezilong

木虫 (著名写手)

应助: 35 (小学生)
金币: 2231.4
散金: 41
帖子: 1123
在线: 278.2小时
虫号: 355130
注册: 2007-04-25
性别: MM
专业: 有机分子功能材料化学

[交流] [求助]Gaussian计算中怎么来求超极化率呀

本人是个菜鸟级的。对一个有机染料在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
中这些数值的含义是什么呀？单位呢？

[ Last edited by zhugezilong on 2008-10-11 at 22:27 ]

回复此楼

» 猜你喜欢

青岛大学化学化工学院分子测量学研究院2026年招收博士研究生已经有0人回复
香港科技大学（广州）诚招电催化方向博士生（2026秋入学）已经有0人回复
物理化学论文润色/翻译怎么收费? 已经有281人回复
2026第二届光电子与半导体器件前沿技术研讨会——光电子赋能·半导体创芯！已经有1人回复
光电子赋能·半导体创芯！已经有1人回复
光电子赋能·半导体创芯！已经有0人回复
求助Cu2+1O的CIF文件(PDF: 05-0667) 已经有1人回复
QE计算电声耦合的时候报错Error in routine lambda (100)wrong or too many modes 已经有2人回复
KAUST（阿卜杜拉国王科技大学）MXene 器件方向博士后招聘已经有0人回复
沙特阿拉伯阿卜杜拉国王科技大学（KAUST）电池方向博士后招聘已经有0人回复
福州大学新能源材料与工程研究院招收2026年入学博士已经有0人回复

1楼 2008-10-11 22:26:09

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

zhugezilong

木虫 (著名写手)

应助: 35 (小学生)
金币: 2231.4
散金: 41
帖子: 1123
在线: 278.2小时
虫号: 355130
注册: 2007-04-25
性别: MM
专业: 有机分子功能材料化学

真是非常感谢。我还想问问
对于计算得到的超极化率共有十个数值，他们分别代表那几个张量呀？
是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

赞一下

回复此楼

3楼2008-10-12 20:42:10

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 5 个回答

yuhuobuku

木虫 (正式写手)

应助: 2 (幼儿园)
贵宾: 0.09
金币: 1661.1
帖子: 332
在线: 6.4小时
虫号: 577652
注册: 2008-06-27
性别: MM
专业: 量子化学

★ ★ ★ ★ ★ ★ ★ ★ ★
zzgyb(金币+4,VIP+0):谢谢你的参与，欢迎再次光临计算模拟版！
zhugezilong(金币+5,VIP+0):谢谢解答

保留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文件里就能很清晰的找到。

[ Last edited by yuhuobuku on 2008-10-12 at 13:44 ]

赞一下(11人)

回复此楼

薄酒可以忘忧，丑妻可以白头，徐行不必驷马，称身不必狐裘平淡才是生活的主旋律