CCL:G: Two questions about ONIOM calculations: Gaussian03 vs Gaussian09
- From: "M. Carmen Ruiz Delgado"
<maycaruiz!^!hotmail.com>
- Subject: CCL:G: Two questions about ONIOM calculations: Gaussian03
vs Gaussian09
- Date: Wed, 18 Nov 2009 22:17:02 -0500
Sent to CCL by: "M. Carmen Ruiz Delgado" [maycaruiz+/-hotmail.com]
Dear all,
I would like to ask you two questions about ONIOM calculations in Gaussian:
(1) I would like to use electronic embedding for a three-layer ONIOM
calculation. In the gaussian manual it looks like this is possible, however, my
calculation crashed with this final message:
Rotational constants (GHZ): 0.0067026 0.0043168 0.0027992
Leave Link 202 at Wed Nov 18 11:26:54 2009, MaxMem= 117964800 cpu: 0.1
(Enter /export/apps/g03/l120.exe)
Electronic embedding only for two layers
Error termination via Lnk1e in /export/apps/g03/l120.exe at Wed Nov 18 11:26:55
2009.
Job cpu time: 0 days 0 hours 0 minutes 20.7 seconds.
Here, it goes the keywords of my input file (in case someone could give me some
advice):
#p oniom(b3lyp/6-31+g*:AM1:uff=qeq)=EmbedCharge scf=(maxcyc=500,tight)
# pop=chelpg nosym
(2) For a two-layer oniom calculation, I have specified a type and a charge for
the atoms in the high level region and I want to assign QEq charges for all
atoms which have charge zero or were untyped. For that purpose, I have used the
UnCharged keyword in the molecule specification input using the EmbedCharged
scheme.I got an error when using gaussian09 while exactly the same input
finishes successfully in Gaussian03.
These are the keywords of my input file:
#p oniom(b3lyp/6-31+g*:uff=uncharged)=embed scf=(maxcyc=500,tight)
# pop=chelpg nosym
This calculation is unsuccessfully completed in Gaussian 09, Revision A.02,
although it was successfully finished using Gaussian03 (when using exactly the
same input).
However, if I do not specified any particular charges for the atoms in the high
level region (using just the keyword Qeq), the calculation is completed.
This is the error that appears in the oputput:
NMagM= 0.0000000 0.0000000 2.7928460 2.7928460 2.7928460 2.7928460
Generating MM parameters.
InvSVX failed for QEq1.
Error termination via Lnk1e in /opt/g09/l101.exe at Mon Nov 2 16:47:03 2009.
Job cpu time: 0 days 0 hours 0 minutes 0.5 seconds.
I would really appreciate if someone could help me and give me some advice.
Thank you so much in advance!
Best regards,
MCarmen.
M. Carmen Ruiz Delgado, PhD
School of Chemistry and Biochemistry
Georgia Institute of Technology
901 Atlantic Dr NW
MSE Building 2100 C
Atlanta, GA 30332-0400
Phone: 404-385-8182