Electrostatics in chemistry: 4. Electrostatic models for weak molecular complexation

This article does not have an abstract

Electrostatics in chemistry. 3. Molecular electrostatic potential: visualization and topography

Topographical features of the electrostatic potential of molecules (MESP) are discussed. Three-dimensional visualization of the MESP of some small molecules is presented. The details are correlated with the reactivity of molecules

Many-body interaction analysis: Algorithm development and application to large molecular clusters

A completely automated algorithm for performing many-body interaction energy analysis of clusters (MBAC) [M. J. Elrodt and R. J. Saykally, Chem. Rev. 94, 1975 (1994); S. S. Xantheas, J. Chem. Phys. 104, 8821 (1996)] at restricted Hartree-Fock (RHF)/MA Plesset 2nd order perturbation theory (MP2)/density functional theory (DFT) level of theory is reported. Use of superior guess density matrices (DM's) for smaller fragments generated from DM of the parent system and elimination of energetically insignificant higher-body combinations, leads to a more efficient performance (speed-up up to 2) compared to the conventional procedure. MBAC approach has been tested out on several large-sized weakly bound molecular clusters such as (H(2)O)(n), n=8, 12, 16, 20 and hydrated clusters of amides and aldehydes. The MBAC results indicate that the amides interact more strongly with water than aldehydes in these clusters. It also reconfirms minimization of the basis set superposition error for large cluster on using superior quality basis set. In case of larger weakly bound clusters, the contributions higher than four body are found to be repulsive in nature and smaller in magnitude. The reason for this may be attributed to the increased random orientations of the interacting molecules separated from each other by large distances.Financial support from the Council of Scientific and Industrial Research (CSIR), New Delhi, India and the Center for Development of Advanced Computing (C-DAC), Pune, India, is gratefully acknowledged

Enabling ab initio Hessian and frequency calculations of large molecules

A linear scaling method, termed as cardinality guided molecular tailoring approach, is applied for the estimation of the Hessian matrix and frequency calculations of spatially extended molecules. The method is put to test on a number of molecular systems largely employing the Hartree-Fock and density functional theory for a variety of basis sets. To demonstrate its ability for correlated methods, we have also performed a few test calculations at the Moller-Plesset second order perturbation theory. A comparison of central processing unit and memory requirements for medium-sized systems with those for the corresponding full ab initio computation reveals substantial gains with negligible loss of accuracy. The technique is further employed for a set of larger molecules, Hessian and frequency calculations of which are not possible on commonly available personal-computer-type hardware.Financial support from the Center for Development of Advanced Computing C-DAC, Pune, Naval Research Board NRB, New Delhi, and Council of Scientific and Industrial Research CSIR, New Delhi is gratefully acknowledged

Signatures of molecular recognition from the topography of electrostatic potential

The recognition of interaction between two molecules is analysed via the topography of their molecular electrostatic potentials (MESP). The point of recognition between two species is proposed to be the geometry at which there is a change in the nature of the set of MESP critical points of one of the molecules vis-a-vis with its MESP topography at infinite separation. These results are presented for certain model systems such as pyridine and benzene dimers, cytosine-guanine and adenine-thymine base pairs in various orientations of approach of the two species

Nonexistence of local maxima in molecular electrostatic potential maps

It has been rigorously established by means of classical electrostatic arguments, that molecular electrostatic potential maps are devoid of local maxima. This forms a generalization of the earlier works of Politzer and co-workers which were restricted to the case of atoms

Molecular tailoring approach for geometry optimization of large molecules: Energy evaluation and parallelization strategies

A linear-scaling scheme for estimating the electronic energy, gradients, and Hessian of a large molecule at ab initio level of theory based on fragment set cardinality is presented. With this proposition, a general, cardinality-guided molecular tailoring approach (CG-MTA) for ab initio geometry optimization of large molecules is implemented. The method employs energy gradients extracted from fragment wave functions, enabling computations otherwise impractical on PC hardware. Further, the method is readily amenable to large scale coarse-grain parallelization with minimal communication among nodes, resulting in a near-linear speedup. CG-MTA is applied for density-functional-theory-based geometry optimization of a variety of molecules including alpha-tocopherol, taxol, gamma-cyclodextrin, and two conformations of polyglycine. In the tests performed, energy and gradient estimates obtained from CG-MTA during optimization runs show an excellent agreement with those obtained from actual computation. Accuracy of the Hessian obtained employing CG-MTA provides good hope for the application of Hessian-based geometry optimization to large molecules.The authors thank C-DAC, UGC (CAS program to the University of Pune) and the Council of Scientific and Industrial Research (CSIR), New Delhi for financial assistance

Molecular tailoring approach for exploring structures, energetics and properties of clusters

Molecular Tailoring Approach (MTA) is a method developed for enabling ab initio calculations on prohibitively large molecules or atomic/molecular clusters. A brief review of MTA, a linear scaling technique based on set inclusion and exclusion principle, is provided. The Molecular Electrostatic Potential (MESP) of smaller clusters is exploited for building initial geometries for the larger ones, followed by MTA geometry optimization. The applications of MTA are illustrated with a few test cases such as (CO2)n and Lin clusters employing Density Functional theory (DFT) and a nanocluster of orthoboric acid at the Hartree-Fock (HF) level. Further, a discussion on the geometries and energetics of benzene tetramers and pentamers, treated at the Moller-Plesset second order (MP2) perturbation theory, is given. MTA model is employed for evaluating some cluster properties viz. adiabatic ionization potential, MESP, polarizability, Hessian matrix and infrared frequencies. These property evaluations are carried out on a series of test cases and are seen to offer quite good agreement with those computed by an actual calculation. These case studies highlight the advantages of MTA model calculations vis-a-vis the actual ones with reference to the CPU-time, memory requirements and accuracy

Rapid topography mapping of scalar fields: Large molecular clusters

The following article appeared in Journal of Chemical Physics 137.7 (2012): 074116 and may be found at http://scitation.aip.org/content/aip/journal/jcp/137/7/10.1063/1.4746243An 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 2O) 25, (C 6H 6) 8 and also to a unit cell of valine crystal at MP26-31G(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 systemsAuthors thank the Center for Development of Advanced Computing (C-DAC), Pune for financial and computational support. S.R.G. is grateful to the Department of Science and Technology (DST) for the award of J. C. Bose National Fellowship. R. López acknowledges partial funding from the CAM (S2009_PPQ-1545, LIQUORGAS) and MICINN (CTQ2010-19232). Authors are also thankful to Dr. Graeme M. Day, University of Cambridge, for providing the coordinates of unit cell of valine crystal and to Dr. V. Subramanian, CLRI, Chennai for providing some test run

A novel approach for the study of intermolecular interactions: molecular deformation densities

A new method to study intermolecular interactions via direct electron density measurements has been proposed. This approach involves the sum of constituent molecular electron densities as a reference rather than the usual sum of spherical atom densities for the evaluation of deformation densities. One may thus be able to focus specifically on the regions involving dominant intermolecular interactions such as hydrogen bonding, charge transfer, etc