| 10 | 1/1 | 返回列表 |
| 查看: 1681 | 回復(fù): 9 | ||
Illusionist銀蟲 (正式寫手)
|
[求助]
該物質(zhì)是對(duì)稱結(jié)構(gòu),計(jì)算結(jié)果卻是HOMO在左,LUMO在右。
|
|
|
安裝g09說明書,計(jì)算熒光光譜的方法算了一遍。 又有直接使用 # opt freq td=(nstates=6,root=1) b3lyp/6-31+g(d,p) scrf=(solvent=n,n-dimethylformamide) geom=connectivity 然后SP算了一遍。 結(jié)果都是一樣的。 該物質(zhì)是對(duì)稱結(jié)構(gòu),在DMF溶劑中計(jì)算,結(jié)果卻是HOMO在左,LUMO在右,不再對(duì)稱了。 是溶劑選擇的模型有問題 還是什么? 新手,求教! 不勝感激。 (不方便上圖,抱歉) |
銀蟲 (正式寫手)
銀蟲 (正式寫手)
|
================================= 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 |
木蟲 (著名寫手)

鐵桿木蟲 (著名寫手)
銀蟲 (正式寫手)
銀蟲 (正式寫手)
|
按照說明書上計(jì)算,到第四步的時(shí)候,沒有使用OPT=RCFC,而直接使用OPT。因?yàn)槭褂肙PT=RCFC無法進(jìn)行。這一步之前,看的HOMO,LUMO 都是對(duì)稱的·。 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 |
銀蟲 (初入文壇)
|
既然是對(duì)稱結(jié)構(gòu),那么在左邊和在右邊應(yīng)該是沒啥分別的吧。 你這種情況是因?yàn)镚aussian 09在有的優(yōu)化過程中會(huì)把整個(gè)分子翻轉(zhuǎn)180度造成的,在最新的C01版本中不再存在這個(gè)問題。具體可以查看C01的Release Notes。 http://www.gaussian.com/g_tech/rel_notes.pdf |
銀蟲 (初入文壇)
鐵桿木蟲 (著名寫手)
|
應(yīng)該不是ls說的那么回事。 覺得激發(fā)態(tài)的幾何或是波函數(shù)沒有優(yōu)化好,激發(fā)態(tài)和基態(tài)點(diǎn)群是否應(yīng)該一致。 看需不需要做一下激發(fā)態(tài)的freq以及相應(yīng)幾何的stable測試。 |
| 10 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 322求調(diào)劑 +4 | 宋明欣 2026-03-27 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 295求調(diào)劑 +4 | 1428151015 2026-03-27 | 5/250 |
|
|
[考研] 考研調(diào)劑 +4 | Sanmu-124 2026-03-26 | 4/200 |
|
|
[考研] 安徽大學(xué)專碩生物與醫(yī)藥專業(yè)(086000)324分,英語已過四六級(jí),六級(jí)521,求調(diào)劑 +4 | 美味可樂雞翅 2026-03-26 | 4/200 |
|
|
[考研] 調(diào)劑 +3 | 李嘉圖·S·路 2026-03-27 | 3/150 |
|
|
[考研] 求調(diào)劑,一志愿 南京航空航天大學(xué)大學(xué) ,080500材料科學(xué)與工程學(xué)碩 +4 | @taotao 2026-03-26 | 5/250 |
|
|
[考研] 324求調(diào)劑 +4 | wysyjs25 2026-03-21 | 4/200 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3 | 崔wj 2026-03-26 | 3/150 |
|
|
[考研] 總分322求生物學(xué)/生化與分子/生物信息學(xué)相關(guān)調(diào)劑 +5 | 星沉uu 2026-03-26 | 6/300 |
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +3 | 邱gl 2026-03-26 | 6/300 |
|
|
[考研] 281求調(diào)劑 +6 | Koxui 2026-03-24 | 7/350 |
|
|
[考研] 一志愿 南京郵電大學(xué) 288分 材料考研 求調(diào)劑 +3 | jl0720 2026-03-26 | 3/150 |
|
|
[考研] 296求調(diào)劑 +4 | 汪!?! 2026-03-25 | 7/350 |
|
|
[考研] 一志愿吉林大學(xué)材料與化工303分求調(diào)劑 +4 | 為學(xué)666 2026-03-24 | 4/200 |
|
|
[考研] B區(qū)考研調(diào)劑 +4 | yqdszhdap- 2026-03-22 | 5/250 |
|
|
[考研] 材料專碩找調(diào)劑 +5 | 哈哈哈吼吼吼哈 2026-03-23 | 5/250 |
|
|
[考研] 化工專碩求調(diào)劑 +3 | question挽風(fēng) 2026-03-24 | 3/150 |
|
|
[考研] 接收2026碩士調(diào)劑(學(xué)碩+專碩) +4 | allen-yin 2026-03-23 | 6/300 |
|
|
[考研] 石河子大學(xué)(211、雙一流)碩博研究生長期招生公告 +3 | 李子目 2026-03-22 | 3/150 |
|
|
[考研] 求助 +5 | 夢里的無言 2026-03-21 | 6/300 |
|