From owner-chemistry@ccl.net Mon Jul 29 10:53:01 2013 From: "Belkhiri, lotfi lotfi.belkhiri{:}umc.edu.dz" To: CCL Subject: CCL:G: Gaussian 09: how to fragment molecules Message-Id: <-49031-130729104755-31999-lbfEz6gawauIZ6g78jPxAA^^server.ccl.net> X-Original-From: "Belkhiri, lotfi" Content-Transfer-Encoding: 8bit Content-Type: text/plain; charset=ISO-8859-1 Date: Mon, 29 Jul 2013 15:47:46 +0100 MIME-Version: 1.0 Sent to CCL by: "Belkhiri, lotfi" [lotfi.belkhiri * umc.edu.dz] Hi, You have to use ADF program (www.scm.com) with EDA (Energy Decomposition Analysis) which allows fragment computations with specific charges. Best greetings Lotfi 2013/7/24, Pratima Satish ps544- -cornell.edu : > > Sent to CCL by: "Pratima Satish" [ps544^cornell.edu] > Hi, > > I am trying to define a transition metal complex in Gaussian 09 with > localised charges - for instance, a hexaaquairon complex [Fe(H2O)6]2+, where > > the iron centre specifically has a +2 charge and all the water ligands have > 0 > charge. However, I am not fully sure as to how to use the Guess=Fragment > feature in Gaussian, which keeps giving me the following error message: > > "This type of calculation cannot be archived." > > I am including my input file: > ------------------------- > > #n UB3LYP/Gen Opt Guess=(Only,Fragment=7) > > Hexaaquairon(2) bond length 2.0 A high spin (multiplicity 5) fragmented into > > Fe2+ and 6 waters > > 2 5, 2 5, 0 1, 0 1, 0 1, 0 1, 0 1, 0 1 > Fe(Fragment=1) 0.000000 0.000007 0.000000 > O(Fragment=2) 0.020930 1.999897 -0.000928 > O(Fragment=3) -1.766147 0.018925 0.938278 > O(Fragment=4) -0.938235 0.009006 -1.766248 > O(Fragment=5) 0.938236 -0.008993 1.766248 > O(Fragment=6) 1.766148 -0.018912 -0.938278 > O(Fragment=7) -0.020929 -1.999884 0.000928 > H(Fragment=2) -0.508951 2.591673 -0.557495 > H(Fragment=2) 0.562411 2.580882 0.555837 > H(Fragment=3) -2.253533 0.798638 1.248259 > H(Fragment=3) -2.270468 -0.750283 1.247407 > H(Fragment=4) -0.542006 0.006388 -2.651473 > H(Fragment=4) -1.898222 0.017369 -1.902694 > H(Fragment=5) 1.898228 -0.016775 1.902692 > H(Fragment=5) 0.542004 -0.007343 2.651475 > H(Fragment=6) 2.270516 0.750215 -1.247533 > H(Fragment=6) 2.253469 -0.798707 -1.248151 > H(Fragment=7) -0.562448 -2.580867 -0.555812 > H(Fragment=7) 0.508968 -2.591675 0.557476 > > Fe 0 > aug-cc-pVTZ > **** > H O 0 > aug-cc-pVDZ > **** > > ------------------------- > > Can someone please let me know how to get around this? Thanks so much! > > - Pratima Satish. > Email: ps544[-]cornell.edu> > > -- Pr. Lotfi Belkhiri Directeur de Recherches - Groupe de Chimie Théorique et Modélisation Unité de Recherche CHEMS - Université Mentouri de Constantine jct9_cne2010 ~ yahoo.fr Fax : +213 (0)31 68 15 45 Homepage : http://www.umc.edu.dz/vf/Labo/FacScienceExact/site_lacmom/crbst_21.html From owner-chemistry@ccl.net Mon Jul 29 12:02:01 2013 From: "Oliver Koch oliver.koch _ tu-dortmund.de" To: CCL Subject: CCL: Research Telegrams at 9th German Conference on Chemoinformatics Message-Id: <-49032-130729105100-32405-+U5AKrJKf0UwziNjq9tGUA**server.ccl.net> X-Original-From: "Oliver Koch" Date: Mon, 29 Jul 2013 10:50:56 -0400 Sent to CCL by: "Oliver Koch" [oliver.koch*|*tu-dortmund.de] The divison Chemistry-Information-Computer (CIC) of the German Chemical Society wants to give PhD student the opportunity to present their work in form of short talks as Research Telegrams at the 9th German Conference on Chemoinformatics The conference will take place November, 10. - 12., 2013 FULDA, Germany. The candidates should be PhD students with excellent results at the end of their thesis work. Submissions of Research Telegrams should be passed by the supervisor directly via Email to Dr. Oliver Koch (oliver.koch ~ tu-dortmund.de). Deadline: September 1, 2013 There will be sessions on: Cheminformatics and Drug Discovery Molecular Modelling Computational Materials Science New Developments The aim of the conference is to reflect and highlight the new role of cheminformatics in the modern information world. It will span a wide range of subjects related to the use of computers in chemistry, pharmacy, materials science and biology; > from chem- and bioinformatics to explicit modelling and from industrial applications to fundamental academic research. The scientific programme will include plenary and contributed lectures, posters and software presentations. The language of the conference is English. Registration is now open at: http://www.gdch.de/gcc2013 From owner-chemistry@ccl.net Mon Jul 29 12:37:00 2013 From: "Bijan Mondal mondal.bijan()gmail.com" To: CCL Subject: CCL: ToPMoD package Message-Id: <-49033-130729122049-13015-oP42rbygZO+bgjQ1bIlZ/A{=}server.ccl.net> X-Original-From: Bijan Mondal Content-Type: multipart/alternative; boundary=001a11c361fcbfa92a04e2a8df9d Date: Mon, 29 Jul 2013 21:50:43 +0530 MIME-Version: 1.0 Sent to CCL by: Bijan Mondal [mondal.bijan+*+gmail.com] --001a11c361fcbfa92a04e2a8df9d Content-Type: text/plain; charset=ISO-8859-1 Dear Mates, I want to calculate ELF (electron localization function) and plot them. Can anyone suggest me any program package which will be helpful? The program should be freely available for windows. I heard a package called ToPMoD can calculate ELF. Does anyone knows whether it is available for windows and it is free? Thanks you very much. Any suggestions will be appreciated. Regards, Bijan -- *"There are no such electron-deficient compounds,* * **only theory-deficient chemists." --* R.E. Rundle Bijan Mondal Research Scholar Mob: 7667156593 Department of Chemistry I.I.T Madras Chennai --001a11c361fcbfa92a04e2a8df9d Content-Type: text/html; charset=ISO-8859-1 Content-Transfer-Encoding: quoted-printable
De= ar Mates,

I want to calculate ELF (electron localization function) and plot them. <= br>
Can anyone sugge= st me any program package which will be helpful? The program should be free= ly available for windows.
I heard a package called ToPMoD can calculate = ELF. Does anyone knows whether it is available for windows and it is free?<= br clear=3D"all">

Thanks you ver= y much.
Any suggestions will be appreciated.

Regards= ,
Bijan<= br>
--
<= /span>
"There are no such electron-deficient compounds,<= /span>
=A0= =A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 =A0only theory-deficient chemist= s." -- R.E. Rundle
=A0
Bijan Mondal
Research Scholar
Mob: 7667156593
Department of Chemistry
I.I.T Madras
Ch= ennai

--001a11c361fcbfa92a04e2a8df9d-- From owner-chemistry@ccl.net Mon Jul 29 13:11:00 2013 From: "Philippe Bourassa philippe.bourassa#%#uqtr.ca" To: CCL Subject: CCL:G: Gaussian 09: how to fragment molecules Message-Id: <-49034-130729125007-433-JOiXl1s335UmY/2LZmLC5g#,#server.ccl.net> X-Original-From: Philippe Bourassa Content-Type: multipart/alternative; boundary=047d7bf0f0046e282904e2a948c6 Date: Mon, 29 Jul 2013 12:49:59 -0400 MIME-Version: 1.0 Sent to CCL by: Philippe Bourassa [philippe.bourassa:_:uqtr.ca] --047d7bf0f0046e282904e2a948c6 Content-Type: text/plain; charset=ISO-8859-1 I don't think Gaussian can do calculations with that many fragments. I would try the calculation using 2 fragments. 1 for iron and the other for all 6 water molecules. 2013/7/24 Pratima Satish ps544- -cornell.edu > > Sent to CCL by: "Pratima Satish" [ps544^cornell.edu] > Hi, > > I am trying to define a transition metal complex in Gaussian 09 with > localised charges - for instance, a hexaaquairon complex [Fe(H2O)6]2+, > where > the iron centre specifically has a +2 charge and all the water ligands > have 0 > charge. However, I am not fully sure as to how to use the Guess=Fragment > feature in Gaussian, which keeps giving me the following error message: > > "This type of calculation cannot be archived." > > I am including my input file: > ------------------------- > > #n UB3LYP/Gen Opt Guess=(Only,Fragment=7) > > Hexaaquairon(2) bond length 2.0 A high spin (multiplicity 5) fragmented > into > Fe2+ and 6 waters > > 2 5, 2 5, 0 1, 0 1, 0 1, 0 1, 0 1, 0 1 > Fe(Fragment=1) 0.000000 0.000007 0.000000 > O(Fragment=2) 0.020930 1.999897 -0.000928 > O(Fragment=3) -1.766147 0.018925 0.938278 > O(Fragment=4) -0.938235 0.009006 -1.766248 > O(Fragment=5) 0.938236 -0.008993 1.766248 > O(Fragment=6) 1.766148 -0.018912 -0.938278 > O(Fragment=7) -0.020929 -1.999884 0.000928 > H(Fragment=2) -0.508951 2.591673 -0.557495 > H(Fragment=2) 0.562411 2.580882 0.555837 > H(Fragment=3) -2.253533 0.798638 1.248259 > H(Fragment=3) -2.270468 -0.750283 1.247407 > H(Fragment=4) -0.542006 0.006388 -2.651473 > H(Fragment=4) -1.898222 0.017369 -1.902694 > H(Fragment=5) 1.898228 -0.016775 1.902692 > H(Fragment=5) 0.542004 -0.007343 2.651475 > H(Fragment=6) 2.270516 0.750215 -1.247533 > H(Fragment=6) 2.253469 -0.798707 -1.248151 > H(Fragment=7) -0.562448 -2.580867 -0.555812 > H(Fragment=7) 0.508968 -2.591675 0.557476 > > Fe 0 > aug-cc-pVTZ > **** > H O 0 > aug-cc-pVDZ > **** > > ------------------------- > > Can someone please let me know how to get around this? Thanks so much! > > - Pratima Satish. > Email: ps544[-]cornell.edu> > > --047d7bf0f0046e282904e2a948c6 Content-Type: text/html; charset=ISO-8859-1 Content-Transfer-Encoding: quoted-printable
I don't think Gaussian can do calculations with that m= any fragments. I would try the calculation using 2 fragments. 1 for iron an= d the other for all 6 water molecules.
=

2013/7/24 Pratima Satish ps544- -cornell.edu <owner-chemistry:+:ccl.net>=

Sent to CCL by: "Pratima =A0Satish" [ps544^cornell.edu]
Hi,

I am trying to define a transition metal complex in Gaussian 09 with
localised charges - for instance, a hexaaquairon complex [Fe(H2O)6]2+, wher= e
the iron centre specifically has a +2 charge and all the water ligands have= 0
charge. However, I am not fully sure as to how to use the Guess=3DFragment<= br> feature in Gaussian, which keeps giving me the following error message:

"This type of calculation cannot be archived."

I am including my input file:
-------------------------

#n UB3LYP/Gen Opt Guess=3D(Only,Fragment=3D7)

Hexaaquairon(2) bond length 2.0 A high spin (multiplicity 5) fragmented int= o
Fe2+ and 6 waters

2 5, 2 5, 0 1, 0 1, 0 1, 0 1, 0 1, 0 1
Fe(Fragment=3D1) =A0 =A0 =A00.000000 =A0 =A00.000007 =A0 =A00.000000
O(Fragment=3D2) =A0 =A0 =A0 0.020930 =A0 =A01.999897 =A0 -0.000928
O(Fragment=3D3) =A0 =A0 =A0-1.766147 =A0 =A00.018925 =A0 =A00.938278
O(Fragment=3D4) =A0 =A0 =A0-0.938235 =A0 =A00.009006 =A0 -1.766248
O(Fragment=3D5) =A0 =A0 =A0 0.938236 =A0 -0.008993 =A0 =A01.766248
O(Fragment=3D6) =A0 =A0 =A0 1.766148 =A0 -0.018912 =A0 -0.938278
O(Fragment=3D7) =A0 =A0 =A0-0.020929 =A0 -1.999884 =A0 =A00.000928
H(Fragment=3D2) =A0 =A0 =A0-0.508951 =A0 =A02.591673 =A0 -0.557495
H(Fragment=3D2) =A0 =A0 =A0 0.562411 =A0 =A02.580882 =A0 =A00.555837
H(Fragment=3D3) =A0 =A0 =A0-2.253533 =A0 =A00.798638 =A0 =A01.248259
H(Fragment=3D3) =A0 =A0 =A0-2.270468 =A0 -0.750283 =A0 =A01.247407
H(Fragment=3D4) =A0 =A0 =A0-0.542006 =A0 =A00.006388 =A0 -2.651473
H(Fragment=3D4) =A0 =A0 =A0-1.898222 =A0 =A00.017369 =A0 -1.902694
H(Fragment=3D5) =A0 =A0 =A0 1.898228 =A0 -0.016775 =A0 =A01.902692
H(Fragment=3D5) =A0 =A0 =A0 0.542004 =A0 -0.007343 =A0 =A02.651475
H(Fragment=3D6) =A0 =A0 =A0 2.270516 =A0 =A00.750215 =A0 -1.247533
H(Fragment=3D6) =A0 =A0 =A0 2.253469 =A0 -0.798707 =A0 -1.248151
H(Fragment=3D7) =A0 =A0 =A0-0.562448 =A0 -2.580867 =A0 -0.555812
H(Fragment=3D7) =A0 =A0 =A0 0.508968 =A0 -2.591675 =A0 =A00.557476

Fe 0
aug-cc-pVTZ
****
H O 0
aug-cc-pVDZ
****

-------------------------

Can someone please let me know how to get around this? Thanks so much!

- Pratima Satish.
Email: ps544[-]cornell.edu=



-=3D This is automatically added to each message by the mailing script =3D-=
E-mail to subscribers: CHEMISTRY:+:ccl.n= et or use:
=A0 =A0 =A0 http://www.ccl.net/cgi-bin/ccl/send_ccl_message

E-mail to administrators: CHEM= ISTRY-REQUEST:+:ccl.net or use
=A0 =A0 =A0 http://www.ccl.net/cgi-bin/ccl/send_ccl_message

Subscribe/Unsubscribe:
=A0 =A0 =A0 http://www.ccl.net/chemistry/sub_unsub.shtml

Before posting, check wait time at: http://www.ccl.net

Job: http://www.ccl.n= et/jobs
Conferences: http://server.ccl.net/chemistry/announcements/co= nferences/

Search Messages: http://www.ccl.net/chemistry/searchccl/index.shtml
=A0 =A0 =A0
h= ttp://www.ccl.net/spammers.txt

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



--047d7bf0f0046e282904e2a948c6-- From owner-chemistry@ccl.net Mon Jul 29 17:33:00 2013 From: "Francois Zielinski francois.zielinski|*|manchester.ac.uk" To: CCL Subject: CCL: ToPMoD package Message-Id: <-49035-130729173141-8329-in8UpaX1aKOTbTmdl7tGDw!A!server.ccl.net> X-Original-From: Francois Zielinski Content-Language: en-US Content-Type: multipart/alternative; boundary="_000_3068BCD5C0404D4E9C44F64F4000B64225C77EMBXP02dsmanacuk_" Date: Mon, 29 Jul 2013 21:31:33 +0000 MIME-Version: 1.0 Sent to CCL by: Francois Zielinski [francois.zielinski|a|manchester.ac.uk] --_000_3068BCD5C0404D4E9C44F64F4000B64225C77EMBXP02dsmanacuk_ Content-Type: text/plain; charset="iso-8859-1" Content-Transfer-Encoding: quoted-printable Dear Bijan, TopMoD can be found here, in its 2009 version : http://www.lct.jussieu.fr/p= agesperso/silvi/topmod.html I personally use DGrid, that is a bit more recent and still maintained : ht= tp://www.cpfs.mpg.de/~kohout/dgrid.html Both are free, but I did not tried these on Windows. Cheers, Fran=E7ois ________________________________ > From: owner-chemistry+francois.zielinski=3D=3Dmanchester.ac.uk],[ccl.net [own= er-chemistry+francois.zielinski=3D=3Dmanchester.ac.uk],[ccl.net] on behalf of= Bijan Mondal mondal.bijan()gmail.com [owner-chemistry],[ccl.net] Sent: Monday, July 29, 2013 5:20 PM To: Francois Zielinski Subject: CCL: ToPMoD package Dear Mates, I want to calculate ELF (electron localization function) and plot them. Can anyone suggest me any program package which will be helpful? The progra= m should be freely available for windows. I heard a package called ToPMoD can calculate ELF. Does anyone knows whethe= r it is available for windows and it is free? Thanks you very much. Any suggestions will be appreciated. Regards, Bijan -- "There are no such electron-deficient compounds, only theory-deficient chemists." -- R.E. Rundle Bijan Mondal Research Scholar Mob: 7667156593 Department of Chemistry I.I.T Madras Chennai --_000_3068BCD5C0404D4E9C44F64F4000B64225C77EMBXP02dsmanacuk_ Content-Type: text/html; charset="iso-8859-1" Content-Transfer-Encoding: quoted-printable
Dear Bijan,

TopMoD can be found here, in its 2009 version : http://www.lct.jussieu.fr/pagesperso/silvi/topmod.html
I personally use DGrid, that is a bit more recent and still maintained : http://www.cpfs.mpg.de/~kohout/dgrid.html

Both are free, but I did not tried these on Windows.

Cheers,
Fran=E7ois






From: owner-chemistry+francois.zielin= ski=3D=3Dmanchester.ac.uk],[ccl.net [owner-chemistry+francois.zielinski= =3D=3Dmanchester.ac.uk],[ccl.net] on behalf of Bijan Mondal mondal.bijan()gma= il.com [owner-chemistry],[ccl.net]
Sent: Monday, July 29, 2013 5:20 PM
To: Francois Zielinski
Subject: CCL: ToPMoD package

Dear Mates,

I want to calculate ELF (elec= tron localization function) and plot them.
Can anyone suggest me any pro= gram package which will be helpful? The program should be freely available = for windows.
I heard a package called ToPMoD can calculate ELF. Does anyone knows whethe= r it is available for windows and it is free?

Thanks you very much.
Any suggestions will be = appreciated.

Regards,
Bijan
--
"There are no such electron-deficient compounds,=
 = ;            only theory-deficient chemists." -- R.E. Rundle<= br>
 
Bijan Mondal
Research Scholar
Mob: 7667156593
Department of Chemistry
I.I.T Madras
Chennai

--_000_3068BCD5C0404D4E9C44F64F4000B64225C77EMBXP02dsmanacuk_-- From owner-chemistry@ccl.net Mon Jul 29 18:08:01 2013 From: "sobereva-$-sina.com" To: CCL Subject: CCL: ToPMoD package Message-Id: <-49036-130729163823-4160-SvKToa1BbSGTxsA/R4cw1g,server.ccl.net> X-Original-From: Content-Type: multipart/alternative; boundary="=-sinamail_alt_9204f17ef5f4a5f0e81c303bc92c0a4d" Date: Tue, 30 Jul 2013 04:01:57 +0800 MIME-Version: 1.0 Sent to CCL by: [sobereva^sina.com] --=-sinamail_alt_9204f17ef5f4a5f0e81c303bc92c0a4d Content-Type: text/plain; charset=GBK Content-Transfer-Encoding: base64 Content-Disposition: inline RGVhciBCaWphbiwNCiANCk11bHRpd2ZuIGlzIG9uZSBvZiB0aGUgYmVzdCBjaG9pY2VzLCB3aGlj aCBpcyBlYXN5LXRvLXVzZSwgZmFzdCBhbmQgZWZmaWNpZW50LCBpdHMgZXhlY3V0YWJsZSBmaWxl IGFuZCBzb3VyY2UgY29kZSBjYW4gYmUgZnJlZWx5IGRvd25sb2FkZWQgYXQgaHR0cDovL011bHRp d2ZuLmNvZGVwbGV4LmNvbS4gV2luZG93cywgTGludXggYW5kIE1hY09TIFggdmVyc2lvbnMgYXJl IGF2YWlsYWJsZS4gDQpNdWx0aXdmbiBpcyBjYXBhYmxlIHRvIHBsb3QgY3VydmUgbWFwLCBwbG90 IHBsYW5lIG1hcCwgZ2VuZXJhdGUgZ3JpZCBkYXRhIGFuZCBkaXNwbGF5IGlzb3N1cmZhY2UsIHBl cmZvcm0gdG9wb2xvZ3kgYW5hbHlzaXMgYW5kIGJhc2luIGFuYWx5c2lzIGZvciBFTEYgKGFzIHdl bGwgYXMgRUxGLXBpLCBFTEYtc2lnbWEgYW5kIGxvY2FsaXplZCBvcmJpdGFsIGxvY2F0b3IgKExP TCkpLiBUaGVyZSBhcmUgbWFueSByZWxhdGVkIGV4YW1wbGVzIGdpdmVuIGluIHRoZSBwcm9ncmFt IG1hbnVhbCwgc2VlIFNlY3Rpb24gNC4yLjIsIDQuNC4yLCA0LjUuMywgNC4xNy4yLg0KIA0KQmVz dCB3aXNoZXMsDQogDQpUaWFuIEx1DQogDQotLS0tLSBPcmlnaW5hbCBNZXNzYWdlIC0tLS0tCkZy b206ICJCaWphbiBNb25kYWwgbW9uZGFsLmJpamFuKClnbWFpbC5jb20iIDxvd25lci1jaGVtaXN0 cnlAY2NsLm5ldD4KVG86ICJMdSwgVGlhbiAtaWQjM3BoLSIgPHNvYmVyZXZhQHNpbmEuY29tPgpT dWJqZWN0OiBDQ0w6IFRvUE1vRCBwYWNrYWdlCkRhdGU6IDIwMTMtMDctMzAgMDA6MjAKDQoKDQoN Cg0KRGVhciBNYXRlcywgCgpJIHdhbnQgdG8gY2FsY3VsYXRlIEVMRiAoZWxlY3Ryb24gbG9jYWxp emF0aW9uIGZ1bmN0aW9uKSBhbmQgcGxvdCB0aGVtLiAKQ2FuIGFueW9uZSBzdWdnZXN0IG1lIGFu eSBwcm9ncmFtIHBhY2thZ2Ugd2hpY2ggd2lsbCBiZSBoZWxwZnVsPyBUaGUgcHJvZ3JhbSBzaG91 bGQgYmUgZnJlZWx5IGF2YWlsYWJsZSBmb3Igd2luZG93cy4KSSBoZWFyZCBhIHBhY2thZ2UgY2Fs bGVkIFRvUE1vRCBjYW4gY2FsY3VsYXRlIEVMRi4gRG9lcyBhbnlvbmUga25vd3Mgd2hldGhlciBp dCBpcyBhdmFpbGFibGUgZm9yIHdpbmRvd3MgYW5kIGl0IGlzIGZyZWU/DQoNCg0KCg0KVGhhbmtz IHlvdSB2ZXJ5IG11Y2guIAoNCkFueSBzdWdnZXN0aW9ucyB3aWxsIGJlIGFwcHJlY2lhdGVkLiAK Cg0KUmVnYXJkcywKDQpCaWphbgoNCi0tIAoNCiJUaGVyZSBhcmUgbm8gc3VjaCBlbGVjdHJvbi1k ZWZpY2llbnQgY29tcG91bmRzLA0KICAgICAgICAgICAgIG9ubHkgdGhlb3J5LWRlZmljaWVudCBj aGVtaXN0cy4iIC0tIFIuRS4gUnVuZGxlCg0KIA0KQmlqYW4gTW9uZGFsClJlc2VhcmNoIFNjaG9s YXINCk1vYjogNzY2NzE1NjU5MwpEZXBhcnRtZW50IG9mIENoZW1pc3RyeQpJLkkuVCBNYWRyYXMK Q2hlbm5haQo= --=-sinamail_alt_9204f17ef5f4a5f0e81c303bc92c0a4d Content-Type: text/html; charset=GBK Content-Transfer-Encoding: base64 Content-Disposition: inline PERJVj5EZWFyIEJpamFuLDwvRElWPg0KPERJVj4mbmJzcDs8L0RJVj4NCjxESVY+TXVsdGl3Zm4g aXMgb25lIG9mIHRoZSBiZXN0IGNob2ljZXMsIHdoaWNoJm5ic3A7aXMgZWFzeS10by11c2UsIGZh c3QgYW5kIGVmZmljaWVudCwgaXRzIGV4ZWN1dGFibGUgZmlsZSBhbmQgc291cmNlIGNvZGUgY2Fu IGJlIGZyZWVseSBkb3dubG9hZGVkIGF0IDxBIGhyZWY9Imh0dHA6Ly9NdWx0aXdmbi5jb2RlcGxl eC5jb20iPmh0dHA6Ly9NdWx0aXdmbi5jb2RlcGxleC5jb208L0E+LiZuYnNwO1dpbmRvd3MsIExp bnV4IGFuZCBNYWNPUyBYIHZlcnNpb25zIGFyZSBhdmFpbGFibGUuJm5ic3A7PC9ESVY+DQo8RElW Pk11bHRpd2ZuIGlzIGNhcGFibGUgdG8mbmJzcDtwbG90IGN1cnZlIG1hcCwgcGxvdCBwbGFuZSBt YXAsIGdlbmVyYXRlIGdyaWQgZGF0YSBhbmQgZGlzcGxheSBpc29zdXJmYWNlLCBwZXJmb3JtIHRv cG9sb2d5IGFuYWx5c2lzIGFuZCBiYXNpbiBhbmFseXNpcyZuYnNwO2ZvciBFTEYgKGFzIHdlbGwg YXMmbmJzcDtFTEYtcGksIEVMRi1zaWdtYSBhbmQgbG9jYWxpemVkIG9yYml0YWwgbG9jYXRvciAo TE9MKSkuIFRoZXJlIGFyZSZuYnNwO21hbnkgcmVsYXRlZCBleGFtcGxlcyBnaXZlbiBpbiB0aGUg cHJvZ3JhbSBtYW51YWwsIHNlZSBTZWN0aW9uIDQuMi4yLCA0LjQuMiwgNC41LjMsIDQuMTcuMi48 L0RJVj4NCjxESVY+Jm5ic3A7PC9ESVY+DQo8RElWPkJlc3Qgd2lzaGVzLDwvRElWPg0KPERJVj4m bmJzcDs8L0RJVj4NCjxESVY+VGlhbiBMdTwvRElWPg0KPERJVj4mbmJzcDs8L0RJVj4NCjxESVY+ LS0tLS0gT3JpZ2luYWwgTWVzc2FnZSAtLS0tLTxCUj5Gcm9tOiAiQmlqYW4gTW9uZGFsIG1vbmRh bC5iaWphbigpZ21haWwuY29tIiAmbHQ7b3duZXItY2hlbWlzdHJ5QGNjbC5uZXQmZ3Q7PEJSPlRv OiAiTHUsIFRpYW4gLWlkIzNwaC0iICZsdDtzb2JlcmV2YUBzaW5hLmNvbSZndDs8QlI+U3ViamVj dDogQ0NMOiBUb1BNb0QgcGFja2FnZTxCUj5EYXRlOiAyMDEzLTA3LTMwIDAwOjIwPEJSPjwvRElW Pg0KPERJViBpZD1vcmlnYm9keT48QlI+DQo8RElWIGRpcj1sdHI+DQo8RElWPg0KPERJVj48U1BB TiBzdHlsZT0iRk9OVC1GQU1JTFk6IHRhaG9tYSxzYW5zLXNlcmlmIj5EZWFyIE1hdGVzLCA8QlI+ PEJSPjwvU1BBTj48L0RJVj48U1BBTiBzdHlsZT0iRk9OVC1GQU1JTFk6IHRhaG9tYSxzYW5zLXNl cmlmIj5JIHdhbnQgdG8gY2FsY3VsYXRlIEVMRiAoZWxlY3Ryb24gbG9jYWxpemF0aW9uIGZ1bmN0 aW9uKSBhbmQgcGxvdCB0aGVtLiA8QlI+PC9TUEFOPjwvRElWPjxTUEFOIHN0eWxlPSJGT05ULUZB TUlMWTogdGFob21hLHNhbnMtc2VyaWYiPkNhbiBhbnlvbmUgc3VnZ2VzdCBtZSBhbnkgcHJvZ3Jh bSBwYWNrYWdlIHdoaWNoIHdpbGwgYmUgaGVscGZ1bD8gVGhlIHByb2dyYW0gc2hvdWxkIGJlIGZy ZWVseSBhdmFpbGFibGUgZm9yIHdpbmRvd3MuPEJSPkkgaGVhcmQgYSBwYWNrYWdlIGNhbGxlZCBU b1BNb0QgY2FuIGNhbGN1bGF0ZSBFTEYuIERvZXMgYW55b25lIGtub3dzIHdoZXRoZXIgaXQgaXMg YXZhaWxhYmxlIGZvciB3aW5kb3dzIGFuZCBpdCBpcyBmcmVlPzxCUiBjbGVhcj1hbGw+PC9TUEFO Pg0KPERJVj4NCjxESVY+DQo8RElWPjxTUEFOIHN0eWxlPSJGT05ULUZBTUlMWTogdGFob21hLHNh bnMtc2VyaWYiPjxCUj48L1NQQU4+PC9ESVY+DQo8RElWPjxTUEFOIHN0eWxlPSJGT05ULUZBTUlM WTogdGFob21hLHNhbnMtc2VyaWYiPlRoYW5rcyB5b3UgdmVyeSBtdWNoLiA8QlI+PC9TUEFOPjwv RElWPg0KPERJVj48U1BBTiBzdHlsZT0iRk9OVC1GQU1JTFk6IHRhaG9tYSxzYW5zLXNlcmlmIj5B bnkgc3VnZ2VzdGlvbnMgd2lsbCBiZSBhcHByZWNpYXRlZC4gPEJSPjxCUj48L1NQQU4+PC9ESVY+ DQo8RElWPjxTUEFOIHN0eWxlPSJGT05ULUZBTUlMWTogdGFob21hLHNhbnMtc2VyaWYiPlJlZ2Fy ZHMsPEJSPjwvU1BBTj48L0RJVj4NCjxESVY+PFNQQU4gc3R5bGU9IkZPTlQtRkFNSUxZOiB0YWhv bWEsc2Fucy1zZXJpZiI+QmlqYW48QlI+PC9TUEFOPjwvRElWPg0KPERJVj48U1BBTiBzdHlsZT0i Rk9OVC1GQU1JTFk6IHRhaG9tYSxzYW5zLXNlcmlmIj4tLSA8QlI+PC9TUEFOPg0KPERJVj48U1BB TiBzdHlsZT0iRk9OVC1GQU1JTFk6IHRhaG9tYSxzYW5zLXNlcmlmIj48U1BBTiBzdHlsZT0iQ09M T1I6IHJnYigwLDAsMTUzKSI+PEk+PEZPTlQgc3R5bGU9IkJBQ0tHUk9VTkQtQ09MT1I6IHJnYigy NTUsMjU1LDI1NSkiPiI8Rk9OVCBzaXplPSswPjxGT05UIHNpemU9KzA+VGhlcmUgYXJlIG5vPC9G T05UPiBzdWNoIGVsZWN0cm9uLWRlZmljaWVudCBjb21wb3VuZHMsPC9GT05UPjwvRk9OVD48L0k+ PC9TUEFOPjwvU1BBTj48L0RJVj4NCjxESVY+PFNQQU4gc3R5bGU9IkZPTlQtRkFNSUxZOiB0YWhv bWEsc2Fucy1zZXJpZiI+PFNQQU4gc3R5bGU9IkNPTE9SOiByZ2IoMCwwLDE1MykiPjxJPjxGT05U IHN0eWxlPSJCQUNLR1JPVU5ELUNPTE9SOiByZ2IoMjU1LDI1NSwyNTUpIj48Rk9OVCBzaXplPSsw PiZuYnNwOyZuYnNwOyZuYnNwOyZuYnNwOyZuYnNwOyZuYnNwOyZuYnNwOyZuYnNwOyZuYnNwOyZu YnNwOyZuYnNwOyA8Qj4mbmJzcDs8L0I+PC9GT05UPjwvRk9OVD48L0k+PC9TUEFOPjxTUEFOIHN0 eWxlPSJDT0xPUjogcmdiKDAsMCwxNTMpIj48Rk9OVCBzaXplPSswPjxGT05UIHNpemU9KzA+PEZP TlQgc2l6ZT0rMD48Rk9OVCBzdHlsZT0iQkFDS0dST1VORC1DT0xPUjogcmdiKDI1NSwyNTUsMjU1 KSI+PEZPTlQgY29sb3I9IzAwMDA5OT48ST5vbmx5IHRoZW9yeS1kZWZpY2llbnQgY2hlbWlzdHMu IiAtLTwvST4gUi5FLiBSdW5kbGU8QlI+PC9GT05UPjwvRk9OVD48L0ZPTlQ+PC9GT05UPjwvRk9O VD48L1NQQU4+PC9TUEFOPjwvRElWPjxTUEFOIHN0eWxlPSJGT05ULUZBTUlMWTogdGFob21hLHNh bnMtc2VyaWYiPjwvU1BBTj4NCjxESVY+PFNQQU4gc3R5bGU9IkZPTlQtRkFNSUxZOiB0YWhvbWEs c2Fucy1zZXJpZiI+PFNQQU4gc3R5bGU9IkNPTE9SOiByZ2IoMCwwLDE1MykiPjwvU1BBTj4mbmJz cDs8L1NQQU4+PC9ESVY+DQo8RElWPjxTUEFOIHN0eWxlPSJGT05ULUZBTUlMWTogdGFob21hLHNh bnMtc2VyaWYiPjxTUEFOIHN0eWxlPSJDT0xPUjogcmdiKDAsMCwxNTMpIj5CaWphbiBNb25kYWw8 QlI+UmVzZWFyY2ggU2Nob2xhcjwvU1BBTj48L1NQQU4+PC9ESVY+DQo8RElWPjxTUEFOIHN0eWxl PSJGT05ULUZBTUlMWTogdGFob21hLHNhbnMtc2VyaWYiPjxTUEFOIHN0eWxlPSJDT0xPUjogcmdi KDAsMCwxNTMpIj5Nb2I6IDc2NjcxNTY1OTM8QlI+RGVwYXJ0bWVudCBvZiBDaGVtaXN0cnk8QlI+ SS5JLlQgTWFkcmFzPEJSPkNoZW5uYWk8L1NQQU4+PEJSPjwvU1BBTj48L0RJVj48L0RJVj48L0RJ Vj48L0RJVj48L0RJVj48L0RJVj4= --=-sinamail_alt_9204f17ef5f4a5f0e81c303bc92c0a4d-- From owner-chemistry@ccl.net Mon Jul 29 23:08:00 2013 From: "Jose S ez josesaez(a)uv.es" To: CCL Subject: CCL: Singlet biradicals in CASSCF Message-Id: <-49037-130729155700-32029-ee9Z4tw7a2KL1DO1P/sxVQ]-[server.ccl.net> X-Original-From: "Jose S ez" Date: Mon, 29 Jul 2013 15:56:58 -0400 Sent to CCL by: "Jose S ez" [josesaez]|[uv.es] Hi all, Actually, I am trying to compute the dissociation energy of the C1-C2 carbon bond of 1,2-diphenylcyclobutane (neutral singlet). I have tried to do so with DFT methods, but all of them fail in representing the final product as an open shell singlet biradical yielding a closed shell singlet with = 0.000 (any atom has spin density on it and testing the stability of the closed shell wavefunction against the open shell one tells me that the closed shell wavefunction is stable under the perturbations considered). Since the dissociation product is a singlet biradical, I need a multideterminantal method such as CASSCF. The problem I have is the size of the active space I must consider in my calculations. If this product would be a pure biradical, a CASSCF(2,2) would be enough, but I think that it is more a biradicaloid, so it needs a bigger active space. To clutter more this decision, both unpaired electrons are on benzylic positions. Therefore, a huge CASSCF(14,14) calculation with a proper selection of orbitals is my unique option? (the 14 orbitals of the active space come from the addition of 6 orbitals per phenyl ring plus 2 orbitals of the benzylic positions). Thank you in advance, Jose Sez