versión On-line ISSN 0717-9707
J. Chil. Chem. Soc. vol.56 no.2 Concepción 2011
J. Chil. Chem. Soc., 56, N° 2 (2011), págs.: 656-662
THEORETICAL STUDY OF THE ISOMERIZATION OF MALEIC ACID INTO FUMARIC ACID
RICARDO UGARTE, CARLOS BUSTOS, IGNACIO MORENO-VILLOSLADA
Instituto de Química. Facultad de Ciencias. Universidad Austral de Chile. Casilla 567, Valdivia-Chile. e-mail address: firstname.lastname@example.org
The isomerization reaction of 2-butenedioic acid in gas phase has been studied using Hartree-Fock (HF), Meller-Plesset to 2nd order (MP2) and BLYP, B3LYP Density Functional Theory (DFT) levels of approximation. Potential energy surfaces (PES) of the singlet and triplet states of the molecule are constructed in order to corroborate the hypothesis that the reaction proceeds by a non-adiabatic path. A total of ten stationary points (eight singlet and two triplet) were characterized as energy minima or saddle points. Thermochemical analysis of the triplet conformers and singlet most stable structures yields a value of the AH0 of activation for isomerization that is in good agreement with the expected results.
Keywords: Potential energy surface; Isomerization Maleic acid; DFT calculations.
As is well known, substituted ethylenes exhibit cistrans isomerism. For instance maleic acid transforms into fumaric acid at high temperatures (melting point); this transformation is often accompanied by partial descomposition. The isomerization also occurs by heating an aqueous solution of maleic acid in a sealed tube above 130° C, by catalyzing with various acids and salts at ~ 100° C and as a result of bromine atoms.19 The first investigation on the temperature dependence of the isomerization kinetics for the maleic acid in the liquid state was performed by Hojendal.10 He found an activation energy of 15.8 kcal mol1.
Experimental results have supported the classification of cis-trans isomerizations in two groups which involve rotation about ethylenic double bond. They are both first order reactions. The first group have an Arrhenius rate constant of ~ 1011 exp (-45,000/RT), with a high frequency factor. 11 Phenyl-substituted ethylenes (stilbene, methyl cinnamate) isomerizations belong to this group, and it is considered that the molecules remain in a singlet state throughout the isomerization process. The second group, to which maleic acid isomerization belongs, have an Arrhenius rate constant of ~ 104 exp (-25,000/ RT). In this case, it is considered that the isomerization mechanism involves changing from singlet state to triplet state, and then to singlet state; the low frequency factor is a result of the difficulty associated with these changes. In chemical kinetics the word "adiabatic" has come to refer to a process in which there is no change of quantum state. According to the so-called "adiabatic hypothesis", it is supposed that the electrons remain in equilibrium with the particular nuclear configuration which occurs at a given moment, due to the large speed of the electronic motions relative to those of the nuclei. There are no electronic transitions (e.g. multiplicity changes) during the course of the reaction and the whole process takes place on a single potential energy surface. In consequence, the maleic acid isomerization will be referred as a non-adiabatic reaction (singlet → triplet → singlet).
The aim of this paper is to justify the assumption that the interconversion proceeds by a non-adiabatic path. That is, we will show that the possible triplet structures of the transition state for the maleic acid → fumaric acid interconversion have a lower energy than the corresponding singlet structures. In order to find the most stable conformations of maleic and fumaric acids, the potential energy surfaces (PES) at HF/6-31G(d) level for both electronic states, singlet and triplet, are constructed. The most suitable conformers of the aforementioned PES are chosen and geometry optimizations of these structures at different levels of theory are performed. Thermodynamic functions of the resulting isomerization reaction in gas phase are calculated and the value of the energy for the intramolecular hydrogen bond in maleic acid is estimated. As a consequence of the previous study, the potential energy surfaces at BLYP/STO-6G level for both singlet and triplet states are calculated; these potential energy surfaces are superimposed and stationary points and low energy reaction paths are localized. Our theoretical results are compared with experimental results.
All ab initio Hartree-Fock (HF), Meller-Plesset (MP2) and Density Functional Theory (DFT) calculations in vacuo were carried out with the Gaussian 03 program package12 running on a Dell PowerEdge R805 server.
Potential energy surfaces (PES) are constructed by varying the dihedral angles Φ1 and Φ2 (see Figure 1). These angles determine the possible conformations of the molecule. Bear in mind that the dihedral angle Φ1 of the bonded atoms C5-C1-C3-C6 is the angle between the planes C5-C1-C3 and C1-C3-C6. Viewed from the direction of C5, Φ1 is positive for clockwise and negative for anticlockwise rotations. The value Φ1 = 0° corresponds to the planar-cis arrangement of the bonds C5-C1 and C3-C6 (maleic acid), while the value Φ1 = 180° corresponds to the planar-trans arrangement (fumaric acid). The PES of the 2-butenedioic acid were calculated using a 19 x 19 grid generated by rotating through Φ1 and Φ2 in 20° increments from -180 to +180°. At each point a complete geometry optimization was performed with Φ1 and Φ2 frozen at their respective grid values. In order to obtain local minima, unconstrained geometry optimizations (full optimizations) were performed with a set of initial conformations chosen to lie close to the minima found on the PES. Conformational energy map (contour diagram) were drawn in order to facilitate the analysis of the data and to suggest possible reaction paths.
In order to construct the preliminary potential energy surfaces, the Restricted Hartree-Fock (RHF) and Unrestricted Hartree-Fock (UHF) ab initio procedure with the standard 6-31G(d) basis set was employed. These calculations are relatively inexpensive in computational terms and yield reasonable structural data. Vibrational frequency calculations were carried out to confirm the nature of all stationary points: local minima had no negative (imaginary) frequency, saddle points had exactly one negative frequency and second-order saddle points had exactly two negative frequencies. With the purpose of having theoretical results of the isomerization energies as close as possible to the experimental ones, the most relevant HF/6-31G(d) optimized structures obtained from the singlet and triplet PES were used as the first estimate or starting point for calculations at different levels of theory: HF/6-31G(d,p), MP2/6-31G(d,p), DFT-B3LYP/6-31G(d,p),13 DFT-BLYP/6-31G(d,p) and DFT-BLYP/STO-6G.
Following the indications provided by the preliminary calculations (trial and error), the BLYP hybrid functional14,15 with the minimal STO-6G basis set produces results which are closer to the experimental values. In DFT calculations, BLYP can be invoked in Gaussian 03 using the Iop(3/76), Iop(3/77), and Iop(3/78) statements. In this work, we utilized a user-defined model that contains the following keywords: # blyp/sto-6g iop(3/76=1000005800) iop(3/77=0500004200) iop(3/78=0420004000). Thus, we have considered 58% Hartree-Fock and 42% DFT exchange energies.
In order to obtain standard termochemical data (298.15 K, 1 atm), the optimized structures at the different levels of theory were submitted to frequency calculations. All thermodynamic quantities are evaluated under ideal gas and rigid rotor-harmonic oscillator approximation. The dominant error in the above approximation may be the failure of the harmonic approximation for low-frequency vibrations which are rotational in nature. In our system, there are some low-vibrational frequencies; in this case, it is often more appropriate to compute the thermal component associated with these modes using a hindered rotor approximation. In this work, we request to the Gaussian 03 program the identification of internal rotation modes during the harmonic vibrational analysis, in order to correct the thermodynamic functions. However, when the calculation failed under this approximation, no hindered rotor approximation was used. Associated errors are not significant and therefore do not invalidate the analysis of results. The rotational and the translational contributions to the thermal quantities are taken from the classical expressions. The vibrational contributions are evaluated with the unscaled normal mode frequencies. The relative energies account for the zero-point-vibrational energy correction.
As a consequence of the analysis of the obtained results, the potential energy surfaces at BLYP/STO-6G level for both singlet and triplet states were constructed. In order to obtain information about possible reaction paths through a non-adiabatic process, we superimpose both PES and obtain another combined PES containing the lowest value of the energy for each pair of dihedral angles.
3. RESULTS AND DISCUSSION
The HF/6-31G(d) conformational energy map (contour diagram) of 2-butenedioic acid is shown in Figure 2 for both the singlet and triplet electronic states. The molecular structures close to the value Φ1 = 0° represent the maleic acid, while those close to the value Φ1 = 180° represent the fumaric acid. The global minimum at singlet potential energy surface (S-PES) are located in [Φ1, Φ2] = [180°, 180°]; at triplet PES (T-PES) in [Φ1, Φ2] = [± 90°, 180°]. Starting from an analysis of these contour diagrams, full optimizations were performed with a set of initial conformations chosen to lie close to the minima found on the PES. Thus, different local minima (stable conformations) were found.
From the standard enthalpy change of formation (AfH0) in gaseous phase of the acids,16 it can be inferred that the enthalpy difference for the isomerization reaction cis → trans is +0.86 kcal mol-1. This energy value indirectly corroborates the stabilizing effect of the intramolecular hydrogen bonding in maleic acid (see Figure 1). In solid phase, the experimental isomerization enthalpy is approximately of -4.9 kcal mol-1, revealing that the preferred form corresponds to the fumaric acid. The above, is a consequence of the larger number of intermolecular hydrogen bonds involved in the trans isomer with respect to the cis isomer in solid phase.
In order to compare with experimental results, we decided to submit some of the previously obtained HF/6-31G(d) conformers, namely SI, SIV, TII and Sv (see Figure 3), to geometry optimizations at different levels of theory. The justification of the choice of these structures is based on the following: the enthalpy difference of the SI →SV process can be compared to the previously mentioned experimental gas phase value; TII conformer represent a minimum on the T-PES, and could be identified as hypothetical "activated complex" and finally, the energy difference between Si and Siv gives us an estimate of the intramolecular hydrogen bond energy. According to this, a structural and thermochemical analysis at 298.15 K and 1 atm was performed.
Inspection of Table 1 reveals that the calculated BLYP/STO-6G enthalpy change (+0.80 kcal mol-1) for the cis →trans isomerization reaction (SI →SV) approximates to the experimental value in gas phase. All methods of calculation, except the BLYP/STO-6G, overestimate the energy gap between SI and SV, and give negative enthalpy differences, which is qualitatively in disagreement with experiments. Apparently, and without regard of other aspects, the intramolecular hydrogen bonding are badly reproduced by these calculations, in particular by Hartree-Fock method.
The ΔH° difference between SI and TII is a measure of the activation energy of the cis → trans reaction in a non-adiabatic process. The lower value for this difference is obtained with HF method (~ +26 kcal mol-1); BLYP/STO-6G yields a value of ~ +29 kcal mol-1, while BLYP/6-31G(d,p), B3LYP/6-31G(d,p) and MP2/6-31G(d,p) overestimate this difference (~ +37, +51 and +61 kcal mol-1). The BLYP/STO-6G value agrees with the expected results, which are found to be ~ +25 kcal mol-1.11
The ΔH° difference between Si and Siv (Si → Siv) is underestimated by HF method (~ +3.3 kcal mol-1), BLYP/6-31G(d,p) (~ +3.6 kcal mol-1) and MP2/6-31G(d,p) (~ +3.8 kcal mol-1). B3LYP/6-31G(d,p) yields a value of ~ +6.1 kcal mol-1, while with BLYP/STO-6G the ΔH° difference is of ~ +6.5 kcal mol-1, a reasonable value for a strong intramolecular hydrogen bond.17 Thus, the stabilization energy of the hydrogen bond allows maleic acid being as stable as fumaric acid in gas phase.
If we compare the structural parameters of chosen conformers obtained with different approaches (see Table 2), one arrives to the conclusion that the hydrogen bond in Si conformer, and the triplet state is better modeled by using BLYP/STO-6G. For example, in regard to the hydrogen bond, the O9-O12 and O12-H10 distances are smaller, the O9-H10 bond length is larger, and the O9-H10-O12 angle is closer to 180°; in relation to the triplet state, C1=C3 bond length is larger, so the spin-spin repulsion is decreased. On the other hand, the intrinsic repulsion between O(9) and O(12) in Siv is minimized as a result of the BLYP/STO-6G calculation and thus, the O9-O12 distance is smaller. Maleic acid have planar ring structure in the crystalline state, with O9-O12 distance (246 pm) which indicates strong hydrogen bond.18,19 The planar structure of Si, Siv and Sv forms and the right Φ1 dihedral angle in Tii is obtained by the various calculation methods.
Trends in the electronic charge distributions may provide some additional insight into the intramolecular hydrogen bond. Table 3 displays the computed Mulliken charges20,21 of the atoms involved in the hydrogen bond interaction, i.e. H(10), O(9) and O(12). Substantial change are observed in going from HF/6-31G(d) to BLYP/STO-6G, where the absolute values of the charges diminishes. Thus, the electrostatic interaction between the two oxygens is expected to be less repulsive with BLYP/STO-6G. Besides, the decreasing of the charge points up the covalent character of the intramolecular hydrogen bond, which is indicative of its strong nature.
In summary, it is shown that the BLYP/STO-6G level of theory is the most reliable in this study since this theoretical description reproduce reasonably well the experimental results. in consequence, the BLYP/STO-6G PES were calculated for both singlet and triplet electronic states of the molecule, with the purpose to show that the cis trans isomerization evolves through a non-adiabatic path.
The BLYP/STO-6G PES of 2-butenedioic acid is shown in Figure 4 for the singlet (S-PES) and triplet (T-PES) electronic states. Qualitatively, both surfaces show distinctive features. For example, at S-PES large variations in energy take place through the Φ1 dihedral angle, irrespectively of Φ2, while the opposite in regard to the angles is generally true at T-PES (except for [Φ1, Φ2] = [180°, 180°]). Briefly, these facts could be explained as follows: we shall assume, for simplicity, that the state of the molecule can be represented by the two electrons forming the p system. At the singlet state, when Φ1 = 0° or 180° maximum overlap of the two p-orbitals occurs and spin pairing of the p electrons has its lowest energy, so that a pair of electrons occupy a bonding molecular orbital. If Φ1 = 90°, the overlap is not possible, and the two electrons will occupy a non-bonding orbital. At the triplet state, with Φ1 = 0° the electrons occupy different orbitals, one of them a bonding molecular orbital the other an antibonding molecular orbital. If Φ1 = 90°, both electrons will occupy non-bonding orbitals and so, low energy differences are found with respect to the situation at which Φ1 = 0°. Therefore, since these dihedral angles denote extreme conditions, the energy should not vary notoriously with Φ1.
The global minimum at the S-PES are located in [Φ1, Φ2] = [0°, 180°], for T-PES at [Φ1, Φ2] = [180°, 0°]. There is a correspondence between the symmetry of the PES and the symmetry of the molecule.
The energy barriers for the interconversion between maleic acid and fumaric acid following an adiabatic path can be observed in the S-PES (see Figure 4, left). in order to obtain information about possible reaction paths through a non-adiabatic process, we superimpose both potential energy surfaces and find another combined PES containing the lowest value of the energy for each pair of dihedral angles (see Figure 5). Thus, reaction mechanisms can be found starting from a singlet state, going through a triplet state and ending in a singlet state, following the lowest energy paths.
Some important conformers obtained of the BLYP/STO-6G PES wereat BLYP/STO-6G level and resembles to Si in the sense that also display an submitted to a thermochemical analysis (298.15 K and 1 atm). BLYP/STO-6Gintramolecular hidrogen bond (O3 ~ 7°, O9-H10-O7 ~ 173°, O9-O7 ~ 253 pm). Ti results of these and of the previously studied structures, are shown in Table 4 and Tii represent minima on the T-PES, and become saddle points in the STand can be compared to the obtained with different levels of theory. The set ofPES, clearly identified as hypothetical "activated complexes". conformers Si - Siv correspond to the maleic acid, while Sv and Svi correspondin general, by using large basis set, Hartree-Fock method overestimate to fumaric acid. Sllla: [+2.0°, +26.00] and Slllb: [-2.0°, -26.0°] constitute a pairthe enthalpy difference between the maleic and fumaric acid, while MP2 and of mirror images or enantiomers (the same can be said for SIIa /Sllb, Tla /Tlb,B3LYP overestimate the enthalpy difference between the triplet state and the T /T , S /S , S /S ). The S conformers could only be calculatedmost stable cis conformer.
Notoriously, the standard B3LYP/6-31G(d,p) procedure is not appropriate for the study of the isomerization process. For instance, the enthalpy change of S1 → Sv reaction is -2.53 kcal mol-1, which is in clear contradiction with the experimental result of +0.86 kcal mol-1; the entropy change is negative and it is overestimated for the isomerization reaction. As a consequence, S1 appears as more stable than SV, because the ΔG° for Si → Sv reaction is 2.20 kcal mol-1. In addition, S1 and SVI are very similar structures and therefore, the large difference in ΔG° values is not justified. In other respects, in virtue of the difference in activation enthalpy between triplet and singlet states, the adiabatic mechanism would be more probable than the non-adiabatic path.
It can be said that isomerization energies closer to experimental are found with BLYP/STO-6G level. The principal driving force for S1 → Sv isomerization reaction is of entropic nature, and as a consequence of the ΔG° values, both trans conformers are equally stable. The Siv conformer (saddle point at BLYP/ STO-6G) is not stabilized by hydrogen bonding and there should also be more repulsion between the lone pair electrons of the two oxygens than in S1. These two effects are expected to result in a fairly high energy for the SIV in agreement with the present calculations. As the BLYP/STO-6G ΔG° of Ti and Tii is lower than the energy barrier in the S-PES (depicted by Sw and SVm) we can conclude that the non-adiabatic mechanism is more probable than the adiabatic path.
Finally, conformational energy map based on ST-PES of 2-butenedioic acid is shown in Figure 6. In this contour diagram some of the possible reaction paths that interconnect cis and trans isomers are shown. The internal rotations through Φ2 (Φ1 ~ 0°) involves, respectively, energy barriers < 5 kcal mol-1 and < 2 kcal mol-1 for the S1 ↔ SIII and SIIIa ↔ SIIIb interconversions; while through Φ2 (Φ1 = 180°) the rotation barriers are of ~ 5 kcal mol-1. Thus, in this context, conformational equilibria are possible at ambient temperature.
Experimental thermodynamic functions concerning isomerization of maleic acid to fumaric acid are well reproduced by theoretical calculation using a user-defined model at a BLYP / STO-6G level of theory. Stabilization of maleic acid by a strong hydrogen bond is well characterized with this method. Other methods, such as HF/6-31G(d,p), MP2/6-31G(d,p) and B3LYP/6-31G(d,p) overestimate the enthalpy gap between maleic and fumaric acid, and qualitatively disagree with experiments in the sign of the enthalpy difference. The results of this study show the possibility that the cis → trans isomerization reaction in gas phase of the 2-butenedioic acid proceeds by a non-adiabatic path, with an enthalpy difference between the triplet state and the cis-planar most stable conformer of approximately 29 kcal mol-1.
The authors wish to acknowledge the DID-UACH (Grant S-2008-26) and Fondecyt (Grant 1090341) for finantial support.
1.E.M. Terry, L. Eichelberger, J. Am. Chem. Soc. 47, 1402, (1925) [ Links ]
2.C. Horrex, Trans. Faraday. Soc. 33, 570, (1937) [ Links ]
3.K. Nozaki, R. Ogg, J. Am. Chem. Soc. 63, 2583, (1941) [ Links ]
4.D.H. Derbyshire, W.A. Waters, Trans. Faraday Soc. 45, 749, (1949) [ Links ]
5.R. Kiyama, S. Minomura, The Review of Physical Chemistry of Japan. 22, 4, (1954) [ Links ]
6.M. Davies, F.P. Evans, Trans. Faraday Soc. 52, 74, (1956) [ Links ]
7.S. Seltzer, J. Am. Chem. Soc. 83, 1861, (1961) [ Links ]
8.J.S Meek, J. Chem. Educ. 52, 541, (1975) [ Links ]
9.A.J. Castro, S.R. Ellenberger, J.P Sluka, J. Chem. Educ. 60, 521, (1983) [ Links ]
10.K. Hojendahl, J. Phys. Chem. 28, 758, (1924) [ Links ]
11.J.L. Magee, W. Shand Jr., H. Eyring, J. Am. Chem. Soc. 63, 677, (1941) [ Links ]
12.Gaussian 03, Revisión-B.03, M.J. Frisch, G.W. Trucks, H.B. Schlegel, G.E. Scuseria, M.A. Robb, J.R. Cheeseman, J.A. Montgomery, Jr., T. Vreven, K.N. Kudin, J.C. Burant, J.M. Millam, S.S. Iyengar, J. Tomasi, V. Barone, B. Mennucci, M. Cossi, G. Scalmani, N. Rega G.A. Petersson, H.Nakatsuji, M. Hada, M. Ehara, K. Toyota, R. Fukuda, J. Hasegawa, M. Ishida, T. Nakajima Y. Honda, O. Kitao, H. Nakai, M. Klene, X. Li, , J.E. Knox, H.P. Hratchian, J.B. Cross, C. Adamo, J. Jaramillo, R. Gomperts, R.E. Stratmann, O. Yazyev, A.J. Austin, R. Cammi, C. Pomelli, J.W. Ochterski, P.Y. Ayala, K. Morokuma G.A. Voth, P. Salvador, J.J. Dannenberg, V.G. Zakrzewski, S. Dapprich, A.D. Daniels, M.C. Strain, O. Farkas, D.K. Malick, A.D. Rabuck, K. Raghavachari, J.B. Foresman, , J.V. Ortiz, Q. Cui, A.G. Baboul, S. Clifford, J. Cioslowski, B.B. Stefanov, G. Liu, A. Liashenko, P. Piskorz, I. Komaromi, R.L. Martin, D.J. Fox, T. Keith, M.A. Al-Laham, C.Y. Peng, A. Nanayakkara, M. Challacombe, P.M.W. Gill, B. Johnson, W. Chen, M.W. Wong, C. Gonzalez, J.A. Pople, Gaussian: Pittsburgh, PA, 2003. [ Links ]
13.E.M.S. Macoas, R. Fausto, J. Lundell, M. Pettersson, L. Khriachtchev, M. Rasanen, J. Phys. Chem. A. 105, 3922, (2001) [ Links ]
14.A.D. Becke, Phys. Rev. A. 38, 3098, (1988) [ Links ]
15.C. Lee, W. Yang, R.G. Parr, Phys. Rev. B. 37, 785, (1988) [ Links ]
16.David R. Lide, ed., CRC Handbook of Chemistry and Physics, 89th Edition (Internet Version 2009), CRC Press/Taylor and Francis, Boca Raton, FL. [ Links ]
17.P. Lipkowski, A. Koll, A. Karpfen, P. Wolschann, Chemical Physics Letters. 360, 256, (2002) [ Links ]
18.M. Shahat, Acta Cryst. 5, 763, (1952) [ Links ]
19.R.E. Dodd, R.E. Miller, W.F.K. Wynne-Jones, J. Chem. Soc. (Resumed), 2790, (1961) [ Links ]
20.R.S. Mulliken, J. Chem. Phys. 23, 1833, (1955) [ Links ]
21.J.S. Gómez-Jeria, J. Chil. Chem. Soc. 54, 482, (2009) [ Links ]
(Received: April 24, 2010 - Accepted: May 9, 2011).