1
|
Ahirwar MB, Gadre SR, Deshmukh MM. Molecular Tailoring Approach for the Direct Estimation of Individual Noncovalent Interaction Energies in Molecular Systems. J Phys Chem A 2024; 128:6099-6115. [PMID: 39037864 DOI: 10.1021/acs.jpca.4c01176] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/24/2024]
Abstract
The noncovalent interactions (NCIs) are omnipresent in chemistry, physics, and biology. The study of such interactions offers insights into various physicochemical phenomena. Some indirect approaches proposed in the literature for exploring the NCIs are briefly reviewed in Section 1 of this Perspective. These include: (i) Shift in the stretching frequency of an X-Y bond involved in X-Y···Z interaction. (ii) Topological analysis of molecular electron density. (iii) Empirical equations derived employing experimental and theoretical quantities. However, a direct method for estimating individual intramolecular/intermolecular interaction energies has been conspicuous by its absence from the literature. We have developed a molecular tailoring approach (MTA)-based method enabling a direct and reliable estimation of the energy of intra- as well as intermolecular interactions. This method offers a direct and reliable estimation of these interactions, in particular of the hydrogen bonds (HB) in molecules/weakly bound clusters along with the respective cooperativity contribution. In Section 2, the basis of our method is discussed, along with some illustrative examples. The application of this method to a variety of molecules and clusters, with a special emphasis on estimating the HB energy along with the energy of other NCIs is presented in Section 3. Section 4 discusses some computational strategies for applying our method to large molecular clusters. The last Section provides a summary and a discussion on future developments.
Collapse
Affiliation(s)
- Mini Bharati Ahirwar
- Department of Chemistry, Dr. Harisingh Gour Vishwavidyalaya (A Central University), Sagar 470003, India
| | - Shridhar R Gadre
- Department of Scientific Computing, Modelling, & Simulation, Savitribai Phule Pune University, Pune 411007, India
- Department of Chemistry, Savitribai Phule Pune University, Pune 411007, India
| | - Milind M Deshmukh
- Department of Chemistry, Dr. Harisingh Gour Vishwavidyalaya (A Central University), Sagar 470003, India
| |
Collapse
|
2
|
Ahirwar MB, Khire SS, Gadre SR, Deshmukh MM. Hydrogen bond energy estimation (H-BEE) in large molecular clusters: A Python program for quantum chemical investigations. J Comput Chem 2024; 45:274-283. [PMID: 37792345 DOI: 10.1002/jcc.27237] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/13/2023] [Accepted: 09/18/2023] [Indexed: 10/05/2023]
Abstract
A procedure, derived from the fragmentation-based molecular tailoring approach (MTA), has been proposed and extensively applied by Deshmukh and Gadre for directly estimating the individual hydrogen bond (HB) energies and cooperativity contributions in molecular clusters. However, the manual fragmentation and high computational cost of correlated quantum chemical methods make the application of this method to large molecular clusters quite formidable. In this article, we report an in-house developed software for automated hydrogen bond energy estimation (H-BEE) in large molecular clusters. This user-friendly software is essentially written in Python and executed on a Linux platform with the Gaussian package at the backend. Two approximations to the MTA-based procedure, viz. the first spherical shell (SS1) and the Fragments-in-Fragments (Frags-in-Frags), enabling cost-effective, automated evaluation of HB energies and cooperativity contributions, are also implemented in this software. The software has been extensively tested on a variety of molecular clusters and is expected to be of immense use, especially in conjunction with correlated methods such as MP2, CCSD(T), and so forth.
Collapse
Affiliation(s)
- Mini Bharati Ahirwar
- Department of Chemistry, Dr. Harisingh Gour Vishwavidyalaya (A Central University), Sagar, India
| | | | - Shridhar R Gadre
- Department of Scientific Computing, Modelling & Simulation, Savitribai Phule Pune University, Pune, India
- Department of Chemistry, Savitribai Phule Pune University, Pune, India
| | - Milind M Deshmukh
- Department of Chemistry, Dr. Harisingh Gour Vishwavidyalaya (A Central University), Sagar, India
| |
Collapse
|
3
|
Ahirwar MB, Deshmukh MM. Fragments-in-fragments method for efficient and reliable estimates of individual hydrogen bond energies in large molecular clusters. J Comput Chem 2023. [PMID: 37191018 DOI: 10.1002/jcc.27133] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/04/2023] [Revised: 05/01/2023] [Accepted: 05/03/2023] [Indexed: 05/17/2023]
Abstract
The knowledge of individual hydrogen bond (HB) strength in molecular clusters is indispensable to get insights into the bulk properties of condensed systems. Recently, we have developed the molecular tailoring approach based (MTA-based) method for the estimation of individual HB energy in molecular clusters. However, the direct use of this MTA-based method to large molecular clusters becomes progressively difficult with the increase in the size of a cluster. To overcome this caveat, herein, we propose the use of linear scaling method (such as the original MTA method) for the estimation of single-point (SP) energies of large-sized parent molecular cluster and their respective fragments. Because the fragments of the MTA-based method, for the estimation of HB energy, are further fragmented, this proposed strategy is called as Fragments-in-Fragments (Frags-in-Frags) method. The SP energies of fragments and parent cluster calculated by the Frags-in-Frags approach were utilized to estimate the individual HB energy. The estimated individual HB energies, in various molecular clusters, by Frags-in-Frags method are found to be in excellent linear agreement with their MTA-based counterparts (R2 = 0.9975 of 348 data points). The difference being less than 0.5 kcal/mol in most of the cases. Furthermore, RMSD is 0.43 kcal/mol, MAE is 0.33 kcal/mol, and the standard deviation is 0.44 kcal/mol. Importantly, the Frags-in-Frags method not only enables the reliable estimation of HB energy in large molecular clusters but also requires less computational time and can be possible even with off-the-shelf hardware.
Collapse
Affiliation(s)
- Mini Bharati Ahirwar
- Department of Chemistry, Dr. Harisingh Gour Vishwavidyalaya (A Central University), Sagar, India
| | - Milind M Deshmukh
- Department of Chemistry, Dr. Harisingh Gour Vishwavidyalaya (A Central University), Sagar, India
| |
Collapse
|
4
|
Ahirwar MB, Gadre SR, Deshmukh MM. On the Short-Range Nature of Cooperativity in Hydrogen-Bonded Large Molecular Clusters. J Phys Chem A 2023; 127:4394-4406. [PMID: 37186960 DOI: 10.1021/acs.jpca.3c00359] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/17/2023]
Abstract
The variation in the hydrogen bond (HB) strength has considerable consequences on the physicochemical properties of molecular clusters. Such a variation mainly arises due to the cooperative/anti-cooperative networking effect of neighboring molecules connected by HBs. In the present work, we systematically study the effect of neighboring molecules on the strength of an individual HB and the respective cooperativity contribution toward each of them in a variety of molecular clusters. For this purpose, we propose a use of a small model of a large molecular cluster called the spherical shell-1 (SS1) model. This SS1 model is constructed by placingg the spheres of an appropriate radius centered on X and Y atoms of the X-H···Y HB under consideration. The molecules falling within these spheres constitute the SS1 model. Utilizing this SS1 model, the individual HB energies are calculated within the molecular tailoring approach-based framework and the results are compared with their actual counterparts. It is found that the SS1 is a reasonably good model of large molecular clusters, providing 81-99% of the total HB energy estimated using the actual molecular clusters. This in turn suggests that the maximum cooperativity contribution toward a particular HB is due to the fewer number of molecules (in the SS1 model) directly interacting with two molecules involved in its formation. We further demonstrate that the remaining part of the energy or cooperativity (∼1 to 19%) is captured by the molecules falling in the second spherical shell (SS2) centered on the hetero-atom of the molecules in the SS1 model. The effect of increasing size of a cluster on the strength of a particular HB, calculated by the SS1 model, is also investigated. The calculated value of the HB energy remains unchanged with the increase in the size of a cluster, emphasizing the short-ranged nature of the HB cooperativity in neutral molecular clusters.
Collapse
Affiliation(s)
- Mini Bharati Ahirwar
- Department of Chemistry, Dr. Harisingh Gour Vishwavidyalaya (A Central University), Sagar 470003, India
| | - Shridhar R Gadre
- Department of Scientific Computing, Modelling, and Simulation, Savitribai Phule Pune University, Pune 411007, India
| | - Milind M Deshmukh
- Department of Chemistry, Dr. Harisingh Gour Vishwavidyalaya (A Central University), Sagar 470003, India
| |
Collapse
|
5
|
Ahirwar MB, Deshmukh MM. Two-Step ONIOM Method for the Accurate Estimation of Individual Hydrogen Bond Energy in Large Molecular Clusters. J Phys Chem A 2023; 127:1219-1232. [PMID: 36705264 DOI: 10.1021/acs.jpca.2c08087] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Abstract
The study of molecular clusters to understand the properties of condensed systems has been the subject of immense interest. To get insight into these properties, the knowledge of various noncovalent interactions present in these molecular clusters is indispensable. Our recently developed molecular tailoring approach-based (MTA-based) method for the estimation of the individual hydrogen bond (HB) energy in molecular clusters is useful for this purpose. However, the direct application of this MTA-based method becomes progressively difficult with the increase in the size of the cluster. This is because of the difficulty in the evaluation of single-point energy at the correlated level of theory. To overcome this caveat, herein, we propose a two-step method within the our own N-layer integrated molecular orbital molecular mechanics (ONIOM) framework. In this method, the HB energy evaluated by the MTA-based method employing the actual molecular cluster at a low Hartree-Fock (HF) level of theory is added to the difference in the HB energies evaluated by the MTA-based method, employing an appropriate small model system, called the shell-1 model, calculated at high (MP2) and low (HF) levels of theory. The shell-1 model of a large molecular cluster is made up of only a few molecules that are in direct contact (by a single HB) with the two molecules involved in the formation of an HB under consideration. We tested this proposed two-step ONIOM method to estimate the individual HB energies in various molecular clusters, viz., water (Wn, n = 10-16, 18 and 20), (H2O2)12, (H2O3)8, (NH3)n and strongly interacting (HF)15 and (HF)m(W)n clusters. Furthermore, these estimated individual HB energies by the ONIOM method are compared with those calculated by the MTA-based method using actual molecular clusters. The estimated individual HB energies by the ONIOM method, in all these clusters, are in excellent linear one-to-one agreement (R2 = 0.9996) with those calculated by the MTA-based method using actual molecular clusters. Furthermore, the small values of root-mean-square deviation (0.06), mean absolute error (0.04), |ΔEmax| (0.21) and Sε (0.06) suggest that this two-step ONIOM method is a pragmatic approach to provide accurate estimates of individual HB energies in large molecular clusters.
Collapse
Affiliation(s)
- Mini Bharati Ahirwar
- Department of Chemistry, Dr. Harisingh Gour Vishwavidyalaya (A Central University), Sagar470003, India
| | - Milind M Deshmukh
- Department of Chemistry, Dr. Harisingh Gour Vishwavidyalaya (A Central University), Sagar470003, India
| |
Collapse
|
6
|
Patkar D, Bharati Ahirwar M, Deshmukh MM. A Tug of War between the Self- and Cross-associating Hydrogen Bonds in Neutral Ammonia-Water Clusters: Energetic Insights by Molecular Tailoring Approach. Chemphyschem 2022; 23:e202200476. [PMID: 36127809 DOI: 10.1002/cphc.202200476] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/05/2022] [Revised: 09/09/2022] [Indexed: 01/05/2023]
Abstract
In the present work, the energies of various types of individual HBs observed in neutral (NH3 )m (H2 O)n , (m+n=2 to 7) clusters were estimated using the molecular tailoring approach (MTA)-based method. The calculated individual HB energies suggest that the O-H…N HBs are the strongest (1.21 to 12.49 kcal mol-1 ). The next ones are the O-H…O (3.97 to 9.30 kcal mol-1 ) HBs. The strengths of N-H…N (1.09 to 5.29 kcal mol-1 ) and N-H…O (2.85 to 5.56 kcal mol-1 ) HBs are the weakest. The HB energies in dimers also follow this rank ordering. However, the HB energies in dimers are much smaller than those obtained by the MTA-based method due to the loss in cooperativity contribution in the dimers. Thus, the calculated cooperativity contributions, for different types of HBs, fall in the range 0.64 to 5.73 kcal mol-1 . We wish to emphasize based on the energetic rank ordering obtained by the MTA-based method that the O-H of water is a better HB donor than the N-H of ammonia. The reasons for the observed energetic rank ordering are two folds: (i) intrinsically stronger O-H…N HBs than the O-H…O ones as revealed by dimer energies and (ii) the higher cooperativity contribution in the former than the later ones. Indeed, the MTA-based method is useful in providing the missing energetic rank ordering of various type of HBs in neutral (NH3 )m (H2 O)n clusters, in the literature.
Collapse
Affiliation(s)
- Deepak Patkar
- Department of Chemistry, Dr. Harisingh Gour Vishwavidyalaya, (A Central University), 470003, Sagar, India
| | - Mini Bharati Ahirwar
- Department of Chemistry, Dr. Harisingh Gour Vishwavidyalaya, (A Central University), 470003, Sagar, India
| | - Milind M Deshmukh
- Department of Chemistry, Dr. Harisingh Gour Vishwavidyalaya, (A Central University), 470003, Sagar, India
| |
Collapse
|
7
|
Patkar D, Ahirwar MB, Deshmukh MM. Energetic Ordering of Hydrogen Bond Strengths in Methanol-Water Clusters: Insights via Molecular Tailoring Approach. Chemphyschem 2022; 23:e202200143. [PMID: 35302702 DOI: 10.1002/cphc.202200143] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/02/2022] [Revised: 03/17/2022] [Indexed: 12/31/2022]
Abstract
In this work, we examine the strength of various types of individual hydrogen bond (HB) in mixed methanol-water Mn Wm , (n+m=2 to 7) clusters, with an aim to understand the relative order of their strength, using our recently proposed molecular tailoring-based approach (MTA). Among all the types of HB, it is observed that the OM -H…OW HBs are the strongest (6.9 to 12.4 kcal mol-1 ). The next ones are OM -H…OM HBs (6.5 to 11.6 kcal mol-1 ). The OW -H…OW (0.2 to 10.9 kcal mol-1 ) and OW -H…OM HBs (0.3 to 10.3 kcal mol-1 ) are the weakest ones. This energetic ordering of HBs is seen to be different from the respective HB energies in the dimer i. e., OM -H…OM (5.0 to 6.0 kcal mol-1 )>OW -H…OM (1.5 to 6.0 kcal mol-1 )>OM -H…OW (3.8 to 5.6 kcal mol-1 )>OW -H…OW (1.2 to 5.0 kcal mol-1 ). The plausible reason for the difference in the HB energy ordering may be attributed to the increase or decrease in HB strengths due to the formation of cooperative or anti-cooperative HB networks. For instance, the cooperativity contribution towards the different types of HB follows: OM -H…OW (2.4 to 8.6 kcal mol-1 )>OM -H…OM (1.3 to 6.3 kcal mol-1 )>OW -H…OW (-1.0 to 6.5 kcal mol-1 )>OW -H…OM (-1.2 to 5.3 kcal mol-1 ). This ordering of cooperativity contribution is similar to the HB energy ordering obtained by the MTA-based method. It is emphasized here that, the interplay between the cooperative and anti-cooperative contributions are indispensable for the correct energetic ordering of these HBs.
Collapse
Affiliation(s)
- Deepak Patkar
- Department of Chemistry, Dr. Harisingh Gour Vishwavidyalaya, A Central University), Sagar, 470003, India
| | - Mini Bharati Ahirwar
- Department of Chemistry, Dr. Harisingh Gour Vishwavidyalaya, A Central University), Sagar, 470003, India
| | - Milind M Deshmukh
- Department of Chemistry, Dr. Harisingh Gour Vishwavidyalaya, A Central University), Sagar, 470003, India
| |
Collapse
|
8
|
Ahirwar MB, Gurav ND, Gadre SR, Deshmukh MM. Hydration Shell Model for Expeditious and Reliable Individual Hydrogen Bond Energies in Large Water Clusters. Phys Chem Chem Phys 2022; 24:15462-15473. [DOI: 10.1039/d2cp01663j] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Recently, we have developed and tested a method, based on the molecular tailoring approach (MTA-based) to directly estimate the individual hydrogen bond (HB) energies in molecular clusters. Application of this...
Collapse
|
9
|
Patkar D, Ahirwar MB, Shrivastava SP, Deshmukh MM. Assessment of hydrogen bond strengths and cooperativity in self- and cross-associating cyclic (HF)m(H2O)n (m + n = 2 to 8) clusters. NEW J CHEM 2022. [DOI: 10.1039/d1nj05431g] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023]
Abstract
In this work, we investigated the strengths of various self- and cross-associating hydrogen bonds (HBs) in mixed hydrogen fluoride–water cyclic (HF)m(H2O)n (m + n = 2 to 8) clusters, employing a molecular tailoring approach (MTA)-based method.
Collapse
Affiliation(s)
- Deepak Patkar
- Department of Chemistry, Dr Harisingh Gour Vishwavidyalaya (A Central University), Sagar, 470003, India
| | - Mini Bharati Ahirwar
- Department of Chemistry, Dr Harisingh Gour Vishwavidyalaya (A Central University), Sagar, 470003, India
| | - Satya Prakash Shrivastava
- Department of Chemistry, Dr Harisingh Gour Vishwavidyalaya (A Central University), Sagar, 470003, India
| | - Milind M. Deshmukh
- Department of Chemistry, Dr Harisingh Gour Vishwavidyalaya (A Central University), Sagar, 470003, India
| |
Collapse
|
10
|
Patkar D, Ahirwar MB, Gadre SR, Deshmukh MM. Unusually Large Hydrogen-Bond Cooperativity in Hydrogen Fluoride Clusters, (HF) n, n = 3 to 8, Revealed by the Molecular Tailoring Approach. J Phys Chem A 2021; 125:8836-8845. [PMID: 34612647 DOI: 10.1021/acs.jpca.1c06478] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
In this work, our recently proposed molecular tailoring approach (MTA)-based method is employed for the evaluation of individual hydrogen-bond (HB) energies in linear (L) and cyclic (C) hydrogen fluoride clusters, (HF)n (n = 3 to 8). The estimated individual HB energies calculated at the MP2(full)/aug-cc-pVTZ level for the L-(HF)n are between 6.2 to 9.5 kcal/mol and those in the C-(HF)n lie between 7.9 to 11.4 kcal/mol. The zero-point energy corrections and basis set superposition corrections are found to be very small (less than 0.6 and 1.2 kcal/mol, respectively). The cooperativity contribution toward individual HBs is seen to fall between 1.0 to 4.8 kcal/mol and 3.2 to 6.9 kcal/mol for linear and cyclic clusters, respectively. Interestingly, the HB energies in dimers, cleaved from these clusters, lie in a narrow range (4.4 to 5.2 kcal/mol) suggesting that the large HB strength in (HF)n clusters is mainly due to the large cooperativity contribution, especially for n ≥ 5 (50 to 62% of the HBs energy). Furthermore, the HB energies in these clusters show a good qualitative correlation with geometrical parameters (H···F distance and F-H···F angles), stretching frequencies of F-H bonds, and electron density values at the (3, -1) bond critical points.
Collapse
Affiliation(s)
- Deepak Patkar
- Department of Chemistry, Dr. Harisingh Gour Vishwavidyalaya, (A Central University), Sagar, 470003, India
| | - Mini Bharati Ahirwar
- Department of Chemistry, Dr. Harisingh Gour Vishwavidyalaya, (A Central University), Sagar, 470003, India
| | - Shridhar R Gadre
- Department of Scientific Computing, Modelling and Simulation, Savitribai Phule Pune University, Pune, 411 007, India
| | - Milind M Deshmukh
- Department of Chemistry, Dr. Harisingh Gour Vishwavidyalaya, (A Central University), Sagar, 470003, India
| |
Collapse
|
11
|
Estimation of resonance assisted hydrogen bond (RAHB) energies using properties of ring critical points in some dihydrogen-bonded complexes. J Mol Struct 2021. [DOI: 10.1016/j.molstruc.2021.130710] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
12
|
Ahirwar MB, Gurav ND, Gadre SR, Deshmukh MM. Molecular Tailoring Approach for Estimating Individual Intermolecular Interaction Energies in Benzene Clusters. J Phys Chem A 2021; 125:6131-6140. [PMID: 34251827 DOI: 10.1021/acs.jpca.1c03907] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
There is no general method available for the estimation of individual intermolecular interaction energies in weakly bound molecular clusters, and such studies are limited only to the dimer. Recently, we proposed a molecular tailoring approach-based method for the estimation of individual O-H···O hydrogen bond energies in water clusters. In the present work, we extend the applicability of this method for estimating the individual intermolecular interaction energies in benzene clusters, which are expected to be small. The basis set superposition error (BSSE)-corrected individual intermolecular interaction energies in linear (LN) benzene clusters, LN-(Bz)n n = 3-7, were calculated to be in the range from -1.75 to -2.33 kcal/mol with the cooperativity contribution falling between 0.05 and 0.20 kcal/mol, calculated at the MP2.5/aug-cc-pVDZ level of theory. In the case of non-linear (NLN) benzene clusters, NLN-(Bz)n n = 3-5, the BSSE-corrected individual intermolecular interaction energies exhibit a wider range from -1.16 to -2.55 kcal/mol with cooperativity contribution in the range from 0.02 to -0.61 kcal/mol. The accuracy of these estimated values was validated by adding the sum of interaction energies to the sum of monomer energies. These estimated molecular energies of clusters were compared with their actual calculated values. The small difference (<0.3 kcal/mol) in these two values suggests that our estimated individual intermolecular interaction energies in benzene clusters are quite reliable.
Collapse
Affiliation(s)
- Mini Bharati Ahirwar
- Department of Chemistry, Dr. Harisingh Gour Vishwavidyalaya (A Central University), Sagar 470003, India
| | - Nalini D Gurav
- Department of Scientific Computing, Modelling and Simulation, Savitribai Phule Pune University, Pune 411 007, India
| | - Shridhar R Gadre
- Department of Scientific Computing, Modelling and Simulation, Savitribai Phule Pune University, Pune 411 007, India
| | - Milind M Deshmukh
- Department of Chemistry, Dr. Harisingh Gour Vishwavidyalaya (A Central University), Sagar 470003, India
| |
Collapse
|
13
|
Molecular Tailoring Approach for the Estimation of Intramolecular Hydrogen Bond Energy. Molecules 2021; 26:molecules26102928. [PMID: 34069140 PMCID: PMC8155843 DOI: 10.3390/molecules26102928] [Citation(s) in RCA: 30] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/17/2021] [Revised: 05/11/2021] [Accepted: 05/12/2021] [Indexed: 11/29/2022] Open
Abstract
Hydrogen bonds (HBs) play a crucial role in many physicochemical and biological processes. Theoretical methods can reliably estimate the intermolecular HB energies. However, the methods for the quantification of intramolecular HB (IHB) energy available in the literature are mostly empirical or indirect and limited only to evaluating the energy of a single HB. During the past decade, the authors have developed a direct procedure for the IHB energy estimation based on the molecular tailoring approach (MTA), a fragmentation method. This MTA-based method can yield a reliable estimate of individual IHB energy in a system containing multiple H-bonds. After explaining and illustrating the methodology of MTA, we present its use for the IHB energy estimation in molecules and clusters. We also discuss the use of this method by other researchers as a standard, state-of-the-art method for estimating IHB energy as well as those of other noncovalent interactions.
Collapse
|
14
|
Yu D, Wu D, Liu JY, Li SY, Li Y. On single-electron magnesium bonding formation and the effect of methyl substitution. RSC Adv 2020; 10:34413-34420. [PMID: 35514394 PMCID: PMC9056782 DOI: 10.1039/d0ra06591a] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2020] [Accepted: 09/08/2020] [Indexed: 12/14/2022] Open
Abstract
The complexes formed between MgX2 (X = F, H) molecules and alkyl radicals Y [Y = CH3, CH2CH3, CH(CH3)2, and C(CH3)3] have been characterized by using quantum chemical methods. The binding distance in all cases is less than the sum of vdW radii of Mg and C, indicating the formation of a non-covalent interaction, namely single-electron magnesium bond. Energy decomposition analysis reveals that electrostatic and polarization contributions are the major components responsible for the stability of the studied complexes. According to interaction energy, atoms in molecules, and independent gradient model analyses, methyl substitution on electron donor Y imposes a positive effect on its complexation with MgX2. When compared with other nonbonded interactions, the single-electron magnesium bond is found to have strength comparable to those of the single-electron beryllium bond and π-magnesium bond. The complexes formed between MgX2 (X = F, H) molecules and alkyl radicals Y [Y = CH3, CH2CH3, CH(CH3)2, and C(CH3)3] have been characterized by using quantum chemical methods.![]()
Collapse
Affiliation(s)
- Dan Yu
- Institute of Theoretical Chemistry, Laboratory of Theoretical and Computational Chemistry, College of Chemistry, Jilin University Changchun 130023 P. R. China
| | - Di Wu
- Institute of Theoretical Chemistry, Laboratory of Theoretical and Computational Chemistry, College of Chemistry, Jilin University Changchun 130023 P. R. China
| | - Jing-Yao Liu
- Institute of Theoretical Chemistry, Laboratory of Theoretical and Computational Chemistry, College of Chemistry, Jilin University Changchun 130023 P. R. China
| | - Si-Yi Li
- Department of Transdisciplinary Science and Engineering, Tokyo Institute of Technology 2-12-1 Ookayama, Meguro Tokyo 152-8551 Japan
| | - Ying Li
- Institute of Theoretical Chemistry, Laboratory of Theoretical and Computational Chemistry, College of Chemistry, Jilin University Changchun 130023 P. R. China
| |
Collapse
|
15
|
Ahirwar MB, Gadre SR, Deshmukh MM. Direct and Reliable Method for Estimating the Hydrogen Bond Energies and Cooperativity in Water Clusters, W n, n = 3 to 8. J Phys Chem A 2020; 124:6699-6706. [PMID: 32786666 DOI: 10.1021/acs.jpca.0c05631] [Citation(s) in RCA: 24] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
No direct method for estimating the individual O-H···O hydrogen bond (H-bond) energies in water clusters (Wn) exists in the literature. In this work, we propose such a direct method based on the molecular tailoring approach, which also enables the estimation of the cooperativity contributions. The calculated H-bond energies at MP2(full)/aug-cc-pVTZ and CCSD(T)/aug-cc-pVDZ levels for Wn, n = 3 to 8, agree well with one another and fall between 0.3 and 11.6 kcal mol-1 with the cooperativity contributions in the range of -1.2 and 7.0 kcal mol-1. For gauging the accuracy of our H-bond energies for a cluster, the H-bond energy sum is added to the sum of monomer energies, and the results are compared with the respective total energy. These two values agree with each other to within 8.3 mH (∼5 kcal mol-1), testifying the accuracy of our estimated H-bond energies. Further, these H-bond strengths show a good correlation with the respective O-H stretching frequencies and the molecular electron density values at the (3, -1) O-H···O H-bond critical point.
Collapse
Affiliation(s)
- Mini Bharati Ahirwar
- Department of Chemistry, Dr. Harisingh Gour Vishwavidyalaya (A Central University), Sagar 470003, India
| | - Shridhar R Gadre
- Interdisciplinary School of Scientific Computing and Department of Chemistry, Savitribai Phule Pune University, Pune 411 007, India
| | - Milind M Deshmukh
- Department of Chemistry, Dr. Harisingh Gour Vishwavidyalaya (A Central University), Sagar 470003, India
| |
Collapse
|
16
|
Zins EL. Microhydration of a Carbonyl Group: How does the Molecular Electrostatic Potential (MESP) Impact the Formation of (H 2O) n:(R 2C═O)Complexes? J Phys Chem A 2020; 124:1720-1734. [PMID: 32049521 DOI: 10.1021/acs.jpca.9b09992] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The presence of a carbonyl group in a molecule usually leads to the identification of a π-hole on the molecular electrostatic potential (MESP) of the species. How does this electrophilic site influence the formation of microhydrated complexes? To address this point, a panel of R2CO solutes with various MESPs was selected, and we identified the structures and properties of several complexes containing one, two, three and six water molecules. The following solutes were considered in the present study: H2CO, F2CO, Cl2CO,(NC)2CO and H2C═CO. Geometry optimizations and frequency calculations were carried out at the LC-ωPBE/6-311++G(d,p) level, with the GD3BJ empirical correction for dispersion. For a number of n water molecules around the R2CO solute, the structure and the features of the most stable (H2O)n:(R2CO) complexes are highly dependent on the MESP of the isolated R2CO solute. The formation of pi-hole bondings appears to play a decisive role in the initiation of a three-dimensional organization of water molecules around the solute.
Collapse
Affiliation(s)
- Emilie-Laure Zins
- De la Molécule aux Nano-Objets: Réactivité, Interactions Spectroscopies, MONARIS, CNRS, Sorbonne Université, 75005, Paris, France
| |
Collapse
|
17
|
Singh V, Ibnusaud I, Gadre SR, Deshmukh MM. Fragmentation method reveals a wide spectrum of intramolecular hydrogen bond energies in antioxidant natural products. NEW J CHEM 2020. [DOI: 10.1039/d0nj00304b] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
Very strong and weak IHBs in curcumin.
Collapse
Affiliation(s)
- Vijay Singh
- Department of Chemistry
- Dr. Harisingh Gour Vishwavidyalaya (A Central University)
- Sagar
- India
| | - Ibrahim Ibnusaud
- Institute for Intensive Research in Basic Sciences
- Mahatma Gandhi University Campus
- P.O. Kottayam
- India
| | - Shridhar R. Gadre
- Interdisciplinary School of Scientific Computing and Department of Chemistry
- Savitribai Phule Pune University
- Pune 411 007
- India
| | - Milind M. Deshmukh
- Department of Chemistry
- Dr. Harisingh Gour Vishwavidyalaya (A Central University)
- Sagar
- India
| |
Collapse
|
18
|
Geiger DK, Geiger HC, Moore SM. Inter-molecular inter-actions in a phenol-substituted benzimidazole. Acta Crystallogr E Crystallogr Commun 2019; 75:272-276. [PMID: 30800466 PMCID: PMC6362671 DOI: 10.1107/s2056989019001270] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/22/2019] [Accepted: 01/23/2019] [Indexed: 11/10/2022]
Abstract
Hydrogen bonding plays an important role in the design of solid-state structures and gels with desirable properties. 1-(4-Hydroxybenzyl)-2-(4-hydroxyphenyl)-5,6-dimethyl-1H-benzimidazole was isolated as the acetone disolvate, C22H20N2O2·2C3H6O. O-H⋯N hydrogen bonding between benz-imidazole mol-ecules results in chains parallel to [010]. One of the acetone solvate mol-ecules participates in O-H⋯O hydrogen bonding with the benzimidazole derivative. C-H⋯π inter-actions are observed in the extended structure. Hirshfeld surface analysis was used to explore the inter-molecular inter-actions and density functional theory was used to estimate the strength of the hydrogen bonds.
Collapse
Affiliation(s)
- David K. Geiger
- Department of Chemistry, SUNY-College at Geneseo, Geneseo, NY 14454, USA
| | - H. Cristina Geiger
- Department of Chemistry, SUNY-College at Geneseo, Geneseo, NY 14454, USA
| | - Shawn M. Moore
- Department of Chemistry, SUNY-College at Geneseo, Geneseo, NY 14454, USA
| |
Collapse
|
19
|
Li C, Guo W, Zhou P, Tang Z. The effects of the heteroatom and position on excited-state intramolecular proton transfer of new hydroxyphenyl benzoxazole derivatives: a time-dependent density functional theory study. Org Chem Front 2019. [DOI: 10.1039/c9qo00295b] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
The effects of the heteroatom and position on excited-state intramolecular proton transfer (ESIPT) of 2-[4′-(N-4,6-dichloro-1,3,5-triazi-n-2-yl)2′hydroxyphenyl]benzoxazole (4THBO) have been investigated via time-dependent density functional theory studies.
Collapse
Affiliation(s)
- Changming Li
- School of Electrical Engineering
- University of South China
- Hengyang 421001
- China
- State Key Lab of Molecular Reaction Dynamics
| | - Wei Guo
- School of Electrical Engineering
- University of South China
- Hengyang 421001
- China
| | - Panwang Zhou
- Institute of Molecular Sciences and Engineering
- Shandong University
- Qingdao 266235
- P. R. China
| | - Zhe Tang
- Institute of Molecular Sciences and Engineering
- Shandong University
- Qingdao 266235
- P. R. China
| |
Collapse
|
20
|
Athar M, Das S, Jha PC, Jha AM. Conformational equilibrium study of calix[4]tetrolarenes using Density Functional Theory (DFT) and Molecular Dynamics simulations. Supramol Chem 2018. [DOI: 10.1080/10610278.2018.1517876] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
Affiliation(s)
- Mohd Athar
- School of Chemical Sciences, Central University of Gujarat, Gandhinagar, India
| | - Soubhik Das
- Centre for Applied Chemistry, Central University of Gujarat, Gandhinagar, India
| | - Prakash C. Jha
- Centre for Applied Chemistry, Central University of Gujarat, Gandhinagar, India
| | - Anand M. Jha
- Department of Chemistry, KSR College Sarairanjan, Samastipur, India
| |
Collapse
|
21
|
Hao J, Yang Y. Exploring the ESIPT dynamical processes of two novel chromophores: symmetrical structure CHC and asymmetric structure CHN. Org Chem Front 2018. [DOI: 10.1039/c7qo01127j] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The detailed ESIPT dynamical processes of CHC (symmetrical structure) and CHN (asymmetric structure) chromophores were revealed and compared using the TDDFT method at the B3LYP/6-31+G(d,p) level.
Collapse
Affiliation(s)
- Jiaojiao Hao
- State Key Laboratory of Molecular Reaction Dynamics
- Dalian Institute of Chemical Physics
- Chinese Academy of Sciences
- Dalian 116023
- China
| | - Yang Yang
- State Key Laboratory of Molecular Reaction Dynamics
- Dalian Institute of Chemical Physics
- Chinese Academy of Sciences
- Dalian 116023
- China
| |
Collapse
|
22
|
Thorave RG, Lande DN, Athare SV, Gejji SP, Gonnade RG, Malkhede DD. X-ray structure, spectral characteristics, thermal and redox behavior of quinoline encapsulated in sulfonatocalix[4]arene. J Mol Liq 2017. [DOI: 10.1016/j.molliq.2017.09.061] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
|
23
|
Zyubina TS, Shmygleva LV, Pisarev RV, Zyubin AS, Pisareva AV, Dobrovolsky YA, Volokhov VM. Crystalline hydrates of calix[4]arene-para-sulfonic acid with n (n = 6–16) water molecules: a structure modeling. Russ Chem Bull 2017. [DOI: 10.1007/s11172-017-1700-4] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
24
|
Investigation of structural and conformational equilibrium of Oxacalix[4]arene: A density functional theory approach. J Mol Liq 2017. [DOI: 10.1016/j.molliq.2017.04.065] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|
25
|
Lande DN, Gejji SP. Cooperative Hydrogen Bonding, Molecular Electrostatic Potentials, and Spectral Characteristics of Partial Thia-Substituted Calix[4]arene Macrocycles. J Phys Chem A 2016; 120:7385-97. [DOI: 10.1021/acs.jpca.6b07568] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Dipali N. Lande
- Department of Chemistry, Savitribai Phule Pune University, 411 007 Pune, India
| | - Shridhar P. Gejji
- Department of Chemistry, Savitribai Phule Pune University, 411 007 Pune, India
| |
Collapse
|
26
|
Siddiqui N, Singh V, Deshmukh MM, Gurunath R. Structures, stability and hydrogen bonding in inositol conformers. Phys Chem Chem Phys 2016; 17:18514-23. [PMID: 26108975 DOI: 10.1039/c5cp02690c] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
Various ab initio calculations using the density-functional (DFT), the second order Möller-Plesset perturbation (MP2) and self-consistent reaction field (SCRF) theories were performed on thirteen theoretically possible inositol stereoisomers. Gas phase calculations reveal that the myo- and neo-isomers of inositol (bearing one and two axial hydroxyl groups, respectively) are marginally more stable (by 0.5 kcal mol(-1)) than the all equatorially substituted scyllo-inositol. The calculations when done in different polar solvents show that the scyllo-inositol becomes the most stable inositol isomer, a fact attributed to weaker intramolecular hydrogen bonds. The individual hydrogen bond energy in all the isomers of inositol was also estimated using the molecular tailoring approach (MTA). The calculated hydrogen bond energies in these isomers are in excellent agreement with reported O-H···O hydrogen bond distances and ν(O-H) stretching frequencies. The estimated H-bond energy values suggest that the order of the intramolecular hydrogen bond strength follows: axial-axial > equatorial-axial > axial-equatorial > equatorial-equatorial hydrogen bonds. The intramolecular hydrogen bonds in the scyllo isomer are much weaker than those in other conformers, thus making this isomer more stable in polar solvents.
Collapse
Affiliation(s)
- Nazia Siddiqui
- Department of Chemistry, Indian Institute of Technology Kanpur, Kanpur, 208016 India.
| | | | | | | |
Collapse
|
27
|
Affiliation(s)
- A. Subha Mahadevi
- Centre for Molecular Modelling, CSIR-Indian Institute of Chemical Technology, Tarnaka, Hyderabad, India 500607
| | - G. Narahari Sastry
- Centre for Molecular Modelling, CSIR-Indian Institute of Chemical Technology, Tarnaka, Hyderabad, India 500607
| |
Collapse
|
28
|
Kumbhar S, Jana S, Anoop A, Waller MP. Cooperativity in bimetallic glutathione complexes. J Mol Graph Model 2015; 62:1-10. [DOI: 10.1016/j.jmgm.2015.05.003] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/18/2015] [Revised: 04/30/2015] [Accepted: 05/05/2015] [Indexed: 01/03/2023]
|
29
|
Zyubina TS, Shmygleva LV, Pisarev RV, Zyubin AS, Pisareva AV, Dobrovolsky YA, Volokhov VM. Calix[4]arene-p-sulfonic acid hydrates [CH2(OH)C6H2SO3H]4•nH2O (n = 6, 8, 9, 10, 12): a quantum chemical study. Russ Chem Bull 2015. [DOI: 10.1007/s11172-015-1201-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
30
|
Biswal HS, Bhattacharyya S, Bhattacherjee A, Wategaonkar S. Nature and strength of sulfur-centred hydrogen bonds: laser spectroscopic investigations in the gas phase and quantum-chemical calculations. INT REV PHYS CHEM 2015. [DOI: 10.1080/0144235x.2015.1022946] [Citation(s) in RCA: 49] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
31
|
Hydrogen Bonds Involving Sulfur: New Insights from ab Initio Calculations and Gas Phase Laser Spectroscopy. CHALLENGES AND ADVANCES IN COMPUTATIONAL CHEMISTRY AND PHYSICS 2015. [DOI: 10.1007/978-3-319-14163-3_2] [Citation(s) in RCA: 37] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
|
32
|
Deshmukh MM, Gadre SR, Cocinero EJ. Stability of conformationally locked free fructose: theoretical and computational insights. NEW J CHEM 2015. [DOI: 10.1039/c5nj02106e] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Total stabilization is governed by three predominant factors viz. the sum of energy of all H-bonds, ring strain and anomeric stabilization.
Collapse
Affiliation(s)
- Milind M. Deshmukh
- Department of Chemistry
- Dr. Harisingh Gour Central University
- Sagar, 470003
- India
| | - Shridhar R. Gadre
- Department of Chemistry
- Indian Institute of Technology Kanpur
- Kanpur 208 016
- India
| | - Emilio J. Cocinero
- Departamento de Química Física
- Facultad de Ciencia y Tecnología
- Universidad del País Vasco (UPV-EHU)
- 48080 Bilbao
- Spain
| |
Collapse
|
33
|
Liu X, Li Q, Cheng J, Li W. Influence of cooperativity on the frequency shift of the Ar–H stretch vibration in HArF complexes. Mol Phys 2012. [DOI: 10.1080/00268976.2012.728639] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
|
34
|
Li QZ, Sun L, Liu XF, Li WZ, Cheng JB, Zeng YL. Enhancement of Iodine-Hydride Interaction by Substitution and Cooperative Effects in NCX-NCI-HMY Complexes. Chemphyschem 2012; 13:3997-4002. [DOI: 10.1002/cphc.201200430] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/29/2012] [Revised: 08/13/2012] [Indexed: 11/05/2022]
|
35
|
Nagata T, Fedorov DG, Kitaura K. Analytic gradient for the embedding potential with approximations in the fragment molecular orbital method. Chem Phys Lett 2012. [DOI: 10.1016/j.cplett.2012.07.004] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
|