Evolution of the atomic valence observed by the reaction fragility spectra on the reaction path

Shannon entropy variation as a global indicator of electron density contraction at interatomic regions in chemical reactions

CONTEXT

The negative of the Shannon entropy derivative is proposed to account for electron density contraction as the chemical bonds are breaking and forming during a chemical reaction. We called this property the electron density contraction index, EDC, which allows identifying stages in a reaction that are dominated by electron contraction or expansion. Four different reactions were analyzed to show how the EDC index changes along the reaction coordinate. The results indicate that the rate of change of Shannon entropy is directly related to the rate of change of the electron density at the bond critical points between all the atomic pairs in the molecular systems. It is expected that EDC will complement the detailed analysis of reaction mechanisms that can be performed with the theoretical tools available to date.

METHODS

Density functional theory calculations at the B3LYP/6-31G(d,p) level of theory were carried out using Gaussian 16 to analyze the reaction mechanisms of the four reactions studied. The reaction paths were obtained via the intrinsic reaction coordinate method, which served as the reaction coordinate to obtain the reaction force and the EDC profiles in each case. Shannon entropy and electron density at the bond critical points were calculated using the Multiwfn 3.7 package.

Hyperhardness and hypersoftness of atoms and their ions

CONTEXT

The theory of reactivity based on cDFT has been supplemented with the new method of calculating the atomic and local indices. With the use of previously derived relationship of the electron density gradient to the softness kernel and to the linear response function, we deliver theoretical analysis to obtain significant reactivity indices-the electron density derivatives: local softness and local hypersoftness together with the global hyperhardness index and the derivative of the global softness with respect to the number of electrons. The local derivatives have been applied in the calculation of responses of atoms to perturbation by an external potential by the alchemical approach. The vital role of the local softness has been confirmed; the potential role of the hypersoftness has been indicated.

METHOD

Our original theoretical scheme has been numerically illustrated with the results obtained with electron density calculations with B3LYP method implemented in Gaussian 16 package. The aug-cc-pvqz basis set has been routinely applied, except for the Ca atom (cc-pvqz). Using the pVTZ basis set recommended by Sadlej was necessary for the potassium atom.

The Connectivity Matrix: A Toolbox for Monitoring Bonded Atoms and Bonds

The concept of a connectivity matrix, essential for the reaction fragility (RF) spectra technique for monitoring electron density evolution in a chemical reaction, has been supported with a novel formulation for the diagonal matrix elements; their direct link to the electron density function ρ(r) has been demonstrated. By combining the concept with the atomization energy of a system, the separation of the potential energy into atomic and/or bond contributions has been achieved. The energy derivative diagrams for atoms and bonds that are variable along a reaction path provide new insight into the reaction mechanism. Diagonalization of the connectivity matrix resulted in the eigenvectors that provide information on a role of individual atoms in the development of structural changes along a reaction path.

Bond Softening Indices Studied by the Fragility Spectra for Proton Migration in Formamide and Related Structures

Computational scheme to obtain bond softening index λ, defined within the conceptual DFT, has been obtained with the use of the reaction fragility (RF) concept. Numerical results were obtained with the RF spectra for the proton transfer reaction in formamide molecule (H₂NCHO) and the water assisted proton migration in H₂NCHO·H₂O complex. Double proton transfer reaction in the formamide dimer, (H₂NCHO)₂, and its analogues, (H₂NCHS)₂ and (H₂NCHO)·(H₂NCHS), have also been studied. The atomic and bond RF spectra clearly describe the density reorganization in the backbone of each molecule, resulting from proton displacement in the systems. The obtained softening indices have been calculated for hydrogen atoms in the reactant state (RS) and product state (PS) configuration. These indices provide fine characteristics for the local sensitivity of the reacting system to a disturbance of the position of a chosen atom.

Bond Fragility Spectra for the Double Proton-Transfer Reaction in the Formic Acid-Type Dimers

The newly developed method of fragility spectra for observation of bond breaking and formation upon a reaction has been applied to the canonical reaction series of the double proton transfer (DPT). Formic acid and its thio-analogues HCXYH (X, Y = O, S) have been chosen for the analysis. Very accurate linear correlations have been determined between the nondiagonal elements of the connectivity matrix, essential for the method, and the Wiberg bond orders for the corresponding bonds. Relation of the slope of this correlation to the global softness and to the atomic numbers of the bonded atoms has been proved, thus corroborating the c-DFT formula describing the fragility spectra. The electron density changes in bonds, as observed by the fragility spectra, are in harmony with the curvature diagrams reported by other authors.