| 10 | 1/1 | 返回列表 |
| 查看: 1687 | 回復(fù): 9 | ||
Illusionist銀蟲 (正式寫手)
|
[求助]
該物質(zhì)是對(duì)稱結(jié)構(gòu),計(jì)算結(jié)果卻是HOMO在左,LUMO在右。
|
|
安裝g09說(shuō)明書,計(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 |
鐵桿木蟲 (著名寫手)
銀蟲 (正式寫手)
銀蟲 (正式寫手)
|
按照說(shuō)明書上計(jì)算,到第四步的時(shí)候,沒有使用OPT=RCFC,而直接使用OPT。因?yàn)槭褂肙PT=RCFC無(wú)法進(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)化過(guò)程中會(huì)把整個(gè)分子翻轉(zhuǎn)180度造成的,在最新的C01版本中不再存在這個(gè)問題。具體可以查看C01的Release Notes。 http://www.gaussian.com/g_tech/rel_notes.pdf |
銀蟲 (初入文壇)
鐵桿木蟲 (著名寫手)
|
應(yīng)該不是ls說(shuō)的那么回事。 覺得激發(fā)態(tài)的幾何或是波函數(shù)沒有優(yōu)化好,激發(fā)態(tài)和基態(tài)點(diǎn)群是否應(yīng)該一致。 看需不需要做一下激發(fā)態(tài)的freq以及相應(yīng)幾何的stable測(cè)試。 |
| 10 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 320分,材料與化工專業(yè),求調(diào)劑 +9 | 一定上岸aaa 2026-03-27 | 13/650 |
|
|---|---|---|---|---|
|
[考研] 本科新能源科學(xué)與工程,一志愿華理能動(dòng)285求調(diào)劑 +3 | AZMK 2026-03-28 | 5/250 |
|
|
[考研] 311求調(diào)劑 +4 | 冬十三 2026-03-24 | 4/200 |
|
|
[考研] 321求調(diào)劑 +6 | 材料cailiao 2026-03-21 | 6/300 |
|
|
[考研] 299求調(diào)劑 +7 | 嗯嗯嗯嗯2 2026-03-27 | 7/350 |
|
|
[考研] 數(shù)一英一271專碩(085401)求調(diào)劑,可跨 +4 | 前行必有光 2026-03-28 | 5/250 |
|
|
[考研] 331環(huán)境科學(xué)與工程求調(diào)劑 +3 | 熠然好運(yùn)氣 2026-03-27 | 3/150 |
|
|
[考研] 330一志愿中國(guó)海洋大學(xué) 化學(xué)工程 085602 有讀博意愿 求調(diào)劑 +3 | wywy.. 2026-03-27 | 4/200 |
|
|
[考研] 0703化學(xué)求調(diào)劑,各位老師看看我!。 +5 | 祁祺祺 2026-03-25 | 5/250 |
|
|
[考研] 265求調(diào)劑11408 +3 | 劉小鹿lu 2026-03-27 | 3/150 |
|
|
[考研] 085600,材料與化工321分,求調(diào)劑 +9 | 大饞小子 2026-03-27 | 9/450 |
|
|
[考研] 085601 材料工程 313分 求調(diào)劑 +5 | Ong3 2026-03-27 | 5/250 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3 | 崔wj 2026-03-26 | 3/150 |
|
|
[考研] 352求調(diào)劑 +4 | 大米飯! 2026-03-22 | 4/200 |
|
|
[考研] 一志愿河工大 081700 276求調(diào)劑 +4 | 地球繞著太陽(yáng)轉(zhuǎn) 2026-03-23 | 4/200 |
|
|
[考研] 07化學(xué)303求調(diào)劑 +5 | 睿08 2026-03-25 | 5/250 |
|
|
[考研] 生物技術(shù)與工程 +3 | 1294608413 2026-03-25 | 4/200 |
|
|
[考研] 0854AI CV方向招收調(diào)劑 +4 | 章小魚567 2026-03-23 | 4/200 |
|
|
[考研] 一志愿武理085500機(jī)械專業(yè)總分300求調(diào)劑 +3 | an10101 2026-03-24 | 7/350 |
|
|
[考研] 305分求調(diào)劑(食品工程) +5 | Sxy112 2026-03-21 | 7/350 |
|