Information theoretic properties from the quantum theory of atoms in molecules

Orbital-Free Quantum Crystallographic View on Noncovalent Bonding: Insights into Hydrogen Bonds, π⋅⋅⋅π and Reverse Electron Lone Pair⋅⋅⋅π Interactions

A detailed analysis of a complete set of the local potentials that appear in the Euler equation for electron density is carried out for noncovalent interactions in the crystal of a uracil derivative using experimental X-ray charge density. The interplay between the quantum theory of atoms in molecules and crystals and the local potentials and corresponding inner-crystal electronic forces of electrostatic and kinetic origin is explored. Partitioning of crystal space into atomic basins and atomic-like potential basins led us to the definite description of interatomic interaction and charge transfer. Novel physically grounded bonding descriptors derived within the orbital-free quantum crystallography provided the detailed examination of π-stacking and intricate C=O⋅⋅⋅π interactions and nonclassical hydrogen bonds present in the crystal. The donor-acceptor character of these interactions is revealed by analysis of Pauli and von Weizsäcker potentials together with well-known functions, e. g., deformation electron density and electron localization function. In this way, our analysis throws light on aspects of these closed-shell interactions hitherto hidden from the description.

Density Functional Theory and Information-Theoretic Approach Study on the Origin of Homochirality in Helical Structures

Homochirality of macromolecules such as proteins and DNA is one of the most striking features in nature; yet, there is still no convincing theory to explain its origin. In a recent work by one of the present authors (J. Phys. Chem. Lett. 2020, 11, 8690-8696), a general proposal from the viewpoint of thermodynamics has been put forward. It proposes that it is the handedness of helices ubiquitous in biological macromolecules that plays the decisive role. It also unveiled that there exist strong cooperativity effects dominated by favorable electrostatic interactions in the homochiral conformer. In this work, making use of analytical tools, we recently developed a density functional theory and an information-theoretic approach and through four sets of helical structures we designed for the present study, we examine these systems to provide new insights about these properties. We found that the 3₁₀-helix and the α-helix are markedly different in cooperativity from the viewpoint of both the total energy and its three components. The electrostatic dominance of homochiral species is manifested by both the electron charge distribution and information gain. At the atomic level, different elements behave significantly differently because they play different roles in the systems. Our results from this work validate that these analytical tools can be applied to homochiral systems, which can be further extended to others with potential interest in asymmetric synthesis and macromolecular assembly where the Principle of Homochirality Hierarchy comes into play.

N-Derivatives of Shannon entropy density as response functions

The exact first and second order partial derivatives of Shannon entropy density with respect to the number of electrons at constant external potential are introduced as new descriptors for prediction of the active sites of a molecule. The derivatives, which are a measure of the inhomogeneity of electron density, are calculated both exactly (from analytical forms) and approximately (using the finite difference method) for some molecular systems. According to the maximum entropy principle, the extreme value of the first order derivative on the surface of a given molecule should determine the active sites of the molecule in electrophilic and nucleophilic attack. The second order derivative indicates where the Shannon entropy is more concentrated or depleted during the electron exchange. Although these derivatives on the surfaces of helium and neon atoms are uniform, the corresponding values for argon, krypton and xenon atoms are not. This could explain the greater tendency of heavy noble gas atoms to form stable compounds. A dual descriptor is also defined as the difference between the left and right first order derivatives of Shannon entropy density, which allows one to simultaneously predict the preferable sites for electrophilic and nucleophilic attack over the system at point r. Therefore, the reactivity of an atom in a molecule requires the non-uniformity of the first and second order derivatives of Shannon entropy density on the surface of that atom.

Is It Possible To Determine Oxidation States for Atoms in Molecules Using Density-Based Quantities? An Information-Theoretic Approach and Conceptual Density Functional Theory Study

The oxidation state, also called oxidation number, of atoms in molecules is a fundamental chemical concept. It is defined as the charge of an atom in a molecule after the ionic approximation of its heteronuclear bonds is applied. Even though for simple molecules the assignment of oxidation states is straightforward, redundancy and ambiguity do exist for others. In this work, we present a density-based framework to determine the oxidation state using the quantities from the information-theoretic approach. As a proof of concept, we examined six elements for a total of 49 molecules. Strong linear correlations were obtained with Shannon entropy, Ghosh-Berkowitz-Parr entropy, information gain, relative Rényi entropy of orders 2 and 3, and Hirshfeld charge. We also discovered that the crystal radius of elements plays the key role in rationalizing the system dependent nature of these strong linear correlations. The validity and effectiveness of our results were demonstrated by the examples of molecules containing elements with two or more oxidation states. Our results should be applicable to more complicated systems in assigning different oxidation states.

Scaling properties of information-theoretic quantities in density functional reactivity theory

A number of strong linear correlations between information-theoretic quantities and electron populations for atoms, molecules, and atoms-in-molecules have been disclosed.

DevelopingWinXPRO: a software for determination of the multipole-model-based properties of crystals

The new release of the computer program packageWinXPRO v.3xfor determination of the crystal properties from parameters of the multipole-modeled experimental electron density and anharmonic atomic displacement coefficients is described. The set of properties is significantly extended by using the density functional and information theories. In addition, a built-in multi-functional viewer and programs to display the output data, including the mapping of the chosen functional bonding descriptors onto surfaces of the other properties, are included.

The structures of heterocyclic complexes ruled by hydrogen bonds and halogen interactions: interaction strength and IR modes

In this work, the existence of multiple interactions in heterocyclic complexes of C2H4O⋯nHCCl3 and C2H4S⋯nHCCl3 with n=2 and 3 was unveiled at the B3LYP/6-311++G(d,p) level of theory. The forward analyses of the vibrational spectra revealed the appearing of red-shifts in the H-C bond. In agreement with this and through the optimized geometries of these systems, an increase in the H-C bond length was also observed. Besides O⋯H and S⋯H, other hydrogen bonds formed between chlorine⋯hydrogen and mainly the halogen interactions formed by chlorine⋯chlorine were identified. Thereby, the vibration spectra of the heterocyclic complexes were reanalyzed with the purpose to locate new red-shifts, although only those characterized in H-C have been detected up to then. In addition to the correlation between the frequencies shifted to downward values followed by increases in the bond lengths, the interpretation of the red-shifts was conducted by means of the Bent rule of the hybridization theory. The interaction strength was examined in several viewpoints, and one of them was the relationship between the H-bond energies and the intermolecular electronic density computed by means of the Quantum Theory of Atoms in Molecules (QTAIM). Moreover, the prediction of the interaction strength was also made through the combination between vibration modes (red-shifts) and variation of topological parameters, such as the electronic density and Laplacian of the proton donor bond (C-H).