From owner-chemistry@ccl.net Wed Aug 17 08:09:00 2016 From: "Norrby, Per-Ola Per-Ola.Norrby^^astrazeneca.com" To: CCL Subject: CCL:G: triplet TDDFT calculation Message-Id: <-52343-160817024824-15270-0bJ1kYxP+TpRkE1Pru4Mow%a%server.ccl.net> X-Original-From: "Norrby, Per-Ola" Content-Language: en-US Content-Type: multipart/alternative; boundary="_000_HE1PR04MB209269C49328872FB7ED378BCA140HE1PR04MB2092eurp_" Date: Wed, 17 Aug 2016 06:48:14 +0000 MIME-Version: 1.0 Sent to CCL by: "Norrby, Per-Ola" [Per-Ola.Norrby---astrazeneca.com] --_000_HE1PR04MB209269C49328872FB7ED378BCA140HE1PR04MB2092eurp_ Content-Type: text/plain; charset="us-ascii" Content-Transfer-Encoding: quoted-printable Dear Michael, You get variable because you use the unrestricted approximation. For= a triplet, the exact value of is 2; any value higher than that indi= cates contamination from higher states. There is no easy or cheap way out o= f this. You could use restricted open shell calculations, RO instead of the= default U (not sure about Gaussian, I've done those in Jaguar), but in my = experience, energies are less reliable even though you're guaranteed pure s= pin states. For ground state calculations, you can project out most of the = higher spins in the unrestricted solution, but again, I'm not sure the meth= od is available for TDDFT. For accurate identification of "real" excited st= ates, I've used CASSCF/CASPT2, but even there, the energies from unrestrict= ed TDDFT can be as good or better. I did find the CASPT2 calculations good = for evaluating which of the DFT excited states were "real", and which were = artifacts in the unrestricted approximation. Not published, just experience= from a long project we never managed to finalize. Best regards, Per-Ola > From: owner-chemistry+per-ola.norrby=3D=3Dastrazeneca.com\a/ccl.net [mailto:o= wner-chemistry+per-ola.norrby=3D=3Dastrazeneca.com\a/ccl.net] On Behalf Of Mi= chael Morgan michaelmorgan937%gmail.com Sent: den 17 augusti 2016 05:42 To: Norrby, Per-Ola Subject: CCL:G: triplet TDDFT calculation Dear All, Standard TDDFT calculations calculate singlet->singlet excitations. In Gaus= sian, I can see that for each excitation, =3D0.000. While I calculate molecules like O2 which the ground state is triplet, I ne= ed calculate triplet->triplet excitations by assigning TD=3Dtriplet . In output, I found all kinds of values (from 2 to 4). I am confused:= for those are not close to 3, are they still spin-allowed therefore experimentally observable? Should I include them if I want to predict an ab= sorption spectrum for the molecule? Thank you very much! Michael ________________________________ Confidentiality Notice: This message is private and may contain confidentia= l and proprietary information. If you have received this message in error, = please notify us and remove it from your system and note that you must not = copy, distribute or take any action in reliance on it. Any unauthorized use= or disclosure of the contents of this message is not permitted and may be = unlawful. --_000_HE1PR04MB209269C49328872FB7ED378BCA140HE1PR04MB2092eurp_ Content-Type: text/html; charset="us-ascii" Content-Transfer-Encoding: quoted-printable

Dear Michael,

 

You get variable <S= **2> because you use the unrestricted approximation. For a triplet, the = exact value of <S**2> is 2; any value higher than that indicates cont= amination from higher states. There is no easy or cheap way out of this. You could use restricted open shell calculations, R= O instead of the default U (not sure about Gaussian, I’ve done those = in Jaguar), but in my experience, energies are less reliable even though yo= u’re guaranteed pure spin states. For ground state calculations, you can project out most of the higher spins in= the unrestricted solution, but again, I’m not sure the method is ava= ilable for TDDFT. For accurate identification of “real” excited= states, I’ve used CASSCF/CASPT2, but even there, the energies from unrestricted TDDFT can be as good or better. I did find = the CASPT2 calculations good for evaluating which of the DFT excited states= were “real”, and which were artifacts in the unrestricted appr= oximation. Not published, just experience from a long project we never managed to finalize.

 

Best regards,

 

Per-Ola

 

From: owner-chemistry+per-ola.norrby=3D= =3Dastrazeneca.com\a/ccl.net [mailto:owner-chemistry+per-ola.norrby=3D=3D= astrazeneca.com\a/ccl.net] On Behalf Of Michael Morgan michaelmorgan937%gmail.com
Sent: den 17 augusti 2016 05:42
To: Norrby, Per-Ola <Per-Ola.Norrby\a/astrazeneca.com>
Subject: CCL:G: triplet TDDFT calculation

 

Dear All,

 

Standard TDDFT calculations calculate singlet->si= nglet excitations. In Gaussian, I can see that for each excitation, <S**= 2>=3D0.000.

 

While I calculate molecules like O2 which the ground= state is triplet, I need calculate triplet->triplet excitations by assi= gning TD=3Dtriplet .
In output, I found all kinds of <S**2> values (from 2 to 4). I am con= fused: for those are not close to 3, are they still spin-allowed therefore
experimentally observable? Should I include them if I want to predict an ab= sorption spectrum for the molecule?

 

Thank you very much!

 

Michael


Confidentiality Notice: This message is private and may contain confidential and proprieta= ry information. If you have received this message in error, please notify u= s and remove it from your system and note that you must not copy, distribut= e or take any action in reliance on it. Any unauthorized use or disclosure of the contents of this message = is not permitted and may be unlawful.

--_000_HE1PR04MB209269C49328872FB7ED378BCA140HE1PR04MB2092eurp_--