DOI QR코드

DOI QR Code

Computational Investigation of Isomeric and Conformeric Structures of Methyl Fluoroperoxide and Fluoromethyl Fluoroperoxides

Methyl fluoroperoxide와 fuoromethyl fluoroperoxides의 conformers와 isomers 구조에 대한 이론연구

  • Lee, Kyoung-Min (Department of Chemical Education, Chungbuk National University) ;
  • Sung, Eun-Mo (Department of Chemical Education, Chungbuk National University)
  • 이경민 (충북대학교 화학교육과) ;
  • 성은모 (충북대학교 화학교육과)
  • Received : 2011.05.04
  • Accepted : 2011.05.26
  • Published : 2011.06.20

Abstract

The ab initio calculations for fluoromethyl fluoroperoxides have been carried out using MP2/6-311G(d,p) and B3LYP/6-311++G(d,p) method. The structural optimizations were performed for several isomers and conformers of methyl fluoroperoxide, $CH_3OOF$ and the vibrational frequencies were calculated. The most stable conformer of $CH_3OOF$ is skew form and has fairly short O-O bond distance. The trans and cis conformers have 8-12 kcal/mol higher energies than skew form and the other isomers are very unstable. The structures of $CH_2FOOF$, $CHF_2OOF$ and $CF_3OOF$ are also optimized and vibrational frequencies were calculated. These molecules also have skew forms as the lowest energy conformers. The O-O bond distances are longer and C-O bond distances are shorter than $CH_3OOF$, but the structural parameters are almost independent of the number of fluorine atoms in methyl group.

Methyl hydro-peroxide의 methyl group과 O-H에 F를 치환 시킨 fluoromethyl fluoroperoxide($CH_nF_{3-n}OOF$)의 안정한 conformers와 isomers에 대하여 MP2/6-311G(d,p) 방법과 B3LYP/6-311++G(d,p) 방법으로 최적화 계산을 수행하였고 vibrational frequencies를 구하였다. $CH_3OOF$$CH_2FOOF$, $CHF_2OOF$, $CF_3OOF$ 분자 모두 가장 안정한 형태는 skew 구조로 나타났으며 $CH_3OOH$와 비교할 때 대단히 짧은 O-O 결합길이를 갖는다. Trans, cis conformers의 경우 skew 형태보다 8-12 kcal/mol 정도 높은 에너지를 갖으며 O-O 결합길이는 훨씬 길어진다. $CH_2FOOF$, $CHF_2OOF$, $CF_3OOF$ 분자들은 methyl group의 F에 의한 induction 효과로 짧은 C-O 결합길이를 갖고 O-F 결합길이는 길어지나 치환된 F 원자의 수에 따른 구조변화는 그리 크지 않다.

Keywords

서 론

Peroxide는 O-O 결합의 특이한 성질로 인해 생화학이나 연소화학 등 여러 분야에서 폭 넓게 연구되어 왔다. 특히 hydrogen peroxide와 halogen peroxide는 환경문제의 중요성 때문에 많은 주의를 끌고 있으며 이중 대표적인 것들이 극지방에서의 오존층 파괴, 산성비 등으로 대기오염문제와 관련되어 다루어지고 있다.1-3 Alkyl hydroperoxide는 대기압 하에서는 비교적 안정하나 대류권에서는 침전에 의해 제거 되는 것으로 알려져 있고 이들 구조에 대한 특성이나 열역학적 특성은 연소화학에서 산화되는 경로를 이해하는 데 도움이 된다.4 또한 대기화학이나 생화학 분야의 여러 산화과정에 대한 연구에서도 그 자료는 많은 도움을 준다. 따라서 간단한 alkyl hydroperoxide의 구조적 특성이나 분광학적, 열역학적 특성이 많은 실험연구와 이론연구에 의해 밝혀졌다.5-12 연구가 이루어진 alkyl hydro-peroxide의 alkyl group으로는 methyl, ethyl, isopropyl, tert-butyl, isobutyl sec-butyl13등이 있으며 이들 aliphatic chain의 길이나 가지달린 구조가 결합 에너지나 열역학적 안정성에 어떤 효과를 주는 가에 초점을 맞추고 있다.

그러나 이들 alkyl group에 halogen이 치환된 halogenated alkyl hydro-peroxide의 치환에 의한 유도효과를 조사한 연구는 별로 많지 않다. 이들 중에서 실험적으로 알려진 자료를 이용한 Knyazev et al.14의 연구가 있고 Sun et al.15은 CH3OOH의 methyl group에 단계적으로 염소를 치환시킨 치환체에 대한 이론적 연구를 했으며 El-Taher16는 fluorinated methyl hydro-peroxides에 대해 연구를 하였다.

최근에 halogenated methyl hydro-peroxide, CHnX3-nOOH와 halomethyl haloperoxide isomer들, CHn+1X2-nOOX (X=F, Cl, Br, I)에 대하여 density functional theory와 ab initio electronic structure를 이용한 이론 연구가 보고되었다.17 이 연구에서 hydro-peroxide의 열역학적 안정성과 이에 따른 안정한 구조를 찾았으며 치환된 isomeric haloperoxides의 구조적 변화를 조사하였다. 불소로 치환된 hydro-peroxide의 경우 ΔHf가 감소하는 경향을 보였으나 좀 더 무거운 halogenated compound에서는 이런 경향이 일관되게 나타나지 않았음을 보여주었다. 그러나 methylic group에 치환되는 halogen atom에 따라 상대적 안정성이나 생성열의 크기에 꽤 큰 영향을 주는 것으로 나타났다. 이 연구에서는 각각의 치환체의 isomeric structure에 대한 연구는 다루지 못했고 이러한 연구로는 methyl chloroperoxide의 isomeric structure와 conformeric structure에 대한 연구가 Kosmas 등18,19 에 의해 보고 되었다. 이 연구에서는 CH3OOCl의 구조적 isomer들, CH3OOCl(skew, trans, cis), CH3OClO, CH3ClO2 등에 대하여 구조적 안정성과 isomerization barriers를 밝히고 이들 각각의 conformational isomer에 대한 에너지를 비교하였다. 이와 비슷한 연구로 methyl bromoperoxide의 isomers에 대한 양자역학적 계산 연구20가 있으며 methyl iodoperoxide의 isomers 연구21도 Kosmas등에 의해 수행되었다. 이 연구에서는 CH3OOI의 isomers에 대한 안정성 연구뿐 아니라 ICH2OOH와 CH2IOOH 그리고 CH3OOI의 에너지를 비교하여 가장 안정한 형태가 ICH2OOH로 CH3OOI보다 25.7 kcal/mol나 더 낮다고 밝혔다.

Methyl chloroperoxide에서도 이와 같은 구조적 안정성을 보여 주었는데 ClCH2OOH와 CH3OOCl의 경우 ClCH2OOH가 훨씬 짧은 C-O 결합길이를 갖으며 CH3OOCl보다 낮은 에너지 상태로 나타났다. 이것은 halomethyl isomer에서 halogenated methylic group과 -OOH group과의 상호작용이 치환되지 않은 methyl group과 -OOX와의 상호작용에 비해 훨씬 강하게 작용하고 있음을 의미한다고 보고하였다.18

본 연구에서는 fluoromethyl fluoroperoxide에 대하여 그 isomers의 구조와 conformational structures를 알아보고 fluorine으로 치환된 효과가 그 구조에 어떤 영향을 주는지 fluorine에 의한 induction 효과가 peroxide O-O bond의 안정성과 torsional frequencies에 어떤 영향을 주는지 살펴보고자 한다.

 

계산 및 결과

모든 계산은 Gaussian 03 program22으로 수행하였으며 MP2/6-311g(d,p) 방법과 B3LYP/6-311++g(d,p) 방법으로 계산하였다.

Methyl fluoroperoxide(CH3OOF)

가장 안정한 isomer로 알려진 CH3OOF에 대하여 구조의 최적화를 시킨 결과 C-O와 O-F가 skew구조를 갖는 conformer가 가장 안정하게 나타났다. 이 C-O와 O-F의 이면각을 변화시키며 최적화시켜 또 다른 conformer인 trans와 cis 형태를 얻었다. 이 두 구조의 conformer는 그리 안정한 형태는 아니나 methyl chloroperoxide나 methyl iodoperoxide 에서 가능한 conformer로 구조 연구를 하였던 형태이다. 역시 이 계산에서도 이 두 구조는 skew구조 보다 꽤 높은 에너지를 갖는 것으로 나타났다. CH3OOF외에 가능한 isomer로 보는 CH3OFO와 CH3FO2의 구조를 최적화시켜 Table 1과 2에 함께 나타내었다. CH3OFO와 CH3FO2 역시 가능한 isomer의 형태로 보나 에너지는 CH3OOF보다 비교적 높다. 이들 isomer와 conformer들을 Fig. 1에 나타내었다.

먼저 CH3OOF conformer들의 구조를 보면 O-O 결합길이의 경우 안정한 skew conformer에서 가장 짧고 이보다 불안정한 trans와 cis에서 훨씬 긴 것을 볼 수 있다. 이 skew conformer의 O-O 결합길이는 대단히 짧아 분자 O2의 O-O 결합길이인 1.208Å에 근접하는 길이인 1.356Å으로 CH3OOCl의 1.405Å이나 CH3OOBr의 1.437Å, CH3OOI의 1.416Å보다 매우 짧은 것으로 나타났다. CH3OOH에서도 O-O 결합길이는 이보다 훨씬 긴 길이로 보고되었다.17 반면 C-O 결합길이는 CH3OOF가 1.432Å이고 CH3OOCl이 1.421Å CH3OOBr이 1.429Å CH3OOI가 1.419Å으로 오히려 CH3OOF가 제일 긴 것을 알 수 있다. O-F결합길이도 C-O결합길이와 같은 경향을 보인다. 이러한 경향은 CH3OOCl에서도 확인된 결과이다.19 B3LYP 계산에서도 이러한 효과는 나타나고 있으며 그 차이가 더 크게 나타남을 볼 수 있다. 이들 결합길이 외에 CH3 group의 구조에는 별로 영향을 주지않음을 알 수 있다. CH3OOF의 conformers와 isomers에서도 O-O 결합길이는 분자 전체의 안정화 에너지와 관련이 있는 것으로 나타났다. 안정한 conformer의 O-O 결합길이가 짧고 에너지가 높아질수록 O-O 결합길이는 길어지고 C-O 결합길이가 짧아진다. CH3OFO isomer의 경우 에너지는 CH3FO2 isomer보다 낮지만 최적화 과정에서 C-O, O-F, F-O' 결합길이의 조정이 어려웠고 O-F 결합길이가 예상보다 길게 나타났다. CH3OOF의 O-F 결합길이보다 0.05-0.10Å 정도가 길게 나타났다. 이렇게 O-F나 F-O' 길이에 민감하게 에너지가 변한다는 것은 이 isomer 자체가 별로 안정하지 못하다는 것을 보여주는 것 같다.

Table 1.The Optimized structural parameters of isomeric and conformeric structures of CH3OOF at MP2/6-311G(d,p) level of theory

Table 2.The optimized structural parameters of isomeric and conformeric structures of CH3OOF at B3LYP/6-311++G(d,p) level of theory

Fig. 1.The stable conformers and isomers of CH3OOF.

이들 conformers와 isomers의 vibrational frequency를 계산하여 Table 3에 나타내었다. 이 값들은 MP2/6-311g(d,p) 계산에 의한 것으로 frequency가 B3LYP/6-311++g(d,p) 계산에 의한 값보다 약간씩 더 높다. Table 4에 이들 conformers와 isomers의 에너지를 비교하였다. 이 에너지는 zero-point 에너지를 포함시킨 값으로 가장 안정한 conformer인 CH3OOF (skew)에 대한 상대적인 ΔE 값을 나타내었다. Isomer 들의 에너지가 생각보다 많이 높아 안정한 isomer로 존재하기 힘들 것 같다. CH3OOCl이나 CH3OOBr, CH3OOI 분자들의 경우 isomer들과의 에너지 차이가 훨씬 작아 10-60 kcal/mol 정도로 나타났는데 CH3OOF isomers가 이렇게 높은 것은 CH3OFO’나 CH3FO2가 존재하기 힘든 불안정한 상태인 것으로 보인다. CH3OOF conformers의 에너지 차이는 다른 methyl haloperoxide 분자들과 비슷하게 나타났고 trans form이 cis form보다 낮게 나타난 것도 일치한다. F 원자가 Cl, Br, I 원자보다 작아서 cis form이 조금은 유리할 것으로 생각되나 실제로는 거의 차이가 없다.

Table 3.aNumbers in parentheses are IR intensities.

Table 4.The relative energies (Kcal·mol-1) including ZPE corrections for isomers and conformers of CH3OOF

Fluoromethyl fluoroperoxide(CH2FOOF, CHF2OOF, CF3OOF)

Methyl group에 불소가 치환된 분자, CH2FOOF, CHF2OOF, CF3OOF의 구조에 대한 최적화 계산을 하여 MP2 계산 결과를 Table 5에 나타내었다. 전체적으로 MP2 계산과 B3LYP 계산이 비슷한 경향성을 보이며 서로 일관성 있는 결과를 나타내었다. 우선 F 원자가 치환되면서 induction effect에 의해 C-F 결합이 강하게 극성화 됨으로서 결과적으로 C-O 결합거리가 줄어든 것을 알 수 있다. 그러나 CH3 group에 F 원자가 2개 3개로 늘어나는데 따른 변화는 크지 않다. 이 C-O 결합길이의 변화는 methyl hydroperoxide, CH3OOH의 경우 CH3 group에 할로겐 원소가 치환될 때 훨씬 크게 나타나는 것으로 보고되었다.15,16 흥미로운 것은 O-O 결합길이의 변화로 C-O 결합길이와는 다른 경향성을 보인다. O-O 결합길이는 대단히 짧은 길이를 가진 CH3OOF보다 일반적으로 증가한다. 그러나 다른 CH3OOH와 같은 분자에 비하면 아직도 꽤 짧은 길이로 나타났고 이 경우도 F 원자수의 증가에 따른 변화는 별로 없다. 그러나 O-F 결합거리는 CH3 group의 F 원자 치환에 따라 줄어든 것으로 나타났으며 치환된 F 원자 수가 증가함에 따라 약간씩 더 줄어듬을 볼 수 있다.

Table 5.Optimized structural parameters for conformers of CH2FOOF, CH2FOOF and CF3OOF at MP2/6-311G(d,p) level of Theory

그 외의 다른 구조는 비슷한 경향을 보이며 그 중 C-F 결합길이는 F 원자의 치환에 따라 F 원자수가 늘어날수록 결합길이가 현저히 줄어든다. C-O-O-F 이면각은 가장 안정한 skew conformer에서 CH3OOF를 기준으로 CH3 group에 F 원자가 한 개씩 더 치환될수록 늘어남을 볼 수 있는데 F 원자가 2개에서 3개로 늘어날 때는 별 변화가 없다. 이것은 Kosmas 등17의 계산에서 나타난 경향성과 일치한다. 일반적으로 MP2 계산에 의한 것보다 B3LYP 계산에 의한 값들이 더 변화의 폭이 커 보이나 C-O-O-F 이면각의 경우 두 계산 결과가 거의 일치한다. 이 계산 결과로 나타난 세 분자들의 구조 중에서 가장 안정한 skew conformers를 Fig. 2에 나타내었다. 이 그림에서 볼 수 있는 것처럼 안정한 conformation에서 O-O의 O 원자가 CH3 group의 치환된 F 원자와 gauche 형태를 갖는 것을 볼 수 있다. 이것은 이면각 F-C-O-O에서도 볼 수 있는데 이 conformation에서 산소와 F 원자의 비결합 전자쌍 사이의 반발력은 최소화되고 분자 내 상호작용이 가능하게 되는 것으로 보인다.

Table 6.Harmonic vibrational frequencies (cm-1) of CH2FOOF, CHF2OOF and CF3OOF at MP2/6-311G(d,p) level of theory

Fig. 2.The stable conformers of fluoromethyl fluoroperoxide.

Table 7.The relative energies(Kcal·mol-1) including ZPE corrections for isomers of CHnF3-nOOF.

Table 6에는 이 세 분자들의 vibrational frequencies를 MP2 방법에 의해 계산하여 나타내었다. 그리고 Table 7에는 이 세 분자들의 conformers에 대하여 안정한 skew 형태에 대한 상대적인 에너지를 나타내었다. Trans 형태가 cis 형태보다 안정한 것은 CH3OOF의 경우와 같고 methyl group에 치환된 F 원자의 수가 증가할수록 상대적 에너지는 조금씩 더 높아지는 것으로 나타났다. Methyl haloperoxide의 에너지를 isomeric hydro-peroxide의 에너지와 비교해 보면 할로겐 원소 중에서 fluoroperoxide가 가장 불안정한 것으로 나타난다. Methyl group에 F 원자가 치환된 CH2FOOH의 경우 이 F 원자에 의해 이 분자는 대단히 안정화 되는 것을 볼 수 있으나 이 F 원자가 O-H의 H원자에 치환된 CH3OOF는 매우 불안정하게 되는 것을 볼 수 있다. 이 두 isomer의 에너지 차이는 69.7 kcal/mol이나 된다.17 이것은 불안정한 O-F 결합에서 기인하는 것으로 보인다.

 

결 론

Fluoromethyl fluoroperoxide의 안정한 구조에 대한 계산 결과 methyl group에 F 원자가 치환되거나 hydro-peroxide, O-O-H의 H 대신 F가 치환되었을 때 전체 분자구조에 큰 영향을 주는 것으로 나타났다. CH3OOF의 isomer와 conformers에 대한 계산에서 skew 형태의 conformer가 가장 안정한 것으로 나타났고 trans, cis conformer는 이보다 높은 에너지를 갖는 것으로 나타났다. 가장 안정한 skew conformer와 trans, cis conformer를 비교했을 때 가장 차이가 나는 것은 O-O 결합길이와 C-O 결합길이로 그 중에서도 O-O 결합길이는 거의 0.1Å 정도나 skew conformer가 짧다. Skew form의 CH3OOF는 다른 methyl haloperoxide나 halomethyl haloperoxide중에서도 가장 짧은 O-O 결합길이를 갖는다. 이 결합길이는 분자의 안정화에 기여하게 되고 따라서 trans나 cis conformer에서는 이 O-O 결합길이가 크게 늘어남을 볼 수 있다. 반면에 O-F 결합은 F 원자의 큰 전기 음성도에 의해 약해지고 결합길이가 조금 길어진다. CH3OOF의 isomer CH3OFO나 CH3FO2는 이에 비해 대단히 불안정하고 에너지차이가 크다. Mehtyl chloroperoxide나 methyl iodoperoxide의 경우는 그 isomer들과의 에너지 차이가 그리 크지 않아 10-50 kcal/mol 정도의 차이를 보이며, 안정한 skew 구조에서도 O-O 결합길이가 CH3OOF의 O-O 결합길이보다 훨씬 길어짐을 볼 수 있다.19,21

Fluoromethyl fluoroperoxides의 경우 methyl group에 치환된 F의 induction 효과로 O-O 결합길이는 조금 길어지고 C-O와 O-F 결합길이는 짧아지는 것으로 나타났으며 가장 안정한 conformer는 C-O-O-F 이면각이 CH3OOF와 비슷한 skew 형태를 갖는 것으로 나타났다. 흥미로운 것은 methyl group에 치환된 F 원자의 수가 증가함에 따라 C-F 결합길이가 조금씩 짧아지고 C-O-O-F 이면각이 약간 증가할 뿐 구조의 변화가 크지 않다는 것이다. 다른 cis, trans conformer들과의 상대적 에너지 차이도 거의 보이지 않았다.

그러나 중요한 것은 구조의 변화에서 이들 결합길이가 길어지고 짧아지는 것이 그대로 분자의 안정화 에너지에 연결되지 않는다는 것이다. 분자 내에서의 각 원자들의 상호 작용은 여러 결합길이의 변화라든지 구조의 변화에 그대로 직접적으로 비례해서 변화하기보다는 좀 더 복잡한 방식으로 분자들의 특성에 영향을 준다는 것이다. 예를들면 CH2FOOH에서 F 원자를 O-H의 H 위치로 보내면, CH3OOF 분자의 O-O 결합길이는 대단히 짧아지는데 ΔHf0는 CH2FOOH가 CH3OOF보다 훨씬 큰 음의 값을 갖는다고 보고되었다.15,17

이 논문에서는 ΔHf0와 같은 직접적으로 비교할 수 있는 값이 없어 직접적인 비교는 어렵고 methyl group에 F가 치환되었을 때 분명히 더 안정화 되리라 예측할 수 있으나 구조에서 나타난 변화, O-O 결합길이나 C-O 결합길이의 변화를 분자의 안정화 에너지와 연관 지어 설명하기 어렵다.

References

  1. Lee, M.; Heikes, B. G.; O’ Sullivon, D. W. Atmos. Environ. 2000, 34, 3475. https://doi.org/10.1016/S1352-2310(99)00432-X
  2. Vogel, B.; Feng, W.; Streibel, M.; Müller, R. Atmos. Chem. Phys. 2006, 6, 3099. https://doi.org/10.5194/acp-6-3099-2006
  3. Messer, B.; Stielstra, D. E.; Cappa, C. D.; Scholtens, K. W.; Elrod, M. J. J. Mass. Spectrom. 2000, 197, 219. https://doi.org/10.1016/S1387-3806(99)00260-2
  4. Tyndall, G. S.; Cox, R. A.; Granier, C.; Lesclaux, R.; Moortggat, G. K.; Pilling, M. J.; Ravishankara, A. P.; Wallington, T. J. J. Geophys. Res. 2001, 106, 12157. https://doi.org/10.1029/2000JD900746
  5. Schneider, W. F.; Wallington, T. J. J. Phys. Chem. 1993, 97, 12783. https://doi.org/10.1021/j100151a025
  6. Catoire, V.; Lesclaux, R.; Schneider, W. F.; Wallington, T. J. J. Phys. Chem. 1996, 100, 14356. https://doi.org/10.1021/jp960572z
  7. Wallington, T. J.; Hurley, M. D.; Schneider, W. F. Chem. Phys. Lett. 1996, 251, 164. https://doi.org/10.1016/0009-2614(96)00080-2
  8. Jungkamp, T. P. W.; Seinfeld, J. H. Chem. Phys. Lett. 1996, 257, 15. https://doi.org/10.1016/0009-2614(96)00520-9
  9. Lay, T. H.; Krasnoperov, L. N.; Venanzi, C. A.; Bozzelli, J. W. J. Phys. Chem. A, 1996, 100, 8240. https://doi.org/10.1021/jp952976h
  10. Lay, T. H.; Bozzelli, J. W. J. Phys. Chem. A, 1997, 101, 9505 https://doi.org/10.1021/jp972103i
  11. Sumathi, R.; Green Jr., W. H.; Phys. Chem. Chem. Phys. 2003, 5, 3402. https://doi.org/10.1039/b307050f
  12. Matthews, J.; Sinha, A.; Francisco, J. S. J. Chem. Phys. 2005, 122, 221101. https://doi.org/10.1063/1.1928228
  13. Simmie, J. M.; Black, G.; Curran, H. J.; Hinde, J. P. Phys. Chem. A, 2008, 112, 5010. https://doi.org/10.1021/jp711360z
  14. Knyazev, V. D.; Slagle, I. R. J. Phys. Chem. A, 1998, 102, 1770. https://doi.org/10.1021/jp9726091
  15. Sun, H.; Chen, C.-J.; Bozzelli, J. W., J. Phys. Chem. A, 2000, 104, 8270. https://doi.org/10.1021/jp0013917
  16. El-Taher, S. J. Fluor. Chem. 2006, 127, 54. https://doi.org/10.1016/j.jfluchem.2005.08.018
  17. Kosmas, A. M.; Mpellos, C.; Salta, Z.; Drougas, E. Chem. Phys. 2010, 371, 36. https://doi.org/10.1016/j.chemphys.2010.03.026
  18. Drougas, E.; Kosmas, A. M. Chem. Phys. letters, 2003, 369, 269. https://doi.org/10.1016/S0009-2614(02)01981-4
  19. Drougas, E.; Kosmas, A. M. Int. J. Quant. Chem. 2004, 98, 335. https://doi.org/10.1002/qua.20012
  20. Papayannis, D. K.; Melissas, V. S.; Kosmas, A. M. Phys. Chem. Chem. Phys, 2003, 5, 2976. https://doi.org/10.1039/b302666n
  21. Drougas, E.; Kosmas, A. M. Can. J. Chem. 2005, 83, 9. https://doi.org/10.1139/v04-156
  22. Frisch, M. J.; Trucks, G. W.; Schegel, H. B.; Scuseria, G.E.; Robb, M. A.;Cheeseman, J. R.; Motgomery, J. A., Jr;Verven, T.; Kudin, K. N.; Burant. C.; Millam, J. M.; Lyengar,S. S.; Tomasi, J.; Barone, V.; Mennucci, B.; Cossi,M.; Scalmani, G.; Rega, N.; Petersson, G. A.; Nakatsuji,H.; Hada, M.; Ehara, M.; Toyota, K.; Fukuda, R.; Hasegawa,J.; Ishida, M.; Nakajima, T.; Honda, Y.; Kitao, O.;Nakai, H.; Klene, M.; Li, X.; Knox, J. E.; Hratchian, H.P.; Cross, J. B.; Bakken, V.; Adamo, C.; Jaramillo, J.;Gomperts, R.; Stratmann, R. E.; Yazyev, O.; Austin, A. J.;Cammi, R.; Pomelli, C.; Ochterski, J.; Ayala, P. Y.; Morokuma,K.; Voth, G. A.; Salvador, P.; Dannenberg J. J.;Zakzewski, A. G.; Dapprich, S.; Daniels, A. D.; Strain, M.C.; Farkas, O.; Malick, D.; Raghavachari, K.; Foresman,J. B.; Ortiz, J. V.; Cui, Q.; Baboul, A. G.; Clifford, S.;Cioslowski, J.; Stefanov, B. B.; Liu, G.; Liashenko, A.;Piskorz, P.; Komaromi, I.; Martin, R. L.; Fox. D. J.; Keith,T.; Al-Laham, M. A.; Peng, C. Y.; Nanayakkara, A.;Challacombe, M.; Gill, P. M. W.; Johnson, B. G.; Chen,W.; Wong, M. W.; Gonzalez, C.; Pople, J. A.; Gaussian 03, revision B.03; Gaussian, Inc.; Pittsburgh PA, 2003.