From owner-chemistry@ccl.net Fri Apr 22 11:26:01 2022 From: "Nikunj Kumar nkk749 a gmail.com" To: CCL Subject: CCL: How to eliminate the spin contamination in the ORCA Message-Id: <-54642-220422053810-7040-caNEr/Scf9ADxwhjyM7LmA|server.ccl.net> X-Original-From: "Nikunj Kumar" Date: Fri, 22 Apr 2022 05:38:05 -0400 Sent to CCL by: "Nikunj Kumar" [nkk749..gmail.com] Dear all, I am working on an interesting catalyst system which has non-innocent ligand. I have optimized the triplet structure using the ORCA 4.2.1 employing the TPPSh functional. I was expecting the S^2 value to be around 2.0 but it came out to be 2.71. Can we consider this as a triplet? Is this value reliable, if not how can we get rid of the spin contamination. Thanks Nikunj Kumar From owner-chemistry@ccl.net Fri Apr 22 13:29:00 2022 From: "Marcel Swart marcel.swart.:.gmail.com" To: CCL Subject: CCL: How to eliminate the spin contamination in the ORCA Message-Id: <-54643-220422132558-23511-JgQ3z0Koi+41aPndmjrwpw]_[server.ccl.net> X-Original-From: Marcel Swart Content-Transfer-Encoding: 7bit Content-Type: multipart/alternative; boundary=Apple-Mail-A51D5EEC-D1A5-4A58-A359-2A4D317925FE Date: Fri, 22 Apr 2022 19:25:51 +0200 Mime-Version: 1.0 (1.0) Sent to CCL by: Marcel Swart [marcel.swart^gmail.com] --Apple-Mail-A51D5EEC-D1A5-4A58-A359-2A4D317925FE Content-Type: text/plain; charset=utf-8 Content-Transfer-Encoding: quoted-printable Partially this can be the result of the basis set used, see e.g. J. Phys. Chem. A 2010, 114, 26, 7191=E2=80=937197 https://pubs.acs.org/doi/10.1021/jp102712z In the same paper you can find the formulas for spin decontamination, which c= an be done at the energy, gradient, Hessian level with QUILD. https://www.marcelswart.eu/quild/index.html MS Marcel Swart FRSC FYAE MAE, Prof. Dr. ICREA Research Professor at University of Girona Director of Institut de Qu=C3=ADmica Computacional i Cat=C3=A0lisi Univ. Girona, Campus Montilivi (Ci=C3=A8ncies) c/ Maria Aur=C3=A8lia Capmany i Farn=C3=A9s, 69 17003 Girona, Spain www.marcelswart.eu marcel.swart*gmail.com vCard addressbook://www.marcelswart.eu/MSwart.vcf > On 22 Apr 2022, at 18:55, Nikunj Kumar nkk749 a gmail.com wrote: >=20 > =EF=BB=BF > Sent to CCL by: "Nikunj Kumar" [nkk749..gmail.com] > Dear all, > I am working on an interesting catalyst system which has non-innocent liga= nd. I=20 > have optimized the triplet structure using the ORCA 4.2.1 employing the TP= PSh=20 > functional. I was expecting the S^2 value to be around 2.0 but it came out= to=20 > be 2.71. Can we consider this as a triplet? Is this value reliable, if not= how=20 > can we get rid of the spin contamination. >=20 > Thanks > Nikunj Kumar >=20 >=20 >=20 > -=3D This is automatically added to each message by the mailing script =3D= ->=20>=20>=20 > Subscribe/Unsubscribe:=20>=20>=20 > Job: http://www.ccl.net/jobs=20>=20>=20>=20>=20 >=20 --Apple-Mail-A51D5EEC-D1A5-4A58-A359-2A4D317925FE Content-Type: text/html; charset=utf-8 Content-Transfer-Encoding: quoted-printable Partially this can be the result of the bas= is set used, see e.g.

J. Phys. Chem. A = 2010, 114, 267191=E2=80=937197
https://pubs.acs.org/doi/= 10.1021/jp102712z

In the same paper you can find the formulas for= spin decontamination, which can be done at the energy, gradient, Hessian le= vel with QUILD.

MS

= Marcel Swart FRSC FYAE MAE, Prof. Dr.
ICREA Research Professor at University of Gir= ona
Director of Institut de Qu=C3=AD= mica Computacional i Cat=C3=A0lisi

Univ. Girona, C= ampus Montilivi (Ci=C3=A8ncies)
c/ M= aria Aur=C3=A8lia Capmany i Farn=C3=A9s, 69


vCard
<= div class=3D"" style=3D"font-family: Arial; font-size: 12px; -webkit-text-si= ze-adjust: auto; font-variant-ligatures: normal; font-variant-position: norm= al; font-variant-numeric: normal; font-variant-alternates: normal; font-vari= ant-east-asian: normal; word-wrap: break-word; -webkit-nbsp-mode: space; lin= e-break: after-white-space;"><= br class=3D"">

On 22 Apr 2022, at 18:55, Nikunj Kumar nkk749 a gmail.com <owne= r-chemistry*ccl.net> wrote:

=EF=BB=BF
Sent to CCL by: "Nik= unj  Kumar" [nkk749..gmail.com]
Dear all,
I am working on an interesting catalyst system which has non-innocent lig= and. I
have optimized the triplet structure using the ORCA 4= .2.1 employing the TPPSh
functional. I was expecting the S^= 2 value to be around 2.0 but it came out to
be 2.71. Can we= consider this as a triplet? Is this value reliable, if not how
<= span>can we get rid of the spin contamination.

Thanks
Nikunj Kumar



-=3D This is automatically added to each mess= age by the mailing script =3D-
To recover the email address o= f the author of the message, please change
the strange chara= cters on the top line to the * sign. You can also
look up th= e X-Original-From: line in the mail header.


 &nb= sp;   http://www.ccl.net/cgi-bin/ccl/send_ccl_message<= br>
E-mail to administrators: CHEMISTRY-REQUEST*ccl.ne= t or use
     http://www.ccl.net/c= gi-bin/ccl/send_ccl_message

Subscribe/Unsub= scribe:
     http://www.ccl.net/c= hemistry/sub_unsub.shtml

Before posting, ch= eck wait time at: http://www.ccl.net

Job: h= ttp://www.ccl.net/jobs
Conferences: http://server.ccl.net/c= hemistry/announcements/conferences/

Search M= essages: http://www.ccl.net/chemistry/searchccl/index.shtml
=

     http://www.ccl.net/spammers.txt

RTFI: http://www.ccl.net/chemistry/aboutccl/= instructions/


= --Apple-Mail-A51D5EEC-D1A5-4A58-A359-2A4D317925FE-- From owner-chemistry@ccl.net Fri Apr 22 14:07:00 2022 From: "Kathrin Helen Hopmann kathrin.hopmann[#]uit.no" To: CCL Subject: CCL: How to eliminate the spin contamination in the ORCA Message-Id: <-54644-220422140231-6460-Cm8e+CG48w7nrJaEuJP3jg(0)server.ccl.net> X-Original-From: Kathrin Helen Hopmann Content-Language: en-US Content-Transfer-Encoding: 8bit Content-Type: text/plain; charset="us-ascii" Date: Fri, 22 Apr 2022 18:02:23 +0000 MIME-Version: 1.0 Sent to CCL by: Kathrin Helen Hopmann [kathrin.hopmann]_[uit.no] Dear Nikunj Kumar, Have you checked where the unpaired spin is located? 2.71 sounds as if you may have a "broken symmetry triplet", which would imply excess beta spin in one place of the molecule (the ligand?) and excess alpha spin in another (for example a metal). If you have 1 unpaired beta spin on the ligand and 3 unpaired alpha spins on the metal you would get a broken symmetry triplet with a higher than expected S^2 value (the actual computed spin populations are probably smaller, for example 0.6 on ligand and 2.6 on metal). In my experience, functionals with less (or no) HF exchange are more likely to give pure spin states (= less spin contamination). So if you are certain that the electronic structure should be a pure triplet, you could go for another (pure) functional and see what solution you get. If you obtain the pure triplet (2 unpaired alpha spins), you could try to also find the broken symmetry triplet with the other functional (this can be done by using the broken symmetry TPSSh electronic structure or the quintet electronic structure as input for the optimization). In that way you could get an idea if the pure or the broken symmetry triplet is lower in energy. It may be noted though that the energy of a broken symmetry state is approximate (you would have to do a spin projection of the energy to get a more correct value). You should probably also check the energy of the singlet as well as the quintet, just to be sure that none of these is lower than the triplet. And of course, as always, different functionals may give you different answers as to what is the preferred electronic state. best regards, Kathrin -----Original Message----- > From: owner-chemistry+kathrin.hopmann==uit.no . ccl.net On Behalf Of Nikunj Kumar nkk749 a gmail.com Sent: fredag 22. april 2022 11:38 To: Kathrin Helen Hopmann Subject: CCL: How to eliminate the spin contamination in the ORCA Sent to CCL by: "Nikunj Kumar" [nkk749..gmail.com] Dear all, I am working on an interesting catalyst system which has non-innocent ligand. I have optimized the triplet structure using the ORCA 4.2.1 employing the TPPSh functional. I was expecting the S^2 value to be around 2.0 but it came out to be 2.71. Can we consider this as a triplet? Is this value reliable, if not how can we get rid of the spin contamination. Thanks Nikunj Kumarhttp://www.ccl.net/cgi-bin/ccl/send_ccl_messagehttp://www.ccl.net/chemistry/sub_unsub.shtmlhttp://www.ccl.net/spammers.txt From owner-chemistry@ccl.net Fri Apr 22 15:55:01 2022 From: "=?utf-8?B?TWFyaXVzeiBSYWRvxYQ=?= mariusz.radon^_^uj.edu.pl" To: CCL Subject: CCL: How to eliminate the spin contamination in the ORCA Message-Id: <-54645-220422155357-23120-Fd0iIj5/vXirgw/6E0RDpQ _ server.ccl.net> X-Original-From: =?utf-8?B?TWFyaXVzeiBSYWRvxYQ=?= Content-ID: <71EF152EE5C712438CD8BC1AEBA91D93 _ eurprd04.prod.outlook.com> Content-Language: en-US Content-Transfer-Encoding: 8bit Content-Type: text/plain; charset="utf-8" Date: Fri, 22 Apr 2022 19:53:48 +0000 MIME-Version: 1.0 Sent to CCL by: =?utf-8?B?TWFyaXVzeiBSYWRvxYQ=?= [mariusz.radon:-:uj.edu.pl] > On 22 Apr 2022, at 11:38, Nikunj Kumar nkk749 a gmail.com wrote: > > > Sent to CCL by: "Nikunj Kumar" [nkk749..gmail.com] > Dear all, > I am working on an interesting catalyst system which has non-innocent ligand. I > have optimized the triplet structure using the ORCA 4.2.1 employing the TPPSh > functional. I was expecting the S^2 value to be around 2.0 but it came out to > be 2.71. Can we consider this as a triplet? Is this value reliable, if not how > can we get rid of the spin contamination. > > Thanks > Nikunj Kumar Dear Nikunj: Assuming that the contamination of your triplet originates from one single quintet state, the value of 2.71 would imply ~18% of the contaminant: = (1-x)*2 + x*6 => x=0.18 To know how important is this effect for the energy, you need to also compute the energy of contaminating quintet state. As Marcel Swart already told you, there are relatively simple formulas to correct the energy for spin-contamination. You can find them in the Marcel’s paper or in many other papers, such as mine: https://dx.doi.org/10.1021/ic503109a. One important remark: you need to compute the the proper quintet state (i.e., the one which contaminates your triplet state). To know how to obtain the proper quintet state, it can be useful to take a look at natural orbitals of your spin-contaminated triplet state. You will most likely see the following natural orbitals: - two orbitals with occupation numbers exactly 1.0 (phi_1, phi_2) - a pair of orbitals with fractional occupation numbers 2-v and v (phi_3, phi_4). For the proper quintet state, you need to put unpaired electrons into these four orbitals (phi_1, phi_2, phi_3, phi_4). You take a look into the above paper from Inorganic Chemistry (https://dx.doi.org/10.1021/ic503109a) for some examples. I don’t know if it is possible to get natural orbitals directly from Orca, but you can always compute them with my program Natorbs. https://tungsten.ch.uj.edu.pl/~mradon/natorbs/ or https://pypi.org/project/natorbs/ Natorbs imports UHF-type molecular orbitals via either Molden format or cclib interface. Best wishes, Mariusz Radon -- Mariusz Radon, Ph.D., D.Sc. Associate Professor Faculty of Chemistry, Jagiellonian University Address: Gronostajowa 2, 30-387 Krakow, Poland Room C1-06, Phone: 48-12-686-24-89 E-mail: mradon : chemia.uj.edu.pl (mariusz.radon : uj.edu.pl) Web: https://tungsten.ch.uj.edu.pl/~mradon ORCID: https://orcid.org/0000-0002-1901-8521 From owner-chemistry@ccl.net Fri Apr 22 18:24:00 2022 From: "Christoph Riplinger riplinger__faccts.de" To: CCL Subject: CCL: How to eliminate the spin contamination in the ORCA Message-Id: <-54646-220422171710-1665-2py5r2DUxG7iFy65Kqjj0g||server.ccl.net> X-Original-From: Christoph Riplinger Content-Language: en-US Content-Type: multipart/alternative; boundary="------------33iVECioqCfex3kQgs6o4CDh" Date: Fri, 22 Apr 2022 23:17:01 +0200 MIME-Version: 1.0 Sent to CCL by: Christoph Riplinger [riplinger#%#faccts.de] This is a multi-part message in MIME format. --------------33iVECioqCfex3kQgs6o4CDh Content-Type: text/plain; charset=UTF-8; format=flowed Content-Transfer-Encoding: base64 RGVhciBOaWt1bmosDQoNCmluIGFkZGl0aW9uIHRvIHRoZSBwcmV2aW91cyBjb21tZW50czog SWYgeW91IGhhdmUgYSBicm9rZW4gc3ltbWV0cnkgc3RhdGUsIHlvdSBjYW4gZ2V0IGFuZCBh bmFseXplIHRoZSAiY29ycmVzcG9uZGluZyAgb3JiaXRhbHMiIHdpdGggT1JDQSBieSBhZGRp bmcgdGhlIGZvbGxvd2luZyBrZXl3b3JkIHRvIHlvdXIgaW5wdXQ6ICEgVUNPDQoNCkZvciBh IGRlc2NyaXB0aW9uIG9uIGhvdyB0byBsb29rIGF0IHRoZSBVQ09zIHNlZSByZWZlcmVuY2U6 IE5lZXNlLCBGLiwgRGVmaW5pdGlvbiBvZmNvcnJlc3BvbmRpbmcgIG9yYml0YWxzIGFuZCB0 aGUgZGlyYWRpY2FsIGNoYXJhY3RlciBpbmJyb2tlbiAgc3ltbWV0cnkgIERGVCBjYWxjdWxh dGlvbnMgb24gc3BpbiBjb3VwbGVkIHN5c3RlbXMuIEpvdXJuYWwgb2YgUGh5c2ljcyBhbmQg Q2hlbWlzdHJ5IG9mIFNvbGlkcyAyMDA0LCA2NSAoNCksIDc4MS03ODUuDQoNCkJlc3Qgd2lz aGVzLA0KQ2hyaXN0b3BoDQoNClBTOiBBbHNvLCBJIGVuY291cmFnZSB5b3UgdG8gdXBncmFk ZSB0byBPUkNBIDUsIHdoaWNoIHNob3dzIHNpZ25pZmljYW50bHkgaW1wcm92ZWQgcm9idXN0 bmVzcyBhbmQgc3BlZWR1cCwgaW4gcGFydGljdWxhciBmb3IgaHlicmlkIGZ1bmN0aW9uYWxz Lg0KDQpPbiAyMi4wNC4yMiAyMDowMiwgS2F0aHJpbiBIZWxlbiBIb3BtYW5uIGthdGhyaW4u aG9wbWFublsjXXVpdC5ubyB3cm90ZToNCj4gU2VudCB0byBDQ0wgYnk6IEthdGhyaW4gSGVs ZW4gSG9wbWFubiBba2F0aHJpbi5ob3BtYW5uXV9bdWl0Lm5vXQ0KPg0KPiBEZWFyIE5pa3Vu aiBLdW1hciwNCj4NCj4gSGF2ZSB5b3UgY2hlY2tlZCB3aGVyZSB0aGUgdW5wYWlyZWQgc3Bp biBpcyBsb2NhdGVkPw0KPg0KPiAyLjcxIHNvdW5kcyBhcyBpZiB5b3UgbWF5IGhhdmUgYSAi YnJva2VuIHN5bW1ldHJ5IHRyaXBsZXQiLCB3aGljaCB3b3VsZCBpbXBseSBleGNlc3MgYmV0 YSBzcGluIGluIG9uZSBwbGFjZSBvZiB0aGUgbW9sZWN1bGUgKHRoZSBsaWdhbmQ/KSBhbmQg ZXhjZXNzIGFscGhhIHNwaW4gaW4gYW5vdGhlciAoZm9yIGV4YW1wbGUgYSBtZXRhbCkuIElm IHlvdSBoYXZlIDEgdW5wYWlyZWQgYmV0YSBzcGluIG9uIHRoZSBsaWdhbmQgYW5kIDMgdW5w YWlyZWQgYWxwaGEgc3BpbnMgb24gdGhlIG1ldGFsIHlvdSB3b3VsZCBnZXQgYSBicm9rZW4g c3ltbWV0cnkgdHJpcGxldCB3aXRoIGEgaGlnaGVyIHRoYW4gZXhwZWN0ZWQgU14yIHZhbHVl ICh0aGUgYWN0dWFsIGNvbXB1dGVkIHNwaW4gcG9wdWxhdGlvbnMgYXJlIHByb2JhYmx5IHNt YWxsZXIsIGZvciBleGFtcGxlIDAuNiBvbiBsaWdhbmQgYW5kIDIuNiBvbiBtZXRhbCkuDQo+ DQo+IEluIG15IGV4cGVyaWVuY2UsIGZ1bmN0aW9uYWxzIHdpdGggbGVzcyAob3Igbm8pIEhG IGV4Y2hhbmdlIGFyZSBtb3JlIGxpa2VseSB0byBnaXZlIHB1cmUgc3BpbiBzdGF0ZXMgKD0g bGVzcyBzcGluIGNvbnRhbWluYXRpb24pLiBTbyBpZiB5b3UgYXJlIGNlcnRhaW4gdGhhdCB0 aGUgZWxlY3Ryb25pYyBzdHJ1Y3R1cmUgc2hvdWxkIGJlIGEgcHVyZSB0cmlwbGV0LCB5b3Ug Y291bGQgZ28gZm9yIGFub3RoZXIgKHB1cmUpIGZ1bmN0aW9uYWwgYW5kIHNlZSB3aGF0IHNv bHV0aW9uIHlvdSBnZXQuIElmIHlvdSBvYnRhaW4gdGhlIHB1cmUgdHJpcGxldCAoMiB1bnBh aXJlZCBhbHBoYSBzcGlucyksIHlvdSBjb3VsZCB0cnkgdG8gYWxzbyBmaW5kIHRoZSBicm9r ZW4gc3ltbWV0cnkgdHJpcGxldCB3aXRoIHRoZSBvdGhlciBmdW5jdGlvbmFsICh0aGlzIGNh biBiZSBkb25lIGJ5IHVzaW5nIHRoZSBicm9rZW4gc3ltbWV0cnkgVFBTU2ggZWxlY3Ryb25p YyBzdHJ1Y3R1cmUgb3IgdGhlIHF1aW50ZXQgZWxlY3Ryb25pYyBzdHJ1Y3R1cmUgYXMgaW5w dXQgZm9yIHRoZSBvcHRpbWl6YXRpb24pLiBJbiB0aGF0IHdheSB5b3UgY291bGQgZ2V0IGFu IGlkZWEgaWYgdGhlIHB1cmUgb3IgdGhlIGJyb2tlbiBzeW1tZXRyeSB0cmlwbGV0IGlzIGxv d2VyIGluIGVuZXJneS4gSXQgbWF5IGJlIG5vdGVkIHRob3VnaCB0aGF0IHRoZSBlbmVyZ3kg b2YgYSBicm9rZW4gc3ltbWV0cnkgc3RhdGUgaXMgYXBwcm94aW1hdGUgKHlvdSB3b3VsZCBo YXZlIHRvIGRvIGEgc3BpbiBwcm9qZWN0aW9uIG9mIHRoZSBlbmVyZ3kgdG8gZ2V0IGEgbW9y ZSBjb3JyZWN0IHZhbHVlKS4gWW91IHNob3VsZCBwcm9iYWJseSBhbHNvIGNoZWNrIHRoZSBl bmVyZ3kgb2YgdGhlIHNpbmdsZXQgYXMgd2VsbCBhcyB0aGUgcXVpbnRldCwganVzdCB0byBi ZSBzdXJlIHRoYXQgbm9uZSBvZiB0aGVzZSBpcyBsb3dlciB0aGFuIHRoZSB0cmlwbGV0LiBB bmQgb2YgY291cnNlLCBhcyBhbHdhIQ0KPiAgIHlzLCBkaWZmZXJlbnQgZnVuY3Rpb25hbHMg bWF5IGdpdmUgeW91IGRpZmZlcmVudCBhbnN3ZXJzIGFzIHRvIHdoYXQgaXMgdGhlIHByZWZl cnJlZCBlbGVjdHJvbmljIHN0YXRlLg0KPg0KPiBiZXN0IHJlZ2FyZHMsDQo+IEthdGhyaW4N Cj4NCj4NCj4gLS0tLS1PcmlnaW5hbCBNZXNzYWdlLS0tLS0NCj4NCj4gU2VudCB0byBDQ0wg Ynk6ICJOaWt1bmogIEt1bWFyIiBbbmtrNzQ5Li5nbWFpbC5jb21dIERlYXIgYWxsLCBJIGFt IHdvcmtpbmcgb24gYW4gaW50ZXJlc3RpbmcgY2F0YWx5c3Qgc3lzdGVtIHdoaWNoIGhhcyBu b24taW5ub2NlbnQgbGlnYW5kLiBJIGhhdmUgb3B0aW1pemVkIHRoZSB0cmlwbGV0IHN0cnVj dHVyZSB1c2luZyB0aGUgT1JDQSA0LjIuMSBlbXBsb3lpbmcgdGhlIFRQUFNoIGZ1bmN0aW9u YWwuIEkgd2FzIGV4cGVjdGluZyB0aGUgU14yIHZhbHVlIHRvIGJlIGFyb3VuZCAyLjAgYnV0 IGl0IGNhbWUgb3V0IHRvIGJlIDIuNzEuIENhbiB3ZSBjb25zaWRlciB0aGlzIGFzIGEgdHJp cGxldD8gSXMgdGhpcyB2YWx1ZSByZWxpYWJsZSwgaWYgbm90IGhvdyBjYW4gd2UgZ2V0IHJp ZCBvZiB0aGUgc3BpbiBjb250YW1pbmF0aW9uLg0KPg0KLS0gDQpGQWNjVHMgR21iSA0KUm9s YW5kc3RyYXNzZSA2NywgNTA2NzcgS8O2bG4NCkFtdHNnZXJpY2h0IEvDtmxuIEhSQiA4ODQw Ng0KR2VzY2jDpGZ0c2bDvGhyZXI6IERyLiBDaHJpc3RvcGggUmlwbGluZ2VyDQpodHRwczov L3d3dy5mYWNjdHMuZGUgIHxodHRwczovL3R3aXR0ZXIuY29tL2ZhY2N0c19vcmNhDQoNCg== --------------33iVECioqCfex3kQgs6o4CDh Content-Type: text/html; charset=UTF-8 Content-Transfer-Encoding: 8bit
Dear Nikunj,

in addition to the previous comments: If you have a broken symmetry state, you can get and analyze the "corresponding orbitals" with ORCA by adding the following keyword to your input: ! UCO

For a description on how to look at the UCOs see reference: Neese, F., Definition of corresponding orbitals and the diradical character in broken symmetry DFT calculations on spin coupled systems. Journal of Physics and Chemistry of Solids 2004, 65 (4), 781-785.

Best wishes,
Christoph

PS: Also, I encourage you to upgrade to ORCA 5, which shows significantly improved robustness and speedup, in particular for hybrid functionals.

On 22.04.22 20:02, Kathrin Helen Hopmann kathrin.hopmann[#]uit.no wrote:
Sent to CCL by: Kathrin Helen Hopmann [kathrin.hopmann]_[uit.no]

Dear Nikunj Kumar, 

Have you checked where the unpaired spin is located? 

2.71 sounds as if you may have a "broken symmetry triplet", which would imply excess beta spin in one place of the molecule (the ligand?) and excess alpha spin in another (for example a metal). If you have 1 unpaired beta spin on the ligand and 3 unpaired alpha spins on the metal you would get a broken symmetry triplet with a higher than expected S^2 value (the actual computed spin populations are probably smaller, for example 0.6 on ligand and 2.6 on metal). 

In my experience, functionals with less (or no) HF exchange are more likely to give pure spin states (= less spin contamination). So if you are certain that the electronic structure should be a pure triplet, you could go for another (pure) functional and see what solution you get. If you obtain the pure triplet (2 unpaired alpha spins), you could try to also find the broken symmetry triplet with the other functional (this can be done by using the broken symmetry TPSSh electronic structure or the quintet electronic structure as input for the optimization). In that way you could get an idea if the pure or the broken symmetry triplet is lower in energy. It may be noted though that the energy of a broken symmetry state is approximate (you would have to do a spin projection of the energy to get a more correct value). You should probably also check the energy of the singlet as well as the quintet, just to be sure that none of these is lower than the triplet. And of course, as alwa!
 ys, different functionals may give you different answers as to what is the preferred electronic state. 

best regards,
Kathrin 


-----Original Message-----

Sent to CCL by: "Nikunj  Kumar" [nkk749..gmail.com] Dear all, I am working on an interesting catalyst system which has non-innocent ligand. I have optimized the triplet structure using the ORCA 4.2.1 employing the TPPSh functional. I was expecting the S^2 value to be around 2.0 but it came out to be 2.71. Can we consider this as a triplet? Is this value reliable, if not how can we get rid of the spin contamination.

-- 
FAccTs GmbH
Rolandstrasse 67, 50677 Köln
Amtsgericht Köln HRB 88406
Geschäftsführer: Dr. Christoph Riplinger
https://www.faccts.de | https://twitter.com/faccts_orca
--------------33iVECioqCfex3kQgs6o4CDh--