From owner-chemistry@ccl.net Fri Jun 10 07:57:00 2011 From: "Krishna Kuben Govender kk.govender]|[gmail.com" To: CCL Subject: CCL: Parameterization of semi-empirical methods Message-Id: <-44874-110610075557-26262-GsHsO4tbpoyhlrEOqGZcxw a server.ccl.net> X-Original-From: "Krishna Kuben Govender" Date: Fri, 10 Jun 2011 07:55:56 -0400 Sent to CCL by: "Krishna Kuben Govender" [kk.govender*_*gmail.com] Good day CCL users, I am currently working on reparameterization of some semi-empirical methods. I wish to find out if it is possible for anyone, who has already worked on parameterization of such methods, can provide me with some scripts that I maybe able to use as a starting point to begin the reparameterization process. I would appreciate any assistance with this matter, because as alot of you might already know parameterization of semi-empirical methods is a considerably difficult process and I would just like an idea of how to go about doing this. Regards, Krishna Kuben Govender University of Cape Town Scientific Computing Research Unit Rondebosch Cape Town 7701 From owner-chemistry@ccl.net Fri Jun 10 08:32:00 2011 From: "=?ISO-8859-1?Q?Juan_Arag=F3?= juan.arago{:}uv.es" To: CCL Subject: CCL:G: Correctly Determining Binding Energies Using Counterpoise Message-Id: <-44875-110610050637-4719-cJ33iu7WORPfWEq0nbDjDQ-,-server.ccl.net> X-Original-From: =?ISO-8859-1?Q?Juan_Arag=F3?= Content-Type: multipart/alternative; boundary="------------000203030601090707080209" Date: Fri, 10 Jun 2011 11:06:22 +0200 MIME-Version: 1.0 Sent to CCL by: =?ISO-8859-1?Q?Juan_Arag=F3?= [juan.arago%a%uv.es] This is a multi-part message in MIME format. --------------000203030601090707080209 Content-Type: text/plain; charset=ISO-8859-1; format=flowed Content-Transfer-Encoding: 8bit Hi Sam, I would like to clarify you some points. For instance, in a complex formed by two fragments (A and B) counterpoise correction in G09 uses five single-point energy calculations: * The energy of the complex AB with the basis set X, E(AB). * The energy of the fragments A and B with the basis set X, E(A) and E(B). * The energy of the fragments A and B, but now with the basis set of the complex, E(A') and E(B'). With these energies you can compute the counterpoise-corrected energy of complex as follows: Counterpoise-corrected energy of complex= E(AB)-E(A')-E(B'); Basis set superposition error (BSSE) energy = E(A)+E(B)-E(A')-E(B'); In this sense, you cannot talk about "counterpoise corrected energy fragment 1 or 2"the only counterpoise-corrected energy is the complex one. You can find each single-point energy in the output after each SCF. For instance, in a dimer computed with wB97X-D the energies follow this order: SCF Done: E(RwB97XD) = -2335.23274970 A.U. after 21 cycles *E(AB)* SCF Done: E(RwB97XD) = -1656.67203385 A.U. after 16 cycles *E(A')* SCF Done: E(RwB97XD) = -678.537009287 A.U. after 15 cycles *E(B')* SCF Done: E(RwB97XD) = -1656.67122073 A.U. after 16 cycles *E(A)* SCF Done: E(RwB97XD) = -678.536266897 A.U. after 15 cycles *E(B)* I hope that has been helpful. Juan El 09/06/2011 20:25, Abrash, Sam sabrash~!~richmond.edu escribió: > > Hi Folks, > > I'm about to do some calculations using the Counterpoise keyword in > G09. I need some help in calculating the binding energy of my complex > from the output. I know, for example, how and where to find the > counterpoise corrected energy, and I know that the binding energy is > > counterpoise corrected energy of complex -- counterpoise corrected > energy fragment 1 -- counterpoise corrected energy fragment 2. > > What I don't know is where to find the counterpoise corrected energy > of the fragments in the Gaussian 09 or Gaussian 03 output. I would > greatly appreciate it if someone could help me with this (the Gaussian > User's Reference isn't much help since it doesn't show the output.) I > would especially help if someone could attach at least a portion of an > output file which includes the counterpoise corrected energies of the > fragments and the complex, so I can see where they are in a real case. > > Thank you very much for your help. > > Sam > > Samuel A. Abrash > Department of Chemistry > University of Richmond > Richmond, VA 23173 > Phone: 804-289-8248 > Fax: 804-287-1897 > E-mail: sabrash:-:richmond.edu > Web-page: http://www.richmond.edu/~sabrash > > > "In 1893 Charles Hinton left Japan to become a mathematics instructor > at Princeton University, where he invented a baseball-pitching machine > that used gunpowder to propel the balls, like a cannon. After several > accidents, the device was abandoned and Hinton lost his job ..." Terry > Pratchett, Ian Steward and Jack Cohen, The Science of Diskworld III > > *From:*owner-chemistry+sabrash==richmond.edu:-:ccl.net > [mailto:owner-chemistry+sabrash==richmond.edu:-:ccl.net] *On Behalf Of > *Miguel Quiliano Meza rifaximina===gmail.com > *Sent:* Tuesday, March 22, 2011 6:15 PM > *To:* Abrash, Sam > *Subject:* CCL: Carbocation calculations, an approach > > Dear colleagues. > > My apologies if I disturb you with a couple of questions. > > I have never worked with carbocations, but now I have a good case to start > > I would like to answer some chemicals questions using computational > chemistry, but unfortunately I obtained rare energy values. > > Could you tell me please, what should I do? or what kind of chemistry > model should I use?. > > I want to answer some questions: > > *Which reaction is the most probable to happen? > *What energy is necessary to reach transitional compounds before > obtained products? > > I have performed some calculations using this: > > # opt=tight freq b3lyp/6-31+g scf=qc with the respect spin > multiplicity, but my calculations told me that compound A is most > stable than B, besides compound C is most stable than A and B. I > stopped here because I think that is impossible. > > A: is a carbocation very similar in structure like B > B is not a carbocation > C is a carbocation product of A+B > > Final energy for A = -195.60813 > > Final energy for B = -195.27260 > > Final energy for C = -390.92315 > > Thank in advance for your time. Maybe if you don´t have time please > recommen me one good paper or good tutorial. > > Best regards. > Miguel. > -- ------------------------------------------------------------- Juan Aragó March Universitat de València Instituto de Ciencia Molecular Catedrático José Beltrán nº2 E-46980 Paterna (Valencia) España Despacho 0.1.6 e-mail: juan.aragó Phone: +34 963544428 fax: +34 963543274 ------------------------------------------------------------- --------------000203030601090707080209 Content-Type: text/html; charset=ISO-8859-1 Content-Transfer-Encoding: 7bit Hi Sam,

I would like to clarify you some points. For instance, in a complex formed by two fragments (A and B) counterpoise correction in G09 uses five single-point energy calculations:
  • The energy of the complex AB with the basis set X, E(AB).
  • The energy of the fragments A and B with the basis set X, E(A) and E(B).
  • The energy of the fragments A and B, but now with the basis set of the complex, E(A') and E(B').
With these energies you can compute the counterpoise-corrected energy of complex as follows:
Counterpoise-corrected energy of complex= E(AB)-E(A')-E(B');
Basis set superposition error (BSSE) energy = E(A)+E(B)-E(A')-E(B');

In this sense, you cannot talk about "counterpoise corrected energy fragment 1 or 2"the only counterpoise-corrected energy is the complex one.
You can find each single-point energy in the output after each SCF. For instance, in a dimer computed with wB97X-D the energies follow this order:
 SCF Done:  E(RwB97XD) =  -2335.23274970     A.U. after   21 cycles    E(AB)
 SCF Done:  E(RwB97XD) =  -1656.67203385     A.U. after   16 cycles    E(A')
 SCF Done:  E(RwB97XD) =  -678.537009287     A.U. after   15 cycles    E(B')
 SCF Done:  E(RwB97XD) =  -1656.67122073     A.U. after   16 cycles    E(A)
 SCF Done:  E(RwB97XD) =  -678.536266897     A.U. after   15 cycles    E(B)

I hope that has been helpful.

Juan

  


El 09/06/2011 20:25, Abrash, Sam sabrash~!~richmond.edu escribió:

Hi Folks,

 

I’m about to do some calculations using the Counterpoise keyword in G09.  I need some help in calculating the binding energy of my complex from the output. I know, for example, how and where to find the counterpoise corrected energy, and I know that the binding energy is

 

counterpoise corrected energy of complex – counterpoise corrected energy fragment 1 – counterpoise corrected energy fragment 2.

 

What I don’t know is where to find the counterpoise corrected energy of the fragments in the Gaussian 09 or Gaussian 03 output.  I would greatly appreciate it if someone could help me with this (the Gaussian User’s Reference isn’t much help since it doesn’t show the output.)  I would especially help if someone could attach at least a portion of an output file which includes the counterpoise corrected energies of the fragments and the complex, so I can see where they are in a real case.

 

Thank you very much for your help.

 

Sam

 

Samuel A. Abrash
Department of Chemistry
University of Richmond
Richmond, VA 23173
Phone:  804-289-8248
Fax:  804-287-1897
E-mail:  sabrash:-:richmond.edu
Web-page:  http://www.richmond.edu/~sabrash

"In 1893 Charles Hinton left Japan to become a mathematics instructor at Princeton University, where he invented a baseball-pitching machine that used gunpowder to propel the balls, like a cannon.  After several accidents, the device was abandoned and Hinton lost his job ..." Terry Pratchett, Ian Steward and Jack Cohen, The Science of Diskworld III

 

From: owner-chemistry+sabrash==richmond.edu:-:ccl.net [mailto:owner-chemistry+sabrash==richmond.edu:-:ccl.net] On Behalf Of Miguel Quiliano Meza rifaximina===gmail.com
Sent: Tuesday, March 22, 2011 6:15 PM
To: Abrash, Sam
Subject: CCL: Carbocation calculations, an approach

 

Dear colleagues.

My apologies if I disturb you with a couple of questions.

I have never worked with carbocations, but now I have a good case to start

I would like to answer some chemicals questions using computational chemistry, but unfortunately I obtained rare energy values.

Could you tell me please, what should I do? or what kind of chemistry model should I use?.

I want to answer some questions:

*Which reaction is the most probable to happen?
*What energy is necessary to reach transitional compounds before obtained products?

I have performed some calculations using this:

# opt=tight freq b3lyp/6-31+g scf=qc  with the respect spin multiplicity, but my calculations told me that compound A is most stable than B, besides compound C is most stable than A and B. I stopped here because I think that is impossible.

A: is a carbocation very similar in structure like B
 B is not a carbocation
C is a carbocation product of A+B

Final energy for A = -195.60813

Final energy for B = -195.27260

Final energy for C = -390.92315

Thank in advance for your time. Maybe if you don´t have time please recommen me one good paper or good tutorial.

Best regards.
Miguel.


-- 
-------------------------------------------------------------
Juan Aragó March

Universitat de València
Instituto de Ciencia Molecular
Catedrático José Beltrán nº2
E-46980 Paterna (Valencia)
España

Despacho 0.1.6               e-mail: juan.aragó
Phone: +34 963544428         fax: +34 963543274
-------------------------------------------------------------
--------------000203030601090707080209-- From owner-chemistry@ccl.net Fri Jun 10 09:07:00 2011 From: "Gregorio Garc a gjgarcia[]ujaen.es" To: CCL Subject: CCL:G: Problems with Band structure in Gaussian use diffuse functions Message-Id: <-44876-110610053823-20324-JSlL0+tx6/pH6ePuolj3SQ-,-server.ccl.net> X-Original-From: "Gregorio Garc a" Date: Fri, 10 Jun 2011 05:38:22 -0400 Sent to CCL by: "Gregorio Garc a" [gjgarcia(0)ujaen.es] Hi all I am trying to obtain band structure (see keywords below) in Gaussian at level B3LYP/6-31+G* // B3LYP/6-31G*. # b3lyp/6-31+g*/Auto PBC=CellRange=5000 pbc=nkpoint=15 iop(5/103=5,5/98=1,5/33=1) If I use the levels B3LYP/6-31G* // B3LYP/6-31G* and B3LYP/6-311G** // B3LYP/6-31G* (basis set without diffuse functions), there are not problems. However, when I use diffuse functions I obtain the next problem: " Possibly not enough cells included in PBC. Numerical problems in PBC. Error termination via Lnk1e in /home/Gaussian/g09/l502.exe at Fri Jun 10 00:20:58 2011. Job cpu time: 0 days 16 hours 21 minutes 59.3 seconds." In such case, the solution is increase the number of cell (N) by means PBC=CellRange=N (see http://scuseria.rice.edu/gau/g09_pbc/pbc_guide.html). However, I have used up to PCB=CellRange=2000, and I have obtained the same results. I think that it is imposible to calculate band structure using basis sets with diffuse functions..... Could someone help me? Someone have calculated bandstructure using diffuse functions? Thanks in advance From owner-chemistry@ccl.net Fri Jun 10 09:41:00 2011 From: "Abrash, Sam sabrash]![richmond.edu" To: CCL Subject: CCL:G: Correctly Determining Binding Energies Using Counterpoise Message-Id: <-44877-110610093226-26960-ulm9QKr0DAmt/onxWVyslQ . server.ccl.net> X-Original-From: "Abrash, Sam" Content-Language: en-US Content-Type: multipart/alternative; boundary="_000_DB1DF332AFC21544BA642BA6B7D4654D2ED57C38A0UREXCHANGESCC_" Date: Fri, 10 Jun 2011 09:32:08 -0400 MIME-Version: 1.0 Sent to CCL by: "Abrash, Sam" [sabrash]-[richmond.edu] --_000_DB1DF332AFC21544BA642BA6B7D4654D2ED57C38A0UREXCHANGESCC_ Content-Type: text/plain; charset="iso-8859-1" Content-Transfer-Encoding: quoted-printable Thank you. Sam Samuel A. Abrash Department of Chemistry University of Richmond Richmond, VA 23173 Phone: 804-289-8248 Fax: 804-287-1897 E-mail: sabrash-x-richmond.edu Web-page: http://www.richmond.edu/~sabrash "In 1893 Charles Hinton left Japan to become a mathematics instructor at Pr= inceton University, where he invented a baseball-pitching machine that used= gunpowder to propel the balls, like a cannon. After several accidents, th= e device was abandoned and Hinton lost his job ..." Terry Pratchett, Ian St= eward and Jack Cohen, The Science of Diskworld III > From: owner-chemistry+sabrash=3D=3Drichmond.edu-x-ccl.net [mailto:owner-chemi= stry+sabrash=3D=3Drichmond.edu-x-ccl.net] On Behalf Of Juan Arag=F3 juan.arag= o{:}uv.es Sent: Friday, June 10, 2011 5:06 AM To: Abrash, Sam Subject: CCL:G: Correctly Determining Binding Energies Using Counterpoise Hi Sam, I would like to clarify you some points. For instance, in a complex formed = by two fragments (A and B) counterpoise correction in G09 uses five single-= point energy calculations: * The energy of the complex AB with the basis set X, E(AB). * The energy of the fragments A and B with the basis set X, E(A) and E(B= ). * The energy of the fragments A and B, but now with the basis set of the= complex, E(A') and E(B'). With these energies you can compute the counterpoise-corrected energy of co= mplex as follows: Counterpoise-corrected energy of complex=3D E(AB)-E(A')-E(B'); Basis set superposition error (BSSE) energy =3D E(A)+E(B)-E(A')-E(B'); In this sense, you cannot talk about "counterpoise corrected energy fragmen= t 1 or 2"the only counterpoise-corrected energy is the complex one. You can find each single-point energy in the output after each SCF. For ins= tance, in a dimer computed with wB97X-D the energies follow this order: SCF Done: E(RwB97XD) =3D -2335.23274970 A.U. after 21 cycles E(= AB) SCF Done: E(RwB97XD) =3D -1656.67203385 A.U. after 16 cycles E(= A') SCF Done: E(RwB97XD) =3D -678.537009287 A.U. after 15 cycles E(= B') SCF Done: E(RwB97XD) =3D -1656.67122073 A.U. after 16 cycles E(= A) SCF Done: E(RwB97XD) =3D -678.536266897 A.U. after 15 cycles E(= B) I hope that has been helpful. Juan El 09/06/2011 20:25, Abrash, Sam sabrash~!~richmond.edu escribi=F3: Hi Folks, I'm about to do some calculations using the Counterpoise keyword in G09. I= need some help in calculating the binding energy of my complex from the ou= tput. I know, for example, how and where to find the counterpoise corrected= energy, and I know that the binding energy is counterpoise corrected energy of complex - counterpoise corrected energy fr= agment 1 - counterpoise corrected energy fragment 2. What I don't know is where to find the counterpoise corrected energy of the= fragments in the Gaussian 09 or Gaussian 03 output. I would greatly appre= ciate it if someone could help me with this (the Gaussian User's Reference = isn't much help since it doesn't show the output.) I would especially help= if someone could attach at least a portion of an output file which include= s the counterpoise corrected energies of the fragments and the complex, so = I can see where they are in a real case. Thank you very much for your help. Sam Samuel A. Abrash Department of Chemistry University of Richmond Richmond, VA 23173 Phone: 804-289-8248 Fax: 804-287-1897 E-mail: sabrash:-:richmond.edu Web-page: http://www.richmond.edu/~sabrash "In 1893 Charles Hinton left Japan to become a mathematics instructor at Pr= inceton University, where he invented a baseball-pitching machine that used= gunpowder to propel the balls, like a cannon. After several accidents, th= e device was abandoned and Hinton lost his job ..." Terry Pratchett, Ian St= eward and Jack Cohen, The Science of Diskworld III > From: owner-chemistry+sabrash=3D=3Drichmond.edu:-:ccl.net [mailto:owner-che= mistry+sabrash=3D=3Drichmond.edu:-:ccl.net] On Behalf Of Miguel Quiliano Me= za rifaximina=3D=3D=3Dgmail.com Sent: Tuesday, March 22, 2011 6:15 PM To: Abrash, Sam Subject: CCL: Carbocation calculations, an approach Dear colleagues. My apologies if I disturb you with a couple of questions. I have never worked with carbocations, but now I have a good case to start I would like to answer some chemicals questions using computational chemist= ry, but unfortunately I obtained rare energy values. Could you tell me please, what should I do? or what kind of chemistry model= should I use?. I want to answer some questions: *Which reaction is the most probable to happen? *What energy is necessary to reach transitional compounds before obtained p= roducts? I have performed some calculations using this: # opt=3Dtight freq b3lyp/6-31+g scf=3Dqc with the respect spin multiplicit= y, but my calculations told me that compound A is most stable than B, besid= es compound C is most stable than A and B. I stopped here because I think t= hat is impossible. A: is a carbocation very similar in structure like B B is not a carbocation C is a carbocation product of A+B Final energy for A =3D -195.60813 Final energy for B =3D -195.27260 Final energy for C =3D -390.92315 Thank in advance for your time. Maybe if you don=B4t have time please recom= men me one good paper or good tutorial. Best regards. Miguel. -- ------------------------------------------------------------- Juan Arag=F3 March Universitat de Val=E8ncia Instituto de Ciencia Molecular Catedr=E1tico Jos=E9 Beltr=E1n n=BA2 E-46980 Paterna (Valencia) Espa=F1a Despacho 0.1.6 e-mail: juan.arag=F3 Phone: +34 963544428 fax: +34 963543274 ------------------------------------------------------------- --_000_DB1DF332AFC21544BA642BA6B7D4654D2ED57C38A0UREXCHANGESCC_ Content-Type: text/html; charset="iso-8859-1" Content-Transfer-Encoding: quoted-printable

= Thank you.

 

Sam

 

Samuel A. Abrash
Department of Chem= istry
University of Richmond
Richmond, VA 23173
Phone:=A0 804-289-= 8248
Fax:=A0 804-287-1897
E-mail:=A0 sabrash-x-richmond.edu
Web-page:=A0 http://www.richmond.edu/~sabrash

"In 1893 Charles Hinton left= Japan to become a mathematics instructor at Princeton University, where he= invented a baseball-pitching machine that used gunpowder to propel th= e balls, like a cannon.  After several accidents, the device was aband= oned and Hinton lost his job ..." Terry Pratchett, Ian Steward and&nbs= p;Jack Cohen, The Science of Diskworld III

 

From: owner-= chemistry+sabrash=3D=3Drichmond.edu-x-ccl.net [mailto:owner-chemistry+sabrash= =3D=3Drichmond.edu-x-ccl.net] On Behalf Of Juan Arag=F3 juan.arago{:}u= v.es
Sent: Friday, June 10, 2011 5:06 AM
To: Abrash, Sa= m
Subject: CCL:G: Correctly Determining Binding Energies Using Co= unterpoise

 = ;

Hi Sam,

I would like to clarify you = some points. For instance, in a complex formed by two fragments (A and B) c= ounterpoise correction in G09 uses five single-point energy calculations: <= o:p>

  • The energy of= the complex AB with the basis set X, E(AB).
  • The energy of the fragments A and B with the basis set X, = E(A) and E(B).
  • The energy of= the fragments A and B, but now with the basis set of the complex, E(A') an= d E(B').

With these energies you c= an compute the counterpoise-corrected energy of complex as follows:
Coun= terpoise-corrected energy of complex=3D E(AB)-E(A')-E(B');
Basis set sup= erposition error (BSSE) energy =3D E(A)+E(B)-E(A')-E(B');

In this se= nse, you cannot talk about "counterpoi= se corrected energy fragment 1 or 2"the only counterpoise-corre= cted energy is the complex one.
You can find each single-point energy i= n the output after each SCF. For instance, in a dimer computed with wB97X-D= the energies follow this order:
 SCF Done:  E(RwB97XD) =3D&n= bsp; -2335.23274970     A.U. after   21 cycle= s    E(AB)
 SCF Done:  E(RwB97XD) =3D&nb= sp; -1656.67203385     A.U. after   16 cycles=     E(A')
 SCF Done:  E(RwB97XD) =3D&nbs= p; -678.537009287     A.U. after   15 cycles&= nbsp;   E(B')
 SCF Done:  E(RwB97XD) =3D = ; -1656.67122073     A.U. after   16 cycles&n= bsp;   E(A)
 SCF Done:  E(RwB97XD) =3D  = -678.536266897     A.U. after   15 cycles&nbs= p;   E(B)

I hope that has been helpful.

Juan=

  


El 09/06/2011 20:25, Abrash, Sam sabrash~!= ~richmond.edu escribi=F3:

Hi Folks,

 

I’m about to do some calculations using the Coun= terpoise keyword in G09.  I need some help in calculating the binding = energy of my complex from the output. I know, for example, how and where to= find the counterpoise corrected energy, and I know that the binding energy= is

&= nbsp;

= counterpoise corrected energy of complex – counterpoise corrected ene= rgy fragment 1 – counterpoise corrected energy fragment 2.

 

What I don̵= 7;t know is where to find the counterpoise corrected energy of the fragment= s in the Gaussian 09 or Gaussian 03 output.  I would greatly appreciat= e it if someone could help me with this (the Gaussian User’s Referenc= e isn’t much help since it doesn’t show the output.)  I wo= uld especially help if someone could attach at least a portion of an output= file which includes the counterpoise corrected energies of the fragments a= nd the complex, so I can see where they are in a real case.

 

Thank you very much = for your help.

 

Sam

 

Samuel A. Abrash
Department of Chemistry
University of Richm= ond
Richmond, VA 23173
Phone:  804-289-8248
Fax:  804-28= 7-1897
E-mail:  sabrash:-= :richmond.edu
Web-page:  http://www.richmond.edu/~sabrash

"In 1893 Charles Hinton left Japan to = become a mathematics instructor at Princeton University, where he invented = a baseball-pitching machine that used gunpowder to propel the balls, l= ike a cannon.  After several accidents, the device was abandoned and H= inton lost his job ..." Terry Pratchett, Ian Steward and Jack Coh= en, The Science of Diskworld III

=  

From: owner-chemistry+sabrash=3D=3Drichmond.edu:= -:ccl.net [mailto:owner-chemistry+sabrash=3D=3Drichmond.edu:-:ccl.net] O= n Behalf Of Miguel Quiliano Meza rifaximina=3D=3D=3Dgmail.com
Sen= t: Tuesday, March 22, 2011 6:15 PM
To: Abrash, Sam
Subj= ect: CCL: Carbocation calculations, an approach

 

Dear colleagues.

My apologies if I disturb yo= u with a couple of questions.

I have never worked with carbocations,= but now I have a good case to start

I would like to answer some che= micals questions using computational chemistry, but unfortunately I obtaine= d rare energy values.

Could you tell me please, what should I do? or= what kind of chemistry model should I use?.

I want to answer some = questions:

*Which reaction is the most probable to happen?
*What = energy is necessary to reach transitional compounds before obtained product= s?

I have performed some calculations using this:

# opt=3Dtig= ht freq b3lyp/6-31+g scf=3Dqc  with the respect spin multiplicity, but= my calculations told me that compound A is most stable than B, besides com= pound C is most stable than A and B. I stopped here because I think that is= impossible.

A: is a carbocation very similar in structure like B B is not a carbocation
C is a carbocation product of A+B

F= inal energy for A =3D -195.60813

Final energy for B =3D -195.27260
Final energy for C =3D -390.92315

Thank in advance for your ti= me. Maybe if you don=B4t have time please recommen me one good paper or goo= d tutorial.

Best regards.
Miguel.




-- 
------=
-------------------------------------------------------
Juan Arag=F3 March
 
Unive=
rsitat de Val=E8ncia
Instituto de Ciencia Molecular
Catedr=E1tico Jos=E9 Beltr=E1n n=BA2
E-46980 Paterna (Valencia)
Espa=F1a
=
 
Despacho 0.1.6=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0 e-mail: juan.arag=F3
Phone: +34 963544=
428=A0=A0=A0=A0=A0=A0=A0=A0 fax: +34 963543274
-------=
------------------------------------------------------
= --_000_DB1DF332AFC21544BA642BA6B7D4654D2ED57C38A0UREXCHANGESCC_-- From owner-chemistry@ccl.net Fri Jun 10 11:38:01 2011 From: "Devendra Mani devdmani10#%#gmail.com" To: CCL Subject: CCL:G: how to calculate distortion constants Message-Id: <-44878-110610112206-17932-eSc6C2yCCSce8pLLLmWjCA::server.ccl.net> X-Original-From: Devendra Mani Content-Type: multipart/alternative; boundary=90e6ba3097147e5c0604a55d2166 Date: Fri, 10 Jun 2011 20:51:59 +0530 MIME-Version: 1.0 Sent to CCL by: Devendra Mani [devdmani10{=}gmail.com] --90e6ba3097147e5c0604a55d2166 Content-Type: text/plain; charset=ISO-8859-1 Dear Abhishek Shahi, You can do these calculations on gaussian. Just use freq=vibrot command in the input line. Hope it will work for you. Devendra Mani On Wed, Jun 8, 2011 at 9:41 AM, ABHISHEK SHAHI shahi.abhishek1984|*| gmail.com wrote: > Dear All, > > > I am looking for a program which can calculate the distortion constants > (Dj, Djk etc...).Thanks in advance. > > > > With regards; > *ABHISHEK SHAHI* > > *Research Scholar > Inorganic and Physical Chemistry > Indian Institute Of Science > Bangalore-12* > INDIA > > > --90e6ba3097147e5c0604a55d2166 Content-Type: text/html; charset=ISO-8859-1 Content-Transfer-Encoding: quoted-printable Dear Abhishek Shahi,
=A0You can do these calculations on gaussian. Just= use freq=3Dvibrot command in the input line.

Hope it will work for= you.

Devendra Mani

On Wed, Jun = 8, 2011 at 9:41 AM, ABHISHEK SHAHI shahi.abhishek1984|*|gmail.com <owner-chemistry^ccl.net> wrote:
Dear All,


=A0=A0 I am looking f= or a program which can calculate the distortion constants (Dj, Djk etc...).= Thanks in advance.

=A0
=
With regards;

=A0ABHISHEK SHAHI
Research Scholar
Inorganic and= Physical Chemistry
In= dian Institute Of Science
Bangalore-12

INDIA



--90e6ba3097147e5c0604a55d2166--