| 查看: 2323 | 回复: 5 | ||
xiaoningli新虫 (小有名气)
|
[求助]
molpro输出文件的对称性问题 已有1人参与
|
|
我想问一下，molpro的输出文件中，怎样看出她的对称性的不可约表示有几个啊？下面是我运行的例子的输出文件， ATOMIC COORDINATES NR ATOM CHARGE X Y Z 1 C 6.00 0.000000000 0.000000000 -1.115737204 2 O 8.00 0.000000000 0.000000000 1.117918915 3 H1 1.00 0.000000000 1.762587230 -2.224791323 4 H2 1.00 0.000000000 -1.762587230 -2.224791323 Bond lengths in Bohr (Angstrom) 1-2 2.233656119 1-3 2.082478039 1-4 2.082478039 (1.182000000) (1.102000000) (1.102000000) Bond angles 2-1-3 122.17890000 2-1-4 122.17890000 3-1-4 115.64220000 BASIS DATA Nr Sym Nuc Type Exponents Contraction coefficients 1.1 A1 1 1s 6665.000000 0.000692 -0.000146 0.000000 2.1 A1 1000.000000 0.005329 -0.001154 0.000000 3.1 A1 228.000000 0.027077 -0.005725 0.000000 64.710000 0.101718 -0.023312 0.000000 21.060000 0.274740 -0.063955 0.000000 7.495000 0.448564 -0.149981 0.000000 2.797000 0.285074 -0.127262 0.000000 0.521500 0.015204 0.544529 0.000000 0.159600 -0.003191 0.580496 1.000000 4.1 A1 1 2pz 9.439000 0.038109 0.000000 5.1 A1 2.002000 0.209480 0.000000 0.545600 0.508557 0.000000 0.151700 0.468842 1.000000 6.1 A1 1 3d0 0.550000 1.000000 7.1 A1 1 3d2+ 0.550000 1.000000 8.1 A1 2 1s 11720.000000 0.000710 -0.000160 0.000000 9.1 A1 1759.000000 0.005470 -0.001263 0.000000 10.1 A1 400.800000 0.027837 -0.006267 0.000000 113.700000 0.104800 -0.025716 0.000000 37.030000 0.283062 -0.070924 0.000000 13.270000 0.448719 -0.165411 0.000000 5.025000 0.270952 -0.116955 0.000000 1.013000 0.015458 0.557368 0.000000 0.302300 -0.002585 0.572759 1.000000 11.1 A1 2 2pz 17.700000 0.043018 0.000000 12.1 A1 3.854000 0.228913 0.000000 1.046000 0.508728 0.000000 0.275300 0.460531 1.000000 13.1 A1 2 3d0 1.185000 1.000000 14.1 A1 2 3d2+ 1.185000 1.000000 15.1 A1 3 1s 13.010000 0.019685 0.000000 16.1 A1 4 1s 1.962000 0.137977 0.000000 0.444600 0.478148 0.000000 0.122000 0.501240 1.000000 17.1 A1 3 2pz 0.727000 1.000000 4 2pz 18.1 A1 3 2py 0.727000 1.000000 4 -2py 1.2 B1 1 2px 9.439000 0.038109 0.000000 2.2 B1 2.002000 0.209480 0.000000 0.545600 0.508557 0.000000 0.151700 0.468842 1.000000 3.2 B1 1 3d1+ 0.550000 1.000000 4.2 B1 2 2px 17.700000 0.043018 0.000000 5.2 B1 3.854000 0.228913 0.000000 1.046000 0.508728 0.000000 0.275300 0.460531 1.000000 6.2 B1 2 3d1+ 1.185000 1.000000 7.2 B1 3 2px 0.727000 1.000000 4 2px 1.3 B2 1 2py 9.439000 0.038109 0.000000 2.3 B2 2.002000 0.209480 0.000000 0.545600 0.508557 0.000000 0.151700 0.468842 1.000000 3.3 B2 1 3d1- 0.550000 1.000000 4.3 B2 2 2py 17.700000 0.043018 0.000000 5.3 B2 3.854000 0.228913 0.000000 1.046000 0.508728 0.000000 0.275300 0.460531 1.000000 6.3 B2 2 3d1- 1.185000 1.000000 7.3 B2 3 1s 13.010000 0.019685 0.000000 8.3 B2 4 1s 1.962000 0.137977 0.000000 0.444600 0.478148 0.000000 0.122000 0.501240 1.000000 9.3 B2 3 2py 0.727000 1.000000 4 2py 10.3 B2 3 2pz 0.727000 1.000000 4 -2pz 1.4 A2 1 3d2- 0.550000 1.000000 2.4 A2 2 3d2- 1.185000 1.000000 3.4 A2 3 2px 0.727000 1.000000 4 -2px NUCLEAR CHARGE: 16 NUMBER OF PRIMITIVE AOS: 68 NUMBER OF SYMMETRY AOS: 66 NUMBER OF CONTRACTIONS: 38 ( 18A1 + 7B1 + 10B2 + 3A2 ) NUMBER OF CORE ORBITALS: 2 ( 2A1 + 0B1 + 0B2 + 0A2 ) NUMBER OF VALENCE ORBITALS: 10 ( 5A1 + 2B1 + 3B2 + 0A2 ) 然后看到手册中解释的是，为c2v点群，这个我知道，5个轨道属于不可约 表示a1，一个属于b1，两个b2，这句话是怎么得到的啊？我看上面那么多a1，b1的，怎么确定哪几个是有用的呢？ 谢谢你们！ |
回复此楼» 猜你喜欢
高效微纳米储能材料与产业化推进实验室 苏州 招收硕士研究生（含调剂）
已经有0人回复
-
温州大学新能源材料微观溯源课题组招收2026级博士研究生（材料、化学方向）
已经有6人回复
-
物理化学论文润色/翻译怎么收费?
已经有97人回复
-
哈尔滨工业大学（深圳）绿色化学团队招收2026年秋季入学博士生
已经有11人回复
-
材料与化工 323 英一+数二+物化，一志愿：哈工大 本人本科双一流
已经有5人回复
-
26申博
已经有2人回复
-
澳大利亚昆士兰科技大学PhD博士奖学金
已经有0人回复
-
巴黎西岱大学（Université Paris Cité）Fayna MAMMERI 导师寻找CSC公派博士生
已经有9人回复
-
中国散裂中子源诚聘“AI材料计算”方向博士后/职工
已经有75人回复
-
【211博士招生】环境化学、地学、毒理方向，擅长污染物迁移转化降解、计算模拟者优先
已经有15人回复
-
求助如何利用vasp计算磁光克尔效应
已经有1人回复
仰视幸福1127
木虫 (正式写手)
- 应助: 7 (幼儿园)
- 金币: 1497.3
- 红花: 3
- 帖子: 598
- 在线: 145.4小时
- 虫号: 2132854
- 注册: 2012-11-17
- 性别: MM
- 专业: 原子和分子物理
2楼2014-11-24 20:16:49
3楼2018-06-07 22:45:23
【答案】应助回帖
|
在 ref：https://www.molpro.net/info/2015.1/doc/manual/node36.html 中得到C2v的不可约molpro表示，A1 B1 B2 A2 你在给出的信息中得到的是基组信息 用 print basis 全局命令激活打印 得到的是该基组在c2v对称性下的不可约归属，你也可以手动推导 会发现结果是一样的。 如何推导呢？ 基组也就是各个原子轨道，那么就依照分别对应的对称性操作来决定，举个例子 C上的1s 轨道，已知c2v的对称性操作有：E C2 sigma（xz） sigma（yz），后两个为经过z轴的镜面操作。 那么该s轨道在这四个对称性操作下得到的不可约表示是： 1 1 1 1。得到的irrep是A1 再举个例子， H上的1s 轨道，按照同样方式得到（定义分子在xz平面上）的不可约表示为：1 -1 1 -1。得到的irrep是B2 大致上就是通过以上推导方式得到原子基组轨道的不可约归属。这个事情molpro已经帮你做了，所以并不用手动操作，你只需要知道这是怎么来的就行了。 至于你给出的信息，对于下面的计算有用的大概就是： NUMBER OF CORE ORBITALS: 2 ( 2A1 + 0B1 + 0B2 + 0A2 ) NUMBER OF VALENCE ORBITALS: 10 ( 5A1 + 2B1 + 3B2 + 0A2 ) 这些了。 这些信息给你molpro默认定义的core轨道以及valence轨道，可以看出，A1下有两个core轨道，很显然是C跟O上的两个1s轨道。valence轨道则是由以上valence原子轨道不可约归属相加得到，原理是一样的。 同样的信息会在计算结束后的分子轨道中读得。你所说的“有用的”信息，应该是指计算后的分子轨道占据，可以用 print orbital全局命令激活打印。 |
4楼2018-06-08 02:41:03
5楼2018-06-08 02:43:16
6楼2018-09-03 22:36:11
