From owner-chemistry@ccl.net Thu Oct 4 02:54:00 2012 From: "Thomas Gkourmpis thomas.gkourmpis|a|borealisgroup.com" To: CCL Subject: CCL: NPT Dynamics Keywords in Tinker Message-Id: <-47710-121004025246-6172-2vasMp3RUHLwpnQKvALTKA,+,server.ccl.net> X-Original-From: "Thomas Gkourmpis" Date: Thu, 4 Oct 2012 02:52:44 -0400 Sent to CCL by: "Thomas Gkourmpis" [thomas.gkourmpis|*|borealisgroup.com] Hello everyone I am trying to perform an NPT molecular dynamics simulation of poly(ethylene oxide) in Tinker and i was wondering about the appropriate parameters of the necessary keywords. In the Tinker manual all the default values are set for water so any help will be most appreciated. I am interested in the appropriate parameters for the tau-temperature, tau-pressure and collision keywords. Thanks a lot in advance Thomas Gkourmpis From owner-chemistry@ccl.net Thu Oct 4 04:03:01 2012 From: "Robert W. Gora robert.gora===pwr.wroc.pl" To: CCL Subject: CCL: Resonance energy Message-Id: <-47711-121004040154-3610-rTEd6iMaa+auc+9fsk6PFw a server.ccl.net> X-Original-From: "Robert W. Gora" Content-type: multipart/mixed; boundary="Boundary_(ID_i6uMHgHihtSxkVNmzZUm/A)" Date: Thu, 04 Oct 2012 10:01:47 +0200 (CEST) MIME-version: 1.0 Sent to CCL by: "Robert W. Gora" [robert.gora-,-pwr.wroc.pl] This message is in MIME format. The first part should be readable text, while the remaining parts are likely unreadable without MIME-aware tools. --Boundary_(ID_i6uMHgHihtSxkVNmzZUm/A) Content-type: TEXT/PLAIN; charset=UTF-8; format=flowed Content-transfer-encoding: 8BIT Dear Amy, Dear All, > Maybe my thinking is too simplistic. It seems like you could calculate the > predicted value for the enthalpy of complete hydrogenation of your system > which would be N times (N=# double bonds) that of one double bond). Then > subtract the actual predicted value for the complete hydrogenation (which > will be less due to the resonance stabilization). That difference should be > the resonance energy of the system. I don't think what you propose to calculate constitutes the resonance stabilization energy even though the enthalpy of hydrogenation is obviously somehow related. Apart from the recent paper by Rashid et al. I have mentioned previously, I could recommend a discussion in the Journal of Chemical Education that originated from Robert C. Kerber's paper "If It’s Resonance, What Is Resonating?" [DOI 10.1021/ed083p223]. Personally, I am particularly fond of Donald G. Truhlar's reply "The Concept of Resonance" ibid [DOI 10.1021/ed084p781], which captures the essence of resonance as it is understood in the quantum chemistry community: "The meaning of resonance energy as it appears in valence bond theory is the lowering of the calculated ground state electronic energy when one improves a quantum mechanical model containing a single valence bond structure in the wave function to include other significant valence bond structures." Where the term structure corresponds to "a single set of spin orbitals with a spin coupling corresponding to a particular covalent or ionic bonding scheme." Best, Robert PS. By the way, I apologize Brian Salter-Duke for forgetting to mention his excellent code ;) -- Robert W. Góra, Theoretical Chemistry Group, ICHFiT, WUT [http://zcht.info] [http://www.researchgate.net/profile/Robert_Gora] [http://www.researcherid.com/rid/B-3919-2010] --Boundary_(ID_i6uMHgHihtSxkVNmzZUm/A)-- From owner-chemistry@ccl.net Thu Oct 4 05:02:00 2012 From: "Sharan sara180681%%gmail.com" To: CCL Subject: CCL:G: Polymer Construction Message-Id: <-47712-121004045919-24241-GKHktYWhdU2/66/WTL6srg a server.ccl.net> X-Original-From: Sharan Content-Type: multipart/alternative; boundary=f46d040838d31c25a404cb37f8b1 Date: Thu, 4 Oct 2012 14:29:13 +0530 MIME-Version: 1.0 Sent to CCL by: Sharan [sara180681#gmail.com] --f46d040838d31c25a404cb37f8b1 Content-Type: text/plain; charset=ISO-8859-1 Hi Can somebody help me in constructing Polymeric unit in Gaussian for optimization. I am not clear on constructing the Translational vector for polymeric unit. I went through the example given in Gaussian 03, I am not clear on few fixed values supplied in the examples in constructing the polymeric unit. Some helping hand will be greateful.. thanks in advance. With best regards ** ************************************************** Dr. I. SARAVANAN, M.Sc., M.Phil., PGDCA., PGDHRM., Ph.D., Consultant, * --f46d040838d31c25a404cb37f8b1 Content-Type: text/html; charset=ISO-8859-1 Content-Transfer-Encoding: quoted-printable
Hi
=A0
Can somebody help me in constructing Polym= eric unit in Gaussian for optimization. I am not clear on constructing the = Translational vector for polymeric unit. I=A0went through the example given= in Gaussian 03, I am not clear on few fixed values supplied in the example= s in constructing the polymeric unit.
=A0
Some helping hand will be greateful.. thanks in advance.=
=A0
With best regards
=A0
*************************************************
Dr. I= . SARAVANAN, M.Sc., M.Phil., PGDCA., PGDHRM., Ph.D.,
Consultant,
--f46d040838d31c25a404cb37f8b1-- From owner-chemistry@ccl.net Thu Oct 4 07:41:00 2012 From: "Sergio Manzetti sergio.manzetti]*[gmx.com" To: CCL Subject: CCL:G: Resonance energy Message-Id: <-47713-121004073514-17002-cpzOklGBIlSoH3KaWmOgiw|-|server.ccl.net> X-Original-From: "Sergio Manzetti" Content-Type: multipart/alternative; boundary="========GMXBoundary93581349350505163550" Date: Thu, 04 Oct 2012 13:35:04 +0200 MIME-Version: 1.0 Sent to CCL by: "Sergio Manzetti" [sergio.manzetti!=!gmx.com] --========GMXBoundary93581349350505163550 Content-Type: text/plain; charset="utf-8" Content-Transfer-Encoding: 8bit Robert, I agree on that Truhlar (and Friedman) are central references on this. However, how can a program (like Gaussian for isntance) be used to apply their model for resonance? Amy's suggestion is quite interesting, and could be a way to get a glance of resonance for certain non-modifed PAHs for instance. Sergio ----- Original Message ----- > From: Robert W. Gora robert.gora===pwr.wroc.pl Sent: 10/04/12 10:01 AM To: Manzetti, Sergio Subject: CCL: Resonance energy Dear Amy, Dear All, > Maybe my thinking is too simplistic. It seems like you could calculate the > predicted value for the enthalpy of complete hydrogenation of your system > which would be N times (N=# double bonds) that of one double bond). Then > subtract the actual predicted value for the complete hydrogenation (which > will be less due to the resonance stabilization). That difference should be > the resonance energy of the system. I don't think what you propose to calculate constitutes the resonance stabilization energy even though the enthalpy of hydrogenation is obviously somehow related. Apart from the recent paper by Rashid et al. I have mentioned previously, I could recommend a discussion in the Journal of Chemical Education that originated from Robert C. Kerber's paper "If It’s Resonance, What Is Resonating?" [DOI 10.1021/ed083p223]. Personally, I am particularly fond of Donald G. Truhlar's reply "The Concept of Resonance" ibid [DOI 10.1021/ed084p781], which captures the essence of resonance as it is understood in the quantum chemistry community: "The meaning of resonance energy as it appears in valence bond theory is the lowering of the calculated ground state electronic energy when one improves a quantum mechanical model containing a single valence bond structure in the wave function to include other significant valence bond structures." Where the term structure corresponds to "a single set of spin orbitals with a spin coupling corresponding to a particular covalent or ionic bonding scheme." Best, Robert PS. By the way, I apologize Brian Salter-Duke for forgetting to mention his excellent code ;) -- Robert W. Góra, Theoretical Chemistry Group, ICHFiT, WUT [http://zcht.info] [http://www.researchgate.net/profile/Robert_Gora] [http://www.researcherid.com/rid/B-3919-2010] --========GMXBoundary93581349350505163550 Content-Type: text/html; charset="utf-8" Content-Transfer-Encoding: quoted-printable Robert, = I agree on that Truhlar (and Friedman) are central references on this. Howe= ver, how can a program (like Gaussian for isntance) be used to apply their = model for resonance?
=20 Amy's suggestion is quite interesting, and could be a way to get a glance o= f resonance for certain non-modifed PAHs for instance.
=20
=20 Sergio
=20
=20

=20 =C2=A0

=20
=20

=20 ----- = Original Message -----

=20

=20 From: = Robert W. Gora robert.gora=3D=3D=3Dpwr.wroc.pl

=20

=20 Sent: = 10/04/12 10:01 AM

=20

=20 To: Ma= nzetti, Sergio

=20

=20 Subjec= t: CCL: Resonance energy

=20
=20
=20
=20
=20
Dear Amy,=20
Dear All,=20

> Maybe my thinking is too simplistic. It seems like you could calculate=
 the=20
> predicted value for the enthalpy of complete hydrogenation of your sys=
tem=20
> which would be N times (N=3D# double bonds) that of one double bond). =
Then=20
> subtract the actual predicted value for the complete hydrogenation (wh=
ich=20
> will be less due to the resonance stabilization). That difference shou=
ld be=20
> the resonance energy of the system.=20

I don't think what you propose to calculate constitutes the resonance=20
stabilization energy even though the enthalpy of hydrogenation is obviously=
=20
somehow related.=20

Apart from the recent paper by Rashid et al. I have mentioned previously, I=
=20
could recommend a discussion in the Journal of Chemical Education that=20
originated from Robert C. Kerber's paper "If It=E2=80=99s Resonance, What I=
s=20
Resonating?" [DOI 10.1021/ed083p223].=20

Personally, I am particularly fond of Donald G. Truhlar's reply "The Concep=
t=20
of Resonance" ibid [DOI 10.1021/ed084p781], which captures the essence of=
=20
resonance as it is understood in the quantum chemistry community:=20

"The meaning of resonance energy as it appears in valence bond theory is th=
e=20
lowering of the calculated ground state electronic energy when one improves=
 a=20
quantum mechanical model containing a single valence bond structure in the=
=20
wave function to include other significant valence bond structures." Where =
the=20
term structure corresponds to "a single set of spin orbitals with a spin=20
coupling corresponding to a particular covalent or ionic bonding scheme."=
=20

Best,=20
Robert=20

PS. By the way, I apologize Brian Salter-Duke for forgetting to mention his=
=20
excellent code ;)=20

--=20
Robert W. G=C3=B3ra, Theoretical Chemistry Group, ICHFiT, WUT [http://zcht.=
info]=20
[http://www.researchgate.net/profile/Robert_Gora]=20
[http://www.researcherid.com/rid/B-3919-2010]
=20
=20
=20
=20

=20 =C2=A0

=20
=20
--========GMXBoundary93581349350505163550-- From owner-chemistry@ccl.net Thu Oct 4 08:16:01 2012 From: "Gkourmpis, Thomas Thomas.Gkourmpis##borealisgroup.com" To: CCL Subject: CCL:G: Polymer Construction Message-Id: <-47714-121004063255-31754-Y9hWNwgDdo39wiKgbKz+LQ a server.ccl.net> X-Original-From: "Gkourmpis, Thomas" Content-Language: en-US Content-Type: multipart/alternative; boundary="_000_1F6EE93266D93848A521EDB3F2B53E303553DFA4CMS02mignetwork_" Date: Thu, 4 Oct 2012 12:32:43 +0200 MIME-Version: 1.0 Sent to CCL by: "Gkourmpis, Thomas" [Thomas.Gkourmpis[A]borealisgroup.com] --_000_1F6EE93266D93848A521EDB3F2B53E303553DFA4CMS02mignetwork_ Content-Type: text/plain; charset="us-ascii" Content-Transfer-Encoding: quoted-printable Hi=0D=0A=0D=0ARegarding the polymer in Gaussian the question is how many re= peating units you need=2E If you need only a few units and each is small yo= u can always construct it by hand in GaussView=2E If on the other hand you = want a "proper" chain with many units then you'll be better off if you use = any of the builders used in Molecular Dynamics packages and construct your = system there=2E Then you can convert it to any file that Gaussian can read = and use it as you see fit=2E=0D=0A=0D=0AHave I answered your question or I = misunderstood something?=0D=0A=0D=0AThomas=0D=0A=0D=0AFrom: owner-chemistry= +thomas=2Egkourmpis=3D=3Dborealisgroup=2Ecom:-:ccl=2Enet [mailto:owner-chemis= try+thomas=2Egkourmpis=3D=3Dborealisgroup=2Ecom:-:ccl=2Enet] On Behalf Of Sha= ran sara180681%%gmail=2Ecom=0D=0ASent: Thursday, October 04, 2012 10:59 AM= =0D=0ATo: Gkourmpis, Thomas=0D=0ASubject: CCL:G: Polymer Construction=0D=0A= =0D=0AHi=0D=0A=0D=0ACan somebody help me in constructing Polymeric unit in = Gaussian for optimization=2E I am not clear on constructing the Translation= al vector for polymeric unit=2E I went through the example given in Gaussia= n 03, I am not clear on few fixed values supplied in the examples in constr= ucting the polymeric unit=2E=0D=0A=0D=0ASome helping hand will be greateful= =2E=2E thanks in advance=2E=0D=0A=0D=0AWith best regards=0D=0A=0D=0A*******= ******************************************=0D=0ADr=2E I=2E SARAVANAN, M=2ES= c=2E, M=2EPhil=2E, PGDCA=2E, PGDHRM=2E, Ph=2ED=2E,=0D=0AConsultant,=0D=0A= =0D=0A=0D=0AThis Email and any files transmitted with it are confidential a= nd=0D=0Aintended solely for the use of the individual or the entity to whom= =0D=0Ait is addressed=2E If you have received this Email by error, please= =0D=0Anotify the sender and delete the material from any storage device=2E= =0D=0ABorealis extends no warranties and makes no representations as to=0D= =0Athe accuracy or completeness of the information provided=2E It is the=0D= =0Acustomer's responsibility to inspect and test our products and=0D=0Atech= nical advice in order to satisfy itself as to the suitability=0D=0Aof the p= roducts and technical advice for the customer's particular=0D=0Apurpose=2E = --_000_1F6EE93266D93848A521EDB3F2B53E303553DFA4CMS02mignetwork_ Content-Type: text/html; charset="us-ascii" Content-Transfer-Encoding: quoted-printable

Hi=

 

Regarding the polymer in Gaussian the question is how many rep= eating units you need=2E If you need only a few units and each is small you= can always construct it by hand in GaussView=2E If on the other hand you w= ant a “proper” chain with many units then you’ll be bette= r off if you use any of the builders used in Molecular Dynamics packages an= d construct your system there=2E Then you can convert it to any file that G= aussian can read and use it as you see fit=2E

 =

Have I answe= red your question or I misunderstood something?

 

Thomas

 

From: owner-= chemistry+thomas=2Egkourmpis=3D=3Dborealisgroup=2Ecom:-:ccl=2Enet [mailto:own= er-chemistry+thomas=2Egkourmpis=3D=3Dborealisgroup=2Ecom:-:ccl=2Enet] On B= ehalf Of Sharan sara180681%%gmail=2Ecom
Sent: Thursday, Octob= er 04, 2012 10:59 AM
To: Gkourmpis, Thomas
Subject: CCL= :G: Polymer Construction

 

Hi

<= p class=3DMsoNormal> 

Ca= n somebody help me in constructing Polymeric unit in Gaussian for optimizat= ion=2E I am not clear on constructing the Translational vector for polymeri= c unit=2E I went through the example given in Gaussian 03, I am not cl= ear on few fixed values supplied in the examples in constructing the polyme= ric unit=2E

 

Some helping hand will be greateful=2E= =2E thanks in advance=2E

&nbs= p;

With best regards

 

*************************************************
Dr=2E I=2E SARAVANAN, M=2ESc=2E, M=2EPhil=2E, PGDCA= =2E, PGDHRM=2E, Ph=2ED=2E,
Consultant,

=0D=0AThis Email and any files transm= itted with it are confidential and intended solely for the use of the indiv= idual or the entity to whom it is addressed=2E If you have received this Em= ail by error, please notify the sender and delete the material from any sto= rage device=2E Borealis extends no warranties and makes no representations = as to the accuracy or completeness of the information provided=2E It is the= customer's responsibility to inspect and test our products and technical a= dvice in order to satisfy itself as to the suitability of the products and = technical advice for the customer's particular purpose=2E --_000_1F6EE93266D93848A521EDB3F2B53E303553DFA4CMS02mignetwork_-- From owner-chemistry@ccl.net Thu Oct 4 09:31:00 2012 From: "Robert W. Gora robert.gora ~~ pwr.wroc.pl" To: CCL Subject: CCL:G: [CCL] Resonance energy Message-Id: <-47715-121004092005-24145-ht0dyurEnJR2VyPJHQunJw]^[server.ccl.net> X-Original-From: "Robert W. Gora" Content-id: Content-type: multipart/mixed; boundary="Boundary_(ID_vAHn+k5dgOsxFW9C+r7DOA)" Date: Thu, 04 Oct 2012 15:19:58 +0200 (CEST) MIME-version: 1.0 Sent to CCL by: "Robert W. Gora" [robert.gora(!)pwr.wroc.pl] This message is in MIME format. The first part should be readable text, while the remaining parts are likely unreadable without MIME-aware tools. --Boundary_(ID_vAHn+k5dgOsxFW9C+r7DOA) Content-id: Content-type: TEXT/PLAIN; CHARSET=utf-8; FORMAT=flowed Content-transfer-encoding: 8BIT Sergio, > Robert, I agree on that Truhlar (and Friedman) are central references on > this. However, how can a program (like Gaussian for isntance) be used to > apply their model for resonance? Well, if you insist on using Gaussian then there is a method called generalized valence-bond GVB but this is really a particular type of MCSCF and I do not think it would be usefull if you are really interested in resonance stabilization. I believe that you need a vb wavefunction for that and Brian has provided references to pretty much every possibilities. My personal preference is CASVB implemented in MOLCAS and MOLPRO, but the resonance stabilization could be extracted from any of the VB codes. In general, and oversimplifying VB is like CI but in AO basis. At some point you need to select a set of structures, which correspond to particular modes of spin coupling. If you subtract the energy of a single valence structure > from that of a full VB (within a chosen approximation), you will end up with a resonance stabilization. This is particularly easy in CASVB or SCVB implemented in VB2000 (the two methods are in fact equivalent). Have a look for instance at Grant Hill's paper in JPCA [doi:10.1021/jp057458d] for the details how this could be done for benzene-like molecules. I can sent you an input file for molcas and/or molpro if you are interested. > Amy's suggestion is quite interesting, and could be a way to get a glance of > resonance for certain non-modifed PAHs for instance. Sure, but also a bit more than that (I mean not only the resonance) :) Cheers, Robert -- Robert W. Góra, Theoretical Chemistry Group, ICHFiT, WUT [http://zcht.info] [http://www.researchgate.net/profile/Robert_Gora] [http://www.researcherid.com/rid/B-3919-2010] --Boundary_(ID_vAHn+k5dgOsxFW9C+r7DOA)-- From owner-chemistry@ccl.net Thu Oct 4 10:48:00 2012 From: "Andrew Voronkov drugdesign~~yandex.ru" To: CCL Subject: CCL: question about software for synthetic accessibility evaluation and design of synthetic routes Message-Id: <-47716-121003101042-14829-G9De3jDuegmxgWKl6f1fqQ[]server.ccl.net> X-Original-From: Andrew Voronkov Content-Transfer-Encoding: base64 Content-Type: text/html; charset=utf-8 Date: Wed, 03 Oct 2012 18:10:28 +0400 MIME-Version: 1.0 Sent to CCL by: Andrew Voronkov [drugdesign*o*yandex.ru] PGRpdj5EZWFyIFdvbGYsIFRIRVJFU0EgaXMgdW5mb3J1bmF0ZWx5IGlzIGRpc2NvbnRpbnVlZCBi eSBNb2xlY3VsYXIgTmV0d29ya3MgYW5kIHRoZXkgZG9uJ3Qgc2F5IGJ5IHdoaWNoIHJlYXNvbiBh bmQgZG9uJ3Qgd2FudCB0byBwcm92aWRlIGl0LjwvZGl2PjxkaXY+RG8geW91IGtub3cgaWYgdGhl cmUgYXJlIGFueSBvdGhlciBhbmFsb2dzIG9mIGl0PzwvZGl2PjxkaXY+SSBkb24ndCBrbm93IHdo YXQgaXMgdGhlIHJlYXNvbiB0byBkaXNjb250aW51ZSBzdWNoIGludGVyZXN0aW5nIGFuZCB1c2Vm dWwgcHJvZ3JhbS4gRXZlbiBpZiBpdCBpcyBub3QgYWNjdXJhdGUgdGhpcyBjYW4gYmUgYWx3YXlz IGNoZWNrZWQgYnkgb3JnYW5pYyBjaGVtaXN0IGFuZCBhbnl3YXkuLi5hbG1vc3QgYWxsIHNvZnR3 YXJlIHByb2dyYW1zIGFyZSBwcm92aWRlZCB3aXRob3V0IGFueSB3YXJyYW50LiBTbyBpZiBhbnlv bmUgZWxzZSBrbm93cyBzaW1pbGFyIHByb2dyYW1zLCBvciBoYXMgb2xkIHZlcnNpb24gb2YgVGhl cmVzYSwgcGxlYXNlIGxldCBtZSBrbm93LjwvZGl2PjxkaXY+VGhleSBoYXZlIG5vdyBzb2Z0d2Fy ZSwgY2FsbGVkIFN5bHZpYSwgd2hpY2ggYWxzbyBsb29rcyB2ZXJ5IGludGVyZXN0aW5nIGFuZCBJ IHdhbnQgdG8gZXZhbHVhdGUgaXQsIGJ1dCBpdCBkb2Vzbid0IHN1Z2dlc3Qgc3ludGhldGljIHJv dXRlcyB1bmZvcnR1bmF0ZWx5LjwvZGl2PjxkaXY+wqA8L2Rpdj48ZGl2PkJlc3QgcmVnYXJkcyw8 L2Rpdj48ZGl2PkFuZHJldzwvZGl2PjxkaXY+wqA8L2Rpdj48ZGl2PjEyLjA5LjIwMTIsIDAwOjQ3 LCAiV29sZiBJaGxlbmZlbGR0IHdkaXwqfHhlbWlzdHJ5LmNvbSIgJmx0O293bmVyLWNoZW1pc3Ry eUBjY2wubmV0Jmd0Ozo8L2Rpdj48YmxvY2txdW90ZSB0eXBlPSJjaXRlIj48YnIgLz5USEVSRVNB LCBieSBNb2xlY3VsYXIgTmV0d29ya3MgKDxhIGhyZWY9Imh0dHA6Ly93d3cubW9sLW5ldC5kZS8i Pnd3dy5tb2wtbmV0LmRlPC9hPik8YnIgLz48YnIgLz48YnIgLz48ZGl2Pk9uIFR1ZSwgU2VwIDEx LCAyMDEyIGF0IDk6MDMgQU0sIEFuZHJldyBWb3JvbmtvdiBkcnVnZGVzaWduLn4hfi48YSBocmVm PSJodHRwOi8veWFuZGV4LnJ1LyI+eWFuZGV4LnJ1PC9hPiA8c3BhbiBkaXI9Imx0ciI+Jmx0Ozxh IGhyZWY9Im1haWx0bzpvd25lci1jaGVtaXN0cnl+IX5jY2wubmV0IiB0YXJnZXQ9Il9ibGFuayI+ b3duZXItY2hlbWlzdHJ5fiF+Y2NsLm5ldDwvYT4mZ3Q7PC9zcGFuPiB3cm90ZTo8YnIgLz48Ymxv Y2txdW90ZSBzdHlsZT0ibWFyZ2luOjAgMCAwIDAuOGV4O2JvcmRlci1sZWZ0OjFweCAjY2NjIHNv bGlkO3BhZGRpbmctbGVmdDoxZXg7Ij48YnIgLz4gU2VudCB0byBDQ0wgYnk6IEFuZHJldyBWb3Jv bmtvdiBbZHJ1Z2Rlc2lnbioqPGEgaHJlZj0iaHR0cDovL3lhbmRleC5ydS8iIHRhcmdldD0iX2Js YW5rIj55YW5kZXgucnU8L2E+XTxiciAvPiBEZWFyIENDTCB1c2Vycyw8YnIgLz4gY2FuIHlvdSBw bGVhc2UsIHJlY29tbWVuZCBzb21lIHNvZnR3YXJlLCB3aGljaCBpcyBhYmxlIHRvIGV2YWx1YXRl IHN5bnRoZXRpYyBhY2Nlc3NpYmlsaXR5IG9mIMKgc3BlY2lmaWVkIG5ld2x5IGRlc2lnbmVkIG1v bGVjdWxlcyBhbmQgd2hpY2ggaXMgYWJsZSB0byBzdWdnZXN0IHN5bnRoZXNpcyByb3V0ZXMgZm9y IHRob3NlLCB3aGljaCBhcmUgc3ludGhlc2l6YWJsZT88YnIgLz4gPGJyIC8+IDxiciAvPiBCZXN0 IHJlZ2FyZHMsPGJyIC8+IEFuZHJldzxiciAvPiA8YnIgLz4gPGJyIC8+IDxiciAvPiAtPSBUaGlz IGlzIGF1dG9tYXRpY2FsbHkgYWRkZWQgdG8gZWFjaCBtZXNzYWdlIGJ5IHRoZSBtYWlsaW5nIHNj cmlwdCA9LSA8YnIgLz4gRS1tYWlsIHRvIHN1YnNjcmliZXJzOiA8YSBocmVmPSJtYWlsdG86Q0hF TUlTVFJZfiF+Y2NsLm5ldCI+Q0hFTUlTVFJZfiF+Y2NsLm5ldDwvYT4gb3IgdXNlOjxiciAvPiDC oCDCoCDCoCA8YSBocmVmPSJodHRwOi8vd3d3LmNjbC5uZXQvY2dpLWJpbi9jY2wvc2VuZF9jY2xf bWVzc2FnZSIgdGFyZ2V0PSJfYmxhbmsiPmh0dHA6Ly93d3cuY2NsLm5ldC9jZ2ktYmluL2NjbC9z ZW5kX2NjbF9tZXNzYWdlPC9hPjxiciAvPiA8YnIgLz4gRS1tYWlsIHRvIGFkbWluaXN0cmF0b3Jz OiA8YSBocmVmPSJtYWlsdG86Q0hFTUlTVFJZLVJFUVVFU1R+IX5jY2wubmV0Ij5DSEVNSVNUUlkt UkVRVUVTVH4hfmNjbC5uZXQ8L2E+IG9yIHVzZTxiciAvPiDCoCDCoCDCoCA8YSBocmVmPSJodHRw Oi8vd3d3LmNjbC5uZXQvY2dpLWJpbi9jY2wvc2VuZF9jY2xfbWVzc2FnZSIgdGFyZ2V0PSJfYmxh bmsiPmh0dHA6Ly93d3cuY2NsLm5ldC9jZ2ktYmluL2NjbC9zZW5kX2NjbF9tZXNzYWdlPC9hPjxi ciAvPiA8YnIgLz4gU3Vic2NyaWJlL1Vuc3Vic2NyaWJlOjxiciAvPiDCoCDCoCDCoCA8YSBocmVm PSJodHRwOi8vd3d3LmNjbC5uZXQvY2hlbWlzdHJ5L3N1Yl91bnN1Yi5zaHRtbCIgdGFyZ2V0PSJf YmxhbmsiPmh0dHA6Ly93d3cuY2NsLm5ldC9jaGVtaXN0cnkvc3ViX3Vuc3ViLnNodG1sPC9hPjxi ciAvPiA8YnIgLz4gQmVmb3JlIHBvc3RpbmcsIGNoZWNrIHdhaXQgdGltZSBhdDogPGEgaHJlZj0i aHR0cDovL3d3dy5jY2wubmV0LyIgdGFyZ2V0PSJfYmxhbmsiPmh0dHA6Ly93d3cuY2NsLm5ldDwv YT48YnIgLz4gPGJyIC8+IEpvYjogPGEgaHJlZj0iaHR0cDovL3d3dy5jY2wubmV0L2pvYnMiIHRh cmdldD0iX2JsYW5rIj5odHRwOi8vd3d3LmNjbC5uZXQvam9iczwvYT48YnIgLz4gQ29uZmVyZW5j ZXM6IDxhIGhyZWY9Imh0dHA6Ly9zZXJ2ZXIuY2NsLm5ldC9jaGVtaXN0cnkvYW5ub3VuY2VtZW50 cy9jb25mZXJlbmNlcy8iIHRhcmdldD0iX2JsYW5rIj5odHRwOi8vc2VydmVyLmNjbC5uZXQvY2hl bWlzdHJ5L2Fubm91bmNlbWVudHMvY29uZmVyZW5jZXMvPC9hPjxiciAvPiA8YnIgLz4gU2VhcmNo IE1lc3NhZ2VzOiA8YSBocmVmPSJodHRwOi8vd3d3LmNjbC5uZXQvY2hlbWlzdHJ5L3NlYXJjaGNj bC9pbmRleC5zaHRtbCIgdGFyZ2V0PSJfYmxhbmsiPmh0dHA6Ly93d3cuY2NsLm5ldC9jaGVtaXN0 cnkvc2VhcmNoY2NsL2luZGV4LnNodG1sPC9hPjxiciAvPiDCoCDCoCDCoCA8YSBocmVmPSJodHRw Oi8vd3d3LmNjbC5uZXQvc3BhbW1lcnMudHh0IiB0YXJnZXQ9Il9ibGFuayI+aHR0cDovL3d3dy5j Y2wubmV0L3NwYW1tZXJzLnR4dDwvYT48YnIgLz4gPGJyIC8+IFJURkk6IDxhIGhyZWY9Imh0dHA6 Ly93d3cuY2NsLm5ldC9jaGVtaXN0cnkvYWJvdXRjY2wvaW5zdHJ1Y3Rpb25zLyIgdGFyZ2V0PSJf YmxhbmsiPmh0dHA6Ly93d3cuY2NsLm5ldC9jaGVtaXN0cnkvYWJvdXRjY2wvaW5zdHJ1Y3Rpb25z LzwvYT48YnIgLz4gPGJyIC8+IDwvYmxvY2txdW90ZT48L2Rpdj48YnIgLz48YnIgY2xlYXI9ImFs bCIgLz48YnIgLz4tLSA8YnIgLz5Xb2xmLUQuIElobGVuZmVsZHQgLcKgIFhlbWlzdHJ5IEdtYkgg LSA8YSBocmVmPSJtYWlsdG86d2RpfiF+eGVtaXN0cnkuY29tIiB0YXJnZXQ9Il9ibGFuayI+d2Rp fiF+eGVtaXN0cnkuY29tPC9hPjxiciAvPlBob25lOiA8c3Bhbj4rNDkgNjE3NCAyMDE0NTU8L3Nw YW4+IC0gRmF4IDxzcGFuPis0OSA2MTc0IDIwOTY2NTwvc3Bhbj48YnIgLz4tLS08YnIgLz54ZW1p c3RyeSBnbWJoIOKAkyBHZXNjaMOkZnRzZsO8aHJlci9NYW5hZ2luZyBEaXJlY3RvcjogRHIuIFcu IEQuIElobGVuZmVsZHQ8YnIgLz4gQWRkcmVzczogSGFpbmhvbHp3ZWcgMTEsIEQtNjE0NjIgS8O2 bmlnc3RlaW4sIEdlcm1hbnk8YnIgLz5IUiBLw7ZuaWdzdGVpbiBCNzUyMiA6IFVzdC9WQVQgSUQg REU8c3Bhbj4yMTUzMTYzMjk8L3NwYW4+IDogRFVOUyA8c3Bhbj4zNC00MDAtMTcxOTwvc3Bhbj48 YnIgLz48YnIgLz48L2Jsb2NrcXVvdGU+ From owner-chemistry@ccl.net Thu Oct 4 12:10:00 2012 From: "Rachelle Bienstock biensto1===niehs.nih.gov" To: CCL Subject: CCL: Programs for dsigning synthesis Message-Id: <-47717-121004120846-904-2Fsauh+3JhF8rsUR7G7rPA(-)server.ccl.net> X-Original-From: "Rachelle Bienstock" Date: Thu, 4 Oct 2012 12:08:45 -0400 Sent to CCL by: "Rachelle Bienstock" [biensto1%a%niehs.nih.gov] I do not work for any of these commercial companies and am not endorsing or advertising their products. This is just informational for the people who asked about software that suggests accessibility and design of chemical synthetic routes. SimBioSys (http://www.simbiosys.com) has products ArChem and CAESA which computationally suggest chemical compound synthetic routes. Elsevier has Reaxys (https://www.reaxys.com/info/about-overview). I hope this information is helpful. I am not aware of free and/or open source products, but they might exist as well. From owner-chemistry@ccl.net Thu Oct 4 14:07:01 2012 From: "Victor Rosas Garcia rosas.victor,,gmail.com" To: CCL Subject: CCL: question about software for synthetic accessibility evaluation and design of synthetic routes Message-Id: <-47718-121004122318-6512-n/neTkqthkc5jxkRZn0YuQ.@.server.ccl.net> X-Original-From: Victor Rosas Garcia Content-Type: multipart/alternative; boundary=0015174a0e5cc01afe04cb3e2b2a Date: Thu, 4 Oct 2012 11:23:09 -0500 MIME-Version: 1.0 Sent to CCL by: Victor Rosas Garcia [rosas.victor###gmail.com] --0015174a0e5cc01afe04cb3e2b2a Content-Type: text/plain; charset=windows-1252 Content-Transfer-Encoding: quoted-printable Dear Andrew, just found this: Route Designer: A Retrosynthetic Analysis Tool Utilizing Automated Retrosynthetic Rule Generation James Law, Zsolt Zsoldos, Aniko Simon, Darryl Reid, Yang Liu, Sing Yoong Khew, A. Peter Johnson, Sarah Major, Robert A. Wade, and Howard Y. Ando J. Chem. Inf. Model., 2009, 49 (3), 593-602=95 DOI: 10.1021/ci800228y it might be of interest Victor 2012/10/3 Andrew Voronkov drugdesign~~yandex.ru > Dear Wolf, THERESA is unforunately is discontinued by Molecular Networks > and they don't say by which reason and don't want to provide it. > Do you know if there are any other analogs of it? > I don't know what is the reason to discontinue such interesting and usefu= l > program. Even if it is not accurate this can be always checked by organic > chemist and anyway...almost all software programs are provided without an= y > warrant. So if anyone else knows similar programs, or has old version of > Theresa, please let me know. > They have now software, called Sylvia, which also looks very interesting > and I want to evaluate it, but it doesn't suggest synthetic routes > unfortunately. > > Best regards, > Andrew > > 12.09.2012, 00:47, "Wolf Ihlenfeldt wdi|*|xemistry.com" < > owner-chemistry]~[ccl.net>: > > > THERESA, by Molecular Networks (www.mol-net.de) > > > On Tue, Sep 11, 2012 at 9:03 AM, Andrew Voronkov drugdesign.~!~.yandex.ru > wrote: > > > Sent to CCL by: Andrew Voronkov [drugdesign**yandex.ru] > Dear CCL users, > can you please, recommend some software, which is able to evaluate > synthetic accessibility of specified newly designed molecules and which = is > able to suggest synthesis routes for those, which are synthesizable? > > > Best regards, > Andrew > > > > -=3D This is automatically added to each message by the mailing script = =3D- > E-mail to subscribers: CHEMISTRY~!~ccl.net or use:> > E-mail to administrators: CHEMISTRY-REQUEST~!~ccl.net or use> > > > > -- > Wolf-D. Ihlenfeldt - Xemistry GmbH - wdi~!~xemistry.com > Phone: +49 6174 201455 - Fax +49 6174 209665 > --- > xemistry gmbh =96 Gesch=E4ftsf=FChrer/Managing Director: Dr. W. D. Ihlenf= eldt > Address: Hainholzweg 11, D-61462 K=F6nigstein, Germany > HR K=F6nigstein B7522 : Ust/VAT ID DE215316329 : DUNS 34-400-1719 > > --0015174a0e5cc01afe04cb3e2b2a Content-Type: text/html; charset=windows-1252 Content-Transfer-Encoding: quoted-printable Dear Andrew,

just found this:

Route Designer: A Retrosyntheti= c Analysis Tool Utilizing Automated Retrosynthetic Rule Generation
James= Law, Zsolt Zsoldos, Aniko Simon, Darryl Reid, Yang Liu, Sing Yoong
Khew= , A. Peter Johnson, Sarah Major, Robert A. Wade, and Howard Y. Ando
J. Chem. Inf. Model., 2009, 49 (3), 593-602=95 DOI: 10.1021/ci800228y
it might be of interest

Victor

2= 012/10/3 Andrew Voronkov drugdesign~~yandex.ru= <owner-chemistry]~[ccl.net>
Dear Wolf, THERESA is unforunately is d= iscontinued by Molecular Networks and they don't say by which reason an= d don't want to provide it.
Do you know if there are any other analogs of it?
I don'= t know what is the reason to discontinue such interesting and useful progra= m. Even if it is not accurate this can be always checked by organic chemist= and anyway...almost all software programs are provided without any warrant= . So if anyone else knows similar programs, or has old version of Theresa, = please let me know.
They have now software, called Sylvia, which also looks very interesti= ng and I want to evaluate it, but it doesn't suggest synthetic routes u= nfortunately.
=A0
Best regards,
Andrew
=A0
12.09.2012, 00:47, "Wolf Ihlenfeldt wdi|*|xemistry.com" <owner-chemistry]~[ccl= .net>:

THERESA, by Molecular Networks (www.mol-net.de)


On Tue, Sep 11, 2012 at 9:03 AM, Andrew Voronkov drugdesign.~!~.yandex.ru <<= a href=3D"mailto:owner-chemistry~!~ccl.net" target=3D"_blank">owner-chemist= ry~!~ccl.net> wrote:

Sent to CCL by: Andrew Voronkov [drugdesign**yandex.ru]
Dear CCL users,
can= you please, recommend some software, which is able to evaluate synthetic a= ccessibility of =A0specified newly designed molecules and which is able to = suggest synthesis routes for those, which are synthesizable?


Best regards,
Andrew



-=3D This is automat= ically added to each message by the mailing script =3D-
E-mail to subs= cribers: CHEMISTRY= ~!~ccl.net or use:
=A0 =A0 =A0 http://www.ccl.net/cgi-bin/ccl/send_ccl_message
E-mail to administrators: CHEMISTRY-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_u= nsub.shtml

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

Job: http://www.ccl.net/jobs
Conference= s: http://server.ccl.net/chemistry/announcements/conferences/=

Search Messages: http://www.ccl.net/chemistry/searchccl/index.s= html
=A0 =A0 =A0 http://www.ccl.net/spammers.txt

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



--
Wolf-D. Ihlenfeldt -=A0 Xemistry GmbH - wdi~!~xemistry.com
Phone: +49 6174= 201455 - Fax +49 6174 209665
---
xemistry gmbh = =96 Gesch=E4ftsf=FChrer/Managing Director: Dr. W. D. Ihlenfeldt
Address: Hainholzweg 11, D-61462 K=F6nigstein, Germany
HR K=F6nigstein = B7522 : Ust/VAT ID DE215316329 : DUNS 34-400-1719=


--0015174a0e5cc01afe04cb3e2b2a-- From owner-chemistry@ccl.net Thu Oct 4 16:21:00 2012 From: "lloyd riggs lloyd.riggs]^[gmx.ch" To: CCL Subject: CCL: question about software for synthetic accessibility evaluation and design of synthetic routes Message-Id: <-47719-121004145230-2371-IvAJbr+Q0jbpCp5LDjYSvA^-^server.ccl.net> X-Original-From: "lloyd riggs" Content-Transfer-Encoding: 8bit Content-Type: text/plain; charset="utf-8" Date: Thu, 04 Oct 2012 20:52:21 +0200 MIME-Version: 1.0 Sent to CCL by: "lloyd riggs" [lloyd.riggs_-_gmx.ch] Dear All I spent a good month trying to track down retrosynth software. Theres two sites with CAOS available, however a chemsketch type additional software, and a library are needed. Other than that the only other free version is RDkit in general, which has all the code to do it. Both require piecing together things and some minimal scripting. In addition a library is needed, although you can find some for free and piece them together, I believe the library itself is for sale from sources such as sigma for around 310 $US, with a few million annotated reactions in smiles type format, so basically, it would be hard to beat this later with freeware. Still I occasionally see folks attempting to gather support for a free annotated smiles or the similes like systems which I would be for. Essentially there just a spread sheet like file (the libraries) with smiles a+b=c or a+b+c=d+f, a reaction type, rxn methods, references (which may be of importance for patenting synthetic routs that work), and possibly an image file for the a , b or , c , etc...while the software is really just a pattern search, which I guess works better if the smiles portion of the search pattern can be broken down into fragments (ie the AI aspect). if your starting a company though I would just forgo this and pay the 300$ for the software and 310$ for a library, although some of these packages charge in the thousands, so you should see if they meet your needs first. If you have a month, id just write one yourself and buy the library, although even that one might be able to track down for free. Hope that helps. Sincerely, Stephan Watkins -------- Original-Nachricht -------- > Datum: Thu, 4 Oct 2012 11:23:09 -0500 > Von: "Victor Rosas Garcia rosas.victor,,gmail.com" > An: "Watkins, Steve " > Betreff: CCL: question about software for synthetic accessibility evaluation and design of synthetic routes > Dear Andrew, > > just found this: > > Route Designer: A Retrosynthetic Analysis Tool Utilizing Automated > Retrosynthetic Rule Generation > James Law, Zsolt Zsoldos, Aniko Simon, Darryl Reid, Yang Liu, Sing Yoong > Khew, A. Peter Johnson, Sarah Major, Robert A. Wade, and Howard Y. Ando > J. Chem. Inf. Model., 2009, 49 (3), 593-602• DOI: 10.1021/ci800228y > > it might be of interest > > Victor > > 2012/10/3 Andrew Voronkov drugdesign~~yandex.ru > > > > Dear Wolf, THERESA is unforunately is discontinued by Molecular Networks > > and they don't say by which reason and don't want to provide it. > > Do you know if there are any other analogs of it? > > I don't know what is the reason to discontinue such interesting and > useful > > program. Even if it is not accurate this can be always checked by > organic > > chemist and anyway...almost all software programs are provided without > any > > warrant. So if anyone else knows similar programs, or has old version of > > Theresa, please let me know. > > They have now software, called Sylvia, which also looks very interesting > > and I want to evaluate it, but it doesn't suggest synthetic routes > > unfortunately. > > > > Best regards, > > Andrew > > > > 12.09.2012, 00:47, "Wolf Ihlenfeldt wdi|*|xemistry.com" < > > owner-chemistry[-]ccl.net>: > > > > > > THERESA, by Molecular Networks (www.mol-net.de) > > > > > > On Tue, Sep 11, 2012 at 9:03 AM, Andrew Voronkov > drugdesign.~!~.yandex.ru > > wrote: > > > > > > Sent to CCL by: Andrew Voronkov [drugdesign**yandex.ru] > > Dear CCL users, > > can you please, recommend some software, which is able to evaluate > > synthetic accessibility of specified newly designed molecules and which > is > > able to suggest synthesis routes for those, which are synthesizable? > > > > > > Best regards, > > Andrew> > E-mail to subscribers: CHEMISTRY~!~ccl.net or use:> > > E-mail to administrators: CHEMISTRY-REQUEST~!~ccl.net or use> > > > > > > > > -- > > Wolf-D. Ihlenfeldt - Xemistry GmbH - wdi~!~xemistry.com > > Phone: +49 6174 201455 - Fax +49 6174 209665 > > --- > > xemistry gmbh – Geschäftsführer/Managing Director: Dr. W. D. > Ihlenfeldt > > Address: Hainholzweg 11, D-61462 Königstein, Germany > > HR Königstein B7522 : Ust/VAT ID DE215316329 : DUNS 34-400-1719 > > > > From owner-chemistry@ccl.net Thu Oct 4 16:56:00 2012 From: "lloyd riggs lloyd.riggs[#]gmx.ch" To: CCL Subject: CCL: question about software for synthetic accessibility evaluation and design of synthetic routes Message-Id: <-47720-121004150028-6724-1/aXtMAc+vw5dBAjgR1+bQ:-:server.ccl.net> X-Original-From: "lloyd riggs" Content-Transfer-Encoding: 8bit Content-Type: text/plain; charset="utf-8" Date: Thu, 04 Oct 2012 21:00:17 +0200 MIME-Version: 1.0 Sent to CCL by: "lloyd riggs" [lloyd.riggs~~gmx.ch] Dear Andrew Voronkov I spent a good month trying to track down retrosynth software. Theres two sites with CAOS and CAOS/OSET available, however a chemsketch type additional software, and a library are needed. Other than that the only other free version is RDkit in general, which has all the code to do it. Both require piecing together things and some minimal scripting. In addition a library is needed, although you can find some for free and piece them together, I believe the library itself is for sale from sources such as sigma for around 310 US, with a few million annotated reactions in smiles type format, so basically, it would be hard to beat this later with freeware. Still I occasionally see folks attempting to gather support for a free annotated smiles or the similes like systems which I would be for. Essentially there just a spread sheet like file (the libraries) with smiles a+b=c or a+b+c=d+f, a reaction type, rxn methods, references (which may be of importance for patenting synthetic routs that work), and possibly an image file for the a , b or , c , etc...while the software is really just a pattern search, which I guess works better if the smiles portion of the search pattern can be broken down into fragments (ie the AI aspect). if your starting a company though I would just forgo this and pay the 300 for the software and 310 for a library, although some of these packages charge in the thousands, so you should see if they meet your needs first. If you have a month, id just write one yourself and buy the library, although even that one might be able to track down for free. Hope that helps. Sincerely, Stephan Watkins -------- Original-Nachricht -------- > Datum: Thu, 4 Oct 2012 11:23:09 -0500 > Von: "Victor Rosas Garcia rosas.victor,,gmail.com" > An: "Watkins, Steve " > Betreff: CCL: question about software for synthetic accessibility evaluation and design of synthetic routes > Dear Andrew, > > just found this: > > Route Designer: A Retrosynthetic Analysis Tool Utilizing Automated > Retrosynthetic Rule Generation > James Law, Zsolt Zsoldos, Aniko Simon, Darryl Reid, Yang Liu, Sing Yoong > Khew, A. Peter Johnson, Sarah Major, Robert A. Wade, and Howard Y. Ando > J. Chem. Inf. Model., 2009, 49 (3), 593-602• DOI: 10.1021/ci800228y > > it might be of interest > > Victor > > 2012/10/3 Andrew Voronkov drugdesign~~yandex.ru > > > > Dear Wolf, THERESA is unforunately is discontinued by Molecular Networks > > and they don't say by which reason and don't want to provide it. > > Do you know if there are any other analogs of it? > > I don't know what is the reason to discontinue such interesting and > useful > > program. Even if it is not accurate this can be always checked by > organic > > chemist and anyway...almost all software programs are provided without > any > > warrant. So if anyone else knows similar programs, or has old version of > > Theresa, please let me know. > > They have now software, called Sylvia, which also looks very interesting > > and I want to evaluate it, but it doesn't suggest synthetic routes > > unfortunately. > > > > Best regards, > > Andrew > > > > 12.09.2012, 00:47, "Wolf Ihlenfeldt wdi|*|xemistry.com" < > > owner-chemistry[-]ccl.net>: > > > > > > THERESA, by Molecular Networks (www.mol-net.de) > > > > > > On Tue, Sep 11, 2012 at 9:03 AM, Andrew Voronkov > drugdesign.~!~.yandex.ru > > wrote: > > > > > > Sent to CCL by: Andrew Voronkov [drugdesign**yandex.ru] > > Dear CCL users, > > can you please, recommend some software, which is able to evaluate > > synthetic accessibility of specified newly designed molecules and which > is > > able to suggest synthesis routes for those, which are synthesizable? > > > > > > Best regards, > > Andrew> > E-mail to subscribers: CHEMISTRY~!~ccl.net or use:> > > E-mail to administrators: CHEMISTRY-REQUEST~!~ccl.net or use> > > > > > > > > -- > > Wolf-D. Ihlenfeldt - Xemistry GmbH - wdi~!~xemistry.com > > Phone: +49 6174 201455 - Fax +49 6174 209665 > > --- > > xemistry gmbh – Geschäftsführer/Managing Director: Dr. W. D. > Ihlenfeldt > > Address: Hainholzweg 11, D-61462 Königstein, Germany > > HR Königstein B7522 : Ust/VAT ID DE215316329 : DUNS 34-400-1719 > > > > From owner-chemistry@ccl.net Thu Oct 4 17:30:00 2012 From: "Salter-Duke, Brian James - brian.james.duke.{:}.gmail.com" To: CCL Subject: CCL:G: [CCL] Resonance energy Message-Id: <-47721-121004165702-3505-T2gvKlgiVSKjKQnhaD9ssA{:}server.ccl.net> X-Original-From: "Salter-Duke, Brian James -" Content-Disposition: inline Content-Transfer-Encoding: 8bit Content-Type: text/plain; charset=iso-8859-1 Date: Fri, 5 Oct 2012 06:56:51 +1000 MIME-Version: 1.0 Sent to CCL by: "Salter-Duke, Brian James -" [brian.james.duke.:.gmail.com] On Thu, Oct 04, 2012 at 03:19:58PM +0200, Robert W. Gora robert.gora ~~ pwr.wroc.pl wrote: > Sergio, > >Robert, I agree on that Truhlar (and Friedman) are central references > >on this. However, how can a program (like Gaussian for isntance) be > >used to apply their model for resonance? > Well, if you insist on using Gaussian then there is a method called > generalized valence-bond GVB but this is really a particular type of > MCSCF and I do not think it would be usefull if you are really > interested in resonance stabilization. I believe that you need a vb > wavefunction for that and Brian has provided references to pretty much > every possibilities. My personal preference is CASVB implemented in > MOLCAS and MOLPRO, but the resonance stabilization could be extracted > from any of the VB codes. > In general, and oversimplifying VB is like CI but in AO basis. At some > point you need to select a set of structures, which correspond to > particular modes of spin coupling. If you subtract the energy of a > single valence structure from that of a full VB (within a chosen > approximation), you will end up with a resonance stabilization. > This is particularly easy in CASVB or SCVB implemented in VB2000 (the > two methods are in fact equivalent). This is not correct. First, however, I have to point out that the term "CASVB" means different things to different people. CASVB implemented in MOLCAS and MOLPRO by David Cooper is indeed identical to SCVB. It is obtained by projecting out a SCVB function from a CASSCF wave function, either using an overlap criteria or an energy minimisation criteria. CASVB as used in VB2000 and by some other groups is identical to the appropriate CASSCF, but using a set of localised VB orbitals rather than MOs. In general there are many more structures in that CASVB than there are in SCVB. The advantage of SCVB is the relatively small number of structures giving a very close approximation to the larger number of structures in CASVB/CASSCF. In some cases, this CASVB (and CASSCF) is equivalent to SCVB. In H2, for example, CASVB(2,2) = CASSCF(2,2) = SCVB(2) = GVB. In benzene, as you can see from the paper by Karadakov et al in my last message, SCVB(6) gives an energy lowering from RHF that is 89.5% of the lowering to CASSCF(6,6). Brian. > Have a look for instance at Grant Hill's paper in JPCA > [doi:10.1021/jp057458d] for the details how this could be done for > benzene-like molecules. I can sent you an input file for molcas and/or > molpro if you are interested. > >Amy's suggestion is quite interesting, and could be a way to get a glance of > >resonance for certain non-modifed PAHs for instance. > > Sure, but also a bit more than that (I mean not only the resonance) :) > > Cheers, > Robert > > -- > Robert W. Góra, Theoretical Chemistry Group, ICHFiT, WUT [http://zcht.info] > [http://www.researchgate.net/profile/Robert_Gora] > [http://www.researcherid.com/rid/B-3919-2010] -- Brian Salter-Duke (Brian Duke) Brian.Salter-Duke[]monash.edu Adjunct Associate Professor Monash Institute of Pharmaceutical Sciences Monash University Parkville Campus, VIC 3052, Australia