Research Article  Open Access
Ana M. Amado, Sónia M. Fiuza, Luis A. E. Batista de Carvalho, Paulo J. A. RibeiroClaro, "On the Relevance of Considering the Intermolecular Interactions on the Prediction of the Vibrational Spectra of Isopropylamine", Journal of Chemistry, vol. 2013, Article ID 682514, 12 pages, 2013. https://doi.org/10.1155/2013/682514
On the Relevance of Considering the Intermolecular Interactions on the Prediction of the Vibrational Spectra of Isopropylamine
Abstract
The effects of implicitly considering the effects of hydrogen bonding on the molecular properties, such as vibrational frequencies, were inferred on the basis of DFT calculations. Several clusters of isopropylamine were assembled and theoretically characterized. The results showed that maximum Hbond cooperativity is achieved when the amine group acts simultaneously as donor and acceptor. The effect of Hbond cooperativity manifests itself in the relative cluster stability and on the structural and vibrational frequency predictions. Referring to the vibrational frequencies it was found that the stretching and torsion vibrational modes are the most affected by the amine involvement in hydrogen bonding. Both stretching modes were found to be significantly redshifted relative to the monomer. The torsional mode, on the other hand, was found to be blueshifted up to 350 . Finally, the comparative study between the theory levels performed allows to conclude that the small 631G* basis set is able to stabilize weak interactions as long as the new dispersion corrected DFT methods are considered. The impairments observed with conventional DFT methods for describing weak interactions may be overcome with the improvement of basis set, but the associated increase of computational costs may turn the calculations unfeasible.
1. Introduction
Intermolecular interactions play an important role in many areas of chemistry and biology, being important in phenomena such as solvation, molecular and macromolecular structure, supramolecular chemistry, and molecular recognition [1–5]. In condensed phase samples (liquid or solid), two types of intermolecular interactions may be of particular relevance, namely, the effects of the interactions between closelying molecular neighbours (short range intermolecular interactions) and the effects of the surrounding medium on the molecular aggregates formed (mediumtolong range intermolecular interactions). The occurrence of these interactions may have a pronounced effect on the experimental spectra, such as vibrational ones, turning their assignment a difficult task.
Quantum chemical calculations are widely recognized as a powerful and helpful methodology for predicting and assigning experimental vibrational spectra. However, most of the theoretical calculations are performed considering the isolated molecule. By neglecting the effect of potential intermolecular interactions, significant deviations from the experimental data are frequently observed. This is particularly important when amine systems are targeted as they are prone to form extensive intermolecular hydrogen bond networks, with each NH_{2} group being able to establish up to three intermolecular interactions with neighbouring groups—two as a donor and one as an acceptor. Moreover the occurrence of multiple hydrogen bonds involving the same amine group opens up the possibility of hydrogen bond cooperativity [6–11].
In this context, the main goal of the present work is to infer the relevance of considering these intermolecular effects on the theoretical predictions of the vibrational spectra of aliphatic amines. Isopropylamine (2aminopropane, iPram) was selected as a case study since, as a smallsized molecular system, it allows the study of different clusters, of increasing dimension and complexity, using different theory levels at a permissive computational cost. In the present work, molecular clusters of up to four iPram units, held by N–HN and/or C–HN intermolecular interactions, were considered.
2. Theoretical Calculations Description
Calculations were performed on a PC computer, using the Gaussian 09W (G09W) program [12]. The allelectron basis set 631 was used in all calculations. The authors are aware of the smallness of the basis set used. Nevertheless, its use is justified as the present work is part of a larger project that aims to study larger aliphatic amines and their palladium and platinum complexes.
Six DFT methods were implemented in G09 under the designations PW91PW, B3LYP, mPW1PW, B971, B97D, and wB97XD. The use of these distinct DFT methods allows to infer on the effect of using either pure (PW91PW) or hybrid (B971, mPW1PW, and B3LYP) density functionals and on considering (or not) the dispersion and longrange effects (B97D and wB97XD). Several studies have shown that the conventional DFT methods, either pure or hybrid, tend to underestimate the interaction energies of hydrogenbonded systems [7, 10, 13–17]. The inclusion of semiempirical corrections for dispersion (as in B97D) coupled with longrange correction (as in wB97XD) has been suggested as viable routes for achieving a better description of the intermolecular interactions, including the interaction energies, without compromising the computational efficiency [13, 15, 17, 18]. Testing hybrid B971 DFT seemed relevant as it is a modification of the hybrid B97 method, the parent DFT of both dispersion corrected DFTs considered [19]. Some additional calculations were performed on specific objectives and, therefore, are presented in the text where appropriate.
As stated previously, the present work is aimed at gauging the cumulative effect of different hydrogen bonds involving a given molecular group (in the present case, a primary amine NH_{2} group) on the theoretical predictions of molecular properties, such as vibrational frequencies. Three main effects were targeted: (i) the relative stabilization relevance of the distinct Hbond donor groups (that is, the relative stabilizing effect of N–HN and C–HN interactions), (ii) the relative preference of the amine group to act as a donor (d) or as an acceptor (a), and finally, (iii) the improvements achieved by fulfilling the hydrogen bonding capabilities of the amine group. In view of these goals, different cluster geometries were assembled on a basis of a progressive increase of the amino group involvement in hydrogen bonds (monomer → dimer → trimer → tetramer). The different cluster geometries studied are presented in Figure 1, along with the correspondent designations used throughout the work and type of intermolecular interactions present. All geometries were fully optimized without any symmetry constraints. All attempts made to assemble clusters including –HN interactions failed as all converged to one of the clusters presented in Figure 1.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
Lines of evidence from the literature show that counterpoise correction (CP) during the gradient optimisation of molecular complexes, such as molecular selfassociated clusters, has a pronounced effect on the calculated results [10, 14, 17]. Thus, to account for the basis set superposition error (BSSE) all cluster geometries were reoptimized within the scheme of Boys and Bernardi [20] as implemented in G09w [12] under the keyword counterpoise.
In all cases, vibrational frequency calculations were performed, at the same theory level, to verify its correspondence to a real minimum (no negative eigenvalues) and to estimate both the zeropoint vibrational energy (zpve) and thermal energy corrections (all energy values considered in the present work are corrected values for both zpve and thermal contributions). The calculated frequency values have not been scaled by any factor, as relative raw effects are aimed for discussion.
Oligomers Overall Stability Evaluation. The relative stabilizing effect of the hydrogen bonding pattern of each oligomer (cluster) was evaluated by comparing two energetic parameters, namely, the cluster formation energy () and the cooperative effect between multiple Hbonds (CEHB). The value of each cluster () was determined as where represents the number of iPram units forming the cluster and represents the energy determined for the monomer. The corresponding CEHB values, on the other hand, were determined using the expression where coefficients and stand, respectively, for the number of N–HN and C–HN interactions present in the cluster.
3. Results and Discussion
3.1. Energetic Analysis
Figure 2 shows the variation of the CPuncorrected values as a function of oligomer size and type of intermolecular interactions—(A) oligomers presenting only N–HN interactions (group I) and (B) oligomers also containing C–HN contacts (group II). The variation of the CPcorrected values is presented in Figure S1 of the Supplementary Material (see Supplementary Material available online at http://dx.doi.org/10.1155/2013/682514).
(a)
(b)
The first evidence that stands from Figure 2 and is the difficulty encountered by the conventional DFT methods in describing the C–HN intermolecular contact. In the case of the CPuncorrected results (Figure 2) this impairment is only observed with B971, which is unable to stabilize the O3c cluster. All attempts made to stabilize the O3c geometry failed with convergence to O3a. The same impairment was encountered when using the MP2 formalism (results not shown). When the CPcorrected calculations are gathered (Figure ), it is found that the inability to stabilize O3c cluster is extended to the remaining conventional DFTs. Moreover, MP2 and all conventional DFT methods, but B971, also become unable to stabilize the dimeric form O2b.
The impairment of stabilizing iPram clusters held by C–HN intermolecular contact affects more strongly the B3LYP calculations. In fact, performing CP correction turns the B3LYP/631 calculations unable to stabilize either of the group II clusters. This impairment of B3LYP calculations has been correlated with the underestimation of the dispersion interactions by this DFT method [21]. In this context, and only to infer on the importance of the basis set size for this B3LYP impairment, calculations were performed at the B3LYP/6311++ theory level (with and without CP correction). The results showed that the improvement given to the basis set allows to overcome the limitations of the calculations B3LYP to stabilize the C–HN intermolecular contact. The relative ordering of the different adducts as a function of the magnitude of the is the same as encountered with the B97D/631 and wB97XD/631 calculations. In fact, the only clear difference encountered between those calculations refers to the magnitude of the values, being the B97D/631 and wB97XD/631 values significantly overestimated relatively to the B3LYP/6311++ results. On the other hand, using the larger 6311++ basis set leads to a dramatic increase of the computational costs, which rapidly turn the calculations prohibitive as amine size and complexity increase.
With no surprise, the magnitude of the interaction energies increases with cluster size and number of N–H⋯N interactions established. In the case of the group II oligomers (Figures 2(b) and S1(b)), the interaction energy enhancement observed is almost linear with respect to cluster size. Regarding group I oligomers (Figures 2(a) and (a)), however, it is found that the magnitude of the depends not only on cluster size and number of N–HN interactions but also on the overall Hbond pattern assembled. It is found that a cumulative da involvement of the amine group as in O3a is more stabilizing than the dd involvement that characterizes O3b by 7–12 kJ mol^{−1} (5–12 kJ mol^{−1} in the case of the CPcorrected results).
Regarding the effect of performing BSSE correction, it is found that the values of group II clusters are generally more affected than the ones related to group I oligomers. The conventional DFT methods are definitely more affected by BSSE than the dispersion corrected methods. This differential effect is probably related to the widely described deficiency of the former methods to cover dispersion energies [14, 16, 22].
An important issue in hydrogen bonding is the manifestation (or not) of Hbond cooperativity as the number of hydrogen bonds increases and the overall spatial arrangement varies. The concept is for long explored in the context of ligandtoreceptor binding but, more recently, has also been applied to the understanding of the behaviour of noncovalently bonded systems, such as molecular clusters [10, 11, 23]. Within the context of hydrogen bonding, positive cooperativity has been defined as the strength enhancement of existent intermolecular Hbond(s) as a new molecular unit is added to the growing aggregate [11]. Conversely, negative cooperativity refers to the weakening of the preexistent Hbond(s) due to cluster growth.
One way of evaluating hydrogen bond cooperativity is by determining the overall stability enhancement achieved when a new intermolecular interaction is added to the preexistent Hbond pattern. This can be done by analysing the variation of the CEHB values with cluster size. Note that by definition (2), a positive CEHB value corresponds to a positive cooperativity.
Figure 3 plots the CPuncorrected and CPcorrected CEHB values as a function of DFT and cluster type. Recall that when BSSE correction was performed, the conventional functionals showed some difficulties in stabilising group II clusters, disabling some comparisons regarding the CP correction of the CEHB values obtained with the conventional DFT. In the cases where comparisons are possible, it is found that, in general, CP correction of the CEHB values, while affecting their absolute magnitude, does not affect their sign. The only exceptions are observed for the B97D and wB97XD CEHB values of O3b, which turn from negative to positive by CP correction, and for the B3LYP CEHB of O4a, which shows the opposite sign switch.
(a)
(b)
Based on the results obtained it is found that the enhancement of the interaction energies achieved by cluster size increase is not always interrelated with the occurrence of Hbond cooperativity. According to the results of Figures 2 and , increase of cluster size from O2a to either O3a or O3b is accomplished by a significant enhancement of . However, the data of Figure 3 clearly show that only the transformation O2aO3a is coupled with a meaningful positive Hbond cooperativity. In the case of the O2aO3b transformation the predicted CEHB are either negative or just slightly positive. This observation suggests that changing from a da Hbond profile (O3a) to a dd pattern (O3b) leads to a decrease of Hbond cooperativity. Reinforcing this observation is the significant Hbond cooperativity predicted for the O3c cluster which, in spite of having a weak C–HN instead of a strong N–HN, presents a high CEHB value. These results are clearly in line with the conclusions of Huyskens [24] for the involvement of OH groups in multiple hydrogen bonds. According to those results, the involvement of an OH group as an Hbond donor strengthens its aptitude to act as an acceptor in a second Hbond but weakens its ability to act as a donor in the second Hbond.
Increasing cluster size from trimer to tetramer by fulfilling the amine group Hbonding capacities (O4a and O4b) increases positive Hbond cooperativity. When the transformations O3aO4a or O3cO4b are regarded the CEHB values are increased, while following either of the transformations O3bO4a or O3bO4b leads to a sign switch of the CEHB from negative to positive.
Just aiming to infer on the maximum CEHB enhancement possible by incrementing cluster size, calculations of the CEHB values were performed for a forth trimericcluster (O3d; Figure 4). This iPram cluster simulates the trimeric water cluster predicted as the most stable [25, 26]. As seen in Figure 4, the overall Hbonding pattern around each iPram unit is the same as in O3a (da). However, in this trimeric cluster each iPram unit is linked to the other two through N–HN interaction, in one as an acceptor and in the other as a donor.
(a)
(b)
The graphical representation shown in Figure 4 compares the CPuncorrected and CPcorrected CEHB values predicted for the O3d cluster with the ones obtained for both O3a and O4a clusters (only the B97D results are presented). As seen, the CEHB predicted for cluster O3d, with and without CP correction, largely overwhelm the corresponding CEHB values predicted for the equally sized O3a cluster. The magnitude of the CEHB values determined for O3d is very close to the ones associated with the largersized O4a cluster. This result, which was found to be independent of the theory level considered (including MP2), may reinforce the view that Hbond cooperativity is mainly correlated with the presence of a da Hbond arrangement around the amine group rather than to a fulfilling of its Hbonding capacities.
3.2. Structural Analysis
Hydrogen bonding and cooperativity effects can manifest in other properties than energetic ones. From a structural point of view, the changes promoted on the length of the covalent bonds, directly involved in the interaction, and on the nonbonded HX distance may be particularly informative. Some studies report a close correlation between those structural parameters and the energetic of the intermolecular interactions as cooperative effects are concerned [10, 11].
Table 1 presents the CPcorrected bond length variations predicted for the –H and –H covalent bonds, relative to the monomer, and for the intermolecular distances between iPram units as a function of cluster type (Figure 1). The corresponding CPuncorrected values are presented in Table of the Supplementary Material.

Comparison of the CPcorrected and CPuncorrected results (Table 1 versus Table ) shows that correction of geometries for BSSE generally leads to a shortening of the X–H distances (X standing for either or ) and to a lengthening of the HN intermolecular distances. This observation suggests that BSSE tends to overestimate the strength of intermolecular interactions, being in previously reported results on other hydrogen bonded molecular clusters [14, 17].
CP correction is found to affect distinctly the different DFT methods. The conventional DFT results are more strongly affected by CP correction than the dispersioncorrected results (B97D and wB97XD). Among the conventional DFT methods, B3LYP and mPW1PW seem to be the most affected ones. The high sensibility of B3LYP to CP correction translates in the abovementioned inability to stabilize group II clusters. The mPW1PW DFT, on the other hand, despite giving evidence of being able to stabilize O4b aggregate, predicted an HN distance (ca. 341 pm) that is too long to be considered a real intermolecular distance.
Analysis of the variation trends of the N–H and C–H covalent bonds evidences the differential behaviour of the two structural parameters to the direct involvement in intermolecular interactions (Tables 1 and ). Considering the monomer as reference, it is seen that while the establishment of N–HN interaction promotes the lengthening of the N–H covalent bond, the formation of C–HN interactions leads to a shortening of the C–H covalent bond. This distinct behaviour of the N–H and C–H covalent bonds has already been observed, both experimentally and theoretically, for other hydrogen bonded molecular clusters [11, 27–33], and has been correlated with a balance between two main effects that act in opposite directions. One is the net electron density gain in the X–H bond region, that promotes a contraction of the covalent X–H bond, and the other is the attractive interaction between the positively charged hydrogen atom and the electron rich acceptor Natom, which leads to lengthening of that bond [28, 29, 34]. The distinct behaviour observed for the N–H and C–H covalent bonds is explained by a distinct balance between those two effects, while in the N–HN closecontacts the preponderance of the HN attraction explains the N–H lengthening relative to the monomer, and in the case of the C–HN interaction, an electron density gain in the C–H bond region overwhelms the effects of the attractive HN interaction, leading to a relative C–H bond shortening.
Opposing to what is observed for the C–H bond distance, the lengthening magnitude predicted for the N–H bonds is clearly dependent on the Hbond pattern assembled. The establishment of the da Hbond pattern by incrementing O2a into O3a promotes a significant lengthening of the N–H covalent bond, involved in the preexistent intermolecular interaction, relative to the monomer. A similar effect, despite being more subtle, is observed on passing from O2a to O3c, which also presents the da Hbond pattern but being the newly added interaction much weaker than in the former transformation. On the contrary, when the cluster incrementing is associated with the assembling of a dd Hbond pattern (O2a → O3b) it is found that the N–H lengthening relative to the monomer is significantly reduced (N–H distance in the trimer is shortened relative to the dimer). The effect of maximising the NH_{2}group involvement in hydrogen bonds on the N–H bond length (O3a → O4a and O3b → O4b) is the balance between the effects due to the simultaneous presence of the two Hbond patterns, da and dd, and the type of intermolecular interaction being added, N–HN or C–HN. On the whole, the lengthening effects of the da pattern overwhelms the strengthening effects of the da pattern (ΔN–H in O4a is smaller than in O3a but larger than in O3b) and adding a N–HN interaction to the pattern has a more significant lengthening effect than adding a C–HN interaction (ΔN–H in O4a is more significant than in O4b).
Referring to the variation of the intermolecular distances (HN), and assuming the cluster interconversion processes O2a →O3a → O4a and O2a →O3b → O4b, it is found that the general trends oppose the ones observed for the variation of the N–H covalent bond. The lengthening of the N–H bonds observed in the formation of O3a and O4b is accomplished by a decrease of the corresponding HN distances. In turn, the strengthening of the N–H bonds observed in the formation of O3b and O4a is coupled with a shortening of the corresponding distances. A similar straightforward correlation is, unfortunately, not observed when the C–HN intermolecular interactions are analysed. In these cases, the direction of the variations is dependent on the type of DFT used and if BSSE correction is performed. However, assuming a higher accuracy of the CPcorrected B97D and wB97XD results, it is expectable that incrementing the involvement of the NH_{2}group in hydrogen bonds (O2b → O3c → O4b) promotes a progressive strengthening of the C–HN intermolecular interaction, as judged by the HN distance reduction pattern predicted.
3.3. Vibrational Analysis
The assignment of the vibrational spectra of pure liquid iPram was recently reviewed on the basis of quantum chemical calculations performed for the isolated molecule (monomer, M) [35], using different DFT approaches. As stated there, some of the predicted vibrational modes, namely, those related to the NH_{2}– and –H groups, were found to be significantly shifted relative to the corresponding experimental value. This observation naturally raised the question to what extent the involvement of those molecular groups in intermolecular interactions, expected to occur in the condensed phase, is responsible for these deviations.
Table 2 presents the CPcorrected deviations, relative to the monomer, predicted for the fundamental vibrational modes ascribed to the NH_{2} and –H groups of iPram in the different clusters considered. The CPuncorrected and CPcorrected values obtained when using all the DFT methods are compared in Tables S2–S7 of the Supplementary Material. The present analysis is limited to the nine vibrational modes related to the NH_{2} and –H groups as they are expectedly the most affected ones, due to their direct involvement in the intermolecular interactions. In the case of the dimeric species, the frequency shifts calculated for both iPram units (A and either B or C) are presented. Referring to the higher species, only the shifts observed for unit A are targeted as gauging the effects of Hbond cooperativity on the predicted vibrational frequencies is aimed.
 
= stretching, = deformation, t = twisting, = wagging, = torsion, = inplane deformation and = outofplane deformation. 
The results show that neither CP correction nor changing DFT (Tables S2–S7) affects the general trends observed for the vibrational frequency shifts promoted by the intermolecular interactions. The most significant deviations are definitely observed for the amine related vibrational modes, when the NH_{2}group participates in one of the Hbonds as a donor. When the amine group acts only as an acceptor (unit B in O2a and both units, A and C, in O2b), the vibrational mode deviations promoted are almost negligible.
Among the six vibrational modes related to the NH_{2}group, the two stretching modes (symmetric, NH_{2}, and antisymmetric, NH_{2}) and the torsional mode (NH_{2}) are by far the most affected ones (Table 2). The three NH_{2} deformation modes (scissoring, NH_{2}, twisting, tNH_{2}, and wagging, NH_{2}) denote blueshiftings that are comparatively very slight. Both NH_{2} and NH_{2} modes are significantly redshifted upon establishment of the N–HN interactions, in line with results reported by others [6, 8, 10]. The largest redshift of the NH_{2} mode is observed when the amine group acts simultaneously as donor and acceptor in the overall Hbonding pattern (O3a, O3c, and O4a). In what concerns the NH_{2} mode, the maximum shifting is observed when the amine group Hbonding capacity is achieved. In both cases, the magnitude of the shifts are attenuated when the acceptor character of the amine group is given in a C–HN interaction.
Referring to the effect of hydrogen bonding cooperativity, the results suggest that the effects promoted on the two amine stretching modes oppose each other. In the case of the NH_{2} mode, Hbond cooperativity tends to enhance the redshifting degree (e.g., O2aO3a and O2aO3c), an effect that is opposed by the increase of the amine group involvement as an Hbond donor (e.g., O2aO3b, O3aO4a, and O3cO4b). The variation trends observed for the NH_{2} mode are, in every way, opposite to those predicted for the NH_{2} mode.
Analysis of the effects of amine involvement in hydrogen bonding on the prediction of the NH_{2} vibrational mode requires a more profound discussion. In the prior publication [35], difficulties were encountered in ascribing the τNH_{2} mode. Durig et al. [36] ascribed the NH_{2} mode of iPram, in the transform, to a weak spectral feature observed around 267 cm^{−1} in the infrared gasphase spectrum. This assignment is supported by the calculations performed for the monomeric form of iPram, which predicted the NH_{2} mode at about 300 cm^{−1} (Tables S2–S7). However, neither Durig et al. [36] nor the more recent work of Padrào et al. [35] was able to ascribe this mode to a particular spectral feature of either infrared or Raman spectra of liquid iPram.
The present predictions for the NH_{2} mode of iPram suggest that the Hbond involvement of the amine group generally leads to a significant blueshifting of the NH_{2} mode. The only requirement is that the amine group participates as a donor in at least one of the intermolecular interactions it is involved in. Incrementing the involvement of the NH_{2}group in N–HN interactions (e.g., O2aO3a, O2aO3b, and O3aO4a) enhances the magnitude of the NH_{2} blueshifting. Similar effects have been reported for other aminecontaining systems such as polyamines. Marques et al. [37] found that the involvement of the amine group in hydrogen bonding shifts the NH_{2} mode upwards by more than 100 cm^{−1}. Even more pronounced NH_{2} shiftings were predicted by Amado et al. [38] in the case of solid ethylenediamine dihydrochloride, where strong N–HCl interactions promoted deviations up to 300 cm^{−1}.
These observations naturally highlight the need to review the assignment of τNH_{2} mode in the spectrum of condensed iPram. Thus, taking into account the magnitude of the shifts predicted to occur upon the establishment of the N–HX interactions, the τNH_{2} mode is most likely shifted to the 450–590 cm^{−1} spectral region. Being so, the τNH_{2} mode becomes overlapped with the inplane and outofplane –NH_{2} deformation modes (–NH_{2} and –NH_{2}, resp.), which are predicted to occur in the 400–500 cm^{−1} spectral region (results not shown).
Finally, regarding the effect of hydrogen bonding on the three vibrational modes related to the –H group, it is found that only the stretching mode (–H) seems to be significantly affected. The direction of the –H shifts is found to be dependent on the kind of Hbond involvement of the amine group in the same iPram unit. In the cases where the amine group acts as an Hbond donor (e.g., unit A in O2a, O3b, and O4b) the –H mode is predicted to be redshifted. But when the amine group is involved as an acceptor (e.g., unit B in O2a and unit C in O2b), the –H denotes a considerable blueshifting. This later effect (C–H blueshifting) has been observed, both theoretically and experimentally, for several other molecular systems, and has attracted great interest over the years [39–46].
The C–H stretching bands ascribed to carbonhydrogen bonds which, as the –H of iPram, are anticoplanar to nitrogen lone pairs have been named as Bohlmann bands [47]. The Bohlmann effect has been correlated with a delocalization of the nitrogen lone pair, Lp(N), onto the αhydrogen bond lying in a relative transorientation. This effect leads to a weakening of the C–H bond and, hence, to a lowering (redshift) of the C–H mode. This effect has been widely reported in the literature [31–33, 47] and explains, at least in part, the lowering of 115 cm^{−1} predicted for the –H mode of iPram on changing from the gauche to the transform [35].
The present results show that formation of molecular clusters of iPram may enhance or attenuate the Bohlmann effect, depending on the kind and number of intermolecular interactions holding up the aggregate. Regarding the dimeric forms, it is found that either the direct Hbond involvement of the –H group (as a donor; unit C in O2b) or the involvement of the amine group as an Hbond acceptor (unit B in O2a and unit A in O2b) leads to an attenuation of the Bohlmann effect. In other words, the –H mode is blueshifted. On the contrary, in the cases where the amine group acts as the Hbond donor (unit A in O2a), the Bohlmann effect is enhanced and the –H mode is redshifted. The overall effect in the higher cluster forms may be seen as a balance of the individual effects observed for dimeric forms. Increasing the donor involvement of the amine group in the Hbonding pattern (e.g., unit A in O2aO2aO4a) clearly enhances the Bohlmann effect with the –H mode being progressively more redshifted.
3.4. Evaluation of the PiMM Methodology
Aiming at evaluating the effect of the shortrange interactions on the theoretical predictions of the vibrational spectra of condensed materials, some of the authors developed the Pairs in Molecular Materials (PiMM) methodology [48]. This approach assumes that the effects of the individual intermolecular interactions, in which a given molecular entity is involved in, are pairwise additive. Application of the method has been used with success in the assignment of the vibrational spectra of several solid pharmaceutical drugs known to present polymorphic forms [48–51]. However, that accuracy evaluation was conditioned by the fact that, according to the Xray data available for those drugs, none of the systems presented a molecular group involved in multiple hydrogen bonds. In view of fulfilling this gap, the present section intends to gauge the errors associated with neglecting the effect of cooperativity between multiple Hbonds centered at the same molecular fragment, in the present case the NH_{2} group of iPram.
Table 3 compares the shifts predicted for the NH_{2} and –H related vibrational modes when applying either the PiMM or the fullcluster approaches to the B97D CPcorrected results. The comparison presented in Table 3 is restricted to the two clusters found to present the highest stability and CEHB, namely, O3a and O4a. Listing of the frequency shifts predicted for the PiMM approach for all clusters (CPuncorrected and CPcorrected) is found in Table S8 of the Supplementary Material. As in Table 2, only the deviations predicted for unit A in each cluster are presented.
 
= stretching, = deformation, t = twisting, = wagging, = torsion, = inplane deformation and = outofplane deformation. Fullcluster approach. 
Recall that, within the PiMM approach, the effects of the different intermolecular interactions on the vibrational modes of a given molecular group are considered independent and additive [48]. Thus, considering cluster O3a, the tabulated deviations correspond to the sum of the deviations predicted for units A and B of cluster O2a. In the case of the O4b cluster, on the other hand, the values are the sum of twice the deviations observed for unit A of O2a and the deviations in unit A of O2b.
Comparison of the deviations predicted for the different vibrational modes by either applying the PiMM methodology or the full cluster geometry optimization method (Table S8 versus Table 2) shows that, in general, the shifting directions are the same. The only exception that is worth to be mentioned refers to δNH_{2} vibrational mode, in all clusters excepting O3b.
The same consistency is, however, not observed with respect to the magnitude of the deviations envisaged. In some cases, the deviations provided by the PiMM method are much higher than those obtained when considering the overall optimized molecular cluster. The deviations envisaged by PiMM for the NH_{2} and NH_{2} vibrational modes in O3b, O4a, and O4b, for example, are well beyond those provided by the method of full cluster optimized. The opposing effect is, however, observed when the NH_{2} and NH_{2} modes in O3a and O4a are regarded. In these cases, the PiMM deviations are significantly underestimated relative to those provided by the full cluster optimization method.
On the whole, these observations suggest that the PiMM method is a quite effective methodology for predicting whose vibrational modes are overdeviated due to the establishment of intermolecular interactions. It is also effective in predicting the direction of the deviations (blue or redshifting). However, as expected, the effectiveness of the method in predicting the magnitude of the deviations is strongly affected by neglecting the potential effects of cooperativity among multiple hydrogen bonds.
4. Conclusions
In the present work, calculations, using different DFT methods, were performed to infer on the effects of intermolecular interactions on the prediction of several molecular properties of isopropylamine. From a practicable point of view, B97D/631 and wB97XD/631 theory levels were found to present the best balance between the computational costs and accuracy. As opposed to the conventional DFT methods, both dispersion corrected DFTs were able to stabilize the weak C–HN despite the smallness of the basis set. The conventional DFT methods were only able to describe that type of intermolecular interactions when a significant improvement of the basis set is performed, which involves a drastic increase of the computational requirements.
With no surprise, the results showed that the overall cluster stability increases with oligomer size. However, the stability enhancement observed is not always coupled with cooperativity effects between multiple hydrogen bonds. The energetic and structural evaluation performed showed that achievement of maximal hydrogen bond cooperativity requires a da hydrogen bonding involvement of the amine group. Involvement of the NH_{2}group in a dd Hbond pattern tends to break Hbond cooperativity.
Regarding the prediction of the vibrational frequencies of iPram, it was found that establishment of intermolecular interactions affects particularly the amine related modes. One of the most affected modes is the NH_{2}. This mode is significantly blueshifted upon the establishment of intermolecular interactions that involve the amine group as an Hbond donor. As a consequence, the τNH_{2} mode of liquid iPram turns overlapped with the with the –NH_{2} and –NH_{2} modes.
The Bohlmann band (–H) of iPram was also found to be sensitive to the hydrogen bonding pattern characterizing the molecular cluster. The results clearly showed that both direct involvement of the –H group and amine group involvement as an acceptor in the hydrogen bonding pattern tend to an attenuation of the Bohlmann effect (i.e., promote a blueshift of the –H mode). The involvement of the amine group as an Hbond donor, on the other hand, tends to enhance the Bohlmann effect (i.e., promote a redshift of the –H mode). The overall effect observed on the –H mode as the iPram hydrogen bonding involvement enhances corresponds to a balance between those two lines of evidence taking into account the type and number of intermolecular interactions holding the cluster geometry.
Finally, for the evaluation of the effects of cooperativity on the reliability of the PiMM method [48] for predicting the vibrational frequencies, it was found that the methodology gives good evidences on the sign of the frequency deviations upon establishment of the intermolecular interactions. However, the magnitudes of the individual shifts are not always well predicted by the PiMM method, particularly when the overall Hbonding pattern allows strong cooperativity effects between multiple hydrogen bonds. Some of the vibrational modes evaluated were found to be more sensitive than others. Among the vibrational modes whose deviations were more poorly predicted by the pairwise PiMM approach are undoubtedly the NH_{2}, NH_{2}, NH_{2}, and NH_{2} modes.
5. Supporting Information
Listing of the CPuncorrected and CPcorrected vibrational frequencies predicted for the NH_{2} and –H modes when using each of the different DFT methods tested can be found in the Supporting Information linked to the online version of the paper.
Acknowledgment
The authors acknowledge the financial support from the Portuguese Foundation for Science and Technology—Unidade de QuímicaFísica Molecular (PEstOE/QUI/UI0070/2011) and Laboratório Associado CICECO.
Supplementary Materials
Figure S1: Variation of the CPcorrected ΔE_{olig} as a function of oligomer size and functional: (A) oligomers presenting just N–H···N interactions (group I) and (B) oligomers also C–H···N contacts (group II). Table S1: Selected CPuncorrected structural parameters characterizing the hydrogen bonding pattern of each molecular cluster, as a function of oligomer type and functional. Atom and cluster nomenclature in accord to Figure 1. Table S2: Vibrational frequency shifts predicted with B97D functional for the NH_{2} and C_{(2)}–H molecular groups. In each case, the second line refers to the CPcorrected value. Cluster nomenclature in accord to Figure 1. Table S3: Vibrational frequency shifts predicted with PW91PW functional for the NH_{2} and C_{(2)}–H molecular groups. In each case, the second line refers to the CPcorrected value. Cluster nomenclature in accord to Figure 1. Table S4: Vibrational frequency shifts predicted with B3LYP functional for the NH_{2} and C_{(2)}–H molecular groups. In each case, the second line refers to the CPcorrected value. Cluster nomenclature in accord to Figure 1. Table S5: Vibrational frequency shifts predicted with mPW1PW functional for the NH_{2} and C_{(2)}–H molecular groups. In each case, the second line refers to the CPcorrected value. Cluster nomenclature in accord to Fig. 1. Table S6: Vibrational frequency shifts predicted with B971 functional for the NH_{2} and C_{(2)}–H molecular groups. In each case, the second line refers to the CPcorrected value. Cluster nomenclature in accord to Figure 1. Table S7: Vibrational frequency shifts predicted with wB97XD functional for the NH_{2} and C_{(2)}–H molecular groups. In each case, the second line refers to the CPcorrected value. Cluster nomenclature in accord to Figure 1. Table S8: Frequency shifts predicted for the NH_{2} and C_{(2)}–H related vibrational modes by applying the PiMM approach to the CPuncorrected and CPcorrected B97D results (second line values are CPcorrected). Cluster nomenclature in accord to Figure 1.
References
 G. A. Jeffrey, An Introduction to Hydrogen Bonding, Oxford University Press, New York, NY, USA, 1997.
 Y. Lu, Y. Wang, Z. Xu et al., “CX⋯H contacts in biomolecular systems: how they contribute to proteinligand binding affinity,” The Journal of Physical Chemistry B, vol. 113, no. 37, pp. 12615–12621, 2009. View at: Publisher Site  Google Scholar
 S. Melandri, “‘Union is strength’: how weak hydrogen bonds become stronger,” Physical Chemistry Chemical Physics, vol. 13, no. 31, pp. 13901–13911, 2011. View at: Publisher Site  Google Scholar
 N. Rawat and P. Biswas, “Shape, flexibility and packing of proteins and nucleic acids in complexes,” Physical Chemistry Chemical Physics, vol. 13, no. 20, pp. 9632–9643, 2011. View at: Publisher Site  Google Scholar
 T. Wyttenbach and M. T. Bowers, “Hydration of biomolecules,” Chemical Physics Letters, vol. 480, no. 13, pp. 1–16, 2009. View at: Publisher Site  Google Scholar
 E. M. CabaleiroLago and M. A. Ríos, “Ab initio study of interactions in methylamine clusters. The significance of cooperative effects,” Journal of Chemical Physics, vol. 112, no. 5, pp. 2155–2163, 2000. View at: Google Scholar
 E. M. CabaleiroLago and M. A. Ríos, “An ab initio study of the interaction in dimethylamine dimer and trimer,” Journal of Chemical Physics, vol. 113, no. 21, pp. 9523–9531, 2000. View at: Google Scholar
 M. F. de Carvalho, R. A. Mosquera, and R. Rivelino, “A density functional theory study of the hydrogen bond interactions in glycine dimers,” Chemical Physics Letters, vol. 445, no. 4–6, pp. 117–124, 2007. View at: Publisher Site  Google Scholar
 Q. Li, X. An, B. Gong, and J. Cheng, “Spectroscopic and theoretical evidence for the cooperativity between redshift hydrogen bond and blueshift hydrogen bond in DMSO aqueous solutions,” Spectrochimica Acta A, vol. 69, no. 1, pp. 211–215, 2008. View at: Publisher Site  Google Scholar
 A. S. Mahadevi, Y. I. Neela, and G. N. Sastry, “A theoretical study on structural, spectroscopic and energetic properties of acetamide clusters [CH_{3}CONH_{2}] ($n=115$),” Physical Chemistry Chemical Physics, vol. 13, no. 33, pp. 15211–15220, 2011. View at: Publisher Site  Google Scholar
 R. D. Parra, M. Furukawa, B. Gong, and X. C. Zeng, “Energetics and cooperativity in threecenter hydrogen bonding interactions. I. DiacetamideX dimers (X = HCN, CH_{3}OH),” Journal of Chemical Physics, vol. 115, no. 13, pp. 6030–6035, 2001. View at: Publisher Site  Google Scholar
 M. J. Frisch, G. W. Trucks, H. B. Schlegel et al., Gaussian 09, Revision A.02, Gaussian 09, Revision A.02, Gaussian, Inc., Wallingford, Conn, USA, 2009.
 V. Barone, M. Biczysko, and M. Pavone, “The role of dispersion correction to DFT for modelling weakly bound molecular complexes in the ground and excited electronic states,” Chemical Physics, vol. 346, no. 1–3, pp. 247–256, 2008. View at: Publisher Site  Google Scholar
 I. Dabkowska, P. Jurečka, and P. Hobza, “On geometries of stacked and Hbonded nucleic acid base pairs determined at various DFT, MP2, and CCSD(T) levels up to the CCSD(T)/complete basis set limit level,” The Journal of Chemical Physics, vol. 122, no. 20, Article ID 204322, 9 pages, 2005. View at: Google Scholar
 S. Grimme, “Semiempirical GGAtype density functional constructed with a longrange dispersion correction,” Journal of Computational Chemistry, vol. 27, no. 15, pp. 1787–1799, 2006. View at: Publisher Site  Google Scholar
 J. Ireta, J. Neugebauer, and M. Scheffler, “On the accuracy of DFT for describing hydrogen bonds: dependence on the bond directionality,” The Journal of Physical Chemistry A, vol. 108, no. 26, pp. 5692–5698, 2004. View at: Publisher Site  Google Scholar
 R. Peverati, M. Macrina, and K. K. Baldridge, “Assessment of DFT and DFTD for potential energy surfaces of rare gas trimersimplementation and analysis of functionals and extrapolation procedures,” Journal of Chemical Theory and Computation, vol. 6, no. 7, pp. 1951–1965, 2010. View at: Publisher Site  Google Scholar
 A. M. Amado, S. M. Fiuza, L. A. E. Batista de Carvalho, and P. J. A. RibeiroClaro, “On the effects of changing Gaussian program version and SCRF defining parameters: Iisopropylamine as a case study,” Bulletin of the Chemical Society of Japan, vol. 85, no. 9, pp. 962–975, 2012. View at: Publisher Site  Google Scholar
 2011, http://www.gaussian.com/g_tech/g_ur/k_dft.htm.
 S. F. Boys and F. Bernardi, “The calculation of small molecular interactions by the differences of separate total energies. Some procedures with reduced errors,” Molecular Physics, vol. 19, no. 4, pp. 553–566, 1970. View at: Publisher Site  Google Scholar
 L. F. Holroyd and T. van Mourik, “Insufficient description of dispersion in B3LYP and large basis set superposition errors in MP2 calculations can hide peptide conformers,” Chemical Physics Letters, vol. 442, no. 1–3, pp. 42–46, 2007. View at: Publisher Site  Google Scholar
 Y.G. Wang, “Examination of DFT and TDDFT methods II,” The Journal of Physical Chemistry A, vol. 113, no. 41, pp. 10873–10879, 2009. View at: Publisher Site  Google Scholar
 I. A. W. Filot, A. R. A. Palmans, P. A. J. Hilbers, R. A. van Santen, E. A. Pidko, and T. F. A. de Greef, “Understanding cooperativity in hydrogenbondinduced supramolecular polymerization: a density functional theory study,” The Journal of Physical Chemistry B, vol. 114, no. 43, pp. 13667–13674, 2010. View at: Publisher Site  Google Scholar
 P. L. Huyskens, “Factors governing the influence of a first hydrogen bond on the formation of a second one by the same molecule or ion,” Journal of the American Chemical Society, vol. 99, no. 8, pp. 2578–2582, 1977. View at: Google Scholar
 M. P. Hodges, A. J. Stone, and S. S. Xantheas, “Contribution of manybody terms to the energy for small water clusters: a comparison of ab initio calculations and accurate model potentials,” The Journal of Physical Chemistry A, vol. 101, no. 48, pp. 9163–9168, 1997. View at: Google Scholar
 H. Kistenmacher, G. C. Lie, H. Popkie, and E. Clementi, “Study of the structure of molecular complexes. VI. Dimers and small clusters of water molecules in the HartreeFock approximation,” The Journal of Chemical Physics, vol. 61, no. 2, pp. 546–561, 1974. View at: Publisher Site  Google Scholar
 A. Ebrahimi, H. Roohi, and S. Mostafa Habibi, “The characterization of stationary points in the potential energy surface of difluoromethane dimer,” Journal of Molecular Structure: THEOCHEM, vol. 684, no. 1–3, pp. 87–93, 2004. View at: Publisher Site  Google Scholar
 J. Joseph and E. D. Jemmis, “Red, blue, or noshift in hydrogen bonds: a unified explanation,” Journal of the American Chemical Society, vol. 129, no. 15, pp. 4620–4632, 2007. View at: Publisher Site  Google Scholar
 A. Kovács, A. Szabó, D. Nemcsok, and I. Hargittai, “Blueshifting CH⋯X, (X = O, halogen) hydrogen bonds in the dimers of formaldehyde derivatives,” The Journal of Physical Chemistry A, vol. 106, no. 23, pp. 5671–5678, 2002. View at: Publisher Site  Google Scholar
 P.P. Zhou and W.Y. Qiu, “Red and blueshifted hydrogen bonds in the cistrans noncyclic formic acid dimer,” ChemPhysChem, vol. 10, no. 11, pp. 1847–1858, 2009. View at: Publisher Site  Google Scholar
 A. K. Chandra, S. Parveen, S. Das, and T. ZeegersHuyskens, “Blue shifts of the CH stretching vibrations in hydrogenbonded and protonated trimethylamine. Effect of hyperconjugation on bond properties,” Journal of Computational Chemistry, vol. 29, no. 9, pp. 1490–1496, 2008. View at: Publisher Site  Google Scholar
 A. K. Chandra, S. Parveen, and T. ZeegersHuyskens, “Anomeric effects in the symmetrical and asymmetrical structures of triethylamine. Blueshifts of the CH stretching vibrations in complexed and protonated triethylamine,” The Journal of Physical Chemistry A, vol. 111, no. 36, pp. 8884–8891, 2007. View at: Publisher Site  Google Scholar
 K.H. Chen, J.H. Lii, Y. Fan, and N. L. Allinger, “Molecular mechanics (MM4) study of amines,” Journal of Computational Chemistry, vol. 28, no. 15, pp. 2391–2412, 2007. View at: Publisher Site  Google Scholar
 X. Li, L. Liu, and H. B. Schlegel, “On the physical origin of blueshifted hydrogen bonds,” Journal of the American Chemical Society, vol. 124, no. 32, pp. 9639–9647, 2002. View at: Publisher Site  Google Scholar
 S. Padrào, S. M. Fiuza, A. M. Amado, A. M. Amorim da Costa, and L. A. E. Batista de Carvalho, “Validation of the mPW1PW quantum chemical calculations for the vibrational study of organic molecules—reassignment of the isopropylamine vibrational spectra,” Journal of Physical Organic Chemistry, vol. 24, no. 2, pp. 110–121, 2011. View at: Publisher Site  Google Scholar
 J. R. Durig, G. A. Guirgis, and D. A. C. Compton, “Analysis of torsional spectra of molecules with two internal C_{3v} rotors. 13. Vibrational assignments, torsional potential functions, and gas phase thermodynamic functions of isopropylamined_{0} and d_{2},” The Journal of Physical Chemistry, vol. 83, no. 10, pp. 1313–1323, 1979. View at: Publisher Site  Google Scholar
 M. P. M. Marques, L. A. E. Batista de Carvalho, and J. Tomkinson, “Study of biogenic and α,ωpolyamines by combined inelastic neutron scattering and Raman spectroscopies and by ab initio molecular orbital calculations,” The Journal of Physical Chemistry A, vol. 106, no. 11, pp. 2473–2482, 2002. View at: Publisher Site  Google Scholar
 A. M. Amado, J. C. Otero, M. P. M. Marques, and L. A. E. Batista de Carvalho, “Spectroscopic and theoretical studies on solid 1,2ethylenediamine dihydrochloride salt,” ChemPhysChem, vol. 5, no. 12, pp. 1837–1847, 2004. View at: Publisher Site  Google Scholar
 E. Cubero, M. Orozco, P. Hobza, and F. J. Luque, “Hydrogen bond versus antihydrogen bond: a comparative analysis based on the electron density topology,” The Journal of Physical Chemistry A, vol. 103, no. 32, pp. 6394–6401, 1999. View at: Google Scholar
 P. Hobza, “The Hindex unambiguously discriminates between hydrogen bonding and improper blueshifting hydrogen bonding,” Physical Chemistry Chemical Physics, vol. 3, no. 13, pp. 2555–2556, 2001. View at: Publisher Site  Google Scholar
 B. J. van der Veken, W. A. Herrebout, R. Szostak, D. N. Shchepkin, Z. Havlas, and P. Hobza, “The nature of improper, blueshifting hydrogen bonding verified experimentally,” Journal of the American Chemical Society, vol. 123, no. 49, pp. 12290–12293, 2001. View at: Publisher Site  Google Scholar
 W. Zierkiewicz, D. Michalska, Z. Havlas, and P. Hobza, “Study of the nature of improper blueshifting hydrogen bonding and standard hydrogen bonding in the X_{3}CH⋯OH_{2} and XH⋯OH_{2} complexes (X = F, Cl, Br, I): a correlated ab initio study,” ChemPhysChem, vol. 3, no. 6, pp. 511–518, 2002. View at: Google Scholar
 P. J. A. RibeiroClaro, L. Batista de Carvalho, and A. M. Amado, “Evidence of dimerization through CH⋯O interactions in liquid 4methoxybenzaldehyde from Raman spectra and Ab Initio calculations,” Journal of Raman Spectroscopy, vol. 28, no. 11, pp. 867–872, 1997. View at: Google Scholar
 P. J. A. RibeiroClaro, M. P. M. Marques, and A. M. Amado, “Experimental and theoretical evidence of CH⋯O hydrogen bonding in liquid 4fluorobenzaldehyde,” ChemPhysChem, vol. 3, no. 7, pp. 599–606, 2002. View at: Google Scholar
 T. Steiner, E. B. Starikov, A. M. Amado, and J. J. C. TeixeiraDias, “Weak hydrogen bonding. Part 2. The hydrogen bonding nature of short CH ⋯π contacts: cristallographic, spectroscopic and quantummechanical studies of some terminal alkynes,” Journal of the Chemical Society, Perkin Transactions 2, no. 7, pp. 1321–1326, 1995. View at: Publisher Site  Google Scholar
 P. D. Vaz, M. Nolasco, N. Fonseca et al., “CH⋯O hydrogen bonding in 4phenylbenzaldehyde: a comprehensive crystallographic, spectroscopic and computational study,” Physical Chemistry Chemical Physics, vol. 7, no. 16, pp. 3027–3034, 2005. View at: Publisher Site  Google Scholar
 S. Wolfe, H. B. Schlegel, M. H. Whangbo, and F. Bernardi, “On the origin of the Bohlmann bands,” Canadian Journal of Chemistry, vol. 52, no. 22, pp. 3787–3792, 1974. View at: Publisher Site  Google Scholar
 M. M. Nolasco, A. M. Amado, and P. J. A. RibeiroClaro, “Computationallyassisted approach to the vibrational spectra of molecular crystals: study of hydrogenbonding and pseudopolymorphism,” ChemPhysChem, vol. 7, no. 10, pp. 2150–2161, 2006. View at: Publisher Site  Google Scholar
 M. Sardo, A. M. Amado, and P. J. A. RibeiroClaro, “Hydrogen bonding in nitrofurantoin polymorphs: a computationassisted spectroscopic study,” Journal of Raman Spectroscopy, vol. 40, no. 12, pp. 1956–1965, 2009. View at: Publisher Site  Google Scholar
 M. M. Nolasco, A. M. Amado, and P. J. A. RibeiroClaro, “Effect of hydrogen bonding in the vibrational spectra of transcinnamic acid,” Journal of Raman Spectroscopy, vol. 40, no. 4, pp. 394–400, 2009. View at: Publisher Site  Google Scholar
 M. M. Nolasco, A. M. Amado, and P. J. A. RibeiroClaro, “Insights into phase stability of anhydrous/hydrate systems: a Ramanbased methodology,” Journal of Raman Spectroscopy, vol. 41, no. 3, pp. 340–349, 2010. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2013 Ana M. Amado et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.