Search

Grid is a free and open-source Python library for constructing numerical grids to integrate, interpolate, and differentiate functions (e.g., molecular properties), with a strong emphasis on facilitating these operations in computational chemistry and conceptual density functional theory. Although designed, maintained, and released as a stand-alone Python library, Grid was originally developed for molecular integration, interpolation, and solving the Poisson equation in the HORTON and ChemTools packages. Grid is designed to be easy to use, extend, and maintain; this is why we use Python and adopt many principles of modern software development, including comprehensive documentation, extensive testing, continuous integration/delivery protocols, and package management. We leverage popular scientific packages, such as NumPy and SciPy, to ensure high efficiency and optimized performance in grid development. This article is the official release note of the Grid library showcasing its unique functionality and scope. [ABSTRACT FROM AUTHOR]

We argue that when one divides a molecular property into atom‐in‐a‐molecule contributions, one should perform the division based on the property density of the quantity being partitioned. This is opposition to the normal approach, where the electron density is given a privileged role in defining the properties of atoms‐in‐a‐molecule. Because partitioning each molecular property based on its own property density is inconvenient, we design a reference‐free approach that does not (directly) refer atomic property densities. Specifically, we propose a stockholder partitioning method based on relative influence of a molecule's atomic nuclei on the electrons at a given point in space. The resulting method does not depend on an “arbitrary” choice of reference atoms and it has some favorable properties, including the fact that all of the electron density at an atomic nucleus is assigned to that nucleus and the fact all the atoms in a molecule decay at a uniform asymptotic rate. Unfortunately, the resulting model is not easily applied to spatially degenerate ground states. Furthermore, the practical realizations of this strategy that we tried here gave disappointing numerical results. © 2017 Wiley Periodicals, Inc. [ABSTRACT FROM AUTHOR]

We establish the physical origins of chemical transferability from the perspective of the nearsightedness of electronic matter. To do this, we explicitly evaluate the response of electron density to a change in the system, at constant chemical potential, by computing the softness kernel, s(r, r'). The softness kernel is nearsighted, indicating that under constant-chemical-potential conditions like dilute solutions changing the composition of the molecule at r has only local effects and does not have any significant impact on the reactivity at positions r' far away from point r. This locality principle elucidates the transferability of functional groups in chemistry. [ABSTRACT FROM AUTHOR]

We show that the statistical divergence measures associated with non-extensive thermodynamic entropy functions-specifically the Tsallis, Réyni, Sharma-Mittal, supraextensive, and H-divergences-are associated with the Hirshfeld atoms-in-molecules partitioning. This extends the treatment of Nalewajski and Parr (J Phys Chem A 109:3957-3959, 2005), (for the extensive Shannon entropy) to non-extensive entropy measures. It also extends the work of Heidar-Zadeh and Ayers (J Chem Phys 142(4):044107, 2015), (for divergence measures that are local density functionals) to non-local functionals. These results dramatically extend the mathematical framework that one can use for similarity-based atoms-in-molecules partitioning. [ABSTRACT FROM AUTHOR]

We consider the problem of defining an appropriate local descriptor corresponding to an arbitrary global descriptor. Although it does not seem easy to rigorously and uniquely define local analogues of derived global descriptors (e.g., the electrophilicity) or the fundamental global descriptors associated with the canonical ensemble (e.g., the hardness), the local response of these global descriptors can be defined unambiguously. We look at the local response of the global electrophilicity and compare it to the conventional, ad hoc, definition of the local electrophilicity. The local response of global nucleofugality and electrofugality is also discussed. [ABSTRACT FROM AUTHOR]

We present explicit formulas for arbitrary-order derivatives of the energy, grand potential, electron density, and higher-order response functions with respect to the number of electrons, and the chemical potential for any smooth and differentiable model of the energy versus the number of electrons. The resulting expressions for global reactivity descriptors (hyperhardnesses and hypersoftnesses), local reactivity descriptors (hyperFukui functions and local hypersoftnesses), and nonlocal response functions are easy to evaluate computationally. Specifically, the explicit formulas for global/local/nonlocal hypersoftnesses of arbitrary order are derived using Bell polynomials. Explicit expressions for global and local hypersoftness indicators up to fifth order are presented. [ABSTRACT FROM AUTHOR]

We provide a general (yet elementary) proof of Pearson's hard/soft acid/base principle. Contrary to previous attempts to justify this rule, we go beyond the venerable parabolic model of Parr and Pearson. This new treatment encompasses more realistic E versus N dependencies, while emphasizing the importance of the correct normalization of these models. [ABSTRACT FROM AUTHOR]

We discuss an extension of the venerable Parr‐Pearson parabolic charge transfer model that explicitly includes 3rd order effects. We then use this new model to provide new arguments supporting the validity of Parr and Yang's "|Δμ| big is good" (DMB) principle. Additionally, we show how to re‐derive this principle from the Parr‐Pearson parabolic model, using simpler arguments than those presented in our previous proof of the DMB rule. [ABSTRACT FROM AUTHOR]

Using results from atomic spectroscopy, we show that there are two types of flat-planes conditions. The first type of flat-planes condition occurs when the energy as a function of the number of electrons of each spin, Nα and Nβ, has a derivative discontinuity on a line segment where the number of electrons, Nα + Nβ, is an integer. The second type of flat-planes condition occurs when the energy has a derivative discontinuity on a line segment where the spin polarization, Nα - Nβ, is an integer, but does not have a discontinuity associated with an integer number of electrons. Type 2 flat planes are rare-we observed just 15 type 2 flat-planes conditions out of the 4884 cases we tested-but their mere existence has implications for the design of exchange-correlation energy density functionals. To facilitate the development of functionals that have the correct behavior with respect to both fractional number of electrons and fractional spin polarization, we present a dataset for the chromium atom and its ions that can be used to test new functionals. [ABSTRACT FROM AUTHOR]