INTRODUCTION Today’s high-performance computers allow for the ground-state DFT geometry optimization of molecules with about 100 atoms in a single day. According to Moore’s law, in 25 years the computing power will increase by a factor of 106 [1], and algorithms will improve, so advances in experimental chemistry will not be comparable to advances in computing power. Quantum chemical computing will become an essential tool for studying chemical problems, even surpassing experimental observations, and by 2050, chemical research methods will be very different. Therefore, it is crucial to focus current research and future efforts on theory development. Many chemical processes and biological systems operate based on weak intermolecular interactions [2], so weak intermolecular interactions have been of particular interest in chemistry [3-20], for example, in the field of gas chromatographic separations using capillary columns, where weak intermolecular interactions play an irreplaceable role [21], We found that hydrogen bonding plays a critical role in the retention time of acetic acid (AA) when using gas chromatography to analyze short-chain fatty acids (SCFAs) dissolved in dimethyl sulfoxide (DMSO) [22, 23]. Hydrogen bonding is a frequent weak intermolecular interaction, similar to chemical bonding in that it is directional and saturated, while the bond energy of hydrogen bonding is much smaller than that of chemical bonding. The bonding mechanism of hydrogen bonds is different from chemical bonds in that hydrogen bonds are a special type of the intermolecular force, dipole-dipole. The formation of hydrogen bonds involves electron donors and electron acceptors, where the electron acceptors Mu R et al., Mong. J. Chem., 23(49), 2022, 19-26 https://doi.org/10.5564/mjc.v23i49.1407 A quantum chemical study of the interaction of carboxylic acids with DMSO Mu Ren1, Ao Rigele2, Na Shun3, Narantsogt Natsagdorj1* 1Department of Chemistry School of Mathematics and Natural Sciences, Mongolian National University of Education, Ulaanbaatar, 14191, Mongolia 2Academic Affairs, Hulunbuir University, Hulunbuir, 021108, China 3Agilent Technologies China Co., Ltd., Beijing, 100102, China * Corresponding author: Narantsogt@msue.edu.mn; ORCID ID: 0000-0002-1040-1276 Received: 23 March 2022; revised: 11 May; accepted: 19 May 2022 ABSTRACT Quantum chemical computational methods, which use quantum mechanics and molecular dynamics theory, have developed rapidly in the past few decades, and quantum chemical computation has penetrated almost all fields of chemistry. Hydrogen bonds are ubiquitously common weak intermolecular interactions. Moreover, the bonding mechanism of hydrogen bonds is considered to be different from that of chemical bonding. Because of the difficulty of experimental studies, a more accurate calculation of hydrogen bonding from theory is a more convenient and direct method to understand hydrogen bonding. Density functional theory (DFT) is the most widely used general function in quantum chemical calculations, giving accurate results for most chemical systems. In this paper, the geometries of the hydrogen-bonded dimer complex of acetic acid and DMSO was structurally optimized and potential energy surface was determined. The geometries of four related hydrogen-bonded dimer complexes were fully optimized using the M06-2X/6-311++G (3d, 2p) exchange- correlation functional with DFT-D3(BJ) empirical dispersion correction. We found that hydrogen bonding is a mixture of electrostatic interactions and covalent bonding, and that hydrogen bonding is a kind of force with different percentages of electrostatic and covalent character, rather than a special force independent of chemical bonding. Thus, more clearly defining our inherent classification of forces between substances provides a new perspective for our future study of weak interactions such as hydrogen bonding. Keywords: Quantum chemistry, hydrogen bonding natures, acetic acid, dimethylsulfoxide Mongolian Journal of Chemistry www.mongoliajol.info/index.php/MJC © The Author(s). 2022 Open access. This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 19 https://doi.org/10.5564/mjc.v23i49.1407 https://orcid.org/0000-0002-1040-1276 Mu R et al., Mong. J. Chem., 23(49), 2022, 19-26 must be positively charged hydrogen atoms [5]. A more accurate calculation of hydrogen bonding from theory is an obvious way to understand hydrogen bonding, since experiments can be a complicating factor. Liu et al. performed AIMD calculations on the structure of hydrogen bonds in liquid water [24], Dereka et al. investigated (HF)2 hydrogen bonding using quantum chemical calculations combined with experimental methods [25], and Alkorta et al. described non- covalent interactions in detail using quantum chemical calculations [26]. However, the definition of hydrogen bonding still seems to be vague. The focus of this paper is to present the research methods of quantum chemistry and its utility for studying hydrogen bonds. The hydrogen-bonded dimer complexes were fully optimized using the M06-2X/6- 311++G (3d, 2p) exchange-correlation function with Basis set superposition error (BSSE) [27] of DFT-D3(BJ) empirical dispersion correction [28]. The hydrogen bonding in the active site was determined using the electrostatic potential (ESP) results of the molecule; a quantum-chemical explanation of hydrogen bonding was derived using molecular orbital theory (MO) [29, 30]; confirmation of the directionality of hydrogen bonding was determined using a potential energy surface scan; estimates and comparisons of the energy of intermolecular hydrogen bonds in the geometries of the related 2 hydrogen-bonded dimer complexes (Formic Acid (FA)-DMSO and Trifluoroacetic acid (TFA)-DMSO) were determined using the M06-2X/6- 311++G (3d, 2p) exchange-correlation function with DFT-D3 (BJ) empirical dispersion correction. Based on the above calculations, we propose a supplement to the traditional understanding of hydrogen bonding. EXPERIMENTAL Methods of quantum chemistry: In classical mechanics, Newton’s second law describes the motion of macroscopic objects. In quantum mechanics, Schrödinger proposed the fluctuation equation to describe the state of motion of microscopic particles, which became the most fundamental equation of quantum mechanics and the theoretical basis of quantum chemical calculations [31, 32]. The time-independent Schrödinger equation can be written as: Where Ĥ is a Hamiltonian operator and includes all possible interactions within that system. Ψ is the wave function of the system containing all measurable information. E is the energy of the system and the eigenvalue of the Schrödinger equation. On this basis, the traditional Hartree-Fock equation [33] or the DFT method [34-36] both seek approximations to solve the Schrödinger equation. Molecular modeling: Molecular modeling is simply the theoretical methods and computational techniques used to model or simulate the behavior of molecules [37]. Each model is an approximate solution of the Schrödinger equation, which includes theoretical approaches and basis sets. The rapid development of computer technology and rapid advances in density general function theory and molecular dynamics theory have made theoretical calculations a powerful tool for studying intermolecular forces. In 1998, the Nobel Prize in Chemistry honored quantum chemists Pople and Kohn for their contributions to quantum chemistry. Pople’s approach to quantum chemical calculations and Kohn’s density functional theory are significant milestones in computational chemistry. Mhamad Chrayteh et al. used second-order perturbation theory (MP2), where aldehyde hydrogen atoms strongly interact with water oxygen atoms, in the case of di- and tri-hydrates [38]. Chen Jun et al. used density flooding theory (DFT) and molecular dynamics (MD) simulations to study the hydration mechanism of kaolinite surfaces [39]. Bilonda and Mammino performed calculations of intramolecular hydrogen bonding in matter using different theoretical approaches and basis sets [40]. Quantum chemical calculations played an essential role in all these studies. Theoretically, the use of higher-level theoretical methods and more sophisticated basis sets can lead to more accurate computational results. However, this can also make the computation costly or even unattainable. Ultimately, a trade-off between computational cost and computational accuracy needs to be made when choosing molecular simulations. In this paper, the trend or qualitative forecast of hydrogen bonding in the AA-DMSO hydrogen-bonded dimer complex was found using the Gaussian16 quantum chemistry package [41] in the DFT calculations because the quantitative structure is not predicted. RESULTS AND DISCUSSION Hydrogen bonding sites - Electrostatic natures: The carboxyl group in AA is a combination of two functional groups connected to a single carbon atom, the hydroxyl group (-OH) and the carbonyl group (=O). This combination imparts unique natures to each functional group, including polarity, high electronegativity, and weak acidity. These functional groups are capable of hydrogen bonding by providing and accepting protons [42] and will produce strong hydrogen bonds with the polar non-protic solvent DMSO. Quantum chemical computations provide a good characterization of the intermolecular interactions and the sites of action. In this paper, structure-optimization and frequency calculations were performed separately for the AA molecule and the DMSO molecule using the M062X/6-311++G(3d,2p) [43, 44] method and the DFT- D3BJ [28] dispersion correction, which includes NBO [45] analysis and characterization of the hydrogen bonding sites between the solvent AA molecule and the DMSO molecule. 20 Ĥ Ψ = E Ψ # (1) Mu R et al., Mong. J. Chem., 23(49), 2022, 19-26 The absence of imaginary frequencies is the end of the structural optimization-stable structure. In Fig. 1, the different colors on the molecular surface correspond to the qualitative size of the surface electrostatic potential (ESP). Regions colored in red indicate the ESP is negative: electrons concentrate more in these locations and positively charged particles have stronger interactions. Regions colored in blue indicate the ESP is positive: negatively charged particles more favorably interact in these locations [46]. The maximum positive ESP value of the AA molecule concentrates in the upper portion of the carboxyl functional group (-OH), while the upper region of the S=O in the DMSO molecule has a negative ESP value; these two regions produce hydrogen bonding [47]. Using the ESP results, the AA molecule and DMSO molecule follow the placement of Fig. 1. ESP has been a popular function for revealing electrostatic interactions between substances [48]. We chose the FA-DMSO hydrogen-bonded dimer complex, which is similar to the AA-DMSO hydrogen-bonded dimer complex. We performed ESP computations, and the results were in concordance with the results of Rohmann’s study [49]. We also performed ESP calculations for hydrogen-bonded dimer complexes of AA with N-methyl-pyrrolidone (NMP) and AA- dimethylformamide (DMF) (Fig. 2), both of which are polar, non-protic solvents with similar dipole moments to DMSO [50]. The results and hydrogen bonding interaction sites are consistent with the results of Mu [51] and Safonova [52] et al. ESP accurately shows the hydrogen bonding action sites and effectively illustrates the electrostatic nature of hydrogen bonding. Covalent nature of hydrogen bonds: HOMO/LUMO Interaction: HF is a typical hydrogen bonding system, and the structure-optimized geometry of the (HF)2 hydrogen bonding system according to previous reports [53-55] and using the computations of section 0 is shown in Fig. 3(b). The angles of F-H...F-H are θ1 = 113° and θ2 = 7.6°, and the four atoms not in a straight line, which matches the action site characterized by the ESP results of HF. However, the angle of action does not follow the results of the ESP (Fig. 3 (a)), and our computations are consistent with the studies of Řezac [55] and Howard [56]. Frontier molecular orbital theory, which incorporates highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) calculations, is an important indicator for chemical reactions [57]. The formation of covalent bonds is the redistribution of electrons in the orbitals of frontier molecules when satisfying certain conditions [58, 59]. Following the arrangement of the AA-DMSO hydrogen-bonded dimer complex according to the hydrogen bonding sites, we used the calculation method in section 0 for structure optimization and frequency calculations. During the mapping of its molecular orbitals, we were able to determine the covalent natures of the AA-DMSO hydrogen bonding: the redistribution of electrons in the AA molecule LUMO+3 (Fig. 4(a)) and the DMSO molecule HOMO-1 (Fig. 4(c)) forms HOMO-4 in the AA- DMSO hydrogen-bonded dimer complex (Fig. 4(b)). According to the above results, we cannot explain the hydrogen bonding energy by van der Waals gravity and 21 Fig. 1. The hydrogen bonding site of AA with DMSO Fig. 2. Hydrogen bonding sites for similar hydrogen-bonded dimer complex of AA and DMSO Fig. 3. Hydrogen bonding of (HF)2 (a) Site and direction of (HF)2 interaction presumed by ESP, (b) Hydrogen bonding angle of (HF)2 Fig. 4. Hydrogen bonding associated molecular orbials of the AA-DMSO hydr gen-bon ed dimer complex (a) LUMO+3 of the AA molecule, (b) HOMO-4 of the AA-DSO hydrogen-bonded dimer complex Mu R et al., Mong. J. Chem., 23(49), 2022, 19-26 chemical bonding gravity alone [60], nor can we explain the hydrogen bonding angle of (HF)2 by electrostatic interactions alone, since the van der Waals force is also an electrostatic interaction, and electrostatic interactions will only place the four atoms of (HF)2 in a straight line. The results of Wolters [61] for hydrogen bonding between compounds containing halogens and our quantum chemical calculations indicate the involvement of covalent bonds in the hydrogen bonding interaction. Directionality of hydrogen bonding: Potential energy surface (PES) determination is one of the central concepts in quantum chemical calculations. According to the Born-Oppenheimer approximation, a PES is a mathematical or graphical relationship between the energy of a molecule (or collection of molecules) and its geometry [62], and we can represent the energy of a system as a function of geometric coordinates. We performed a PES scan of the AA-DMSO hydrogen- bonded dimer complex using the calculation method in Section 0. In order to obtain the surface in three- dimensional space to identify the directionality of the hydrogen bonds, the quantum chemical calculation selects only two variables: the hydrogen bond length L and its angle θ (Fig. 5(a)), so the energy E of the system is a function of L and θ. Fig. 5(b) shows the 3D PES of the system energy variation with hydrogen bond length L and angle θ. When the angle θ tends to 180° and the bond length L tends to 1.6 Å, the energy of the AA-DMSO hydrogen- bonded dimer complex decreases, indicating that the AA-DMSO hydrogen-bonded dimer complex is more stable and the hydrogen bond strength is stronger [63- 66], further indicating the hydrogen bond is directional. Directionality is an essential indicator of covalent bonds. The strength of covalent bonds between atoms is strongly dependent on an angle [67]. The three- dimensional PES of the AA-DMSO hydrogen-bonded dimer complex demonstrate the directionality of hydrogen bonds and also indicate that hydrogen bonds have covalent nature. Estimates of the energy of intermolecular hydrogen bonds: Estimates of the energy of intermolecular hydrogen bonds can be characterized in several ways, of which the binding energy (BE) is one of the most important quantities. The BE, as used herein, is equivalent to the “interaction energy” between two monomers in complex geometry. Quantum chemistry has become a routine and reliable method for estimating the BE of various hydrogen bonds [68]. For example, a study evaluated the performance of DFT generalized functions in calculating the binding energy of furan clusters. The results showed that the functionals M05-2X and M06 are recommended for further affordable investigations of the furan clusters [69]. In this section, the geometries of all three hydrogen- bonded dimer complexes (Table 1) were fully optimized using the M06-2X exchange-correlation functional with DFT-D3(BJ) empirical dispersion correction. 22 Fig. 5. Hydrogen bonding and three-dimensional PES of the AA-DMSO hydrogen-bonded dimer complex (a) Hydrogen bond length L and angle θ of the AA-DMSO hydrogen-bonded dimer complex (b) Hydrogen bonding PES of the AA-DMSO hy- drogen-bonded dimer complex Fig. 6. ESP for FA, AA, and TFA and the NPA charge of the hydrogen atoms of FA, AA, and TFA involved in hydrogen bonding Fig. 7. Hydrogen bond lengths of FA-DMAO, AA-DMSO, and TFA-DMSO Project FA-DMSO AA-DMSO TFA-DMSO NPA charge (H atom)* 0.490 0.492 0.502 NPA charge (O atom)* -1.015 -1.015 -1.015 Hydrogen bond length (Å) 1.587 1.618 1.490 Hydrogen bond BE (kJ/mol) -61.743 -59.079 -73.955 Table 1. NPA charges of hydrogen bond atoms, hydrogen bond length, and energy Notice: * NPA charge (H atom) is the NPA charge of the hydrogen atoms of FA, AA and TFA involved in hydrogen bonding. * NPA charge (O atom) is the NPA charge of the oxygen atoms of DMSO involved in hydrogen bonding. Mu R et al., Mong. J. Chem., 23(49), 2022, 19-26 For fully characterizing hydrogen bond interactions, we have utilized many other quantum chemical calculations in addition to BEs, such as the hydrogen bond length, the NPA charges of the atoms at the ends of the hydrogen bonds, and ESP. Using correlations between various well-known hydrogen bond-related quantum chemical calculations and estimated BE values, we have shown that each of these quantum chemical calculations has its own ability to characterize the strength of hydrogen bond interactions. Our calculations have revealed that the higher the NPA charge difference of the hydrogen bond terminal atoms, the shorter the bond length. Additionally, larger hydrogen bond energy is observed with lower electron density of the electron acceptor as determined by ESP. We found from the strong hydrogen-bonded dimeric complexes that the electron acceptor of the hydrogen atom all contain atoms or groups of atoms with high electronegativity. CONCLUSIONS In the work presented here, we have performed DFT calculations for the AA-DMSO hydrogen-bonded dimer complex using ESP, HOMO/LUMO, PES, and estimated and compared the energy of hydrogen bonds of four related hydrogen-bonded dimer complexes. The quantum chemical calculations have revealed the natures of hydrogen bonding in the AA-DMSO hydrogen- bonded dimer complex: electrostatic interaction is the main factor in the generation of hydrogen bonds, and the ESP can characterize the site of action of hydrogen bonds; the redistribution of electrons in molecular orbitals provides hydrogen bonding with a significant degree of chemical bonding character, which is one of the reasons why the hydrogen bonding energy is between chemical bonding and van der Waals forces, and the directionality of hydrogen bonding revealed by PES, and estimates of the energy of intermolecular hydrogen bonds is also consistent with this finding. Different interactions link the atoms in a compound, and these links do not have a clear demarcation line, but can be considered a mixture of multiple interactions. Therefore, on the premise that hydrogen bonding is a mixture of electrostatic interactions and covalent character, further in-depth study of the percentage of each interaction force is of more practical research significance than classifying hydrogen bonds. REFERENCES 1. Grimme S., Schreiner P.R. (2018) Computational chemistry: The fate of current methods and future challenges. Angewandte Chemie International Edition, 57(16), 4170-6. https://doi.org/10.1002/anie.201709943 2. Tien H.T., Ottova-Leitmannova A. (2003) Planar lipid bilayers (BLM’s) and their applications. Elsevier. 3. Latimer W.M., Rodebush W.H. (1920) Polarity and ionization from the standpoint of the Lewis theory of valence. Journal of the American Chemical Society, 42(7), 1419-33. https://doi.org/10.1021/ja01452a015 4. Van Doren J.B. (1967) Valence and the structure of atoms and molecules (Lewis, GN). ACS Publications. https://doi.org/10.1021/ed044pA82 5. Arunan E., Desiraju G.R., Klein R.A., Sadlej J, Scheiner S., et al. (2011) Definition of the hydrogen bond. Pure and Applied Chemistry, 83(8),1637- 41. https://doi.org/10.1351/PAC-REC-10-01-02 6. Etter M.C., MacDonald J.C., Bernstein J. (1990) Graph-set analysis of hydrogen-bond patterns in organic crystals. Acta Crystallographica Section B: Structural Science, 46(2), 256-62. https://doi.org/10.1107/S0108768189012929 7. Taylor R., Kennard O. (1984) Hydrogen- bond geometry in organic crystals. Accounts of Chemical Research, 17(9), 320-6. https://doi.org/10.1021/ar00105a004 8. Stillinger F.H. (1980) Water revisited. Science, 209(4455), 451-7. https://doi.org/10.1126/science.209.4455.451 9. Diercksen G.H., Kraemer W.P,. Roos B.O. (1975) SCF-CI studies of correlation effects on hydrogen bonding and ion hydration. Theoretica Chimica Acta, 36(4), 249-74. https://doi.org/10.1007/BF00549690 10. Carroll M.T., Bader R.F. (1988) An analysis of the hydrogen bond in BASE-HF complexes using the theory of atoms in molecules. Molecular Physics, 65(3), 695-722. https://doi.org/10.1080/00268978800101351 11. Huggins M.L. (1971) 50 Years of hydrogen bond theory. Angewandte Chemie International Edition in English, 10(3), 147-52. https://doi.org/10.1002/anie.197101471 12. Pauling L., Brockway L.O. The structure of the carboxyl group: I. The investigation of formic acid by the diffraction of electrons. Proceedings of the National Academy of Sciences of the U.S.A., 20(6), 336. https://doi.org/10.1073/pnas.20.6.336 13. Wu S. (2017) Polymer interface and adhesion. Routledge. https://doi.org/10.1201/9780203742860 14. Riel A.M.S., Rowe R.K., Ho E.N., Carlsson A.C.C., Rappé A.K., Berryman O.B., et al. (2019) Hydrogen bond enhanced halogen bonds: A synergistic interaction in chemistry and biochemistry. Accounts of Chemical Research, 52(10), 2870- 80. https://doi.org/10.1021/acs.accounts.9b00189 15. Karas L.J., Wu C.H., Das R., Wu J.I.C. (2020) Hydrogen bond design principles. Wiley Interdisciplinary Reviews. Computational Molecular Science, 10(6), e1477. https://doi.org/10.1002/wcms.1477 16. Gibb B.C. (2020) The centenary (maybe) of the hydrogen bond. Nature Publishing Group, 23 https://doi.org/10.1002/anie.201709943 https://doi.org/10.1021/ja01452a015 https://doi.org/10.1021/ed044pA82 https://doi.org/10.1351/PAC-REC-10-01-02 https://doi.org/10.1107/S0108768189012929 https://doi.org/10.1021/ar00105a004 https://doi.org/10.1126/science.209.4455.451 https://doi.org/10.1007/BF00549690 https://doi.org/10.1080/00268978800101351 https://doi.org/10.1002/anie.197101471 https://doi.org/10.1073/pnas.20.6.336 https://doi.org/10.1201/9780203742860 https://doi.org/10.1021/acs.accounts.9b00189 https://doi.org/10.1002/wcms.1477 Mu R et al., Mong. J. Chem., 23(49), 2022, 19-26 https://doi.org/10.1038/s41557-020-0524-2 17. Niemann T., Strate A., Ludwig R., Zeng H.J., Menges F.S., Johnson M.A. (2019) Cooperatively enhanced hydrogen bonds in ionic liquids: Closing the loop with molecular mimics of hydroxy- functionalized cations. Physical Chemistry Chemical Physics, 21(33), 18092-8. https://doi.org/10.1039/C9CP03300A 18. Mishra K.K., Borish K., Singh G., Panwaria P., Metya S., Madhusudhan M., et al. Observation of an unusually large IR red-shift in an unconventional S-H••• S hydrogen-bond. The Journal of Physical Chemistry Letters, 12(4), 1228-35. https://doi.org/10.1021/acs.jpclett.0c03183 19. Mishra K.K., Singh S.K., Ghosh P., Ghosh D., Das A. (2017) The nature of selenium hydrogen bonding: Gas phase spectroscopy and quantum chemistry calculations. Physical Chemistry Chemical Physics, 19(35), 24179-87. https://doi.org/10.1039/C7CP05265K 20. Scheiner S. (2020) The hydrogen bond: A hundred years and counting. Journal of the Indian Institute of Science, 100(1), 61-76. https://doi.org/10.1007/s41745-019-00142-8 21. Dai Z.Y., Chen W.H., Zhou H., Yang H.Y., (2000) Studies on quantitative chromatographic retention-structure relationships. Chinese Journal of Chromatography, 18(2), 125-127. PMID: 12541586. 22. Ren M., Natsagdorj N., Shun N. (2022) Influence and mechanism of polar solvents on the retention time of short-chain fatty acids in gas chromatography. Separations, 9(5), 124. https://doi.org/10.3390/separations9050124 23. Mu R., Na S., Narantsogt N., LI C.J. (2022) Effect of dimethyl sulfoxide on the retention time of acetic acid in gas chromatography and Its mechanism. Chinese Journal of Applied Chemistry, 39(5), 852- 854. https://doi.org/10.19894/j.issn.1000-0518.210175 24. Liu J., He X., Zhang J.Z., Qi L.W. (2018) Hydrogen- bond structure dynamics in bulk water: Insights from ab initio simulations with coupled cluster theory. Chemical Science, 9(8), 2065-73. https://doi.org/10.1039/C7SC04205A 25. Dereka B., Yu Q, Lewis N.H., Carpenter W.B., Bowman J.M., Tokmakoff A. (2021) Crossover from hydrogen to chemical bonding. Science, 371(6525), 160-4. https://doi.org/10.1126/science.abe1951 26. Alkorta I., Elguero J., Frontera A. (2020) Not only hydrogen bonds: Other noncovalent interactions. Crystals, 10(3), 180. https://doi.org/10.3390/cryst10030180 27. Boys S.F., Bernardi F. (1970) The calculation of small molecular interactions by the differences of separate total energies. Some procedures with reduced errors. Molecular Physics, 19(4), 553- 66. https://doi.org/10.1080/00268977000101561 28. Grimme S., Ehrlich S., Goerigk L. (2011) Effect of the damping function in dispersion corrected density functional theory. Journal of Computational Chemistry, 32(7), 1456-65. https://doi.org/10.1002/jcc.21759 29. Mulliken R.S. (1932) The interpretation of band spectra part III. Electron quantum numbers and states of molecules and their atoms. Reviews of modern physics, 4(1), 1. https://doi.org/10.1103/RevModPhys.4.1 30. Hund F. (1933) Allgemeine quantenmechanik des atom-und molekelbaues. Quantentheorie. Springer, 561-694. https://doi.org/10.1007/978-3-642-52619-0_4 31. Bartlett R.J., Musiał M. (2007) Coupled-cluster theory in quantum chemistry. Reviews of Modern Physics. 79(1), 291. https://doi.org/10.1103/RevModPhys.79.291 32. Lowdin P.O. (1970) Advances in quantum chemistry. Academic Press. 33. Slater J. (1951) Magnetic effects and the Hartree- Fock equation. Physical Review, 82(4), 538. https://doi.org/10.1103/PhysRev.82.538 34. Hohenberg P., Kohn W. (1964) Inhomogeneous electron gas. Physical Review, 136(3B), B864. https://doi.org/10.1103/PhysRev.136.B864 35. Kohn W., Sham L.J. (1965) Self-consistent equations including exchange and correlation effects. Physical Review, 140(4A), A1133. https://doi.org/10.1103/PhysRev.140.A1133 36. Pauli W. (1933) Die allgemeinen prinzipien der wellenmechanik. Quantentheorie. Springer, 83- 272. https://doi.org/10.1007/978-3-642-52619-0_2 37. Leach A.R., Leach A.R. (2001) Molecular modelling: Principles and applications. Pearson education. 38. Chrayteh M., Huet T.R., Dréan P. (2020) Microsolvation of myrtenal studied by microwave spectroscopy highlights the role of quasi-hydrogen bonds in the stabilization of its hydrates. The Journal of Chemical Physics, 153(10), 104304. https://doi.org/10.1063/5.0019957 39. Chen J., Min F.f., Liu L., Liu C., (2019) Mechanism research on surface hydration of kaolinite, insights from DFT and MD simulations. Applied Surface Science, 476, 6-15. https://doi.org/10.1016/j.apsusc.2019.01.081 40. Bilonda M.K., Mammino L. (2017) Intramolecular hydrogen bonds in conformers of quinine and quinidine: An HF, MP2 and DFT study. Molecules, 22(2), 245. https://doi.org/10.3390/molecules22020245 41. Frisch M.J., Trucks G.W., Schlegel H.B., Scuseria G.E., Robb M.A., et al. (2016) Gaussian 16 Rev. C.01. Wallingford, CT2016. 42. Klecker C., Nair L.S. (2017) Matrix chemistry 24 https://doi.org/10.1038/s41557-020-0524-2 https://doi.org/10.1039/C9CP03300A https://doi.org/10.1021/acs.jpclett.0c03183 https://doi.org/10.1039/C7CP05265K https://doi.org/10.1007/s41745-019-00142-8 https://doi.org/10.3390/separations9050124 https://doi.org/10.19894/j.issn.1000-0518.210175 https://doi.org/10.1039/C7SC04205A https://doi.org/10.1126/science.abe1951 https://doi.org/10.3390/cryst10030180 https://doi.org/10.1080/00268977000101561 https://doi.org/10.1002/jcc.21759 https://doi.org/10.1103/RevModPhys.4.1 https://doi.org/10.1007/978-3-642-52619-0_4 https://doi.org/10.1103/RevModPhys.79.291 https://doi.org/10.1103/PhysRev.82.538 https://doi.org/10.1103/PhysRev.136.B864 https://doi.org/10.1103/PhysRev.140.A1133 https://doi.org/10.1007/978-3-642-52619-0_2 https://doi.org/10.1063/5.0019957 https://doi.org/10.1016/j.apsusc.2019.01.081 https://doi.org/10.3390/molecules22020245 Mu R et al., Mong. J. Chem., 23(49), 2022, 19-26 controlling stem cell behavior. Biology and Engineering of Stem Cell Niches, 195-213. https://doi.org/10.1016/B978-0-12-802734- 9.00013-5 43. Frisch M.J., Pople J.A., Binkley J.S. (1984) Self consistent molecular orbital methods 25. Supplementary functions for Gaussian basis sets. The Journal of Chemical Physics, 80(7), 3265-9. https://doi.org/10.1063/1.447079 44. Hariharan P.C., Pople J.A. (1973) The influence of polarization functions on molecular orbital hydrogenation energies. Theoretica Chimica Acta, 28(3), 213-22. https://doi.org/10.1007/BF00533485 45. Glendening E., Reed A., Carpenter J., Weinhold F. (1998) NBO Version 3.1. 46. Singh D.K., Jagannathan R., Khandelwal P., Abraham P.M., Poddar P. (2013) In situ synthesis and surface functionalization of gold nanoparticles with curcumin and their antioxidant properties: An experimental and density functional theory investigation. Nanoscale, 5(5), 1882-93. https://doi.org/10.1039/c2nr33776b 47. Orozco M., Alhambra C., Barril X., López J.M., Busquets M.A., Luque F.J. (1996) Theoretical methods for the representation of solvent. Molecular Modeling Annual, 2(1), 1-15. https://doi.org/10.1007/s0089460020001 48. Lu T., Chen Q. (2020) Van der Waals potential: An important complement to molecular electrostatic potential in studying intermolecular interactions. Journal of Molecular Modeling, 26(11), 1-9. https://doi.org/10.1007/s00894-020-04577-0 49. Rohmann K., Kothe J., Haenel M.W., Englert U., Hölscher M., Leitner W. (2016) Hydrogenation of CO2 to formic acid with a highly active ruthenium acriphos complex in DMSO and DMSO/water. Angewandte Chemie International Edition, 55(31), 8966-9. https://doi.org/10.1002/anie.201603878 50. Leonard J., Lygo B., Procter G. (2013) Advanced practical organic chemistry. CRC press. https://doi.org/10.1201/b13708 51. Mu L., Shi Y., Chen L., Ji T., Yuan R., Wang H., et al. [N-Methyl-2-pyrrolidone][C1-C4 carboxylic acid]: A novel solvent system with exceptional lignin solubility. Chemical Communications, 51(70) 13554-7. https://doi.org/10.1039/C5CC04191K 52. Safonova L., Pryakhin A., Shmukler L., Fadeeva Y.A. (2009) NMR studies of N, N-dimethylformamide mixtures with acetic acid and ethanol. Russian Journal of General Chemistry, 79(2), 246-51. https://doi.org/10.1134/S1070363209020133 53. Van der Lubbe S.C., Fonseca Guerra C. (2019) The nature of hydrogen bonds: A delineation of the role of different energy components on hydrogen bond strengths and lengths. Chemistry-An Asian Journal, 14(16), 2760-9. https://doi.org/10.1002/asia.201900717 54. Herschlag D., Pinney M.M. (2018) Hydrogen bonds: Simple after all?. Biochemistry, 57(24), 3338-3352. https://doi.org/10.1021/acs.biochem.8b00217 55. Řezáč J., Hobza P. (2014) Ab initio quantum mechanical description of noncovalent interactions at its limits: Approaching the experimental dissociation energy of the HF dimer. Journal of Chemical Theory and Computation, 10(8), 3066- 3073. https://doi.org/10.1021/ct500047x 56. Howard J.C., Gray J.L., Hardwick A.J., Nguyen L.T., Tschumper G.S. (2014) Getting down to the fundamentals of hydrogen bonding: Anharmonic vibrational frequencies of (HF) 2 and (H2O) 2 from ab initio electronic structure computations. Journal of Chemical Theory and Computation, 10(12), 5426-35. https://doi.org/10.1021/ct500860v 57. Fukui K. (1971) Recognition of stereochemical paths by orbital interaction. Accounts of Chemical Research, 4(2), 57-64. https://doi.org/10.1021/ar50038a003 58. Frey P.A., Whitt S.A., Tobin J.B. (1994) A low- barrier hydrogen bond in the catalytic triad of serine proteases. Science, 264(5167),1927-30. https://doi.org/10.1126/science.7661899 59. Cleland W., Kreevoy M.M. (1994) Low-barrier hydrogen bonds and enzymic catalysis. Science, 264(5167), 1887-1890. https://doi.org/10.1126/science.8009219 60. Hermansson K. (2017) Hydrogen bonds- on the move. HBOND2017, 10-14. https://doi.org/10.1016/B978-0-12-820042- 1.00015-8 61. Wolters L.P., Bickelhaupt F.M. (2012) Halogen bonding versus hydrogen bonding: A molecular orbital perspective. Chemistry Open, 1(2), 96-105. https://doi.org/10.1002/open.201100015 62. Lewars E.G. (2016) The concept of the potential energy surface. Computational Chemistry, 9-49. https://doi.org/10.1007/978-3-319-30916-3_2 63. Ma J., Liu Y., Xie M.X. (2012) Interaction between three isoflavones and different isomers of human serum albumin. Spectroscopy and Spectral Analysis, 32(1), 1-6. h t t p : / / d x . c h i n a d o i . c n / 1 0 . 3 9 6 4 / j . i s s n . 1 0 0 0 - 0593(2012)01-0001-06 64. Fang R., Leng X.J., Wu X., LI Q., Hao R.F., Ren F.Z., et al. (2012) Intermolecular hydrogen bond between protein and flavonoid and its contribution to the stability of the flavonoids. Spectroscopy and Spectral Analysis, 32(1), 108-112. https://doi.org/10.3964/j.issn.1000-0593(2012)01- 0108-05 65. Zhang L., Wei L., Zhai S., Zhao D., Sun J., An Q. (2020) Hydrogen bond promoted thermal stability enhancement of acetate based ionic liquid. Chinese Journal of Chemical Engineering, 28(5), 1293-301. 25 https://doi.org/10.1016/B978-0-12-802734-9.00013-5 https://doi.org/10.1016/B978-0-12-802734-9.00013-5 https://doi.org/10.1063/1.447079 https://doi.org/10.1007/BF00533485 https://doi.org/10.1039/c2nr33776b https://doi.org/10.1007/s0089460020001 https://doi.org/10.1007/s00894-020-04577-0 https://doi.org/10.1002/anie.201603878 https://doi.org/10.1201/b13708 https://doi.org/10.1039/C5CC04191K https://doi.org/10.1134/S1070363209020133 https://doi.org/10.1002/asia.201900717 https://doi.org/10.1021/acs.biochem.8b00217 https://doi.org/10.1021/ct500047x https://doi.org/10.1021/ct500860v https://doi.org/10.1021/ar50038a003 https://doi.org/10.1126/science.7661899 https://doi.org/10.1126/science.8009219 https://doi.org/10.1016/B978-0-12-820042-1.00015-8 https://doi.org/10.1016/B978-0-12-820042-1.00015-8 https://doi.org/10.1002/open.201100015 https://doi.org/10.1007/978-3-319-30916-3_2 http://dx.chinadoi.cn/10.3964/j.issn.1000-0593(2012)01-0001-06 http://dx.chinadoi.cn/10.3964/j.issn.1000-0593(2012)01-0001-06 https://doi.org/10.3964/j.issn.1000-0593(2012)01-0108-05 https://doi.org/10.3964/j.issn.1000-0593(2012)01-0108-05 Mu R et al., Mong. J. Chem., 23(49), 2022, 19-26 https://doi.org/10.1016/j.cjche.2020.02.019 66. Badbedast M., Izadyar M., Gholizadeh M. (2019) The investigation of the structure and stability of catechol complexes with nitrate and hydrogen sulphate ions through theoretical study of hydrogen bond. Nashrieh Shimi va Mohandesi Shimi Iran, 37(4), 287-300. 67. Gridneva I., Milman Y.V., Trefilov V. (1969) On the Mechanical Properties of Crystals with Covalent Bond. Physica Status Solidi (b), 36(1), 59-67. https://doi.org/10.1002/pssb.19690360104 68. Stasyuk O.A., Sedlak R., Guerra C.F., Hobza P. (2018) Comparison of the DFT-SAPT and canonical EDA schemes for the energy decomposition of various types of noncovalent interactions. Journal of Chemical Theory and Computation, 14(7). https://doi.org/10.1021/acs.jctc.8b00034 69. Malloum A., Conradie J. (2021) Structures, binding energies and non-covalent interactions of furan clusters. Journal of Molecular Graphics and Modelling, 111, 108102. https://doi.org/10.1016/j.jmgm.2021.108102 26 https://doi.org/10.1016/j.cjche.2020.02.019 https://doi.org/10.1002/pssb.19690360104 https://doi.org/10.1021/acs.jctc.8b00034 https://doi.org/10.1016/j.jmgm.2021.108102