CCL:G: calcfc Gaussian 09



Dear Raphael, 

In complement to Jean-Jules's answer I'd add this: 

calcfc is an option that calculate the full hessian matrix (second derivatives of the energy with respect to nuclear coordinates) for your system. As previously mentioned, It can be very cpu-time/memory consuming.

Here are some answers to your questions:

1- No. 
You only rarely need to start by calculating force constants, at least when optimizing a system toward a minimum of the PES. 
By default, Gaussian estimates the force constants (eq. to  opt=newestmfc) that will be used in the minimization algorithm at the beginning (with a valence force field, which is not time consuming). Furthermore, along the optimization, the hessian matrix is updated. 2nd derivatives are used in methods such as ones based on Newton-Raphson approach, e.g.

2- Same answer as above but with a precision: for a minimum search, things happen generally well, in most of the cases, calcfc is not needed.
When you are looking for a transition state, things are different: you have to provide a proper hessian matrix (see Gaussian User's manual for TS search, which can be a tricky task, sometimes), which means that either you provide it from a previous calculation or calculate it with calcfc.

3- Almost a linear combination of the two previous answers:
When looking for a minimum, the use of calcfc is rarely (even almost never) needed. But when looking for a TS, this is needed if you don't provide the hessian from a previous calculation.

To finish, I absolutely want to "modulate" my answers: every chemical system is particular and will behave a particular way. Hence, the calculation of the hessian matrix before a minimum search could be needed, as already written by Jean-Jules.

I hope this clarified things. 

Anyway,I suggest you to  have a look at Gaussian09 User's reference, especially: http://www.gaussian.com/g_tech/g_ur/k_opt.htm which describes very precisely the use of different keywords related to optimization methods. I would also suggest some further readings about optimization techniques and algorithms like in chapter 12 (p380, 2nd edition) of F. Jensen's "Introduction to Computational Chemistry"  or section 4 of chapter 2 (pp 40-46) in C. Cramer's Book: "Essentials of Computational Chemistry".

Best Regards, 

Yohann



-- 

Yohann Moreau

Maître de Conférence, Université Joseph Fourier
iRTSV/CBM/MCT, CEA Grenoble
17 avenue des Martyrs
38 054 Grenoble Cedex 09
Tel. : (33) 4 38 78 29 62
Fax : (33) 4 38 78 54 87

De : "Jean Jules Fifen julesfifen .. gmail.com" <owner-chemistry,ccl.net>
Répondre à : CCL Subscribers <chemistry,ccl.net>
Date : mardi 26 février 2013 08:04
À : MOREAU Yohann 224609 <yohann.moreau,cea.fr>
Objet : CCL:G: calcfc Gaussian 09

What I can say is that,
You can use this option when simple optimizations do not converge. The main drawback of using this keyword is the CPU time consuming.


On 25 February 2013 14:49, Raphael Martinez raphaelmartinez1983###gmail.com<owner-chemistry[*]ccl.net> wrote:

Sent to CCL by: "Raphael  Martinez" [raphaelmartinez1983**gmail.com]
Hello everyone:
I have a very basic/quick question. I am using Gaussian 09.
I know that using the option "calcfc" in the optimization is a bit expensive. So, I have few questions
about it:
1. I was wondering if this option "calcfc" must be always use in the optimizations, or I can avoid it?
2. If I avoid the "calcfc", am I going to obtain a bad optimization?
3. are there certain cases in which is better to use it and cases in which I can avoid it?

Thank you all for the help

Raphael



-= This is automatically added to each message by the mailing script =-
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

Subscribe/Unsubscribe:
      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
      http://www.ccl.net/spammers.txt

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





--
J. Jules Fifen.