From owner-chemistry@ccl.net Wed Jan 2 11:17:00 2013 From: "Gary Breton gbreton:_:berry.edu" To: CCL Subject: CCL:G: Defining Becke's Half and Half hybrid functional Message-Id: <-48037-130102111449-16565-Pkb3HRgqDhgQe/qH7l4W4Q * server.ccl.net> X-Original-From: "Gary Breton" Date: Wed, 2 Jan 2013 11:14:48 -0500 Sent to CCL by: "Gary Breton" [gbreton^^berry.edu] Hello everyone and Happy New Year. I am conducting some studies on some relative simple charge transfer complexes using DFT methods within Gaussian 98 and 03. Several recent publications (e.g., Steinmann, S. N.; Piemontesi, C.; Delachat, A.; Corminbeouf, C. J. Chem. Theory. Comput. 2012, 8, 1629) suggest that Becke's so-called "half and half functional" (Becke, A. D. J. Chem. Phys. 1993, 98 (2), 1372) are best employed. The Gaussian 98 and 03 manuals discuss some hybrid functionals (BHandH and BHandHLPY) but the manual explicitly states that these are NOT the same as the half-and-half functionals proposed by Becke (the manual references the paper above). On the following page: http://www.cup.uni-muenchen.de/ch/compchem/energy/dft1.html I found a clear descriptor of this "half-and-half" method which seems pretty straightforward: Becke-Half-and-Half-LYP (BHandHLYP) uses a 1:1 mixture of DFT and exact exchange energies: EXC = 0.5*EX(HF) + 0.5*EX(B88) + EC(LYP) The same page describes the well-known B3LYP functional as: Becke-3-LYP (B3LYP) uses a different mixing scheme involving three mixing parameters: EXC = 0.2*EX(HF) + 0.8*EX(LSDA) + 0.72*DEX(B88) + 0.81*EC(LYP) + 0.19*EC(VWN) So it is easy to implement the B3LYP functional from within Gaussian since it is already integrated into the Gaussian program. In order for me to use the "half-and-half" model I need to generate a "user defined model". In the Gaussian 09 users manual, user defined models can be created according to the following format (the process is pretty much the same as for the G03 and G98 programs but I could easily copy/paste from the G09 website): P2EX[HF] + P1(P4EX[Slater] + P3Ex[non-local]) + P6EC[local] + P5EC[non-local] the six parameters (P1-P6) may be defined using: IOp(3/76=mmmmmnnnnn) sets P1 to mmmmm/10000 and P2 to nnnnn/10000. P1 is usually set to either 1.0 or 0.0, depending on whether an exchange functional is desired or not, and any scaling is accomplished using P3 and P4. IOp(3/77=mmmmmnnnnn) sets P3 to mmmmm/10000 and P4 to nnnnn/10000 IOp(3/78=mmmmmnnnnn) sets P5 to mmmmm/10000 and P6 to nnnnn/10000 The manual provides an example for defining the B3LYP functional (see above for reference to the parameter values): #P BLYP IOp(3/76=1000002000) IOp(3/77=0720008000) IOp(3/78=0810010000) Finally, now my two questions: First, for B3LYP example provided in the manual why is the last IOp defined as IOp(3/78=0810010000) and not IOp(3/78=0810001900)? Is this just a silly mistake in the manual? (which I believe it is) Secondly, if I wanted to implement Becke's half and half functional according to the definition provided above , would the following be correct? #P BLYP IOp(3/76=1000005000) IOp(3/77=0500000000) IOp(3/78=1000000000) My test on G98 (with appropriate modifications since the IOp assignments are different) was successful (i.e., went to completion). Unfortunately, I really don't have any good tests to make sure I am defining things correctly and want to double check with the much-more knowledgeable members of CCL to make sure I'm not doing something incredibly stupid. Thanks, Gary Breton Berry College From owner-chemistry@ccl.net Wed Jan 2 13:12:00 2013 From: "Frank Jensen frj]*[chem.au.dk" To: CCL Subject: CCL:G: Defining Becke's Half and Half hybrid functional Message-Id: <-48038-130102131008-19920-/sJF2B0kYRkqBwnBe+4mmQ%server.ccl.net> X-Original-From: Frank Jensen Content-Language: en-US Content-Transfer-Encoding: 8bit Content-Type: text/plain; charset="us-ascii" Date: Wed, 2 Jan 2013 18:09:56 +0000 MIME-Version: 1.0 Sent to CCL by: Frank Jensen [frj^^chem.au.dk] The short version: Don't use regular DFT (BLYP, B3LYP, BHHLYP, ...) for charge transfer, they will get both the charge separation and energetics wrong. The range-separated (also called long-range corrected) versions, e.g. CAM-B3LYP, have much better performance. Frank Frank Jensen Assoc. Prof. Dept. of Chemistry Aarhus University http://old.chem.au.dk/~frj -----Original Message----- > From: owner-chemistry+frj==chem.au.dk:-:ccl.net [mailto:owner-chemistry+frj==chem.au.dk:-:ccl.net] On Behalf Of Gary Breton gbreton:_:berry.edu Sent: 2. januar 2013 17:15 To: Frank Jensen Subject: CCL:G: Defining Becke's Half and Half hybrid functional Sent to CCL by: "Gary Breton" [gbreton^^berry.edu] Hello everyone and Happy New Year. I am conducting some studies on some relative simple charge transfer complexes using DFT methods within Gaussian 98 and 03. Several recent publications (e.g., Steinmann, S. N.; Piemontesi, C.; Delachat, A.; Corminbeouf, C. J. Chem. Theory. Comput. 2012, 8, 1629) suggest that Becke's so-called "half and half functional" (Becke, A. D. J. Chem. Phys. 1993, 98 (2), 1372) are best employed. The Gaussian 98 and 03 manuals discuss some hybrid functionals (BHandH and BHandHLPY) but the manual explicitly states that these are NOT the same as the half-and-half functionals proposed by Becke (the manual references the paper above). On the following page: http://www.cup.uni-muenchen.de/ch/compchem/energy/dft1.html I found a clear descriptor of this "half-and-half" method which seems pretty straightforward: Becke-Half-and-Half-LYP (BHandHLYP) uses a 1:1 mixture of DFT and exact exchange energies: EXC = 0.5*EX(HF) + 0.5*EX(B88) + EC(LYP) The same page describes the well-known B3LYP functional as: Becke-3-LYP (B3LYP) uses a different mixing scheme involving three mixing parameters: EXC = 0.2*EX(HF) + 0.8*EX(LSDA) + 0.72*DEX(B88) + 0.81*EC(LYP) + 0.19*EC(VWN) So it is easy to implement the B3LYP functional from within Gaussian since it is already integrated into the Gaussian program. In order for me to use the "half-and-half" model I need to generate a "user defined model". In the Gaussian 09 users manual, user defined models can be created according to the following format (the process is pretty much the same as for the G03 and G98 programs but I could easily copy/paste from the G09 website): P2EX[HF] + P1(P4EX[Slater] + P3Ex[non-local]) + P6EC[local] + P5EC[non-local] the six parameters (P1-P6) may be defined using: IOp(3/76=mmmmmnnnnn) sets P1 to mmmmm/10000 and P2 to nnnnn/10000. P1 is usually set to either 1.0 or 0.0, depending on whether an exchange functional is desired or not, and any scaling is accomplished using P3 and P4. IOp(3/77=mmmmmnnnnn) sets P3 to mmmmm/10000 and P4 to nnnnn/10000 IOp(3/78=mmmmmnnnnn) sets P5 to mmmmm/10000 and P6 to nnnnn/10000 The manual provides an example for defining the B3LYP functional (see above for reference to the parameter values): #P BLYP IOp(3/76=1000002000) IOp(3/77=0720008000) IOp(3/78=0810010000) Finally, now my two questions: First, for B3LYP example provided in the manual why is the last IOp defined as IOp(3/78=0810010000) and not IOp(3/78=0810001900)? Is this just a silly mistake in the manual? (which I believe it is) Secondly, if I wanted to implement Becke's half and half functional according to the definition provided above , would the following be correct? #P BLYP IOp(3/76=1000005000) IOp(3/77=0500000000) IOp(3/78=1000000000) My test on G98 (with appropriate modifications since the IOp assignments are different) was successful (i.e., went to completion). Unfortunately, I really don't have any good tests to make sure I am defining things correctly and want to double check with the much-more knowledgeable members of CCL to make sure I'm not doing something incredibly stupid. Thanks, Gary Breton Berry Collegehttp://www.ccl.net/cgi-bin/ccl/send_ccl_messagehttp://www.ccl.net/chemistry/sub_unsub.shtmlhttp://www.ccl.net/spammers.txt From owner-chemistry@ccl.net Wed Jan 2 13:47:00 2013 From: "nico green nicogreen6*|*gmail.com" To: CCL Subject: CCL: Convergence Message-Id: <-48039-130102134451-1048-+gEolFXpKDC4Q7iroNOiTg.:.server.ccl.net> X-Original-From: nico green Content-Type: multipart/alternative; boundary=f46d0445170dc321ab04d252a38a Date: Wed, 2 Jan 2013 15:44:43 -0300 MIME-Version: 1.0 Sent to CCL by: nico green [nicogreen6() gmail.com] --f46d0445170dc321ab04d252a38a Content-Type: text/plain; charset=ISO-8859-1 Hello to everyone happy 2013 to all of you I'm performing optimizations looking for a TS, but in some cases i found a complex (which is part of the reaction mechanism). I would like to know which keyword do you think is the best in cost/benefit relation to avoid this complex and obtain the TS. Thanks to all of you for being so helpfull --f46d0445170dc321ab04d252a38a Content-Type: text/html; charset=ISO-8859-1
Hello to everyone
happy 2013 to all of you

I'm performing optimizations looking for a TS, but in some cases i found a complex (which is part of the reaction mechanism).
I would like to know which keyword do you think is the best in cost/benefit relation to avoid this complex and obtain the TS.

Thanks to all of you for being so helpfull


--f46d0445170dc321ab04d252a38a--