- About this Journal
- Abstracting and Indexing
- Aims and Scope
- Article Processing Charges
- Articles in Press
- Author Guidelines
- Bibliographic Information
- Citations to this Journal
- Contact Information
- Editorial Board
- Editorial Workflow
- Free eTOC Alerts
- Publication Ethics
- Reviewers Acknowledgment
- Submit a Manuscript
- Subscription Information
- Table of Contents
Advances in Physical Chemistry
Volume 2011 (2011), Article ID 347072, 10 pages
Vibrational Study and Force Field of the Citric Acid Dimer Based on the SQM Methodology
1Cátedra de Química General, Instituto de Química Inorgánica, Facultad de Bioquímica, Química y Farmacia, Universidad Nacional de Tucumán, Ayacucho 471, 4000 S. M. de Tucumán, Argentina
2Cátedra de Fisicoquímica I, Instituto de Química Física, Facultad de Bioquímica, Química y Farmacia, Universidad Nacional de Tucumán, San Lorenzo 456, T4000CAN S, M. de Tucumán., Argentina
3Centro de Química Inorgánica, CEQUINOR/CCT, Departamento de Química, Facultad de Ciencias Exactas, Universidad Nacional de La Plata, 47 y 115. CC. 962 (B1900AVV), 1900 La Plata, Argentina
4Laboratorio de Fisicoquímica, Departamento de Física, Facultad de Ciencias Exactas y Tecnología, Universidad Nacional de Tucumán, Avenida Independencia 1800, T4000CAN S. M. de Tucumán., Argentina
Received 11 February 2011; Accepted 14 April 2011
Academic Editor: Joel Bowman
Copyright © 2011 Laura Cecilia Bichara 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.
We have carried out a structural and vibrational theoretical study for the citric acid dimer. The Density Functional Theory (DFT) method with the B3LYP/6-31G∗ and B3LYP/6-311++ methods have been used to study its structure and vibrational properties. Then, in order to get a good assignment of the IR and Raman spectra in solid phase of dimer, the best fit possible between the calculated and recorded frequencies was carry out and the force fields were scaled using the Scaled Quantum Mechanic Force Field (SQMFF) methodology. An assignment of the observed spectral features is proposed. A band of medium intensity at 1242 together with a group of weak bands, previously not assigned to the monomer, was in this case assigned to the dimer. Furthermore, the analysis of the Natural Bond Orbitals (NBOs) and the topological properties of electronic charge density by employing Bader's Atoms in Molecules theory (AIM) for the dimer were carried out to study the charge transference interactions of the compound.
The structural and vibrational studies of citric acid are of chemical, biochemical, and technological interest [1–8]. The crystal and molecular structure of anhydrous citric acid was determined by Nordman et al.  and Glusker et al.  by means of X-ray crystal analyses. Both authors demonstrated that two molecules of the acid are linked through hydrogen bonds between their carboxylic groups originating the centre-symmetrical cyclic structure of the dimeric species. Many previous studies in relation to the chemical and spectroscopic properties of this acid as monomer have been reported [9–14] while an ab-initio/HF force field study for the vibrations of citric acid and citrate trianion was published by Tarakeshwar and Manogaran . In the latter study, the calculated geometries and vibrational frequencies only for the monomer citric acid were used to assign the vibrational spectra. A complete vibrational characterization of this acid is very important to obtain information about chemical and physical components in food, biological materials, and pharmaceutical products at the molecular level [14, 16]. So far, there are no publications about experimental or high-level theoretical studies on the geometries and force field for the dimer citric acid. Hence, obtaining reliable parameters by DFT theoretical methods is an appealing alternative to have a precise knowledge of the normal vibration modes. The aim of this paper is to report a DFT theoretical and experimental study of dimer citric acid by infrared spectroscopy in order to effect a complete assignment of every observed bands in the vibrational spectra by taking into account the presence of the dimeric species in the solid phase. For this purpose, the optimized geometries and frequencies for the normal vibration modes, considering the dimer with two units of the free acid, were calculated. For a complete and reliable assignment of the compound, the DFT calculations were combined with the SQMFF methodology [17–19] in order to fit the theoretical frequencies values to the experimental ones. In addition, the force constants were calculated and compared with those reported for the monomer by using the HF/4-21G method . Here, we demonstrated that a DFT molecular force field for the dimer calculated by using the B3LYP/6-31G* combination is well represented. DFT normal mode assignments, in terms of the potential energy distribution, are in agreement with those obtained from the normal coordinate analysis. The total energy calculated with the 6-31G* basis set for citric acid dimer was corrected for basis set superposition error (BSSE) by the standard Boys—Bernardi counterpoise method . Besides, the dimer electronic properties were evaluated by NBO [21–24] and Atoms in Molecules (AIMs) [25, 26] studies in order to analyze the nature and magnitude of the intermolecular interactions.
2. Experimental and Theoretical Calculations
A pure anhydrous Mallinckrodt commercial sample of citric acid (C6H8O7) (CA) was used. The FTIR spectrum of the solid substance in KBr pellets was recorded with a Bruker IFS 66 FTIR-spectrophotometer from 4000 to 400 cm−1. FT-Raman spectrum was also collected on a Bruker IFS 66 FT-IR spectrophotometer provided with a NIR Raman attachment equipped with an Nd:YAG laser at 1064 nm laser.
The starting point for the dimer geometry optimization was a structure with bond lengths taken from the X-ray crystal structure [6, 7]. These initial geometries were optimized by using the hybrid B3LYP method [27, 28] with the 6-31G* and 6-311++G** basis sets, as implemented in the GAUSSIAN 03 program . The atom structure and labelling can be seen in Figure 1. The electronic charge density topological analysis was performed by using the Atoms in Molecules (AIMs) methodology , by the AIM2000 program package  while the NBO calculations were performed by using the NBO 3.1 program , as implemented in the GAUSSIAN 03 package . The harmonic wavenumbers and the valence force fields in Cartesian coordinates were calculated at the B3LYP/6-31G* and 6-311++G** levels of approximation. The resulting force fields were transformed into “natural” internal coordinates by using the MOLVIB program [30, 31]. The natural internal coordinates are presented in Table S1 of the Supplementary Material at doi/0.1155/2011/347072 and have been defined as proposed by Fogarasi et al. . In this case, some coordinates are slightly different from those used in the previous study . The six inter-monomer coordinates for this dimer are similar to those defined by other authors for other dimeric acids [33–36], which can be seen in Supplementary Table S2. Following the SQMFF procedure [37, 38] the harmonic force fields were evaluated at B3LYP/6-31G* level. In the scaling procedure the diagonal force constants associated with a particular symmetry coordinate are multiplied by a small number of scale factors (,, etc.). As a result, each off-diagonal constant is automatically scaled by the geometrical mean of factors () 1/2 for the corresponding diagonal constants. Finally, the scale factors are adjusted by an iterative method until they reproduce the experimental wavenumbers as possible. In this case, we use the Rauhut and Pulay scale factors [17, 18]. The potential energy distribution (PEDs) components higher than or equal to 10% are subsequently calculated with the resulting SQM. The nature of every vibration modes was carried out by means of the GaussView program .
3. Results and Discussion
3.1. Geometries Optimization
Supplementary Table S3 shows the total energies and the corresponding dipole moment values for the dimer structures with the B3LYP method by using 6-31G* and 6-311++G** basis sets compared with the ones corresponding to two monomer units. The energy difference between the cyclic dimer and the two citric acid isolated monomers with the B3LYP/6–31G* method is 82.86 kJ mol−1 while by using the 6-311++G** basis set the value at 68.64 kJmol−1 decreases. Reva and Stepanian  predict an energy difference of 86.00 kJ mol−1 between the cyclic dimer and two isolated monomers of benzoic acid (ab initio/4-31G) while the calculated energy difference between their monomers is 40.12 kJ mol−1. In this case, the total energy for the citric acid dimer is smaller than the sum of the energies of each monomer (Supplementary Table S3) by using both basis sets. Such observation might mean that the dimer presence even in the gas phase is preferable to the isolated monomer, as observed in the 4-hydroxybenzoic acid dimer . This result is very important to carry out the assignment of the vibrational spectra of this acid. In general, our theoretical values are in agreement with the experimental values for the citric acid [6, 7]. The dimer stability was investigated by means of NBO [21–23] and AIM calculations . In the NBO analysis, the second-order perturbation energies (donor → acceptor) involving the most important delocalization for the dimer clearly reveal the energies and occupancies of the main contributions: LP(2)O19 →O22-H21 (115.2 kJ mol−1), LP(2)O27 →O19-H20 (38.4 kJ mol−1) and LP(2)O27 →C3-H6 (2.4 kJ mol−1) to the delocalization energy due to the OH–O bonds between both structures. In Bader’s topological analysis  the localization of the Bond Critical Point (BCP) and Ring Critical Point (RCP) in the ρ(r) and the Laplacian values at these points is important to characterize molecular electronic structure. The critical point has the typical properties of the closed-shell interaction when the value of ρ(r) is relatively low, the relationship is , and is positive indicating that the interaction is dominated by the charge contraction away from the interatomic surface toward each nucleus. The topological properties for the citric acid dimer can be seen in Supplementary Table S4. The values obtained clearly reveal the formation of hydrogen bonds between two COOH groups of each monomer. In addition, the seven RCPs calculated by means AIM analysis for the dimer corroborate those contributions found in the NBO study.
4. Vibrational Analysis
The dimer’s structure has C1 symmetry and 120 vibrational normal modes where all vibrations are IR and Raman active. The recorded infrared spectra for the compound in the solid phase can be seen in Figure 2. The assignment of the experimental bands to the 120 expected normal vibration modes was made on the basis of the Potential Energy Distribution (PED) in terms of symmetry coordinates by using the B3LYP/6-31G* data and taking into account the previous assignment for the monomer  and for related molecules [33–36, 40]. Table 1 shows the observed frequencies and the assignment of the vibrational normal modes while the PED based on the 6-31G* basis set is shown in Supplementary Table S5. The theoretical calculations reproduce the normal wavenumbers by using the 6-31G* basis set, with an initial rmsd value of 11.5 cm−1 while when the SQMFF method is applied by using Pulay’s scaling factors, the final value is 1.5 cm−1. The best results are obtained with a B3LYP/6-31G* calculation because the used transferable scale factors are only defined for this method. The SQM force fields for the compound dimer can be obtained at request. The broad band observed between 3400 and 2200 cm−1 in the experimental IR spectra of the compound in the solid phase can be attributed to the O–H hydrogen bond formed by the spatial arrangement of molecules in the lattice crystal [41, 42]. Note that the theoretical infrared spectra for the monomer structure demonstrate a little agreement with the experimental spectrum (Figure 2), especially in the high wavenumbers cm−1 region. On the other hand, in the theoretical Raman spectrum shown in Figure 3, in relation to the experimental one, a higher intensity in the bands attributed to the OH stretching is observed. The existence of molecules in the solid phase as dimer justifies such a difference because this species restricts the motion of the COOH groups and is strongly affected by the double hydrogen bond, as was also observed in the 4-hydroxybenzoic acid in solid state . Thus, the wavenumbers in the dimer are shifted with respect to the monomer and the splitting of the antisymmetric and symmetric modes is particularly large. In this compound, both monomer and dimer structures are present in the solid phase and, for this reason, a comparison between the average calculated infrared spectra from B3LYP/6-31G* level for the monomer and dimer by using average wavenumbers and intensities with the corresponding experimental ones demonstrates a good correlation, especially in the high wavenumbers region. The discussion of the assignment of the most important groups for the compound (see Table 1 and Supplementary Table ) is presented as follows.
4.1. Band Assignments
In the previous assignment for the monomer the theoretical HF/4-21G calculations  predict that the behaviour of the central COOH group is different from the terminal groups, as in our case by using the B3LYP/6-31G* level and, as a consequence, the central OH stretching mode belonging to that central group is calculated at higher wavenumbers than the other ones. In the dimer, those modes are assigned according to the SQM calculations; for this reason, the OH stretching modes are associated with the Raman band at 3535 cm−1, the IR bands at 3498, 3446, 3291, and 2933 cm−1, and the shoulder located at 3350 cm−1. The C–OH scissoring mode belonging to the COOH groups involved in the H-bonding formation is calculated at higher wavenumbers than the other ones (1483 and 1402 cm−1) probably because they are coupled with the corresponding C–O stretching modes; for this reason, both modes are assigned to the bands at 1493 and 1389 cm−1. The remaining C–OH scissoring modes are assigned to the shoulder and IR bands at 1340, 1308, and 1081 cm−1, respectively, as seen in Table 1. In the region of lower wavenumbers, the OH torsion modes are expected and these modes are also affected by the presence of H bonding. Previously, those modes in the monomer were assigned, according to HF/4-21G calculations , to the bands at 645, 598, 556, and 383 cm−1 while by using the B3LYP/6-31G* level they are calculated between 655 and 387 cm−1; for this reason, the OH torsion modes are associated with the Raman bands at 642, 627, 492, and 415 cm−1.
These modes in the monomer were previously observed in the 3000–2900 cm−1 region . In this case, the out-of-phase and in-phase antisymmetric stretching modes of these groups can be assigned to the shoulder and IR band at 3035 and 2994 cm−1, as indicated in Table 1, while the out-of-phase and in-phase symmetric modes are assigned to the weak IR band at 2950 cm−1 and to the medium intensity Raman band at 2982 cm−1. In the monomer, the scissoring modes were assigned at 1427 cm−1 while for the dimer those modes were also clearly calculated in this region; thus, they are assigned to the bands at 1430, 1424, 1420, and 1413 cm−1. In the above-mentioned study , the wagging modes were associated with the very weak IR band at 1325 cm−1. Here, as predicted by calculation, the bands at 1389, 1365, 1308, 1292, and 1036 cm−1 are assigned to the wagging modes. Previously, in the monomer the twisting and rocking modes were assigned at 1350 and 966 cm−1, respectively . In the dimer, the IR bands located at 1358, 1292, and 1214 cm−1 and the Raman band at 1285 cm−1 are associated with the rocking modes, while the expected twisting modes are associated with the IR bands at 904 and 881 cm−1 and with the shoulder in the Raman spectrum at 1036, 656 and 514 cm−1.
In the previous studies for the monomer  the band at 1706 cm−1 and the shoulder at 1695 cm−1 were assigned to the C=O stretching modes. Note that the band at 1706 cm−1 (Table 1) was related to the strong Raman band at 1735 cm−1 and assigned to the same mode in spite of the difference between them (19 cm−1). In the dimer, the six expected C=O stretching modes are clearly predicted by SQM/B3LYP/6-31G* calculation at 1786, 1783, 1775, 1771, 1726, and 1677 cm−1. On the other hand, the IR spectrum clearly shows the presence of more than one band in this region; thus, two very strong bands at 1756 and 1708 cm−1 can be seen, each band presenting a shoulder at 1745 and 1698 cm−1, respectively. Therefore, the IR and Raman bands between 1756 and 1698 cm−1 are clearly assigned to those modes. In this way, those shoulders should be also associated with other C=O stretching modes, and so this fact would confirm the presence of the dimer in the solid state. In both HF/4-21G  and B3LYP/6-31G* calculations, the C–O stretching modes are coupled with the C–OH scissoring modes. In the monomer , those modes were assigned to the IR bands at 1195, 1141, and 1108 cm−1, while in the dimer, the IR bands at 1242, 1174, and 1140 cm−1 and the medium intensity Raman band at 904 cm−1 are assigned to these stretching modes. Previously, the medium intensity band in the IR spectrum at 1242 cm−1 was not observed and therefore not assigned. Taking into consideration their relative position, intensities, and the splitting predicted by calculation, the expected out-of-plane deformation modes (COO), in-phase and out-of-phase modes, were predicted as pure modes between 792 and 485 cm−1 and coupled with other modes at 655 and 499 cm−1; for this reason, they were assigned to the bands at 807, 795, 686, 574, and 492 cm−1, as observed in Table 1. In the previous study , only one deformation mode (δCOO) was assigned to the band at 420 cm−1. Here, as predicted by calculation, those expected deformation modes for the dimer are assigned to the IR band and shoulder at 781 and 789 cm−1, respectively, and to the Raman bands at 700, 627, 574, and 541 cm−1. The rocking modes (ρCOO) for the monomer were not previously assigned and, in this case, the ρCOO modes were easily assigned to the IR bands at 599, 550, and 438 cm−1 and to the Raman bands at 367 and 348 cm−1. Finally, the twisting (τwCOO) modes for the monomer were only predicted by the HF/4-21G level at 88, 68, 44, and 39 cm−1 and, in our case, for the dimer they are predicted strongly coupled between them and in the low wavenumbers region; this way, they could not be assigned because there are not observed bands in the vibrational spectra in those regions.
In the previous assignment , the bands in the infrared spectrum at 1090, 1053, 900, 850, and 786 cm−1 were attributed to the C–C stretching modes, while the CCC deformation modes were not assigned because they appear at very low frequencies. Here, the bands in the IR spectrum at 1430, 1051 and 945 cm−1 and in the Raman spectrum at 1469, 1053, 966, 850, and 842 cm−1 are assigned to the C–C stretching modes. The calculations predicted the CCC deformation modes in the lower wavenumbers region; for this reason, they are assigned to the Raman bands between 323 and 212 cm−1 while the two modes predicted at 108 and 92 cm−1 were not assigned.
The intermonomer vibrational modes of the dimer usually appear in the 270–10 cm−1 region [33–36, 40] and they are related to restricted translations or rotations of one molecule against the other (see Supplementary Table S2). Here, the six intermonomer modes appear strongly coupled with other modes and they are observed in the Raman spectrum with the exception of the (OH–O) mode that was predicted by calculations with a higher contribution PED (42%) at 827 cm−1; for this reason, it is assigned to the very weak Raman band at 842 cm−1. Table 1 shows the theoretical assignments for the dimer intermonomer vibrational modes.
5. Force Field
The force constants for the dimer were estimated by using Pulay et al. [37, 38] scaling procedure as mentioned above. The force constants expressed in terms of simple valence internal coordinates were calculated from the corresponding scaled force fields by using the MOLVIB program [30, 31]. It is interesting to compare the principal force constants, which were collected in Table 2 and calculated at the B3LYP/6-31G* level for some vibrations, with those reported for the monomer of this acid by using the HF/4-21G method . The calculated force constants values for the monomer by using the ab-initio method are slightly different from the values for the dimer by using the B3LYP/6-31G* level, as expected. The HF method overestimates the calculated force constants values due to the fact that it neglects the electron correlation and unharmony. On the other hand, the little variation in the f(C=O) force constant value in the monomer with reference to the dimer is justified because in the latter species two molecules of the acid are linked by hydrogen bonds between their carboxylic groups; therefore, the involved C=O bonds decrease the double bonds character. For this reason, a lower value in the f(O–H) force constant can be also seen in the dimer in relation to the monomer. Obviously, the force constants related to the CH2 groups in the monomer have lower values than the dimer because in the latter species the number of those groups is higher.
We conclude the following.(i)The assignment previously made for the monomer  was corrected and completed in accordance with the present theoretical B3LYP/6-31G* results. Also, the assignment of the intermonomer vibrational modes for the dimer is carried out and the assignments of the 120 normal vibration modes corresponding to citric acid dimer are reported.(ii)Six intense bands in the infrared spectrum at 3498, 3291, 1756, 1708, 1174, and 1140 cm−1 are reported to characterize the compound dimer.(iii)Numerous bands of different intensities observed in the vibracional spectra not previously assigned, in this work were assigned to the citric acid dimer.(iv)The NBO analysis reveals that in the citric acid dimer there are strong intramolecular interactions of charge transfer from O lone pairs to the (O–H) antibonds.(v)The AIM studies of the charge density for the dimer satisfy the criteria of hydrogen bond interactions.(vi)The complete B3LYP/6-31G* force field for the citric acid dimer was determined as well as the principal force constants for stretching and deformation modes.
This work was subsidized with grants from CIUNT (Consejo de Investigaciones, Universidad Nacional de Tucumán) and CONICET (Consejo Nacional de Investigaciones Científicas y Técnicas, R. Argentina). The authors thank Professor Tom Sundius for his permission to use MOLVIB.
- H. Eggerer, U. Remberger, and C. Grünewälder, “On the mechanism of biological transformation of citric acid. v. citrate synthase, a hydrolase for malyl coenzyme a.zum mechanismus der biologischen umwandlung von citronensaeure. v. citrat-synthase, eine hydrolase fuer malyl-coenzym A,” Biochemische Zeitschrift, vol. 339, pp. 436–453, 1964.
- P. Wunderwald and H. Eggerer, “18-O-studies with citrate synthase and malate synthase,” European Journal of Biochemistry, vol. 11, no. 1, pp. 97–105, 1969.
- T. S. Raman and E. R. B. Shanmugasundaram, “The effects of two quaternary ammonium compounds on citric acid and sterol synthesis in Aspergillus niger,” Journal of General Microbiology, vol. 31, pp. 23–29, 1963.
- S. Schlücker, B. Küstner, A. Punge, R. Bonfig, A. Marx, and P. Ströbel, “Immuno-Raman microspectroscopy: in situ detection of antigens in tissue specimens by surface-enhanced Raman scattering,” Journal of Raman Spectroscopy, vol. 37, no. 7, pp. 719–721, 2006.
- J. De Gelder, K. De Gussem, P. Vandenabeele, P. De Vos, and L. Moens, “Methods for extracting biochemical information from bacterial Raman spectra: an explorative study on Cupriavidus metallidurans,” Analytica Chimica Acta, vol. 585, no. 2, pp. 234–240, 2007.
- C. E. Nordman, A. S. Weldon, and A. L. Patterson, “X-ray crystal analysis of the substrates of aconitase—I. Rubidium dihydrogen citrate,” Acta Crystallographica, vol. 13, pp. 414–417, 1960.
- J. P. Glusker, J. A. Minkin, and A. L. Patterson, “X-ray crystal analysis of the substrates of aconitase—9. A refinement of the structure of anhydrous citric acid,” Acta Crystallographica Section B, vol. 25, no. 6, pp. 1066–1072, 1969.
- D. E. Zacharias, J. P. Glusker, R. Guthrie, A. C. Sullivan, and P. A. Srere, “The structure of a citric anhydride derivative, C8H6O7,” Acta Crystallographica Section C, vol. 40, pp. 2100–2103, 1984.
- R. Lemor, M. B. Kruger, D. M. Wieliczka, P. Spencer, and T. May, “Dentin etch chemistry investigated by Raman and infrared spectroscopy,” Journal of Raman Spectroscopy, vol. 31, no. 3, pp. 171–176, 2000.
- J. L. Moreira and L. Santos, “Analysis of organic acids in wines by Fourier-transform infrared spectroscopy,” Analytical and Bioanalytical Chemistry, vol. 382, no. 2, pp. 421–425, 2005.
- A. Nose, M. Myojin, M. Hojo, T. Ueda, and T. Okuda, “Proton nuclear magnetic resonance and Raman spectroscopic studies of Japanese sake, an alcoholic beverage,” Journal of Bioscience and Bioengineering, vol. 99, no. 5, pp. 493–501, 2005.
- J. D. Tsay and T. T. Fang, “Effect of pH on the chemistry of the barium titanium citrate gel and its thermal decomposition behavior,” Journal of the American Ceramic Society, vol. 84, no. 11, pp. 2475–2478, 2004.
- A. Caillet, N. Sheibat-Othman, and G. Fevotte, “Crystallization of monohydrate citric acid. 2. Modeling through population balance equations,” Crystal Growth and Design, vol. 7, no. 10, pp. 2088–2095, 2007.
- Y. Xiong, J. M. McLellan, Y. Yin, and Y. Xia, “Synthesis of palladium icosahedra with twinned structure by blocking oxidative etching with citric acid or citrate ions,” Angewandte Chemie, vol. 46, no. 5, pp. 790–794, 2007.
- P. Tarakeshwar and S. Manogaran, “Ground state vibrations of citric acid and the citrate trianion-an ab initio study,” Spectrochimica Acta Part A, vol. 50, no. 14, pp. 2327–2343, 1994.
- P. K. Ghosh and D. S. Jayas, “Use of spectroscopic data for automation in food processing industry,” in Proceedings of the American Society of Agricultural and Biological Engineers (ASABE '08), Providence, Ri, USA, June 2008.
- G. Rauhut and P. Pulay, “Transferable scaling factors for density functional derived vibrational force fields,” Journal of Physical Chemistry, vol. 99, no. 10, pp. 3093–3100, 1995.
- G. Rauhut and P. Pulay, “Erratumml: transferable scaling factors for density functional derived vibrational force fields,” Journal of Physical Chemistry, vol. 99, no. 39, Article ID 14572, 1995.
- F. Kalincsák and G. Pongor, “Extension of the density functional derived scaled quantum mechanical force field procedure,” Spectrochimica Acta Part A, vol. 58, no. 5, pp. 999–1011, 2002.
- 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.
- A. E. Reed, L. A. Curtiss, and F. Weinhold, “Intermolecular interactions from a natural bond orbital, donor-acceptor viewpoint,” Chemical Reviews, vol. 88, no. 6, pp. 899–926, 1988.
- J. P. Foster and F. Weinhold, “Natural hybrid orbitals,” Journal of the American Chemical Society, vol. 102, no. 24, pp. 7211–7218, 1980.
- A. E. Reed and F. Weinhold, “Natural localized molecular orbitals,” The Journal of Chemical Physics, vol. 83, no. 4, pp. 1736–1740, 1985.
- E. D. Gledening, J. K. Badenhoop, A. D. Reed, J. E. Carpenter, and F. Weinhold, NBO 3.1, Theoretical Chemistry Institute, University of Wisconsin, Madison, Wis, USA, 1996.
- R. F. W. Bader, Atoms in Molecules, A Quantum Theory, Oxford University Press, Oxford, UK, 1990.
- F. B. Köning, J. Schönbohm, and D. Bayles, “AIM2000: a program to analyze and visualize atoms in molecules,” Journal of Computational Chemistry, vol. 22, no. 5, pp. 545–559, 2001.
- A. D. Becke, “Density-functional thermochemistry. III. The role of exact exchange,” The Journal of Chemical Physics, vol. 98, no. 7, pp. 5648–5652, 1993.
- C. Lee, W. Yang, and R. G. Parr, “Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density,” Physical Reiew B, vol. 37, no. 2, pp. 785–789, 1988.
- M. J. Frisch, G. W. Trucks, H. B. Schlegel, et al., Gaussian 03, Revision B.01, Gaussian inc., Pittsburgh, Pa, USA, 2003.
- T. Sundius, “Molvib—a flexible program for force field calculations,” Journal of Molecular Structure, vol. 218, pp. 321–326, 1990.
- T. Sundius, “Scaling of ab initio force fields by MOLVIB,” Vibrational Spectroscopy, vol. 29, no. 1-2, pp. 89–95, 2002.
- G. Fogarasi, X. Zhou, P. W. Taylor, and P. Pulay, “The calculation of ab initio molecular geometries: efficient optimization by natural internal coordinates and empirical correction by offset forces,” Journal of the American Chemical Society, vol. 114, no. 21, pp. 8191–8201, 1992.
- G. M. Florio, T. S. Zwier, E. M. Myshakin, K. D. Jordan, and E. L. Sibert, “Theoretical modeling of the OH stretch infrared spectrum of carboxylic acid dimers based on first-principles anharmonic couplings,” Journal of Chemical Physics, vol. 118, no. 4, pp. 1735–1746, 2003.
- M. A. Palafox, J. L. Núñez, and M. Gil, “Theoretical quantum chemical study of benzoic acid: geometrical parameters and vibrational wavenumbers,” International Journal of Quantum Chemistry, vol. 89, no. 1, pp. 1–24, 2002.
- M. Plazanet, N. Fukushima, M. R. Johnson, A. J. Horsewill, and H. P. Trommsdorff, “The vibrational spectrum of crystalline benzoic acid: inelastic neutron scattering and density functional theory calculations,” Journal of Chemical Physics, vol. 115, no. 7, pp. 3241–3248, 2001.
- S. A. Brandán, F. M. López, M. Montejo, J. J. López-González, and A. Ben Altabef, “Theoretical and experimental vibrational spectrum study of 4-hydroxybenzoic acid as monomer and dimer,” Spectrochimica Acta—Part A, vol. 75, no. 5, pp. 1422–1434, 2010.
- P. Pulay, G. Fogarasi, F. Pang, and J. E. Boggs, “Systematic ab initio gradient calculation of molecular geometries, force constants, and dipole moment derivatives,” Journal of the American Chemical Society, vol. 101, no. 10, pp. 2550–2560, 1979.
- P. Pulay, G. Fogarasi, G. Pongor, J. E. Boggs, and A. Vargha, “Combination of theoretical ab initio and experimental information to obtain reliable harmonic force constants. Scaled Quantum Mechanical (SQM) force fields for glyoxal, acrolein, butadiene, formaldehyde, and ethylene,” Journal of the American Chemical Society, vol. 105, no. 24, pp. 7037–7047, 1983.
- A. B. Nielsen and A. J. Holder, Gauss View 3.0, User’s Reference, Gaussian inc., Pittsburgh, Pa, USA, 2003.
- I. D. Reva and S. G. Stepanian, “An infrared study on matrix-isolated benzoic acid,” Journal of Molecular Spectroscopy and Molecular Structure, vol. 349, pp. 337–340, 1995.
- J. Zinczuk, A. E. Ledesma, S. A. Brandán, O. E. Piro, J. J. López-González, and A. Ben Altabef, “Structural and vibrational study of 2-(-furyl)-4,5-1H-dihydroimidazole,” Journal of Physical Organic Chemistry, vol. 22, no. 12, pp. 1166–1177, 2009.
- A. E. Ledesma, S. A. Brandán, J. Zinczuk, O. E. Piro, J. J. López-González, and A. Ben Altabef, “Structural and vibrational study of 2-(2'-furyl)-1H-imidazole,” Journal of Physical Organic Chemistry, vol. 21, no. 12, pp. 1086–1097, 2008.