Conceptual density functional theory based electronic structure principles.

In this review article, we intend to highlight the basic electronic structure principles and various reactivity descriptors as defined within the premise of conceptual density functional theory (CDFT). Over the past several decades, CDFT has proven its worth in providing valuable insights into various static as well as time-dependent physicochemical problems. Herein, having briefly outlined the basics of CDFT, we describe various situations where CDFT based reactivity theory could be employed in order to gain insights into the underlying mechanism of several chemical processes. This journal is © The Royal Society of Chemistry.

Introduction

Ever since the proclamation by Dirac in 1929 that “The underlying physical laws necessary for the mathematical theory of a large part of physics and the whole of chemistry are thus completely known, and the difficulty is only that the exact application of these laws leads to equations much too complicated to be soluble,”[1] a great deal of efforts have been devoted by researchers over the last several decades in order to develop suitable numerical methods to solve the Schrodinger equation for atoms and molecules which in turn can help to shed light on various physicochemical processes. To this end, the formulation of density functional theory (DFT) is widely accepted to be a major development.[2-9] DFT provides a suitable theoretical alternative to ab initio wave-function based electronic structure methods albeit at a much cheaper computational cost. It is a well-established fact that for a system comprising N electrons and bound by an external potential v(r), the Hamiltonian (Ĥ) is completely defined by v(r) and N. Upon solving the Schrodinger equation for a known H, the many-electron wave function ψ(r1, r2…r) could be obtained which can provide all the information about the system under consideration. The single-particle density ρ(r) could be obtained upon performing integration over the coordinates of (N − 1) electrons,where

Therefore, it is evident that there exists a mapping in between v(r) and ρ(r). In the presence of an inverse mapping between ρ(r) and v(r), ρ(r) can determine v(r) which in turn can help to determine H and thus the many-electron wave function. Therefore, all the properties of the system could be evaluated. Hohenberg and Kohn[4] proved that there indeed exists an inverse mapping between ρ(r) and v(r), thus paving the way to develop and interpret chemical processes based on DFT. Since ρ(r) is a 3 dimensional quantity as opposed to the many-electron wave-function, it is convenient to visualise and provide chemically intuitive interpretations of various processes under consideration.[9] On the other hand, ρ(r) can help to provide classical interpretation of quantum mechanical processes. It is also possible to develop various theoretical models by making use of ρ(r). In this regard, conceptual density functional theory (CDFT), originally developed by Parr and collaborators,[9-16] constitutes a prominent theoretical framework, where ρ(r) is utilised to define several reactivity descriptors, based on which chemical processes are interpreted and understood. The reactivity of a system is basically its potential to react to an action. It measures the response of the system when acted upon by a perturbation caused by another reactant, a reagent, catalyst, solvent, external electric/magnetic field, confinement, change in temperature/pressure, etc. As the Hamiltonian of the system is totally characterized by the number of electrons and the external potential, the variation in energy or density by changing N and/or v(r) is utilised in understanding reactivity. In order to better understand the conceptual genesis of Parr's theoretical formulations, one needs to consider the second Hohenberg–Kohn theorem. According to the second Hohenberg–Kohn theorem,[4] the energy functional E[ρ] attains the minimum value for the true N-electron density for a given N and v(r) at the electronic ground state. Therefore, the density ρ(r) could be obtained via the variational optimization of the following Euler–Lagrange equation:

Herein, μ is the Lagrange multiplier associated with the normalization constraint (eqn (2)) and it is classified as the electronic chemical potential. Parr and co-workers proved that there exists a connection between μ and the thermodynamic chemical potential within a grand canonical ensemble at the zero temperature limit.[17,18] Furthermore, it was proposed that μ is the negative of electronegativity (χ), which helped to establish a direct connection with chemical reactivity.[17] One can consider the atoms in a molecule as a part of a grand canonical ensemble. In such a framework, the change in electron density could be analysed via DFT. At the electronic ground state, E[ρ] attains the minimum value. At a finite temperature for such a set of densities, the equilibrium state is characterized via the minimum Helmholtz free energy A[ρ] and grand potential functionals Ω[ρ] within the framework of canonical and grand canonical ensembles, respectively. The ground electronic state could be conceived as the zero temperature limit of the equilibrium state.[9] These concepts lie at the heart of the formulation of CDFT. Since chemical reactions accompany the change and reorientation of electron densities, several global and local reactivity descriptors could be defined within the premise of CDFT, which in turn can help to understand various physicochemical processes.[20,21] The backbone of CDFT is formed by these reactivity descriptors and associated electronic structure principles. The purpose of this review article is to introduce the general scientific community to these principles. Therefore, we do not intend to reiterate mathematical derivations and subtle technical issues within this article as exhaustive accounts are already available on these issues. Rather, the focus of this article is to highlight where and how one can employ these concepts to better understand a given chemical process. In the following sections of this article, we hope to outline the basic principles of CDFT and highlight various applications in static as well as time-dependent situations.

Global reactivity descriptors

As mentioned before, the Hamiltonian for a given system gets fixed upon specifying N and v(r). If a given system transforms from the ground electronic state to another, the energy change (up to the first order) associated with that process could be expressed as follows:[9,19]

Making use of the Euler–Lagrange equation for a fixed v(r), one can arrive at the following expression:

Herein, the Lagrange multiplier μ (electronic chemical potential) could be correlated with the corresponding thermodynamic quantity measuring the escaping tendency of an electron.[9] Upon considering a finite difference approximation to eqn (5), one can arrive at the following expression:

Here, χ, I and A represent the Mulliken electronegativity, ionization potential and electron affinity, respectively. Therefore, a direct connection with chemical reactivity could be established via this formulation. Upon employing further approximation by making use of Koopmans' theorem, one can arrive at the following expression:[22,23]

Therefore, it is possible to express χ in terms of the energies of the lowest unoccupied (ELUMO) and highest occupied (EHOMO) molecular orbitals. These quantities, i.e., χ and μ, are termed global reactivity descriptors as they can help to shed light on the reactivity of atoms, molecules and ions. Pearson introduced the concept of Hardness (η) within the premise of the hard–soft acid–base principle[24,25] and it could be expressed as follows:[26]

Softness (S)[27] could be expressed as the reciprocal of hardness as follows:

In general, the hardness of a system is associated with low polarizability as well as magnetizability. Parr et al. defined the electrophilicity index (ω)[28,29] as follows:

Here, ω is the measure of electrophilic power (bearing resemblance to classical electrostatics (=V2/R), where μ and η act as the potential (V) and resistance (R) respectively) of a system. Due to the discontinuity in energy versus N curves,[18,30] different definitions for the aforementioned quantities exist for the cases of charge acceptance and depletion.[31,32] Despite the limitations within the definitions of these quantities owing to the problem of this discontinuity, one can invoke the concept of atoms-in-a-molecule within the premise of a grand canonical ensemble at the zero temperature limit[33-39] in order to alleviate any inconsistencies.[37-40] Furthermore, one can also consider the isolated chemical entity as part of a larger system so that the aforementioned limitations in the definitions of the global reactivity descriptors could be eradicated.[37-40] It should be noted that these global reactivity descriptors could be used to describe the reactivity as well as the stability of a given chemical entity. In order to understand the site selectivity, several local reactivity descriptors are defined, which will be mentioned in the following section.

Local reactivity descriptors

In order to understand the activity of a specific atomic site within a molecule, several local reactivity descriptors were developed by the researchers. Electron density itself is the most fundamental local reactivity descriptor and we can denote it as follows:

Drawing inspiration from the spirit of Fukui's frontier orbital theory,[41,42] the Fukui function was developed and defined as follows:[43,44]

The Fukui function follows the normalization condition like the electron density. It also follows a cusp condition.[45-47] Owing to the discontinuity in the electron density versus N curves, several Fukui functions are defined for different chemical reactivity situations (electrophilic, nucleophilic and radical attacks, respectively).[43,44] One can employ the frozen core approximation in conjunction with a finite difference approximation, in order to arrive at the following expressions for the Fukui function for different chemical reactivity situations:

It could be seen from the aforementioned expressions that a direct connection could be drawn with Fukui's frontier orbital theory by defining the Fukui functions in this manner as these quantities are functions of frontier orbital density. The Fukui function assumes a rather large value in the cases where electron transfer is energetically favourable. In such situations, a concomitant large variation in the values of μ is observed thereby signifying greater chemical reactivity. In order to make these aforementioned quantities more convenient for chemical interpretations, the concept of atoms-in-a-molecule could be invoked and upon employing these formulations, one can arrive at the following condensed-to-atom variants of these quantities:[48]

Here, qK represents the electron population of the atom within a molecule. In order to compute qK, one could carry out Mulliken population analysis or Hirshfeld population analysis. It was shown that while using the Mulliken population analysis, the values of the Fukui functions could turn out to be negative.[49-51] So in general, Hirshfeld population analysis[49-51] is usually performed in order to calculate qK, which is positive in most cases and thus the Fukui functions. There are, however, situations where negative Fukui functions are qualitatively important to decipher the nature of some chemical processes.[52-55]

We note, however, that these local reactivity descriptors are poor intermolecular reactivity descriptors even though they could be useful in order to decipher the intramolecular reactivity. When two molecules interact with each other in the long range, the concerned molecules cannot ‘sense’ the local reactivity variations in each other. At long intermolecular separations, primarily electrostatic interactions dictate the intermolecular interactions as compared to the orbital interactions. In these situations, instead of the Fukui functions, one could utilize other reactivity descriptors such as local softness s(r)[56] as well as philicity ω(r).[57] By making use of the resolution of the identity associated with the normalization conditions for the Fukui functions, researchers have arrived at the following expressions for these two quantities:[56,57]

Condensed-to-atom variants for the aforementioned quantities were also developed and they are expressed as follows:

It has been shown that philicity could be utilized to understand nucleophile–electrophile interactions, whereas local softness can describe soft–soft interactions quite well. It was shown that one can arrive at the group quantities by summing over all condensed-to-atom quantities over the group of corresponding atoms.[58,59] The various local reactivity descriptors such as s(r), ω(r), ρ(r) and f(r) could be evaluated by following appropriate variational methods.

Attempts have been made to define a local variant of hardness under certain conditions.[60-62] Because of the inter-dependence of v(r) and ρ(r) within the premise of DFT, ambiguities arise[63-65] for the simultaneous definition of local hardness and local softness when local hardness is described as follows:[66,67]η(r) could be evaluated from the hardness kernel η(r,r′) also in the following way:

The hardness kernel could be expressed in the following way:

Here, F[ρ(r)] is the Hohenberg–Kohn universal functional. It should be noted that η(r) provides a measure of nuclear reactivity, whereas s(r) gives an estimate for the electronic reactivity.[68,69]

In addition to the aforementioned quantities, several other local reactivity descriptors are also utilized. Prominent examples among them include the gradient (∇ρ(r)) and Laplacian (∇2ρ(r)) of electron density,[70] quantum potential,[71,72] electron localization function,[73] molecular electrostatic potential,[74-77] multiphilic descriptor (Δω(r)),[78] dual descriptors (Δf(r)),[79-81]etc. A local temperature Θ(r) has been developed[82-86] by making use of the density and kinetic energy density within the realm of an ideal gas kinetic energy density approximation. This quantity has been also defined in the cases of time dependent situations as well as in the cases of the excited states having non-vanishing current density.[87-91] Even though a condensed-to-atom variant of Θ(r) has been developed, such a condensation suffers from a mathematical drawback.[92-98] As Θ(r) is non-linearly dependent on ρ(r), its condensation using the associated electron population becomes ad hoc. Attempts to define a local variant of hardness also suffer from this mathematical inexact formulation.[99] Since all the aforementioned quantities ultimately depend on the employed density partitioning schemes, these quantities need to be employed with careful scrutiny in order to alleviate any unphysical behaviour. In general, during the course of hard–hard interactions, atomic charge dictates[100] the nature of interaction and thus reactivity. Therefore, in such situations, local reactivity descriptors based on atomic charges can describe the reactivity quite well.[101,102] On the other hand, when the nature of interaction is dominated by frontier orbitals which happens to be the case in the instances of soft–soft interactions, one could employ descriptors such as the Fukui function[48] or the inter-molecular variants like s(r) and ω(r).

The aforementioned reactivity descriptors often obey some electronic structure principles. Thus, these reactivity descriptors are utilized in conjunction with these electronic structure principles. In the following sections, we intend to describe these principles, which in concomitance with the reactivity descriptors form the backbone of CDFT.

Electronic structure principles

The concept of electronegativity was introduced by Pauling.[103] It measures the relative propensity of the electron donation or acceptance capacity of two interacting moieties. It has been argued that during the course of flow of electrons between two interacting moieties, the process continues until their electronegativity values get equalized.[104-107] This is stated as Sanderson's electronegativity equalization principle. According to this principle, the final molecular electronegativity tends to equal the geometric mean of the electronegativity values of constituent atoms at the isolated state. During the course of formation of a complex such as X : Y, from an acid X and base Y, the extent of charge transfer (ΔN) and energy change (ΔE) could be expressed as follows:[26]

The aforementioned equations were originally proposed by Parr and Pearson. Since then, these equations have been invoked in various different contexts. These equations were utilized while seeking an analytical proof for the hard–soft acid–base principle (HSAB).[108-116] They were also used while defining an electrophilicity index.[117-119] However, this approach fails to take into account important factors such as entropy and solvent effects, electrostatic interaction, etc. and relies heavily on the effect of charge transfer.[120] This approach also suffers from the derivative discontinuity problem.[46-49] In order to alleviate these issues, efforts have been made by recasting this problem within the framework of a grand canonical ensemble.[33-39] Statistical mechanical apparatus may be properly exploited in defining reactivity descriptors at a finite temperature.[9] For example, the softness can alternatively be defined in terms of number fluctuations.[56] It is to be noted that in order to overcome the derivative discontinuity problem in the energy functional (E[N, v]), the incorporation of temperature within the definition of E[N, v] was considered by Franco-Perez, Ayers, Gázquez and Vela.[121-123] Within this approach, the premise of CDFT moves from the realm of canonical to a grand-canonical ensemble.[33-39] Within the premise of temperature dependent CDFT, the average electronic energy and its derivatives become the crucial component. Due to the incorporation of finite temperature effects, the deviations in the various response functions are negligible from their zero temperature counterparts (within the temperature independent formulation). In addition, within the premise of temperature dependent CDFT, quantities like heat capacity (and its local variant) could be defined.[122] Thus, energy transfer processes could be studied by making use of temperature dependent CDFT. The thermodynamic hardness and dual descriptor could also be defined by following temperature dependent CDFT.[123] These constitute important formal developments in CDFT.

In addition to the difference in electronegativity, the hardness sum is also an important quantity which can help to understand electron transfer processes. The hardness measures the reluctance of a chemical entity for electron transfer.[26] As the value of ΔE is always negative irrespective of the direction of electron transfer within a molecule, one needs to consider the ΔN values in order to understand the direction of the electron flow. In the cases of local reactivity descriptors, the condensed-to-atom variant of the Fukui function could be utilized in order to understand the direction of electron transfer.

In order to better understand the reactivity of a system, one needs to consider the concept of hardness as well, as electronegativity alone cannot describe the reactivity of chemical moieties. The concept of hardness was proposed by Pearson[24,25] and the utility of this concept was demonstrated while trying to explain the nature of acid–base reactions through his HSAB principle as follows: “Hard acids prefer to coordinate with hard bases and soft acids prefer to coordinate with soft bases for both their kinetic and thermodynamic properties.” By following the definition of hardness viaeqn (8), the hard–soft nature of a large number of acids and bases could be corroborated nicely with available experimental trends. An important electronic structure principle concerning the hardness is termed the maximum hardness principle (MHP).[124-130] According to the MHP, the following statement could be stated: “There seems to be a rule of nature that molecules arrange themselves so as to be as hard as possible.” There exist a number of processes where the validity of the MHP has been computationally proven. A few representative examples include the cases of various chemical reactions,[131] molecular vibrations as well as internal rotations,[132-137] Woodward–Hoffmann rules,[138,139] time dependent situations,[140] aromaticity,[141] chaotic ionization from Rydberg states,[142]etc.

In view of the less magnetizability and polarizability of hard systems, two additional electronic structure principles namely the minimum magnetizability principle (MMP) and minimum polarizability principle (MPP) are formulated. According to the MPP,[143] the following could be stated: “The natural direction of evolution of any system is towards a state of minimum polarizability (α)” whereas the MMP states that[144,145] “A stable configuration/conformation of a molecule or a favourable chemical process is associated with a minimum value of the magnetizability (ξ).”

By making use of the η and μ values, the extremum values of electrophilicity (ω) could be calculated.[146-148] It has been proposed that a minimum electrophilicity principle (MEP) might be operative in several cases.[149-152] By drawing comparison with classical electrostatics, it has been argued[153] that the MEP is inherently defined within the Maynard-Parr definition of ω. During the course of chemical reactions, extrema of ω occur at points where the following condition is satisfied:[146-148]

Herein, λ could be either a reaction coordinate for a chemical reaction or some internal degree of freedom for the cases of vibration or internal rotational processes. We note that η is always positive and μ is always negative because of the convex nature of energy. Therefore, electrophilicity achieves the extremum value, when the dependencies of η and μ on λ are opposed.

In addition to the aforementioned principles, it might be prudent to mention how the concept of entropy is defined and utilized within the realm of DFT. Considering an N electron system and assuming that it consists of N non-interacting particles under the influence of an effective potential veff(r,t), entropy density in a time-dependent situation is formulated within an average density framework as follows:[87,154]

Herein, k is the Boltzmann constant, whereas θ represents a space-time dependent ‘temperature’. In order to define θ, the kinetic energy density is utilized as follows:

Here, j represents the current density. We note that by making use of information theory, one can arrive at eqn (29) as well. To this end, the Shannon entropy for the system could be expressed in terms of the density of the system and one could maximize it under certain constraints in order to arrive at eqn (29). The global entropy could be obtained by integrating eqn (29) over the whole space as follows:

In several time-dependent situations, it was shown that in the cases where the processes are favourable, the entropy gets maximized and a maximum entropy principle remains operative in such situations.[142,143] It should be noted that within the realm of CDFT, concepts derived from information theory[155] could be utilized and suitable descriptors could be defined to provide several important insights. To this end, the utilization of information theory has provided physical insights into atoms-in-molecule theory.[156] It has been shown that Shannon entropy, when expressed as a functional of ρ(r), can on its own describe coulombic moieties.[157] The implication of this observation is profound. It could be argued that the Shannon entropy could turn out to be as important as ρ(r) and it could be utilized to explain several chemical processes. In this context, several other density functionals are also defined such as Onicescu information energy,[158] Fisher information[159] and Renyi entropy.[160] The Kullback–Leibler information measure also contains meaningful information.[89] Several important insights have been provided into many physicochemical processes (such as the steric effect, aromaticity, etc.) by employing these concepts.[161-165] The information conservation principle could be utilized to quantify important reactivity descriptors such as electrophilicity, nucleophilicity and regioselectivity.

There are several reactions associated with a change in spin multiplicity. Spin dependent reactivity descriptors have been derived.[166-168] The above-mentioned quantities are separately defined for up and down spins. Related chemical potential, spin potential and different types of hardnesses are made use of.

Reactivity descriptors within time dependent situations

Runge and Gross showed[169] that it is possible to extend DFT as developed by Hohenberg and Kohn, to time-dependent situations. They proved that there exists a one-to-one correspondence between the time-dependent external potential vext(r,t) and time-dependent density ρ(r,t) in the cases of a many-electron moiety evolving from a given initial state. All the time dependent properties of a given system are unique functionals[170] of time-dependent density ρ(r,t) and the current density j(r,t). The set of time dependent Kohn–Sham equations could be expressed as follows:

Here,and,

Inspired by the work of Madelung on quantum fluid dynamics,[171] a different route could be considered in order to evaluate ρ(r,t) and j(r,t). This approach is termed quantum fluid density functional theory (QFDFT).[87-89] Within this framework, a generalized nonlinear Schrodinger equation (GNLSE) is solved to study the temporal evolution of a system and this GNLSE could be expressed as follows:

Here,andHerein, ξ denotes the velocity potential. Other alternatives like the use of dynamical response functions or direct molecular dynamics simulations may also be envisaged.

It should be noted that the Hohenberg–Kohn theorems have been further developed in order to incorporate the effects of external electric and magnetic fields.[172] In this formulation, the kinetic energy operator, within the time dependent Kohn–Sham equations, would include a term arising from the vector potential. Furthermore, j(r,t) and veff[ρ(r,t)] would include terms arising from both vector and scalar potentials. These additional terms arise in order to incorporate the effects of external time-dependent electric and magnetic fields. These additional terms are also taken care of within the QFDFT framework.

The time dependent variants of both global and local reactivity descriptors could be defined so that these quantities could be evaluated within the time dependent scenario.[143,173] To this end, the chemical potential μ(t) could be defined by making use of j(r,t) and ρ(r,t) within eqn (3). One, however, needs to define a suitable time-dependent energy functional. Time-dependent hardness η(t) could be calculated via the time-dependent variants of eqn (23)–(25). Therefore, it becomes possible to evaluate the time-dependent electrophilicity ω(t) by making use of η(t) and μ(t) viaeqn (10).

Having outlined the basics of CDFT, in the following sections we intend to provide an overview of processes where CDFT turns out to be very useful in understanding various physicochemical processes.

Electronic structure principles under special conditions

Before we describe the utility of CDFT in understanding several special physicochemical processes, it might be prudent to briefly state a few general trends vis-à-vis chemical reactivity as observed within the purview of CDFT. Firstly, we consider the case of reactivity of atoms across the periodic table. In general, along a given period, hardness increases whereas it decreases along a given group.[174] This fact is in conformity with the MHP. Similarly, electrophilicity,[175] polarizability[175] and magnetizability[144,145] increase along a group whereas the values of these quantities decrease along a period. These facts are in accordance with principles such as the MEP, MMP and MPP. Therefore, the utility of global reactivity descriptors in explaining chemical reactivity becomes quite evident from the aforementioned facts. One can utilize the radial distribution function of philicity,[176]i.e., 4πr2ω(r), in order to retrieve the characteristic features of the electronic shell structure much like the plot of 4πr2ρ(r). One can define the following cusp condition for ω(r) from that of f(r). We note that for all r in a molecule, ω remains a constant and thus the following equation could be derived:[45,46]

Utilizing eqn (38), one could construct a map of ω(r), which in turn can help to find the nuclei within a molecule as well as the loci of atoms and thus help to shed light on the atoms-in-molecule picture.

Similarly, global and local reactivity descriptors could be utilized to understand the change in reactivity and stability associated with the changes in molecular vibration and internal rotation. We note that due to the changes in molecular vibration and internal rotation, the electron density distribution within a molecule also changes. In general,[132-136] within a minimum energy configuration and/or conformation, the hardness value becomes maximum, whereas polarizability tends to attain the minimum value. In a number of cases, electrophilicity attains the minimum value in the corresponding situations. In order to understand the nature of the extremals in ω, one may consider the characteristics of the first and second order derivatives of this quantity with respect to a generalized reaction coordinate (λ).[146-148] One can differentiate both sides of eqn (10) with respect to λ as follows:and

One could infer from the expressions stated above that ω would be an extremum when η and μ attain the respective extrema. Even though this criterion is a sufficient condition, it does not always have to be a necessary condition. During the course of molecular vibrations and internal rotations, the local reactivity descriptors also exhibit characteristic oscillations.[133-136]

In order to analyse the reaction between an electrophile and nucleophile, one can utilise local and global electrophilicities. An electrophile would have a larger ω value as compared to a nucleophile and as a result of this difference, they will be drawn towards each other. Once these moieties come in close proximity to one another, several other local reactivity descriptors could be utilized to understand the nature of interaction between them. It should be noted that during the course of the interaction between the electrophile and nucleophile, the sites having the largest local electrophilicity and nucleophilicity values would primarily dictate the process. Since , local philicity values at the concerned atomic sites govern the corresponding global reactivity trend.

In order to analyse the nature of cycloaddition reactions, CDFT based reactivity descriptors have been utilized by several research groups. To this end, the difference in electrophilicity (Δω) between a diene and a dienophile could be employed. It has been suggested[177-179] that in the cases when the Δω is small between a diene–dienophile pair, the reaction mechanism follows a concerted pathway, and that the reaction could be primarily understood by considering the properties of the concerned frontier orbitals. Therefore, these types of reactions occur as a result of soft–soft interactions. In the cases when the Δω is large for a diene–dienophile pair, the reactions proceed via a polar and step-wise mechanism. These reactions are primarily charge-controlled and thus could be classified as an example of hard–hard interactions. These trends could be understood by employing philicity also.[180] Electrophilicity has been utilized in order to understand the salient features of acidity/basicity concepts as defined by the Lewis and Brønsted–Lowry models.[181,182] Similarly, electrophilicity has been employed in order to understand several concepts prevalent in organic chemistry, details of which could be found elsewhere.[20] We state another prominent example where the utility of the electrophilicity concept could help to elucidate many facets of organic reactions. Maynard and co-workers showed[28] that there exists a connection between electrophilicity and kinetic behaviour. Thus, this concept could be employed to understand chemical kinetics. Parr et al. developed[29] the concept of electrophilicity by employing an energy minimization criterion. Therefore, one could gain insights into the thermodynamic stability of a given system under consideration. In this context, the variations in electrophilicity along the intrinsic reaction coordinate (IRC) of several reactions have been analysed in great detail.[146-148] In these cases, the minimum electrophilicity principle remains operative. In order to understand this observation, one could make use of eqn (39) and (40). For the minimum energy configurations, the following conditions hold: , whereas for the transition states the corresponding conditions could be stated as follows: .

Local reactivity descriptors demonstrate opposite trends for bond breaking and bond making processes. They tend to intersect at some point along the IRC, in the cases of several reactions.[140,141] In the cases of thermoneutral reactions, these intersection points coincide with the saddle points within the potential energy surfaces. Therefore, by employing the outlined method, one can locate the transition states. In accordance with a reactivity based Hammond's postulate, the intersection points succeed the transition states in the cases of exothermic reactions, whereas these points precede the transition states in the cases of endothermic reactions.[183,184] In a number of exothermic reactions, it has been proven[151-153] that the least electrophilic moiety lies in the product side. The average electrophilicity value generally decreases in these cases as we move from the reactant to the product's end. This observation could be classified as an outcome of the MEP. In general, for a favourable reaction, the ΔηProduct−Reactant value should be a positive quantity, whereas Δα,E,ωProduct−Reactant should be a negative quantity.[185-187]

Change in chemical reactivity indices in the presence of solvents

Since the reactivity of systems varies as we move from a gas phase to a solution phase, it is extremely important to know how the CDFT based reactivity descriptors behave in such situations. Lipinski and co-workers considered a homogeneous polar medium using a virtual charge based model in order to understand the changes in the hardness and electronegativity of some bonded moieties.[188] They have showed that upon increasing the solvent polarizability, the hardness of ions decreases, whereas the electronegativity increases and decreases, respectively, in the cases of anions and cations. However, the molecular level hardness and electronegativity exhibited only minor variations as a function of the solvent effects. According the models considered by Pearson[189] in the cases of hydration, neutral molecules become less hard upon solvation. Anions and cations, on the other hand, become poorer electron donors and acceptors (so more and less electronegative, respectively) respectively upon hydration. Sivanesan et al.[190] analysed the changes in the Fukui function in going from a gas to a water solvent phase in the cases of several organic molecules and they showed that in general the propensity of nucleophilic and electrophilic attacks increases in moving from a gas to a solvent phase. The change in electrophilicity upon solvation was analysed in several studies.[191-194] It has been suggested that solvation enhances the electrophilic power of neutral electrophilic moieties, whereas the opposite effect seemingly operates in the cases of charged and ionic electrophiles. Upon solvation, the insertion energy could be approximated by the changes in Δη and Δμ. While the former was shown to be times the solvation energy, the latter quantity was shown to be twice the solvation energy. Therefore, the changes in electrophilicity are linearly correlated with the solvation energy. The global electrophilicity is much more sensitive to solvation effects as compared to its local counterpart and thus it could be utilized to understand changes in chemical reactivity.

Change in chemical reactivity indices in the presence of an external field

In order to understand how a many-electron system behaves in the presence of an external perturbation, both perturbative and variational methods were adopted.[195-198] We note that the external perturbation might originate from the presence of another molecule or one could utilize electromagnetic radiation in order to perturb the system under consideration. Numerical results[199] have shown that ω and μ are much more prone to getting affected by some external field as compared to η. Local reactivity descriptors are much more susceptible to changes as a function of the strength of the external fields as compared to their global counterparts.

QFDFT[87,88] has been utilized extensively in order to understand the changes in CDFT reactivity descriptors as a function of external perturbation. The temporal evolution of reactivity descriptors has been analysed therein. In the presence of a low-intensity field, in-phase oscillations in the time evolution of μ have been observed, whereas one needs to employ a much higher intensity field in order to observe analogous behaviour in the case of η.[199] The competition between the axial external field and spherical nuclear Coulomb field could be nicely demonstrated with this approach. Furthermore, the chaotic dynamics of Rydberg atoms in the cases of ionization was analysed in these studies.[142]

In addition to external fields, ion-atom collision dynamics was also analysed within the QFDFT framework. We highlight only the salient features herein.[200] The time dependent μ-profile clearly demonstrates the three scattering regions during the course of the reactive collisions. It was shown that the hardness is maximized and the polarizability is minimized in the interaction region. These results are in accordance with the time dependent variants of the MHP and MPP. The changes in ω resembled the corresponding case of μ. Herein, the MEP was not observed to be operative in the sense that we have discussed in the preceding sections in the cases of static situations. The time dependent variant of the HSAB principle was proven in these studies. Insight was provided into the regioselectivity of some reactions as well.

Change in chemical reactivity indices in a confined environment

Confined quantum systems exhibit fascinating changes in their physicochemical properties as compared to their counterparts in the free state, details of which could be found elsewhere.[201-212] In order to understand the effect of confinement on the reactivity of atoms and ions, the Hartree–Fock–Slater equation has been considered in order to obtain a self-consistent field (SCF) electronic wave function. In order to impose the effect of confinement, Dirichlet boundary conditions have been imposed on the electronic wave function. The CDFT based reactivity descriptors have been analysed from this SCF wave function.[206,207] As one imposes the effect of confinement, systems become harder. The polarizability decreases monotonically as the confinement radius is reduced. It is suggested that electronegativity is not very sensitive to the effect of confinement. In the case of ω, a similar behaviour has been noted. It is suggested that as we keep increasing the effect of confinement, systems tend to become less polarizable and harder. Under these conditions, systems also become more difficult to excite.

Utilizing the theoretical framework of QFDFT, simultaneous effects of confinement and electronic excitations have been analyzed.[208-212] Upon increasing the degree of confinement, hardness increases both in the ground and the excited electronic states. Polarizability decreases as compared to the free atom due to the effect of confinement. As expected from the MPP, in both confined as well as unconfined states, the excited state polarizability is higher as compared to that in the ground state. Chattaraj et al. proved the validity of the HSAB principle in the cases of atoms and hydrogen molecules exposed to an external magnetic field in the confined state.

In order to show some representative plots for the time evolution of some reactivity descriptors discussed above during the course of atom–field interaction as well as during collisional processes, we have presented Fig. 1–5 for the purpose of illustration.

Several important insights have been provided by Morell and co-workers on the excited state properties of several chemical systems as well.[213,214]

Change in chemical reactivity indices upon electronic excitation

As any system becomes more reactive upon excitation, it is expected from the MHP and MPP that it will be softer and more polarizable in its excited state than in the ground state. It has been explicitly demonstrated to be the case in the case of atoms, molecules and ions.[215,216] In a two-state ensemble also systems become softer and more polarizable when the contribution from the excited state increases in the complexion of this ensemble.[215,216] A favorable collision occurs when a hard projectile hits a target in its ground state rather than in the excited state where it is softer and more polarizable, as is expected from the HSAB principle.[200] These features are valid in a confined environment as well and a suitable choice of simultaneous application of confinement and excitation during ion–atom collision and atom–field interaction can bring back the ground state behavior of a corresponding free state.[208,209,217] As electrophilicity is a composite reactivity descriptor, it shows the validity of the MEP when both chemical potential and hardness behave properly. It has been argued that the excited state reactivity may be understood by analyzing the change in the electron density distribution upon electronic excitation[218,219]

Quantitative structure–toxicity relationship

Electrophilicity and its local counterpart have been shown to be important descriptors of cytotoxicity.[220-223] A small number of descriptors can provide robust quantitative structure–toxicity relationship (QSTR) models. A multiple linear regression analysis has also highlighted that when they are used in conjunction with other popular descriptors like hydrophobicity, useful QSTR models are generated which in turn will help screening the initial drug targets saving time and money. The robustness of the models is analyzed by employing multilayer perceptron neural networks.[224,225]

Aromaticity

The relative and absolute hardness could be utilized in order to understand the stability of aromatic molecules. In view of the fact that both hardness and aromaticity provide means to understand the stability and reactivity of a system, the correlation between them could be intriguing. Nucleus independent chemical shift (NICS)[226] values could be utilized in order to understand the aromatic or anti-aromatic properties of a molecule. An aromaticity descriptor could be defined as follows:[227]Yaromatic = Ymolecule − Yreference; Y = E, α, ω, η. The reference molecule could be either the corresponding localized molecule or it could be a corresponding open chain system. Aromatic molecules possess a positive[228,229]ηaromatic value, whereas the other descriptors mentioned above possess negative values, according to the CDFT based electronic structure principles as discussed above. In the cases of anti-aromatic molecules, these trends are reversed. A number of all-metal and non-metal molecules and their aromatic behaviour have been analyzed by following these criteria.[230-233] Aromaticity can drive a system towards its stable structure as can be understood through the dynamical variants of the CDFT based electronic structure principles.[234] They are also helpful in understanding the aromaticity of compounds of multivalent superatoms.[235] Recently, a similar approach has been adopted in explaining the stability of metal clusters.[236]

Conclusion

Based on the discussions stated in this article, we can say that CDFT could be utilized in several physicochemical contexts. Global and local reactivity descriptors can help to shed light on stability, reactivity, dynamics, etc. Global reactivity descriptors are connected with several electronic structure principles and thus are very important. Local reactivity descriptors can help to provide insights into site selectivity. Despite a lot of computational evidence in support of CDFT, it is necessary, formally, to understand why and how CDFT works. To this end, more mathematical analysis is needed especially to highlight the domain of applicability of these electronic structure principles.[237] This in turn can increase the predictive power of CDFT going beyond its usual interpretive value.

Author contributions

PKC supervised all of the relevant research work emanating from his research group . He worked out the overall plan of this review article as well. DC wrote the first draft of this article.

Conflicts of interest

The authors declare no competing financial interest.