| 查看: 1757 | 回复: 9 | |||
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在右，不再对称了。 是溶剂选择的模型有问题 还是什么？ 新手，求教！ 不胜感激。 （不方便上图，抱歉） |
回复此楼» 猜你喜欢
085404（计算机技术）293分 希望调剂到材料相关专业
已经有0人回复
-
ELISA试验中不容忽视的细节盘点(二)
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有113人回复
-
细胞培养，这22个细节一定要注意！（一）
已经有0人回复
-
CSC访学博后项目获批，外方学校暂停合作该怎么办？
已经有51人回复
-
七嗪类物质合成求助
已经有0人回复
-
分子动力学模拟合作需求
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 大家好，有好心人帮我计算四个分子的homo和lumo轨道能级图吗？急需，可以支付适当费用
已经有4人回复
- 关于HOMO 和LUMO轨道的计算
已经有4人回复
- 怎么算lumo和homo能级啊？用MS
已经有14人回复
- 循环伏安曲线基础问题--氧化峰，还原峰，氧化电位，还原电位及计算HOMO，LUMO的方法
已经有6人回复
- PWscf计算homo lumo或者valence conduction band，问题已解决
已经有14人回复
- 【求助】如何用高斯计算一个分子中的某一部分HOMO和LUMO？
已经有7人回复
- 【求助】Gaussian计算HOMO-LUMO问题
已经有17人回复
- 【求助】计算homo和lumo轨道出错
已经有14人回复
- 【求助】请问，如何在MS中分析HOMO与LUMO，谢谢您的帮助
已经有12人回复
tiaozaoququ
木虫 (著名写手)
- 应助: 5 (幼儿园)
- 金币: 7213.3
- 散金: 63
- 红花: 1
- 帖子: 1202
- 在线: 175.6小时
- 虫号: 1011016
- 注册: 2010-05-04
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
2楼2011-12-04 09:12:09
Illusionist
银虫 (正式写手)
- 应助: 20 (小学生)
- 金币: 232.8
- 散金: 1001
- 红花: 15
- 帖子: 431
- 在线: 329.1小时
- 虫号: 1123329
- 注册: 2010-10-15
- 专业: 半导体晶体与薄膜材料
3楼2011-12-04 10:52:20
Illusionist
银虫 (正式写手)
- 应助: 20 (小学生)
- 金币: 232.8
- 散金: 1001
- 红花: 15
- 帖子: 431
- 在线: 329.1小时
- 虫号: 1123329
- 注册: 2010-10-15
- 专业: 半导体晶体与薄膜材料
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
Jasminer
铁杆木虫 (著名写手)
- 应助: 52 (初中生)
- 金币: 7407.8
- 散金: 1845
- 红花: 27
- 帖子: 1002
- 在线: 728.8小时
- 虫号: 620156
- 注册: 2008-10-08
- 性别: GG
- 专业: 应用有机化学
5楼2011-12-04 12:48:07
Illusionist
银虫 (正式写手)
- 应助: 20 (小学生)
- 金币: 232.8
- 散金: 1001
- 红花: 15
- 帖子: 431
- 在线: 329.1小时
- 虫号: 1123329
- 注册: 2010-10-15
- 专业: 半导体晶体与薄膜材料
6楼2011-12-04 13:51:47
Illusionist
银虫 (正式写手)
- 应助: 20 (小学生)
- 金币: 232.8
- 散金: 1001
- 红花: 15
- 帖子: 431
- 在线: 329.1小时
- 虫号: 1123329
- 注册: 2010-10-15
- 专业: 半导体晶体与薄膜材料
|
按照说明书上计算，到第四步的时候，没有使用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
【答案】应助回帖
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
9楼2011-12-04 14:31:58
Jasminer
铁杆木虫 (著名写手)
- 应助: 52 (初中生)
- 金币: 7407.8
- 散金: 1845
- 红花: 27
- 帖子: 1002
- 在线: 728.8小时
- 虫号: 620156
- 注册: 2008-10-08
- 性别: GG
- 专业: 应用有机化学
10楼2011-12-04 15:28:51
