Nuclear reactivity indices in the context of spin polarized density functional theory

On the link between the reaction force constant and conceptual DFT

CONTEXT

The reaction force constant ( κ ), introduced by Professor Alejandro Toro-Labbé, plays a pivotal role in characterizing the reaction pathway by assessing the curvature of the potential energy profile along the intrinsic reaction coordinate. This study establishes a novel link between κ and the reactivity descriptors of conceptual density functional theory (c-DFT). Specifically, we derive expressions that relate the reaction force constant to nuclear softness and variations in chemical potential. Our findings indicate that regions of the reaction pathway where κ is negative match with significant electronic structure rearrangements, while positive κ regions match mostly with geometric rearrangements. This correlation between κ and c-DFT reactivity descriptors enhances our understanding of the underlying forces driving chemical reactions and offers new perspectives for analyzing reaction mechanisms.

METHODS

The internal reaction path for the proton transfer in SNOH, chemical potential, and nuclear softness were computed using DFT with B3LYP exchange-correlation functional and 6-311++G(d,2p) basis set.

Molecular Interactions From the Density Functional Theory for Chemical Reactivity: The Interaction Energy Between Two-Reagents

Reactivity descriptors indicate where a reagent is most reactive and how it is most likely to react. However, a reaction will only occur when the reagent encounters a suitable reaction partner. Determining whether a pair of reagents is well-matched requires developing reactivity rules that depend on both reagents. This can be achieved using the expression for the minimum-interaction-energy obtained from the density functional reactivity theory. Different terms in this expression will be dominant in different circumstances; depending on which terms control the reactivity, different reactivity indicators will be preferred.

How predictive could alchemical derivatives be?

The chemical space contains all possible compounds that can be imagined. Its size easily equals the number of fundamental particles in the observable universe. Rational design of compounds aims to find those sectors of the chemical space where compounds optimize a set of desired properties. Then, rational design demands tools to efficiently navigate the chemical space. Ab initio alchemical derivatives offer the possibility to navigate, without empiricism, the energy landscape through alchemical transformations. An alchemical transformation is any process, physical or fictitious, that connects to points in the chemical space. In this work, those transformations are constructed as a perturbative expansion of the energy with respect to perturbations in the stoichiometry. The response functions of that expansion are what is called alchemical derivatives. In this work we assess how effective alchemical derivatives are in predicting energy changes associated to changes in the composition. We do this by including in the expansion, for the first time, electrostatic, polarization and electron-transfer effects. The system we chose is one that challenges alchemical derivatives because none of these effects dominates its behavior. The transmutations studied here correspond to substitutional doping of Al₁₃ with up to four atoms of Si, Al_13-nSi_n. Two types of transformations are considered, those in which the number of electrons remains constant and those in which the number of electrons also changes. It is found that contrary to what has been reported before, polarization cannot be neglected. If polarization is not included, alchemical derivatives fail to predict the change of energy and the relative energy between isomers. For isoelectronic substitution of four or more atoms, the perturbative approach collapses because the strength of the perturbation becomes too strong to guarantee convergence of the series. It is shown, however, that if only one atom is mutated at a time, alchemical derivatives rank pretty well the isomers of Al_13-nSi_n according to their energy. In the case of non-isoelectronic transformations, it is observed that the series rapidly diverges with increasing number of electrons. In this situation, it becomes more important to keep the degree of transmutation of the parent system small.

Conceptual DFT: the chemical relevance of higher response functions

In recent years conceptual density functional theory offered a perspective for the interpretation/prediction of experimental/theoretical reactivity data on the basis of a series of response functions to perturbations in the number of electrons and/or external potential. This approach has enabled the sharp definition and computation, from first principles, of a series of well-known but sometimes vaguely defined chemical concepts such as electronegativity and hardness. In this contribution, a short overview of the shortcomings of the simplest, first order response functions is illustrated leading to a description of chemical bonding in a covalent interaction in terms of interacting atoms or groups, governed by electrostatics with the tendency to polarize bonds on the basis of electronegativity differences. The second order approach, well known until now, introduces the hardness/softness and Fukui function concepts related to polarizability and frontier MO theory, respectively. The introduction of polarizability/softness is also considered in a historical perspective in which polarizability was, with some exceptions, mainly put forward in non covalent interactions. A particular series of response functions, arising when the changes in the external potential are solely provoked by changes in nuclear configurations (the "R-analogues") are also systematically considered. The main part of the contribution is devoted to third order response functions which, at first sight, may be expected not to yield chemically significant information, as turns out to be for the hyperhardness. A counterexample is the dual descriptor and its R analogue, the initial hardness response, which turns out to yield a firm basis to regain the Woodward-Hoffmann rules for pericyclic reactions based on a density-only basis, i.e. without involving the phase, sign, symmetry of the wavefunction. Even the second order nonlinear response functions are shown possibly to bear interesting information, e.g. on the local and global polarizability. Its derivatives may govern the influence of charge on the polarizability, the R-analogues being the nuclear Fukui function and the quadratic and cubic force constants. Although some of the higher order derivatives may be difficult to evaluate a comparison with the energy expansion used in spectroscopy in terms of nuclear displacements, nuclear magnetic moments, electric and magnetic fields leads to the conjecture that, certainly cross terms may contain new, intricate information for understanding chemical reactivity.