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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 该物质是对称结构，计算结果却是HOMO在左，LUMO在右。
101/1返回列表
查看: 1436  |  回复: 9
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

Illusionist

银虫 (正式写手)

[求助] 该物质是对称结构，计算结果却是HOMO在左，LUMO在右。

安装g09说明书，计算荧光光谱的方法算了一遍。

又有直接使用
# opt freq td=(nstates=6,root=1) b3lyp/6-31+g(d,p)
scrf=(solvent=n,n-dimethylformamide) geom=connectivity
然后SP算了一遍。

结果都是一样的。
该物质是对称结构，在DMF溶剂中计算，结果却是HOMO在左，LUMO在右，不再对称了。


是溶剂选择的模型有问题
还是什么？

新手，求教！
不胜感激。

（不方便上图，抱歉）
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2011-12-03 20:56:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖置顶 ( 共有2个 )

Illusionist

银虫 (正式写手)
引用回帖:
2楼: Originally posted by tiaozaoququ at 2011-12-04 09:12:09:
本人碰到同样的问题，关注一下

是不是模型的问题？
把DMF溶剂的极性指向的随机分布变成了单一方向造成的啊··
求牛哥哥解答啊···
赞一下 回复此楼
3楼2011-12-04 10:52:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Illusionist

银虫 (正式写手)
Illusionist: 回帖置顶 2011-12-04 14:04:30
=================================
Model                : PCM.
Atomic radii         : UFF (Universal Force Field).
Polarization charges : Total charges.
Charge compensation  : None.
Solution method      : Matrix inversion.
Cavity type          : Scaled VdW (van der Waals Surface) (Alpha=1.100).
Cavity algorithm     : GePol (No added spheres)
                        Default sphere list used, NSphG=   56.
                        Lebedev-Laikov grids with approx.  5.0 points / Ang**2.
                        Smoothing algorithm: Karplus/York (Gamma=1.0000).
                        Polarization charges: spherical gaussians, with
                                              point-specific exponents (IZeta= 3).
                        Self-potential: point-specific (ISelfS= 7).
                        Self-field    : sphere-specific E.n sum rule (ISelfD= 2).
1st derivatives      : Analytical E(r).r(x)/FMM algorithm (CHGder, D1EAlg=3).
                        Cavity 1st derivative terms included.
Solvent              : n,n-DiMethylFormamide, Eps=  37.219000 Eps(inf)=   2.046330
赞一下 回复此楼
4楼2011-12-04 11:53:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

tiaozaoququ

木虫 (著名写手)

【答案】应助回帖

Illusionist(金币+1): 2011-12-04 13:52:08
本人碰到同样的问题，关注一下
赞一下 回复此楼
再牛逼的选择也需要更傻逼的坚持
2楼2011-12-04 09:12:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Jasminer

铁杆木虫 (著名写手)

【答案】应助回帖

Illusionist(金币+2): 谢谢 2011-12-04 13:51:59
1、不加pcm是不是这样？
2、对称点群是否被识别？
只是提示一下，没遇到过这种情况。
赞一下 回复此楼
5楼2011-12-04 12:48:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Illusionist

银虫 (正式写手)
内容已删除
回复此楼
6楼2011-12-04 13:51:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Illusionist

银虫 (正式写手)
按照说明书上计算，到第四步的时候，没有使用OPT=RCFC，而直接使用OPT。因为使用OPT=RCFC无法进行。这一步之前，看的HOMO，LUMO
都是对称的·。
Step 1: Ground state geometry optimization and frequencies (equilibrium solvation). This is a standard Opt Freq calculation on the ground state including PCM equilibrium solvation.

%chk=01-ac
# B3LYP/6-31+G(d,p) Opt Freq SCRF=(Solvent=Ethanol)

Acetaldehyde ground state

0 1
C
C,1,RA
X,2,1.,1,A
O,2,RB,3,A,1,180.,0
X,1,1.,2,90.,3,0.,0
H,1,R1,2,A1,5,0.,0
H,1,R23,2,A23,5,B23,0
H,1,R23,2,A23,5,-B23,0
H,2,R4,1,A4,3,180.,0

RA=1.53643
RB=1.21718
R1=1.08516
R23=1.08688
R4=1.10433
A=62.1511
A1=110.51212
A23=109.88119
A4=114.26114
B23=120.56468
Step 2: Vertical excitation with linear response solvation. This is a TD-DFT calculation of the vertical excitation, therefore at the ground state equilibrium geometry, with the default solvation: linear response, non-equilibrium. We perform a single-point TD-DFT calculation, which defaults to non-equilibrium solvation. The results of this job will be used to identify which state or states are of interest and their ordering. These results give a reasonable description of the solvation of the excited state, but not quite as good as that from a state-specific solvation calculation. In this case, we see that the n->π* state is the first excited state. Next, we will use the state-specific method to produce a better description of the vertical excitation step.

%chk=02-ac
# B3LYP/6-31+G(d,p) TD=NStates=6 SCRF=(Solvent=Ethanol)
  Geom=Check Guess=Read

Acetaldehyde: linear response vertical excited states

0 1
Step 3: State-specific solvation of the vertical excitation. This will require two job steps: first the ground state calculation is done, specifying NonEq=write in the PCM input section, in order to store the information about non-equilibrium solvation based on the ground state. Second, the actual state-specific calculation is done, reading in the necessary information for non-equilibrium solvation using NonEq=read.

%chk=03-ac
# B3LYP/6-31+G(d,p) SCRF=(Solvent=Ethanol,Read)
  Geom=Check Guess=Read

Acetaldehyde: prepare for state-specific non-eq solvation
by saving the solvent reaction field from the ground state

0 1

NonEq=write

--link1--
%chk=03-ac
# B3LYP/6-31+G(d,p) TD(NStates=6,Root=1)
  SCRF=(Solvent=Ethanol,ExternalIteration,Read)
  Geom=Check Guess=Read

Acetaldehyde: read non-eq solvation from ground state and
compute energy of the first excited with the state-specific method

0 1

NonEq=read
Step 4: Relaxation of the excited state geometry. Next, we perform a TD-DFT geometry optimization, with equilibrium, linear response solvation, in order to find the minimum energy point on the excited state potential energy surface. Since this is a TD-DFT optimization, the program defaults to equilibrium solvation. As is typical of such cases, the molecule has a plane of symmetry in the ground state but the symmetry is broken in the excited state, so the ground state geometry is perturbed slightly to break symmetry at the start of the optimization.

%chk=04-ac
# B3LYP/6-31+G(d,p) TD=(Read,NStates=6,Root=1) SCRF=(Solvent=Ethanol)
  Geom=Modify Guess=Read Opt=RCFC

Acetaldehyde: excited state opt
Modify geometry to break Cs symmetry
since first excited state is A"

0 1

4 1 2 3 10.0
5 1 2 7 -50.0
赞一下 回复此楼
7楼2011-12-04 14:09:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sjsznh

银虫 (初入文坛)

【答案】应助回帖

Illusionist(金币+5): 本来是左右对称的结构，但是现在HOMO和LUMO都不对称。ＨＯＭＯ只在左边，ＬＵＭＯ只在右边。即使是反转的问题，ＨＯＭＯ，ＬＵＭＯ本身也应该对称才对啊· 2011-12-04 15:06:33
既然是对称结构，那么在左边和在右边应该是没啥分别的吧。
你这种情况是因为Gaussian 09在有的优化过程中会把整个分子翻转180度造成的，在最新的C01版本中不再存在这个问题。具体可以查看C01的Release Notes。
http://www.gaussian.com/g_tech/rel_notes.pdf
赞一下 回复此楼
8楼2011-12-04 14:25:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sjsznh

银虫 (初入文坛)

【答案】应助回帖

引用回帖:
8楼: Originally posted by sjsznh at 2011-12-04 14:25:56:
既然是对称结构，那么在左边和在右边应该是没啥分别的吧。
你这种情况是因为Gaussian 09在有的优化过程中会把整个分子翻转180度造成的，在最新的C01版本中不再存在这个问题。具体可以查看C01的Release Notes。
...

Rel_Notes中关于避免翻转180度的说明：
赞一下 回复此楼
9楼2011-12-04 14:31:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Jasminer

铁杆木虫 (著名写手)

【答案】应助回帖

★ ★
Illusionist(金币+5): 我freq了两天了，没出来。如果是结构没有优化好，那该怎么做，在优化一遍，会不会又是这结构。··· 2011-12-04 15:54:31
ben_ladeng(金币+2): 专家考核存档 2011-12-04 16:45:31
应该不是ls说的那么回事。
觉得激发态的几何或是波函数没有优化好，激发态和基态点群是否应该一致。
看需不需要做一下激发态的freq以及相应几何的stable测试。
赞一下(1人) 回复此楼
10楼2011-12-04 15:28:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Illusionist 的主题更新
101/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭