From owner-chemistry@ccl.net Wed Sep 3 00:33:00 2014 From: "Joaquin Barroso Flores joaco_barroso_-_yahoo.com" To: CCL Subject: CCL:G: Error in Gaussian 09 optimization Message-Id: <-50469-140903003216-21209-7qYS4o1cQqV8XrgTJXNjBA _ server.ccl.net> X-Original-From: Joaquin Barroso Flores Content-Type: multipart/alternative; boundary="850476201-980350817-1409718728=:1857" Date: Tue, 2 Sep 2014 21:32:08 -0700 MIME-Version: 1.0 Sent to CCL by: Joaquin Barroso Flores [joaco_barroso/./yahoo.com] --850476201-980350817-1409718728=:1857 Content-Type: text/plain; charset=iso-8859-1 Content-Transfer-Encoding: quoted-printable Hello Ahmad,=0A=0AThis usually gets fixed by either modifying your structur= e (around those torsion angles marked in your outpu) or by re-defining your= starting coordinates so Gaussian wont try to relate those atoms during the= optimization process (yielding unreasonable changes that ultimately end up= in the error you just got)=0A=0AI hope this helps.=0A=0AHave a nice day.= =0A=0APS Hope to be back on CCL for good this time=0A =0A=3D=3D=3D=3D=3D=3D= =3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D= =3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=0A= Joaquin Barroso-Flores, Ph. D.=0ACentro Conjunto de Investigacion en Quimic= a Sustentable=0AInstituto de Quimica=0AUNAM=0A=0A=0A-> http://joaquinbarros= o.wordpress.com=0A=0A=0Ajoaquin.barroso''a''gmail.com =0A=0A=0A"Blogastrono= m=EDa": http://joaquinbarroso.blogspot.com =0A=3D=3D=3D=3D=3D=3D=3D=3D=3D= =3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D= =3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=0A=0A=0AEl = Jueves, 28 de agosto, 2014 4:39:26, Ahmad Ziaee ahmad.ziaee(~)ul.ie escribi=F3:=0A =0A=0A=0A=0ASent to CCL by: "Ahmad Ziaee= " [ahmad.ziaee*_*ul.ie]=0AHi all,=0AI want to do the below optimization pro= cess:=0A%nprocshared=3D24=0A# opt=3D(calcfc,tight,maxcycle=3D1024) ub97d/3-= 21g scf=3D(qc,tight,maxcycle=3D1024)=0Ageom=3Dconnectivity integral=3Dultra= finegrid=0A=0ATitle Card Required=0A=0A0 2=0AZn -0.20065400 = 0.19049100 2.22860600=0AF 1.34990800 0.67092900 = 1.24805800=0AF 1.67521800 2.64763000 -0.18304500=0A= N -1.44138700 1.16109900 0.85702200=0AC = -1.05320700 1.88469700 -0.23211400=0AH 0.01164500 = 2.01628700 -0.42555100=0ASi 2.41313000 1.09862700 = -0.15375400=0AZn 5.34478000 0.88303200 -1.42285000=0A= Zn -4.67082600 3.08025100 -1.95029800=0AZn = -1.26590900 -5.38007900 -1.49424200=0AF 3.43154500 = 1.46458500 -1.56276800=0AF -1.86112300 -0.44482400 = 3.02037200=0AF 3.01239700 -0.51807100 -0.12208700=0A= F 1.13976800 0.53798600 -1.18046400=0AF = 3.63262500 1.61772900 0.89209800=0AN -3.36905000 = 2.29240900 -0.81848000=0AN -0.34777200 -1.50530300 = 0.99222800=0AN -0.89450400 -3.81306700 -0.48466800=0A= N 0.82204500 -0.87752800 3.70169800=0AN = -0.08118200 1.82297900 3.48669800=0AC -2.02693300 = 2.44745200 -1.06396400=0AC -3.73045400 1.55416500 = 0.28744200=0AC -2.77088100 0.99278100 1.12042900=0A= C 0.49719000 -1.88319000 -0.01024000=0AC = 0.21257200 -3.04265600 -0.74290500=0AC -1.72729800 = -3.40658300 0.53459500=0AC -1.45300700 -2.25843300 = 1.26897300=0AH -4.79971700 1.42631500 0.48680500=0A= H -2.98838300 0.38794800 2.00847400=0AH = -1.73152700 3.03396900 -1.94035400=0AH 1.38022400 = -1.27770900 -0.21937800=0AH -2.61075000 -4.01951400 = 0.74582900=0AH -2.07652500 -1.87323000 2.08532200=0A= H 0.88149700 -3.36289100 -1.54933900=0A=0A1 2 1.0 20 1= .0=0A2 7 1.0=0A3 7 1.0=0A4 5 1.5 23 1.5=0A5 6 1.0 21 1.5=0A6=0A7 11 1.0 13 = 1.0 14 1.0 15 1.0=0A8=0A9 16 1.0=0A10 18 1.0=0A11=0A12=0A13=0A14=0A15=0A16 = 21 1.5 22 1.5=0A17 24 1.5 27 1.5=0A"ub97d-1.gjf" [dos] 77L, 2577C =0A= =0Aalthough the SCF is done and I have the final energy I get this error in= the =0Aoptimization process:=0A=0ABerny optimization.=0AUsing GEDIIS/GDIIS= optimizer.=0ABend failed for angle 3 - 6 - 5=0ATors failed for= dihedral 7 - 3 - 6 - 5=0ATors failed for dihedral 4 - = 5 - 6 - 3=0ATors failed for dihedral 21 - 5 - 6 - = 3=0AFormBX had a problem.=0AError termination via Lnk1e in /ichec/package= s/gaussian/09d01/l103.exe at =0AWed Aug 27 23:51:19 2014.=0AJob cpu time: = 5 days 21 hours 58 minutes 55.8 seconds.=0A=0ACan any body help me thr= ough this. the structure is a metal-organic molecule. =0A=0ABest regards,= =0AAhmad=0A=0A=0A=0A-=3D This is automatically added to each message by the= mailing script =3D-=0ATo recover the email address of the author of the me= ssage, please change=0Athe strange characters on the top line to the ~ sign= . You can also=0A=0A= =0A=0A http://www.ccl.= net/cgi-bin/ccl/send_ccl_message=0A=0AE-mail to administrators: CHEMISTRY-R= EQUEST ~ ccl.net or use=0A http://www.ccl.net/cgi-bin/ccl/send_ccl_messa= ge=0A=0A=0A http://www.ccl.net/chemistry/sub_un= sub.shtml=0A=0A=0A=0A==0AConferences: http://server.ccl.net/chemistr= y/announcements/conferences/=0A=0ASearch Messages: http://www.ccl.net/chemi= stry/searchccl/index.shtml=0A=0AIf your mail bounces from CCL with 5.7.1 er= ror, check:=0A=0A=0ARTFI: http://www.c= cl.net/chemistry/aboutccl/instructions/ --850476201-980350817-1409718728=:1857 Content-Type: text/html; charset=iso-8859-1 Content-Transfer-Encoding: quoted-printable
Hello Ahmad,

This usually gets fixed by either modifying your str= ucture (around those torsion angles marked in your outpu) or by re-defining= your starting coordinates so Gaussian wont try to relate those atoms durin= g the optimization process (yielding unreasonable changes that ultimately e= nd up in the error you just got)

I hope this helps.

Have a nice day.

PS Hope to be back on CCL for good this time
 
=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D= =3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D= =3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D
Joaquin Barros= o-Flores, Ph. D.
Centro Conjunto de Investigacion en Quimica Sustentable=
Instituto de Quimica
UNAM


-> http:/= /joaquinbarroso.wordpress.com


joaqu= in.barroso''a''gmail.com


"Blogastronom=EDa": = http://joaquinbarroso.blogspot.com
=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D= =3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D= =3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D


El Jueves, 28 de agosto, 2014 4:39:26, Ahmad= Ziaee ahmad.ziaee(~)ul.ie <owner-chemistry ~ ccl.net> escribi=F3:
=



Sent to CCL by:= "Ahmad  Ziaee" [ahmad.ziaee*_*ul.ie]
Hi all,
I want to do the b= elow optimization process:
%nprocshared=3D24
# opt=3D(calcfc,tight,ma= xcycle=3D1024) ub97d/3-21g scf=3D(qc,tight,maxcycle=3D1024)
geom=3Dconne= ctivity integral=3Dultrafinegrid

Title Card Required

0 2
= Zn                -0.20065400 =   0.19049100    2.22860600
F       =           1.34990800    0.67092900    1.24805800 F                  1.675218= 00    2.64763000  -0.18304500
N      &nb= sp;         -1.44138700    1.16109900  =   0.85702200
C              &nb= sp; -1.05320700    1.88469700  -0.23211400
H  &nb= sp;               0.01164500   = ; 2.01628700  -0.42555100
Si          &n= bsp;     2.41313000    1.09862700  -0.15375400 Zn                5.34478000&nb= sp;   0.88303200  -1.42285000
Zn        =         -4.67082600    3.08025100  -1.9= 5029800
Zn                -1.26590900  -5.38007900  -1.49424200
F   = ;               3.43154500    = 1.46458500  -1.56276800
F           = ;     -1.86112300  -0.44482400    3.02037200
= F                  3.01239700=   -0.51807100  -0.12208700
F        &nb= sp;         1.13976800    0.53798600  -= 1.18046400
F                &nb= sp; 3.63262500    1.61772900    0.89209800
N  =               -3.36905000    = 2.29240900  -0.81848000
N           = ;     -0.34777200  -1.50530300    0.99222800
= N               =20 -0.89450400  -3.81306700  -0.48466800
N     = ;             0.82204500  -0.87752800&n= bsp;   3.70169800
N             = ;   -0.08118200    1.82297900    3.48669800
C=                 -2.02693300  =   2.44745200  -1.06396400
C         = ;       -3.73045400    1.55416500    0.= 28744200
C                -2.7= 7088100    0.99278100    1.12042900
C    =               0.49719000  -1.88319= 000  -0.01024000
C             = ;     0.21257200  -3.04265600  -0.74290500
C = ;              =20 -1.72729800  -3.40658300    0.53459500
C    =             -1.45300700  -2.25843300&n= bsp;   1.26897300
H             = ;   -4.79971700    1.42631500    0.48680500
H=                 -2.98838300  =   0.38794800    2.00847400
H        =         -1.73152700    3.03396900  -1.= 94035400
H                 = ; 1.38022400  -1.27770900  -0.21937800
H     = ;           -2.61075000  -4.01951400  =   0.74582900
H              &nb= sp; -2.07652500  -1.87323000    2.08532200
H  &nb= sp;               0.88149700  -3.36289100  -1.54933900

1 2 1.0 20 1.0 2 7 1.0
3 7 1.0
4 5 1.5 23 1.5
5 6 1.0 21 1.5
6
7 11 = 1.0 13 1.0 14 1.0 15 1.0
8
9 16 1.0
10 18 1.0
11
12
= 13
14
15
16 21 1.5 22 1.5
17 24 1.5 27 1.5
"ub97d-1.gjf= " [dos] 77L, 2577C     

although the SCF is done an= d I have the final energy I get this error in the
optimization process:=

Berny optimization.
Using GEDIIS/GDIIS optimizer.
Bend fai= led for angle    3 -    6 -    5
Tors = failed for dihedral    7 -    3 -    6 -&n= bsp;   5
Tors failed for dihedral    4 -    = 5 -    6 -    3
Tors failed for dihedral  &= nbsp; 21 -    5 -    6 -    3
FormBX h= ad a problem.
Error termination via Lnk1e in /ichec/packages/gaussian/09d01/l103.exe at
Wed Aug 27 23:51:19 2014. Job cpu time:      5 days 21 hours 58 minutes 55.8 secon= ds.

Can any body help me through this. the structure is a metal-orga= nic molecule.

Best regards,
Ahmad



-=3D This is au= tomatically added to each message by the mailing script =3D-
To recover = the email address of the author of the message, please change
the strang= e characters on the top line to the ~ sign. You can also
look up the X-O= riginal-From: line in the mail header.

E-mail to subscribers: CHEMIST= RY ~ ccl.net or use:
      http://www.ccl.net/cgi-bi= n/ccl/send_ccl_message

E-mail to administrators: CHEMISTRY-REQUEST ~ ccl.net or= use
      http://www.ccl.net/cgi-bin/ccl/send_ccl_mes= sage
      http://www.= ccl.net/chemistry/sub_unsub.shtml

Before posting, check wait tim= e at: http://www.ccl.net<= /a>

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

Search Messages: h= ttp://www.ccl.net/chemistry/searchccl/index.shtml

If your mail b= ounces from CCL with 5.7.1 error, check:
      http://www.ccl.net/spammers.txt

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



--850476201-980350817-1409718728=:1857-- From owner-chemistry@ccl.net Wed Sep 3 06:23:00 2014 From: "Fco. Javier Modrego modrego^unizar.es" To: CCL Subject: CCL:G: Error in Gaussian 09 optimization Message-Id: <-50470-140903055116-9150-lh9FDXCUFRGoU9H75lzgPA(a)server.ccl.net> X-Original-From: "Fco. Javier Modrego" Content-Transfer-Encoding: 8bit Content-Type: text/plain; charset=iso-8859-1 Date: Wed, 3 Sep 2014 11:51:06 +0200 Mime-Version: 1.0 (Apple Message framework v1085) Sent to CCL by: "Fco. Javier Modrego" [modrego*o*unizar.es] El 28/08/2014, a las 11:06, Ahmad Ziaee ahmad.ziaee(~)ul.ie escribió: It seems you have some linear bends around the 3-5-6 angle. You can try several options: You can kill all the dihedrals related to the 3-6-5 angle or make an explicit definition of the linear bends involved. A simpler option is to restart the optimization at point it has stopped reading the current geometry and forcing Gaussian to rebuild the set of red. internals using geom=(checkpoint,newdefinition) I hope this can help F.J. Modrego > > Sent to CCL by: "Ahmad Ziaee" [ahmad.ziaee*_*ul.ie] > Hi all, > I want to do the below optimization process: .... > Berny optimization. > Using GEDIIS/GDIIS optimizer. > Bend failed for angle 3 - 6 - 5 > Tors failed for dihedral 7 - 3 - 6 - 5 > Tors failed for dihedral 4 - 5 - 6 - 3 > Tors failed for dihedral 21 - 5 - 6 - 3 > FormBX had a problem. > Error termination via Lnk1e in /ichec/packages/gaussian/09d01/l103.exe at > Wed Aug 27 23:51:19 2014. > Job cpu time: 5 days 21 hours 58 minutes 55.8 seconds. > Dr. F.J. Modrego Department of Inorganic Chemistry Facultad de Ciencias University of Zaragoza 50009 ZARAGOZA SPAIN Tel <34>-976-762288 Fax <34>-976-761187 E-mail: modrego]_[unizar.es From owner-chemistry@ccl.net Wed Sep 3 07:15:01 2014 From: "Arash Edrisi arashedr(0)student.uef.fi" To: CCL Subject: CCL: Convergence failure Message-Id: <-50471-140903050608-24152-kL9rNYBxAPqsiWoHq2Pxww]-[server.ccl.net> X-Original-From: "Arash Edrisi" Date: Wed, 3 Sep 2014 05:06:06 -0400 Sent to CCL by: "Arash Edrisi" [arashedr()student.uef.fi] To whom it may concern, I am working with organometallic complexes and I want to optimize the structure. But I have a major problem with the optimization of the structures. I used the keywords below to launch my calculations; #p opt=(cartesian,maxcyc=200) gen gfoldprint pbe1pbe pop=full pseudo=read scfcyc=2000 However all of the calculations turned to fail and did not converge. I would be thankful if you will help me. BR, Arash From owner-chemistry@ccl.net Wed Sep 3 07:50:01 2014 From: "Igors Mihailovs igors.mihailovs0-.-gmail.com" To: CCL Subject: CCL:G: about electron affinity Message-Id: <-50472-140903053055-7353-n0BLtUmqxnVGfUWzOjrO4g(0)server.ccl.net> X-Original-From: Igors Mihailovs Content-Type: multipart/alternative; boundary=047d7b10cacfefe9c4050225e30c Date: Wed, 3 Sep 2014 12:30:24 +0300 MIME-Version: 1.0 Sent to CCL by: Igors Mihailovs [igors.mihailovs0.:.gmail.com] --047d7b10cacfefe9c4050225e30c Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: quoted-printable Also, it may be useful to note here that the electron affinity is usually negative for either rather small aromatic molecules and those containing electron-donating groups (sp2-N, too, in this case, because it has lone pair), whereas for aromatics with many fused rings and those having electron accepting groups usually have positive EAs. The experimental results also may give controversial suggestions. For example, naphthalene has EA of =E2=80=930.20 eV, a value from ETS experimen= t, according to [D. M. A. Vera and A. B. Pierini, Phys. Chem. Chem. Phys., 2004, 6, 2899=E2=80=932903.], cited in the second paper suggested by Mr. Co= ok; but in [O. Dolgounitcheva,V. G. Zakrzewski and J. V. Ortiz, J. Phys. Chem. A 1997 101 (45), 8554-8564] value for naphthalene from photoelectron spectra is +0.14 eV. With best regards, Igors Mihailovs Institute of Solid State Physics, University of Latvia 2014-09-03 4:47 GMT+03:00 Ronald Cook cookrl/atda.com < owner-chemistry[-]ccl.net>: > There are four papers that address positive and negative electron > affinities and the latter two address the computation of negative electro= n > affinities > > 1) Tozer, David J., and Frank De Proft. "Computation of the hardness and > the problem of negative electron affinities in density functional theory.= " *The > Journal of Physical Chemistry A* 109.39 (2005): 8923-8929. > > 2)De Proft, Frank, et al. "Calculation of negative electron affinity and > aqueous anion hardness using Kohn=E2=80=93Sham HOMO and LUMO energies." *= Faraday > discussions* 135 (2007): 151-159. > > 3)Puiatti, Marcelo, D. Mariano A. Vera, and Adriana B. Pierini. "Species > with negative electron affinity and standard DFT methods. Finding the > valence anions." *Physical Chemistry Chemical Physics* 10.10 (2008): > 1394-1399. > > 4)Puiatti, Marcelo, D. Mariano A. Vera, and Adriana B. Pierini. "In searc= h > for an optimal methodology to calculate the valence electron affinities o= f > temporary anions." *Physical Chemistry Chemical Physics* 11.40 (2009): > 9013-9024. > > Ronald Cook > Principal Scientist > TDA Research, Inc. > > > On Tue, Sep 2, 2014 at 11:26 AM, Acioli, Paulo p-acioli|,|neiu.edu < > owner-chemistry]_[ccl.net> wrote: > >> You will always get negative electron binding energy at the HF level. >> Negative ions require correlation. Once you co to a post-HF method, if y= our >> molecule can bind an electron, you should be able to obtain a positive >> electron affinity. >> >> Paulo Acioli, Chair Earth Science and Physics >> >> Associate Professor of Physics >> >> Department of Physics and Astronomy >> >> Northeastern Illinois University >> >> 5500 North St. Louis Avenue, Chicago, IL 60625 >> >> Phone: (773) 442-4733 >> >> p-acioli]^[neiu.edu >> >> >> *http://www.neiu.edu/academics/college-of-arts-and-sciences/departments/= physics >> * >> >> >> >> On Tue, Sep 2, 2014 at 11:15 AM, Pierre Archirel pierre.archirel++ >> u-psud.fr wrote: >> >>> >>> Sent to CCL by: "Pierre Archirel" [pierre.archirel~~u-psud.fr] >>> Dear colleagues, >>> I have a closed shell molecule (rather large: 12 first row atoms) and I >>> want to know if this molecule can bind an electron. >>> DFT tells me "yes" with positive values of the electron affinity, in th= e >>> range 0.5-1.0 eV according to the functional. I suspect that these valu= es >>> are false, or at least largely overestimated. >>> I now turn to ab initio methods, but I first see that at the HF level >>> the electron is not bound and the SOMO built of the most diffuse gaussi= ans >>> of the basis. >>> The CCSD(T) and SAC-CI methods (which are affordable in my case) give a >>> negative electron affinity... but are these methods reliable? These met= hods >>> use the HF orbitals, can they bind an electron if it is not bound at th= e HF >>> level? >>> Or shall I use CAS methods (very expansive in my case) which optimize >>> orbitals at the CI level? Is it obligatory in this case? >>> Many thanks in advance, >>> Pierre Archirel >>> LCP, Universite Paris-Sud, Orsay, France >>> >>> >>> >>> -=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>>> >>> >>> >> > --047d7b10cacfefe9c4050225e30c Content-Type: text/html; charset=UTF-8 Content-Transfer-Encoding: quoted-printable
Also, it may be useful to note here that the electron affi= nity is usually negative for either rather small aromatic molecules and tho= se containing electron-donating groups (sp2-N, too, in this case, because i= t has lone pair), whereas for aromatics with many fused rings and those hav= ing electron accepting groups usually have positive EAs.
The experimental results also may give controversial suggestions. For examp= le, naphthalene has EA of =E2=80=930.20 eV, a value from ETS experiment, ac= cording to [D. M. A. Vera and A. B. Pierini, Phys. Chem. Chem. Phys., 2004,= 6, 2899=E2=80=932903.], cited in the second paper suggested by Mr. Cook; b= ut in [O. Dolgounitcheva,V. G. Zakrzewski and J. V. Ortiz,
J. Phys. Chem. A 1997 101 (45), 8554-8564] value for naphthalene from = photoelectron spectra is +0.14 eV.

With best regards,
Igors Mihailovs
Institute of So= lid State Physics,
University of Latvia



2014-09-03 4:47 GMT+03:00 Ronald Cook co= okrl/atda.com <owner-chemistry[-]ccl.n= et>:
There are four papers that address positi= ve and negative electron affinities and the latter two address the computat= ion of negative electron affinities

1) Tozer, David J., and=20 Frank De Proft. "Computation of the hardness and the problem of negati= ve electron affinities in density functional theory." The Journal of = Physical Chemistry A 109.39 (2005): 8923-8929.

2)De Proft, Frank= , et al.=20 "Calculation of negative electron affinity and aqueous anion hardness= =20 using Kohn=E2=80=93Sham HOMO and LUMO energies." Faraday discussion= s 135 (2007): 151-159.

3)Puiatti, Marcelo, D.=20 Mariano A. Vera, and Adriana B. Pierini. "Species with negative electr= on affinity and standard DFT methods. Finding the valence anions." Ph= ysical Chemistry Chemical Physics 10.10 (2008): 1394-1399.

4)Pui= atti, Marcelo, D.=20 Mariano A. Vera, and Adriana B. Pierini. "In search for an optimal=20 methodology to calculate the valence electron affinities of temporary=20 anions." Physical Chemistry Chemical Physics 11.40 (2009): 9013= -9024.

Ronald Cook
Principal Scientist
TDA Research, Inc.


On Tue, Sep 2, 2014 at 11:26 AM, Aciol= i, Paulo p-acioli|,|neiu.edu<= /a> <owner-chemistry]_[ccl.net> wrote:
You will always get negativ= e electron binding energy at the HF level. Negative ions require correlatio= n. Once you co to a post-HF method, if your molecule can bind an electron, = you should be able to obtain a positive electron affinity.=C2=A0

Paulo Acioli, Chair Earth Science and Physics=C2=A0

Associate Professor of Physics=C2=A0

Departm= ent of Physics and Astronomy

Northeastern Illinois University

= 5500 North St. Louis Avenue, Chicago, IL 60625

Phone: (773) 442-4733<= /a>

p-acioli]^[neiu.edu


http://www.neiu.edu/academics/college-of-arts-and-sciences/departm= ents/physics

3D""



On Tue, Sep 2, 2014= at 11:15 AM, Pierre Archirel pierre.archirel++u-psud.fr <owner-chemistry]^[ccl.net&= gt; wrote:

Sent to CCL by: "Pierre=C2=A0 Archirel" [pierre.archirel~~u-psud.fr]
Dear colleagues,
I have a closed shell molecule (rather large: 12 first row atoms) and I wan= t to know if this molecule can bind an electron.
DFT tells me "yes" with positive values of the electron affinity,= in the range 0.5-1.0 eV according to the functional. I suspect that these = values are false, or at least largely overestimated.
I now turn to ab initio methods, but I first see that at the HF level the e= lectron is not bound and the SOMO built of the most diffuse gaussians of th= e basis.
The CCSD(T) and SAC-CI methods (which are affordable in my case) give a neg= ative electron affinity... but are these methods reliable? These methods us= e the HF orbitals, can they bind an electron if it is not bound at the HF l= evel?
Or shall I use CAS methods (very expansive in my case) which optimize orbit= als at the CI level? Is it obligatory in this case?
Many thanks in advance,
Pierre Archirel
LCP, Universite Paris-Sud, Orsay, France



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

E-mail to administrators: CHEMISTRY-REQUEST]^[ccl.net or use



--047d7b10cacfefe9c4050225e30c-- From owner-chemistry@ccl.net Wed Sep 3 09:06:00 2014 From: "Mozhdeh Mohammadpour mozhdehmohammadpour-*-gmail.com" To: CCL Subject: CCL: Convergence failure Message-Id: <-50473-140903090434-4468-KDoO7DoMXoegMyrFiTjnBw-x-server.ccl.net> X-Original-From: Mozhdeh Mohammadpour Content-Type: text/plain; charset=ISO-8859-1 Date: Wed, 3 Sep 2014 17:34:26 +0430 MIME-Version: 1.0 Sent to CCL by: Mozhdeh Mohammadpour [mozhdehmohammadpour^_^gmail.com] lots of work you can do to overcome this problem, 1-Try computing a Hessian from a less expensive method (MM,or HF), and reading that Hessian in for the optimization with your desire method. 2- or in this input you can use SCF=QC 3- using UltraFineGrid 4-if your stucture is openshell, you would choose the nearest optimized closed shell structure bests On 9/3/14, Arash Edrisi arashedr(0)student.uef.fi wrote: > > Sent to CCL by: "Arash Edrisi" [arashedr()student.uef.fi] > To whom it may concern, > > I am working with organometallic complexes and I want to optimize the > structure. But I have a major problem with the optimization of the > structures. > I used the keywords below to launch my calculations; > > #p opt=(cartesian,maxcyc=200) gen gfoldprint pbe1pbe pop=full pseudo=read > scfcyc=2000 > > However all of the calculations turned to fail and did not converge. > I would be thankful if you will help me. > > BR, > Arash> > > From owner-chemistry@ccl.net Wed Sep 3 09:40:00 2014 From: "Ahmad.Ziaee Ahmad.Ziaee(0)ul.ie" To: CCL Subject: CCL:G: Convergence failure Message-Id: <-50474-140903092352-14174-Pw2N6q4E+lZQRme0ftyf7w]![server.ccl.net> X-Original-From: "Ahmad.Ziaee" Content-Language: en-US Content-Transfer-Encoding: 8bit Content-Type: text/plain; charset="us-ascii" Date: Wed, 3 Sep 2014 13:23:44 +0000 MIME-Version: 1.0 Sent to CCL by: "Ahmad.Ziaee" [Ahmad.Ziaee * ul.ie] Hi Arash, 1- You should use maxcycle=N instead of maxcyc=N keyword 2- You should write the scf cycles like this: scf=maxcycle=N 3- As far as I know 'gen' keyword is useful when you want to introduce your own basis set in the Gaussian but here I see that you have use a Gaussian integrated functional and basis set. Moreover, is you are going to use your own basis set you should be much more careful about its input I recommend you to use this link for more information on input style: http://www.gaussian.com/g_tech/g_ur/k_gen.htm Hope it would be helpful for you. Best regards, Ahmad -----Original Message----- > From: owner-chemistry+ahmad.ziaee==ul.ie__ccl.net [mailto:owner-chemistry+ahmad.ziaee==ul.ie__ccl.net] On Behalf Of Arash Edrisi arashedr(0)student.uef.fi Sent: 03 September 2014 10:06 To: Ahmad.Ziaee Subject: CCL: Convergence failure Sent to CCL by: "Arash Edrisi" [arashedr()student.uef.fi] To whom it may concern, I am working with organometallic complexes and I want to optimize the structure. But I have a major problem with the optimization of the structures. I used the keywords below to launch my calculations; #p opt=(cartesian,maxcyc=200) gen gfoldprint pbe1pbe pop=full pseudo=read scfcyc=2000 However all of the calculations turned to fail and did not converge. I would be thankful if you will help me. BR, Arashhttp://www.ccl.net/cgi-bin/ccl/send_ccl_messagehttp://www.ccl.net/chemistry/sub_unsub.shtmlhttp://www.ccl.net/spammers.txt