1
|
Sahu N, Khire SS, Gadre SR. Combining fragmentation method and high-performance computing: Geometry optimization and vibrational spectra of proteins. J Chem Phys 2023; 159:044309. [PMID: 37522406 DOI: 10.1063/5.0149572] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/07/2023] [Accepted: 07/12/2023] [Indexed: 08/01/2023] Open
Abstract
Exploring the structures and spectral features of proteins with advanced quantum chemical methods is an uphill task. In this work, a fragment-based molecular tailoring approach (MTA) is appraised for the CAM-B3LYP/aug-cc-pVDZ-level geometry optimization and vibrational infrared (IR) spectra calculation of ten real proteins containing up to 407 atoms and 6617 basis functions. The use of MTA and the inherently parallel nature of the fragment calculations enables a rapid and accurate calculation of the IR spectrum. The applicability of MTA to optimize the protein geometry and evaluate its IR spectrum employing a polarizable continuum model with water as a solvent is also showcased. The typical errors in the total energy and IR frequencies computed by MTA vis-à-vis their full calculation (FC) counterparts for the studied protein are 5-10 millihartrees and 5 cm-1, respectively. Moreover, due to the independent execution of the fragments, large-scale parallelization can also be achieved. With increasing size and level of theory, MTA shows an appreciable advantage in computer time as well as memory and disk space requirement over the corresponding FCs. The present study suggests that the geometry optimization and IR computations on the biomolecules containing ∼1000 atoms and/or ∼15 000 basis functions using MTA and HPC facility can be clearly envisioned in the near future.
Collapse
Affiliation(s)
- Nityananda Sahu
- Theoretische Chemie, Philipps-Universität Marburg, 35032 Marburg, Germany
| | - Subodh S Khire
- RIKEN Center for Computational Science, Kobe 650-0047, Japan
| | - Shridhar R Gadre
- Departments of Scientific Computing, Modelling & Simulation and Chemistry, Savitribai Phule Pune University, Pune 411007, India
| |
Collapse
|
2
|
Khire SS, Gattadahalli N, Gurav ND, Kumar A, Gadre SR. Constructing Potential Energy Surface with Correlated Theory for Dipeptides Using Molecular Tailoring Approach. Chemphyschem 2023; 24:e202200784. [PMID: 36735449 DOI: 10.1002/cphc.202200784] [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: 10/19/2022] [Revised: 01/31/2023] [Accepted: 02/01/2023] [Indexed: 02/04/2023]
Abstract
We demonstrate a cost-effective alternative employing the fragment-based molecular tailoring approach (MTA) for building the potential energy surface (PES) for two dipeptides viz. alanine-alanine and alanine-proline employing correlated theory, with augmented Dunning basis sets. About 1369 geometries are generated for each test dipeptide by systematically varying the dihedral angles Φ ${{\rm{\Phi }}}$ and Ψ ${{{\Psi }}}$ . These conformational geometries are partially optimized by relaxing all the other Z-matrix parameters, fixing the values of Φ ${{\rm{\Phi }}}$ and Ψ ${{{\Psi }}}$ . The MP2 level PES is constructed from the MTA-energies of chemically intact geometries using minimal hardware. The fidelity of MP2/aug-cc-pVDZ level PES is brought out by comparing it with its full calculation counterpart. Further, we bring out the power of the method by reporting the MTA-based CCSD/aug-cc-pVDZ level PES for these two dipeptides containing 498 and 562 basis functions respectively.
Collapse
Affiliation(s)
- Subodh S Khire
- RIKEN Center for Computational Science, Kobe, 650-0047, Japan.,Department of Scientific Computing Modelling and Simulation, Savitribai Phule Pune University, Pune, 411 007, India
| | - Nandini Gattadahalli
- Department of Scientific Computing Modelling and Simulation, Savitribai Phule Pune University, Pune, 411 007, India
| | - Nalini D Gurav
- Department of Scientific Computing Modelling and Simulation, Savitribai Phule Pune University, Pune, 411 007, India.,Organisch-Chemisches Institut and Center for Multiscale Theory and Computation (CMTC), Westfälische Wilhelms-Universität Münster, Corrensstrasse 36, 48149, Münster, Germany
| | - Anmol Kumar
- School of Pharmacy, University of Maryland, Baltimore, 20 Penn Street, HSFII, Baltimore, Maryland, 21201, U.S.A
| | - Shridhar R Gadre
- Department of Scientific Computing Modelling and Simulation, Savitribai Phule Pune University, Pune, 411 007, India
| |
Collapse
|
3
|
G Prasanna K, Sunil R, Gupta K, Lee SC. DJMol: An open-source modeling platform for computational chemistry and materials science with a Python interpreter. J Comput Chem 2021; 42:2116-2129. [PMID: 34406662 DOI: 10.1002/jcc.26740] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/13/2021] [Revised: 06/29/2021] [Accepted: 07/27/2021] [Indexed: 11/06/2022]
Abstract
We present a modular and extendable software suite, DJMol, for performing molecular simulations and it is demonstrated with DFTB+, Siesta, Atomic Simulation Environment, and OpenMD codes. It supports many of the standard features of an integrated development environment and consists of a structure builder and viewer, which could be connected with these electronic structure codes along with a set of data analyzers. This program comprises Java and Python modules and its libraries to carry out a different set of modeling tasks in materials science and chemistry. By adopting a Python interpreter into the software, a range of scriptable Python codes, such as Pymatgen can be incorporated into this programmable modeling platform. DJMol, through its common application programming interface (API), supports multiple modeling codes in the backend and several post-processing tools. It benefits an experienced user by increasing efficiency, while a nonexpert user by easy to use API.
Collapse
Affiliation(s)
- Krishnamohan G Prasanna
- Department of Science and Humanities, Mar Baselios College of Engineering and Technology, Trivandrum, Kerala, India
| | - Rahul Sunil
- Department of Computer Science, Mar Baselios College of Engineering and Technology, Trivandrum, Kerala, India
| | - Kapil Gupta
- Indo-Korea Science and Technology Center, Bengaluru, Karnataka, India
| | - Seung-Cheol Lee
- Indo-Korea Science and Technology Center, Bengaluru, Karnataka, India
| |
Collapse
|
4
|
Electrostatic Potential Topology for Probing Molecular Structure, Bonding and Reactivity. Molecules 2021; 26:molecules26113289. [PMID: 34072507 PMCID: PMC8198923 DOI: 10.3390/molecules26113289] [Citation(s) in RCA: 52] [Impact Index Per Article: 17.3] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/10/2021] [Revised: 05/16/2021] [Accepted: 05/25/2021] [Indexed: 11/18/2022] Open
Abstract
Following the pioneering investigations of Bader on the topology of molecular electron density, the topology analysis of its sister field viz. molecular electrostatic potential (MESP) was taken up by the authors’ groups. Through these studies, MESP topology emerged as a powerful tool for exploring molecular bonding and reactivity patterns. The MESP topology features are mapped in terms of its critical points (CPs), such as bond critical points (BCPs), while the minima identify electron-rich locations, such as lone pairs and π-bonds. The gradient paths of MESP vividly bring out the atoms-in-molecule picture of neutral molecules and anions. The MESP-based characterization of a molecule in terms of electron-rich and -deficient regions provides a robust prediction about its interaction with other molecules. This leads to a clear picture of molecular aggregation, hydrogen bonding, lone pair–π interactions, π-conjugation, aromaticity and reaction mechanisms. This review summarizes the contributions of the authors’ groups over the last three decades and those of the other active groups towards understanding chemical bonding, molecular recognition, and reactivity through topology analysis of MESP.
Collapse
|
5
|
Kulkarni AD. Molecular Hydration of Carbonic Acid: Ab Initio Quantum Chemical and Density Functional Theory Investigation. J Phys Chem A 2019; 123:5504-5516. [PMID: 31244117 DOI: 10.1021/acs.jpca.9b01122] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Molecular hydration of carbonic acid (H2CO3) is investigated in terms of bonding patterns in H2CO3···(H2O) n [ n = 1-4] hydrogen-bonded clusters within ab initio quantum chemical and density functional theory (DFT) frameworks. Successive addition of water molecules to H2CO3···H2O entails elongation of O-H (hydroxyl) bond as well as contraction of specific intermolecular hydrogen bonds signifying hydration of carbonic acid; these structural features get markedly enhanced under the continuum solvation framework. A comparison between the structurally similar clusters H2CO3···(H2O) n and HCOOH···(H2O) n [ n = 1-3] brings out the structural stability of the former. The present investigation in conjunction with the binding energy behavior of approaching water molecule(s) should serve as a precursor for pathways exploring aqueous dissociation of H2CO3 for larger clusters, as well as development of force-field potentials for acid dissociation process.
Collapse
Affiliation(s)
- Anant D Kulkarni
- Solid State and Structural Chemistry Unit , Indian Institute of Science , Bangalore 560012 , India
| |
Collapse
|
6
|
Harnessing desktop computers for ab initio calculation of vibrational IR/Raman spectra of large molecules. J CHEM SCI 2018. [DOI: 10.1007/s12039-018-1568-3] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|
7
|
Khire SS, Bartolotti LJ, Gadre SR. Harmonizing accuracy and efficiency: A pragmatic approach to fragmentation of large molecules. J Chem Phys 2018; 149:064112. [PMID: 30111143 DOI: 10.1063/1.5036595] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
Abstract
Fragmentation methods offer an attractive alternative for ab initio treatment of large molecules and molecular clusters. However, balancing the accuracy and efficiency of these methods is a tight-rope-act. With this in view, we present an algorithm for automatic molecular fragmentation within Molecular Tailoring Approach (MTA) achieving this delicate balance. The automated code is tested out on a variety of molecules and clusters at the Hartree-Fock (HF)- and Møller-Plesset second order perturbation theory as well as density functional theory employing augmented Dunning basis sets. The results show remarkable accuracy and efficiency vis-à-vis the respective full calculations. Thus the present work forms an important step toward the development of an MTA-based black box code for implementation of HF as well as correlated quantum chemical calculations on large molecular systems.
Collapse
Affiliation(s)
- Subodh S Khire
- Interdisciplinary School of Scientific Computing, Savitribai Phule Pune University, Pune 411007, India
| | - Libero J Bartolotti
- Department of Physical and Computational Chemistry, East Carolina University, Greenville, North Carolina 27858, USA
| | - Shridhar R Gadre
- Interdisciplinary School of Scientific Computing, Savitribai Phule Pune University, Pune 411007, India
| |
Collapse
|
8
|
Sahu N, Gadre SR. Vibrational infrared and Raman spectra of polypeptides: Fragments-in-fragments within molecular tailoring approach. J Chem Phys 2017; 144:114113. [PMID: 27004868 DOI: 10.1063/1.4943966] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022] Open
Abstract
The present work reports the calculation of vibrational infrared (IR) and Raman spectra of large molecular systems employing molecular tailoring approach (MTA). Further, it extends the grafting procedure for the accurate evaluation of IR and Raman spectra of large molecular systems, employing a new methodology termed as Fragments-in-Fragments (FIF), within MTA. Unlike the previous MTA-based studies, the accurate estimation of the requisite molecular properties is achieved without performing any full calculations (FC). The basic idea of the grafting procedure is implemented by invoking the nearly basis-set-independent nature of the MTA-based error vis-à-vis the respective FCs. FIF has been tested out for the estimation of the above molecular properties for three isomers, viz., β-strand, 310- and α-helix of acetyl(alanine)nNH2 (n = 10, 15) polypeptides, three conformers of doubly protonated gramicidin S decapeptide and trpzip2 protein (PDB id: 1LE1), respectively, employing BP86/TZVP, M06/6-311G**, and M05-2X/6-31G** levels of theory. For most of the cases, a maximum difference of 3 cm(-1) is achieved between the grafted-MTA frequencies and the corresponding FC values. Further, a comparison of the BP86/TZVP level IR and Raman spectra of α-helical (alanine)20 and its N-deuterated derivative shows an excellent agreement with the existing experimental spectra. In view of the requirement of only MTA-based calculations and the ability of FIF to work at any level of theory, the current methodology provides a cost-effective solution for obtaining accurate spectra of large molecular systems.
Collapse
Affiliation(s)
- Nityananda Sahu
- Department of Chemistry, Indian Institute of Technology Kanpur, Kanpur 208 016, India
| | - Shridhar R Gadre
- Department of Chemistry, Indian Institute of Technology Kanpur, Kanpur 208 016, India
| |
Collapse
|
9
|
Kumar A, Yeole SD, Gadre SR, López R, Rico JF, Ramírez G, Ema I, Zorrilla D. DAMQT 2.1.0: A new version of the DAMQT package enabled with the topographical analysis of electron density and electrostatic potential in molecules. J Comput Chem 2016; 36:2350-9. [PMID: 26505259 DOI: 10.1002/jcc.24212] [Citation(s) in RCA: 39] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/08/2015] [Revised: 09/04/2015] [Accepted: 09/09/2015] [Indexed: 01/22/2023]
Abstract
DAMQT-2.1.0 is a new version of DAMQT package which includes topographical analysis of molecular electron density (MED) and molecular electrostatic potential (MESP), such as mapping of critical points (CPs), creating molecular graphs, and atomic basins. Mapping of CPs is assisted with algorithmic determination of Euler characteristic in order to provide a necessary condition for locating all possible CPs. Apart from the mapping of CPs and determination of molecular graphs, the construction of MESP-based atomic basin is a new and exclusive feature introduced in DAMQT-2.1.0. The GUI in DAMQT provides a user-friendly interface to run the code and visualize the final outputs. MPI libraries have been implemented for all the tasks to develop the parallel version of the software. Almost linear scaling of computational time is achieved with the increasing number of processors while performing various aspects of topography. A brief discussion of molecular graph and atomic basin is provided in the current article highlighting their chemical importance. Appropriate example sets have been presented for demonstrating the functions and efficiency of the code.
Collapse
Affiliation(s)
- Anmol Kumar
- Department of Chemistry, Indian Institute of Technology Kanpur, Kanpur, 208016, Uttar Pradesh, India
| | - Sachin D Yeole
- Department of Chemistry, Indian Institute of Technology Kanpur, Kanpur, 208016, Uttar Pradesh, India
| | - Shridhar R Gadre
- Department of Chemistry, Indian Institute of Technology Kanpur, Kanpur, 208016, Uttar Pradesh, India
| | - Rafael López
- Departamento De Química Física Aplicada, Facultad De Ciencias, Universidad Autónoma De Madrid, Madrid, E-28049, Spain
| | - Jaime F Rico
- Departamento De Química Física Aplicada, Facultad De Ciencias, Universidad Autónoma De Madrid, Madrid, E-28049, Spain
| | - Guillermo Ramírez
- Departamento De Química Física Aplicada, Facultad De Ciencias, Universidad Autónoma De Madrid, Madrid, E-28049, Spain
| | - Ignacio Ema
- Departamento De Química Física Aplicada, Facultad De Ciencias, Universidad Autónoma De Madrid, Madrid, E-28049, Spain
| | - David Zorrilla
- Departamento De Química Física, Facultad De Ciencias, Universidad De Cádiz, Cádiz, E-11501, Spain
| |
Collapse
|
10
|
Sahu N, Gadre SR. Accurate vibrational spectra via molecular tailoring approach: a case study of water clusters at MP2 level. J Chem Phys 2016; 142:014107. [PMID: 25573553 DOI: 10.1063/1.4905004] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In spite of the recent advents in parallel algorithms and computer hardware, high-level calculation of vibrational spectra of large molecules is still an uphill task. To overcome this, significant effort has been devoted to the development of new algorithms based on fragmentation methods. The present work provides the details of an efficient and accurate procedure for computing the vibrational spectra of large clusters employing molecular tailoring approach (MTA). The errors in the Hessian matrix elements and dipole derivatives arising due to the approximation nature of MTA are reduced by grafting the corrections from a smaller basis set. The algorithm has been tested out for obtaining vibrational spectra of neutral and charged water clusters at Møller-Plesset second order level of theory, and benchmarking them against the respective full calculation (FC) and/or experimental results. For (H2O)16 clusters, the estimated vibrational frequencies are found to differ by a maximum of 2 cm(-1) with reference to the corresponding FC values. Unlike the FC, the MTA-based calculations including grafting procedure can be performed on a limited hardware, yet take a fraction of the FC time. The present methodology, thus, opens a possibility of the accurate estimation of the vibrational spectra of large molecular systems, which is otherwise impossible or formidable.
Collapse
Affiliation(s)
- Nityananda Sahu
- Department of Chemistry, Indian Institute of Technology Kanpur, Kanpur 208 016, India
| | - Shridhar R Gadre
- Department of Chemistry, Indian Institute of Technology Kanpur, Kanpur 208 016, India
| |
Collapse
|
11
|
Gadre SR, Yeole SD, Sahu N. Quantum chemical investigations on molecular clusters. Chem Rev 2014; 114:12132-73. [PMID: 25341561 DOI: 10.1021/cr4006632] [Citation(s) in RCA: 147] [Impact Index Per Article: 14.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Affiliation(s)
- Shridhar R Gadre
- Department of Chemistry, Indian Institute of Technology Kanpur , Kanpur 208 016, India
| | | | | |
Collapse
|
12
|
Sahu N, Singh G, Gadre SR. Exploring Structures and Energetics of Large OCS Clusters by Correlated Methods. J Phys Chem A 2013; 117:10964-72. [DOI: 10.1021/jp408311c] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Nityananda Sahu
- Department of Chemistry, Indian Institute of Technology Kanpur, Kanpur, 208016 India
| | - Gurmeet Singh
- Department of Chemistry, Indian Institute of Technology Kanpur, Kanpur, 208016 India
| | - Shridhar R. Gadre
- Department of Chemistry, Indian Institute of Technology Kanpur, Kanpur, 208016 India
| |
Collapse
|
13
|
Sahu N, Yeole SD, Gadre SR. Appraisal of molecular tailoring approach for large clusters. J Chem Phys 2013; 138:104101. [PMID: 23514459 DOI: 10.1063/1.4793706] [Citation(s) in RCA: 51] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023] Open
Abstract
High level ab initio investigations on molecular clusters are generally restricted to those of small size essentially due to the nonlinear scaling of corresponding computational cost. Molecular tailoring approach (MTA) is a fragmentation-based method, which offers an economical and efficient route for studying larger clusters. However, due to its approximate nature, the MTA-energies carry some errors vis-à-vis their full calculation counterparts. These errors in the MTA-energies are reduced by grafting the correction at a lower basis set (e.g., 6-31+G(d)) onto a higher basis set (e.g., aug-cc-pvdz or aug-cc-pvtz) calculation at MP2 level of theory. Further, better estimates of energies are obtained by making use of many-body interaction analysis. For this purpose, R-goodness (Rg) parameters for the three- and four-body interactions in a fragmentation scheme are proposed. The procedure employing grafting and many-body analysis has been tested out on molecular clusters of water, benzene, acetylene and carbon dioxide. It is found that for the fragmentation scheme having higher three- and four-body Rg-values, the errors in MTA-grafted energies are reduced typically to ~0.2 mH at MP2 level calculation. Coupled with the advantage in terms of computational resources and CPU time, the present method opens a possibility of accurate treatment of large molecular clusters.
Collapse
|
14
|
Yeole SD, Sahu N, Gadre SR. High-Level ab Initio Investigations on Structures and Energetics of N2O Clusters. J Phys Chem A 2013; 117:8591-8. [DOI: 10.1021/jp402649y] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Affiliation(s)
- Sachin D. Yeole
- Department of Chemistry, Indian Institute of Technology, Kanpur,
Kanpur 208016, India
- Department of Chemistry, University of Pune, Pune 411007, India
| | - Nityananda Sahu
- Department of Chemistry, Indian Institute of Technology, Kanpur,
Kanpur 208016, India
| | - Shridhar R. Gadre
- Department of Chemistry, Indian Institute of Technology, Kanpur,
Kanpur 208016, India
| |
Collapse
|
15
|
Yeole SD, López R, Gadre SR. Rapid topography mapping of scalar fields: large molecular clusters. J Chem Phys 2013; 137:074116. [PMID: 22920112 DOI: 10.1063/1.4746243] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
An efficient and rapid algorithm for topography mapping of scalar fields, molecular electron density (MED) and molecular electrostatic potential (MESP) is presented. The highlight of the work is the use of fast function evaluation by Deformed-atoms-in-molecules (DAM) method. The DAM method provides very rapid as well as sufficiently accurate function and gradient evaluation. For mapping the topography of large systems, the molecular tailoring approach (MTA) is invoked. This new code is tested out for mapping the MED and MESP critical points (CP's) of small systems. It is further applied to large molecular clusters viz. (H(2)O)(25), (C(6)H(6))(8) and also to a unit cell of valine crystal at MP2/6-31+G(d) level of theory. The completeness of the topography is checked by extensive search as well as applying the Poincaré-Hopf relation. The results obtained show that the DAM method in combination with MTA provides a rapid and efficient route for mapping the topography of large molecular systems.
Collapse
Affiliation(s)
- Sachin D Yeole
- Department of Chemistry, IIT Kanpur, Kanpur 208016, India
| | | | | |
Collapse
|
16
|
Furtado JP, Rahalkar AP, Shanker S, Bandyopadhyay P, Gadre SR. Facilitating Minima Search for Large Water Clusters at the MP2 Level via Molecular Tailoring. J Phys Chem Lett 2012; 3:2253-2258. [PMID: 26295779 DOI: 10.1021/jz300663u] [Citation(s) in RCA: 61] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
Water clusters (H2O)20 and (H2O)25 are explored at the Møller-Plesset second-order perturbation (MP2) level of theory. Geometry optimization is carried out on favorable structures, initially generated by the temperature basin paving (TBP) method, utilizing the fragment-based molecular tailoring approach (MTA). MTA-based stabilization energies at the complete basis set limit are accurately estimated by grafting the energy correction using a smaller basis set. For prototypical cases, the minima are established via MTA-based vibrational frequency calculations at the MP2/aug-cc-pVDZ level. The potential of MTA in tackling large clusters is further demonstrated by performing geometry optimization at MP2/aug-cc-pVDZ starting with the global minimum of (H2O)30 reported by Monte Carlo (MC) and molecular dynamics (MD) investigations. The present study brings out the efficacy of MTA in performing computationally expensive ab initio calculations with minimal off-the-shelf hardware without significant loss of accuracy.
Collapse
Affiliation(s)
- Jonathan P Furtado
- †Department of Chemistry, Indian Institute of Technology, Kanpur - 208016, India
| | - Anuja P Rahalkar
- †Department of Chemistry, Indian Institute of Technology, Kanpur - 208016, India
| | - Sudhanshu Shanker
- ‡School of Computational and Integrative Sciences, Jawaharlal Nehru University, New Delhi - 110067, India
| | - Pradipta Bandyopadhyay
- ‡School of Computational and Integrative Sciences, Jawaharlal Nehru University, New Delhi - 110067, India
| | - Shridhar R Gadre
- †Department of Chemistry, Indian Institute of Technology, Kanpur - 208016, India
| |
Collapse
|
17
|
|
18
|
|
19
|
Yeole SD, Sahu N, Gadre SR. Structures, energetics and vibrational spectra of CO2 clusters through molecular tailoring and cluster building algorithm. Phys Chem Chem Phys 2012; 14:7718-23. [DOI: 10.1039/c2cp23761j] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
20
|
Yeole SD, Gadre SR. Topography of scalar fields: molecular clusters and π-conjugated systems. J Phys Chem A 2011; 115:12769-79. [PMID: 21932861 DOI: 10.1021/jp2038976] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The pioneering works due to Bader and co-workers have generated widespread interest in the study of the topography of molecular scalar fields, the first step of which is the identification and characterization of the corresponding critical points (CPs). The topography of a molecular system becomes successively richer in going from the bare nuclear potential (BNP) to the molecular electrostatic potential (MESP) through the molecular electron density (MED). The present work clearly demonstrates, through the study of some π-conjugated test molecules as well as molecular clusters, that the CPs could be economically located by following this path within ab initio level theory. Further, the topography mapping of large molecules, especially at a higher level of theory, is known to be a demanding task. However, it is rendered possible by following the above sequential mapping assisted by a divide-and-conquer-type method termed as the molecular tailoring approach (MTA). This is demonstrated with the topography mapping of β-carotene and benzene nonamer at MP2 and a (H(2)O)(32) cluster at the HF level of theory, which are rather challenging problems with contemporary off-the-shelf computer hardware.
Collapse
Affiliation(s)
- Sachin D Yeole
- Department of Chemistry, Indian Institute of Technology Kanpur, Kanpur, India
| | | |
Collapse
|
21
|
“Gold standard” coupled-cluster study of acetylene pentamers and hexamers via molecular tailoring approach. Theor Chem Acc 2011. [DOI: 10.1007/s00214-011-1029-2] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/17/2022]
|
22
|
Yeole SD, Gadre SR. Molecular cluster building algorithm: electrostatic guidelines and molecular tailoring approach. J Chem Phys 2011; 134:084111. [PMID: 21361531 DOI: 10.1063/1.3556819] [Citation(s) in RCA: 43] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Nano-sized clusters of various materials are recent experimental targets, since they exhibit size-dependent physico-chemical properties. A vast amount of literature is available on the study of molecular clusters but general methods for systematic evolution of their growth are rather scarce. The present work reports a molecular cluster building algorithm based on the electrostatic guidelines, followed by ab initio investigations, enabled by the application of molecular tailoring approach. Applications of the algorithm for generating geometries and interaction energies of large molecular clusters of zinc sulfide, benzene, and water are presented.
Collapse
|
23
|
Rahalkar AP, Yeole SD, Ganesh V, Gadre SR. Molecular Tailoring: An Art of the Possible for Ab Initio Treatment of Large Molecules and Molecular Clusters. CHALLENGES AND ADVANCES IN COMPUTATIONAL CHEMISTRY AND PHYSICS 2011. [DOI: 10.1007/978-90-481-2853-2_10] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
|
24
|
Mahadevi AS, Rahalkar AP, Gadre SR, Sastry GN. Ab initio investigation of benzene clusters: Molecular tailoring approach. J Chem Phys 2010; 133:164308. [DOI: 10.1063/1.3494536] [Citation(s) in RCA: 68] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023] Open
|
25
|
Rahalkar AP, Katouda M, Gadre SR, Nagase S. Molecular tailoring approach in conjunction with MP2 and Ri-MP2 codes: A comparison with fragment molecular orbital method. J Comput Chem 2010; 31:2405-18. [DOI: 10.1002/jcc.21533] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
|
26
|
Deshmukh MM, Gadre SR. Estimation of N-H...O=C intramolecular hydrogen bond energy in polypeptides. J Phys Chem A 2009; 113:7927-32. [PMID: 19496581 DOI: 10.1021/jp9031207] [Citation(s) in RCA: 72] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Abstract
The previously proposed molecular tailoring approached (MTA) [Deshmukh, M. M.; Gadre, S. R.; Bartolotti, L. J. J. Phys. Chem. A 2006, 110, 12519] for the estimation of intramolecular O-H...O hydrogen bond energy is extended to that for the N-H...O=C bond within polypeptides. The methodology is initially tested on a tetrapeptide containing two types of N-H...O=C hydrogen bonds and is found to distinguish between them. The estimated values are in good agreement with the trends predicted by the geometrical parameters. Furthermore, this methodology is applied to partially as well as fully substituted, capped polyglycines that contain five glycine residues (acetyl-(gly)(5)-NH(2)) to check the effect of substituents on the energetics of hydrogen bonds. The estimated N-H...O=C bond energy values lie in the range of 4-6 kcal/mol. These estimated values are not only in concurrence with the geometric parameters but also able to reflect the subtle effects of substituents for the substituted polypeptides studied in the present work.
Collapse
Affiliation(s)
- Milind M Deshmukh
- Homi Bhabha Centre for Science Education, Tata Institute of Fundamental Research, V. N. Purav Marg, Mankhurd, Mumbai-400 088, India.
| | | |
Collapse
|