## Firefly and PC GAMESS-related discussion club

Learn how to ask questions correctly

Re^3: Geometry optimization of First Excited state.....

Amir Nasser Shamkhali
amir_n_shamkhali@yahoo.com

Dear Siddheshwar
UV-Vis spectra in most cases is a collection of vertical excitations. Please note that when you measure UV-Vis spectra by a spectrophotometer, the most important interaction of light with molecules is its electric filed. Due to the movement of molecules, the direction of interaction varies by time, also. Therefore the Hamiltonians of the molecule in various times within a period between t=0 and t=T(1/freq)do not commute and expectation values of Hamiltonians (energies) vary with time. However, this time is very short ( in order of 10^-17 sec) in which geometry of the molecule have not enough opportunity to change.Thus your material is stable. This means that molecules are in chemical equilibrium. So, in order to calculate UV-Vis spectra, you should only optimize the ground state, by your DFT method. Then, taking the optimized structure, you should run a single point TDDFT calculation with some higher vertical excitations (NSTATE keyword). So, for example a simple input is required as follows:
\$CONTRL scftyp=RHF CITYP=TDDFT icharg=0 mult=1
NZVAR=??? Coord=ZMT MAXIT=800 DFTTYP=??? UNITS=ANGS \$END
\$TDDFT NSTATE=7 ISTATE=1 MULT=1 \$END
Please note that when you are interested to phosphorescence, it means that the multiplicity of excited state is different that of the ground state. So, if the ground state has 2S+1=1, then the MULT keyword in \$CONTROL should be 1 and in \$TDDFT should be 3.
Electron-phonon interaction is a solid state effect which is responsible for superconductivity. Please note that Firefly is based on the Gaussian-type basis sets which are localized and are appropriate for molecules. But electron-phonon interaction is delocalized in nature and you should use a code based on plane wave basis sets (such as ABINIT, espresso, WIEN2k, exciting, and etc.).

On Tue Apr 1 '14 11:15pm, Siddheshwar Chopra wrote
--------------------------------------------------
>Dear Sir,
>That means if I have to study UV-Vis Spectra, then studying vertical excitation only, is correct?? If yes then running a CITYP=TDDFT is correct? Is that right? In vertical excitation if the geometry of the molecule (in first excited state???) doesnt change, then optimization still is needed? How time consuming is this job? Does the optimization and oscillator strengths get calculated in the same run RUNTYP=optimize?

>Also how do I decide upon studying the adiabatic excitation for my sample?? Say it is a new sample/molecule.. Do you mean phosphorescence transitions by this? Thanks again for explaining so well Sir.

>Also is there any way I can comment on radiative or non-radiative nature of transitions?? or can I find electron-phonon interactions using Firefly?

>Regards,

>On Tue Apr 1 '14 3:00pm, Amir Nasser Shamkhali wrote
>----------------------------------------------------
>>Dear Sidheshwar
>>There are two types of excited states: vertical and adiabatic excited states. Vertical excitation is more important for UV-Vis spectra in which lifetime of excitation is too short in such a way that Franck-Condon effect is important. Thus in vertical excitation, the geometry of the molecule is not changed considerably. TDDFT is appropriate method for vertical excitation. However, in adiabatic excitation, lifetime of excited state is longer and geometry of the molecule is changed. For this excitation you should use configuration interaction (CI) or CASSCF methods.
>>
>>
>>On Tue Apr 1 '14 12:06pm, Siddheshwar Chopra wrote
>>--------------------------------------------------
>>>Dear Sir,
>>>I went through the forum regarding this but couldn't get my answer exactly. I am running TDDFT calculations on the ground state configuration. Steps were:
>>>1) Copied the lowest energy config. coordinates to TDDFT i/p file.
>>>2) Using RUNTYP=optimize and CITYP=TDDFT in \$CONTRL, to optimize for the first excited state only.
>>>Am I performing it in the correct manner? Will I get "Equilibrium geometry located" and the coordinates of Lowest first excited state? Will the same run give me UV-Vis. spectra too?
>>>I wish to know how similar is this optimization to the Hessian run, say in terms of the memory and CPU time requirement? How to speed up these calcs?
>>>And secondly, I read on the forum about CIS.. But do I need that here?

>>>Regards,

Wed Apr 2 '14 1:41am