CCL:G: User defined functional problem



You used the "opt freq" option.

When Gaussian launches the second job (the frequency job) it does not necessarily use all the same keywords in that job. You can check by simply looking at the output, which will show what was used. My guess is that your iops failed to be propagated forward. It will be trivial, though, to fix things as your geometry, etc., is still available. Just run the frequency as a separate job with the iops included in the keywords and you should have no trouble. (The other diagnostic, by the way, will be that the energy in the frequency will not agree with that for the last geometry of the optimization, since it is likely using pure BLYP instead of your iop specifications.)

Best regards,

Chris

On Apr 10, 2009, at 6:26 AM, Tobias Kraemer tobiask__chem.gla.ac.uk wrote:


Sent to CCL by: Tobias Kraemer [tobiask##chem.gla.ac.uk]
Dear CCLs

I came accross a result obtained with Gaussian03 that I can't explain. In order to investigate the dependence on HF exchange, I optimised a copper complex I'm interested in,
using the "user-defined DFT model" keywords as described in the manual (setting IOps 3/76 to 3/78 with BLYP keyword). Apart from different values for the HF exchange parameter,
the parameters for the remaining exchange and correlation functionals were kept constant (B3LYP values were used).
Comparing the result for the user defined functional (20% HF) calculation, with the result for the same calculation but using the B3LYP keyword, the former yields
exactly the same optimisation steps with same energies, and exactly the same final geometry and energy. However, the same calculation gives 6 negative frequencies (-197i, 196i, 195i, 127i, 93i, 29i cm-1) for this structure.
In the 'pure' B3LYP optimisation the same structure corresponds to a true minimum. How can the frequencies be so different although the calculations are in fact the same (this is actually the example from the manual).
Am I missing something here??

Here are my input files (the defined .chk-files are the same in both cases, so the initial conditions are the same)

For user-defined model:

%NProc=4
%mem=6000MB
%chk=cu2_N_no2tmemach_BLYP_hf-20.chk
#P UBLYP/gen IOp(3/76=1000002000) IOp(3/77=0720008000) IOp(3/78=0810010000) opt freq gfprint gfinput pop=full guess=read geom=check

HF=20%

1 2

Cu N O       0
6-311+G(2df)
**************
C H          0
6-31G*
**************


For B3LYP keyword


%NProc=4
%mem=6000MB
%chk=cu2_N_no2tmemach_B3LYP.chk
#P UB3LYP/gen opt freq gfprint gfinput pop=full geom=check

HF=20% with B3LYP

1 2

Cu N O       0
6-311+G(2df)
**************
C H          0
6-31G*
**************

Any ideas will be very much appreciated. Thanks.

Tobias

--
____________________________________

dipl.-chem. tobias kraemer
westCHEM, department of chemistry
university of glasgow
joseph black building
university avenue
glasgow G12 8QQ
scotland, uk
 
|phone| +44 (0)141 330 8121
|email| tobiask%x%chem.gla.ac.uk
____________________________________



-= This is automatically added to each message by the mailing script =-
To recover the email address of the author of the message, please change
the strange characters on the top line to the ^ sign. You can also
look up the X-Original-From: line in the mail header.

E-mail to subscribers: CHEMISTRY^ccl.net or use:
    http://www.ccl.net/cgi-bin/ccl/send_ccl_message

E-mail to administrators: CHEMISTRY-REQUEST^ccl.net or use
    http://www.ccl.net/cgi-bin/ccl/send_ccl_message
http://www.ccl.net/chemistry/sub_unsub.shtml

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

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

Search Messages: http://www.ccl.net/chemistry/searchccl/index.shtml

If your mail bounces > from CCL with 5.7.1 error, check:
    http://www.ccl.net/spammers.txt

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



--


Christopher J. Cramer

University of Minnesota

Department of Chemistry

207 Pleasant St. SE

Minneapolis, MN 55455-0431

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

Phone:  (612) 624-0859 || FAX:  (612) 626-2006

Mobile: (952) 297-2575

cramer^umn.edu

http://pollux.chem.umn.edu/~cramer

(website includes information about the textbook "Essentials

    of Computational Chemistry:  Theories and Models, 2nd Edition")