Quantum Mechanical Studies on the Mode of Hydrogen Bonding, Stability and Molecular Conformations of Phenol-Water Complexes ()
1. Introduction
Non-covalent and hydrogen bonding interactions are very important and essential for modern chemistry, material sciences, molecular biology, etc. These interactions are quite common in aromatic-pi systems and these kinds of forces can exist as pi-pi stacking, pi-cations, pi-anions, etc. Among these non-covalent interactions, the hydrogen bonding interactions are the most important in different molecular system as they play a very crucial role in the formation of fundamental units of bulky molecules [1] [2] [3] [4] . It is very common that phenol and diphenol molecule shows hydrogen bonding with water molecule through its functional group (−OH). But, interestingly the mode of hydrogen bonding between these systems is also a very important factor. Therefore, the formation of intermolecular hydrogen bonding with other chemical species or solvent molecules plays a crucial role in its molecular properties such as bond lengths, bond angles and electron density on atoms which further influence on the chemical reactivity of phenolic systems [5] [6] . The determination and quantification of such types of hydrogen bonding or non-covalent interaction are highly desirable in order to predict and understand the probable molecular interactions or concerned reaction mechanism of molecular association [7] . Moreover, the non-covalent interactions are often low in energy (1 - 5 kcal/mol) as compare to the covalent interactions (80 - 110 kcal/mol) [8] .
Nowadays, people have been investigating to identify the hydrogen bond motifs, one of the main aims is to expand the available list of hydrogen bond donor and acceptor in the periodic table for a wider range of application [9] [10] [11] [12] [13] . A number of molecular complexes have extensively been studied by scientists exhibiting D-H---A interactions, where D-H is the “H-bond” Donor and A is the Acceptor both experimentally and computationally [14] [15] [16] [17] . Computational calculations also show that there exists a significant interaction between a hydrogen bond donor and the centre of an aromatic ring, which acts as a hydrogen bond acceptor. This interaction is as strong as a normal hydrogen bond, which contributes approximately 3 kcal/mol (1 cal = 4.184 J) of stabilizing energy and it plays a significant role in molecular associations [18] . Studies on hydrogen bonding interaction with common aromatic systems are well known, for example, H-bonding between electron donor and aromatic ring acceptor was studied by Ray L. Brinkley et al. [19] . The interaction energy of intermolecular phenol-water complex is in the range of about 15 - 28 kcal/mol, which is about three times higher compared to its analogue i.e. O-H---Se [20] . Even O-H---S hydrogen bond interaction energy in phenol-SH2 complex is higher than that of O-H---Se in Phenol-SeH2 complex. The stabilization energy of O-H---O hydrogen bond was estimated to be similar to its analogue i.e. O-H---Se and O-H---S [21] . It is well known for classical Hydrogen bond interaction that the introduction of a positive charge into the proton donor or negative charge on the acceptor enhanced the stabilization energy [22] . Luminol system was also well studied for their extended hydrogen bonding network in water solution [23] . In some cases, it has been observed that the hydrogen bonding interaction can also affected by infrared spectra at high pressure and temperature, e.g. π-hydrogen bonding interaction between water and aromatic hydrocarbons (like benzene, toluene, ethylbenzene, xylene, and mesitylene) [24] . In this current work, we have studied hydrogen bonding interaction between simple aromatic phenol and diphenol rings with water for different phenol-water conformations and also study the stability of such conformers by using computational method. We also try to study more detailed investigation of this mentioned phenomenon where we could investigate on the minute details like different path of interaction, orientation of the molecules by calculating single point energy, Mullikan charge and Gibbs free energy by MPn methods.
2. Methodology
In our investigations all the models were constructed by using Gaussian 5.0. We have studied the hydrogen bonding interactions between simple phenolic system and water molecules by computational methods. Several phenol-water complexes were constructed in such a way that they have different orientation with respect to each other and optimized these complexes in gas phase at room temperature. The molecular structures were optimized by using Møller-Plesset perturbation methodology i.e. MP2 level of theory and double zeta basis set 6-311G/++(d,p) at 298 K and 1atmospheric pressure. After optimization it gives a constant intermolecular hydrogen bond distances and then interaction energies were calculated by changing the intermolecular hydrogen bond distances for every phenol-water complexes. We have calculated the single point energy by MP4 level of theories with basis set 6-311++G(d,p). For our studies we also calculated Gibbs free energy, Mullikan charges etc. for each phenol-water complexes [25] .
The interaction energies (EIE) are calculated by the following equation;
EIE = EPhenol-water − EPhenol − E water
Here, EPhenol-water = Energy for phenol-water complex, EPhenol = Energy of Phenol, E water = Energy of water.
3. Result and Discussions
In this current study, we have prepared different phenol-water models to analyze the actual hydrogen bonding interactions, stability and effect in phenol-water complexes. Phenol can directly interact with water in many ways, but we have taken only the most prominent and chemically favoured interactions between the water and phenol molecules. As we know, phenol undergoes strong hydrogen bonding interactions with its neighbouring water molecule when it is in aqueous medium and this is already experimentally established. Phenol interacts with water molecule through the −OH group of phenol, depending on the electron denoting and accepting abilities −OH group of phenol can interact with water either through H atom (donor) or through O atom (acceptor). It can also show weak π-dipole interaction i.e., π-electrons of phenol interacts with dipoles of water and many other mode of interactions. To observe the proper phenol-H2O complexes, we have constructed four different chemically favoured phenol-H2O complexes (Figure 1(a) & Figure 1(b)). These proposed constructed phenol-H2O complexes are universally accepted and shows better interaction between phenol and water molecule. Our investigation is based on these four proposed models and to determine the proper phenol-H2O conformation with
![]()
Figure 1. (a) Proposed models for phenol-water complexes; (b) Proposed models for phenol-water complexes.
appropriate phenol-H2O hydrogen bonding and dipole-π interactions. We also carried out computational calculations to determine the stability of these phenol-H2O complexes. For all four proposed phenol-H2O complexes, the optimised geometry of model-1, 2 and 3 shows strong hydrogen bonding and gives a constant hydrogen bond distance between phenol and water molecule (Figures 2(a)(i)-(iii)). On the other hand, only model-4 shows dipole-π interaction, which is a very weak force of interaction (Figure 2(b)(iv)). It has been observed that the calculated hydrogen bond distances between the four proposed models of phenol-water complexes are found to be 1.8 and 2.0 Å respectively (Table 1). However, the crystal structure of phenol-water complexes reveal that the actual phenol-water hydrogen bond distance for any linear model is found to be ~1.77 Å [26] , whereas for our proposed models the optimum hydrogen bond distance for model 1, 2 and 3, lies within a range from 1.8 - 2.0 Å (Table 1). Again, we are trying to change the hydrogen bond distances of four proposed phenol-H2O optimized models within a range from 1.2 to 3.0 Å, to calculate the actual phenol-water interaction energies, by this way we can easily predict the strength of hydrogen bond, stability, bond dissociation energies as well as the chemical behaviour of these models. Generally, if we plot interaction energies (kcal/mol) for
(a)
(b)
Figure 2. (a) Minimized models with hydrogen bonding for phenol-water complexes; (b) Minimized models with hydrogen bonding for diphenol-water complexes.
![]()
Table 1. Computed hydrogen bond distances for minimized models of phenol/diphenol-water complexes (MP4 method).
each model with respect to the change in hydrogen bond distances then we can get an exponentially decreasing curve. The minimum of the curve results the more negative interaction energy value which gives the most favoured point with a stable phenol-water complex (Figure 3(a)). During the process of increasing and decreasing in hydrogen bond lengths for each model, it has been observed that on increasing the hydrogen bond length the interaction energy becomes steady and they behave as two individual molecules rather than a phenol-water complex. Moreover, on decreasing the hydrogen bond distances at a certain point the interaction between phenol-water complexes become stronger, which gives the more negative interaction energy value and it may also form covalent bond if they get closer than the minimum intermolecular distances i.e. 1.2 Å. Since, we have to determine the minimized interaction energy between phenol-water complexes; therefore we must take only one minimized point with more negative interaction energy value. All the minimized structures of proposed phenol-water complexes are shown in the Figures 2(a)(i)-(iv). To check the effect of solvation on the selected phenol-water complexes, we put the model in aqueous medium and we also calculated the interaction energy values of phenol-H2O complexes, and also compare it with gas phase (Figure 3(b)). Here, we can easily distinguish the change in interaction in gas phase and solvent phase i.e. in aqueous phase. The minimized interaction energies for phenol-water complexes in gas phase gives more negative value than that of the aqueous phase (Table 2). It happens because in gas phase only one type of hydrogen bonding interaction exists in the individual phenol-water complexes. But, in aqueous phase we have observed two different types of hydrogen bonding interactions viz. 1) between the individual phenol-water complexes and 2) phenol-water complexes with water molecules of the aqueous phase. In other words, it shows that the polarity of the solvent also plays an important role in the stability of the phenol-water complexes, i.e. in aqueous phase water can enhance the stability of the phenol-water complexes (Table 2).
(a)
(b)
Figure 3. (a) Interaction energy plot for phenol-water complexes in gas phase; (b) Interaction energy plot for Phenol-water complexes in aqueous phase.
![]()
Table 2. Computed interaction energies for minimized models of Phenol-water complexes (MP4 method).
Model-1, where the hydrogen bond occurs through −OH group of phenol with water molecule (−OH-----OH2), gives the most negative interaction energy value (more negative values) for both gas and solvent phase. It gives the most stable and stronger hydrogen bond, where the hydrogen bonds distance is found as 1.8 Å. The calculated ΔGº for the minimized geometry for all four models reveals that the model-1 gives the more negative value, i.e. −138.6503 kcal/mol (Table 3). Therefore, model-1 is the most stable and favoured conformation among the other phenol-H2O complexes. On the other hand, the study of phenol-water complexes for model-2 and 3 are also very important, as they show similar kind of hydrogen bonding interaction (HO-----H2O), with hydrogen bond distance ~2.0 Å, but they have different conformations with respect to the water molecule as shown in Figure 2(a)(ii) & Figure 2(a)(iii). Phenol-water complexes for model-2 and 3 gives almost similar hydrogen bonding interaction energies for the minimized models in gas phase, which are found as −8.3149 and −8.4207 kcal/mol respectively (Table 2). Irrespective of the above studies, in phenol-H2O model-4, there exists very weak dipole-π interaction at a separation of 2.4Å, as the water molecule is placed just above the centre of π-electrons of phenol ring as shown in Figure 2(a)(iv), the minimized interaction energy is found as −7.7345 kcal/mol. Although, it may also form hydrogen bond but that will be very weak as the −OH group of phenol and H2O placed far apart from each other. For all the phenol-water complexes, the hydrogen bonding interaction also shows similar trend in both gas and solvent phase. Again, the chemical behaviour of all the phenol-water complexes show negative ΔHº values, therefore all the systems are exothermic in nature.
Mulliken charge density calculation is also a very important tool to investigate the change in individual charge densities of the molecule before and after interaction. Greater is the change in charge densities (Mulliken charges) of the interacting individual atoms which are responsible for hydrogen bonding in phenol-water complexes, more will be the interaction in between phenol and water in that complex. The actual change in the electron charge densities of the minimized phenol-water complexes were obtained by calculating the Mulliken charges of the complex which shows whether there exist strong or weak interactions between phenol and water molecule (Table 4). In our calculations, we have obtained that the model-1 shows great change in the mulliken charge density value for −OH group of phenol and H2O interacting sites than that of the individual molecules. On the other hand model-2 and 3 shows weak hydrogen bonding interaction as there observed a small change in the mulliken charge
![]()
Table 3. Computed free energy change (ΔGº) and enthalpy change (ΔHº) for minimized models of phenol-water complexes (MP2 method).
![]()
Table 4. Mulliken charge calculation for phenol-water complexes (MP4 method).
density values for the same interacting sites i.e. −OH group of phenol and H2O molecule (Figure 2(a)(ii) & Figure 2(a)(iii)). In case of model-4, which shows a minute change in the Mulliken charge densities as it does not show any hydrogen bonding interaction but the existing interaction shows due to the dipole-π interaction (Table 4).
We also carried out similar types of investigation for diphenol-H2O complexes, prepared the four models for diphenol-water complexes as constructed for phenol-H2O, and also calculated the interaction energy values (Figure 1(b) and Table 5). Interestingly, the change the interaction energies and Mulliken change densities for phenol-H2O system and diphenol-H2O systems give almost similar trend (Table 5 & Table 6). The minimized models for diphenol-water complexes are shown in Figure 2(b). The minima of the interaction energy plots for all four diphenol-water systems, both in gas and aqueous phase shows almost similar types of interactions as it has shown in phenol-water systems (Figure 4(a) & Figure 4(b)). But, these kinds of interactions might be changed on changing the substituent on phenolic systems.
4. Conclusion
From the above calculations, we can conclude that in model-1 hydrogen bonding is taking place through −OH-----OH2, and it shows the strongest interactions in phenol-water system with hydrogen bond distance 1.8 Å, in both gas and water phase. It also shows the most negative interaction energy values i.e. −11.1783 kcal/mol, which reveals the most stable and favoured conformation. Model 2 and 3 also show almost similar types of hydrogen bonding but the hydrogen bond is formed through HO-----H2O, with a hydrogen bond distance 2.0 Å. On
(a)
(b)
Figure 4. (a) Interaction energy plot for diphenol-water complexes in gas phase; (b) Interaction energy plot for diphenol-water complexes in aqueous phase.
![]()
Table 5. Interaction energies for minimized models of diphenol-water complexes (MP4 method).
the other hand, model-4 favourably shows dipole-π interaction as there exists direct interaction between π-electron of phenol and dipole of water molecule at a separation of 2.4 Å, it gives the least negative interaction energy than the other phenol-water models. Although, we have done calculations for both gas and solvent phase systems, it has been observed that the hydrogen bond distances remain
![]()
Table 6. Mulliken charge calculation for diphenol-water complexes (MP4 method).
constant for both the phases, but the solvent phase models show more stability and negative interaction energy values, because of the polarity of the water molecule. It is also found that the diphenol-water systems show the same trends as it has shown by phenol-water system.
Acknowledgements
Authors are highly grateful to the AICTE-TEQIP-3 fund and Ministry of Human Resource Development (MHRD), New Delhi, for providing research assistance.