The Fukui matrix: a simple approach to the analysis of the Fukui function and its positive character

The Fukui matrix is introduced as the derivative of the one-electron reduced density matrix with respect to a change in the number of electrons under constant external potential. The Fukui matrix extends the Fukui function concept: the diagonal of the Fukui matrix is the Fukui function. Diagonalizing the Fukui matrix gives a set of eigenvectors, the Fukui orbitals, and accompanying eigenvalues. At the level of theory used, there is always one dominant eigenvector, with an eigenvalue equal to 1. The remaining eigenvalues are either zero or come in pairs with eigenvalues of the same magnitude but opposite sign. Analysis of the frontier molecular orbital coefficient in the eigenvector with eigenvalue 1 gives information on the quality of the frontier molecular orbital picture. The occurrence of negative Fukui functions can be easily interpreted in terms of the nodal character of the dominant eigenvector versus the characteristics of the remaining eigenvectors and eigenvalues.     

Atom and Bond Fukui Functions and Matrices: A Hirshfeld‐I Atoms‐in‐Molecule Approach

The Fukui function is often used in its atom‐condensed form by isolating it from the molecular Fukui function using a chosen weight function for the atom in the molecule. Recently, Fukui functions and matrices for both atoms and bonds separately were introduced for semiempirical and ab initio levels of theory using Hückel and Mulliken atoms‐in‐molecule models. In this work, a double partitioning method of the Fukui matrix is proposed within the Hirshfeld‐I atoms‐in‐molecule framework. Diagonalizing the resulting atomic and bond matrices gives eigenvalues and eigenvectors (Fukui orbitals) describing the reactivity of atoms and bonds. The Fukui function is the diagonal element of the Fukui matrix and may be resolved in atom and bond contributions. The extra information contained in the atom and bond resolution of the Fukui matrices and functions is highlighted. The effect of the choice of weight function arising from the Hirshfeld‐I approach to obtain atom‐ and bond‐condensed Fukui functions is studied. A comparison of the results with those generated by using the Mulliken atoms‐in‐molecule approach shows low correlation between the two partitioning schemes.     

Tuning electronic eigenvalues of benzene via doping

Using variable atomic numbers within molecular grand-canonical ensemble theory, the highest occupied Kohn-Sham eigenvalue of isoelectronic benzene derivatives is tuned. The performed transmutational changes correspond to the iterative doping with boron and nitrogen. The molecular Fukui function proves to be a reliable index in order to predict the changes in the highest occupied molecular orbital eigenvalue due to doping.     

Removing electrons can increase the electron density: a computational study of negative Fukui functions

Ab initio and density-functional theory calculations for a family of substituted acetylenes show that removing electrons from these molecules causes the electron density along the C-C bond to increase. This result contradicts the predictions of simple frontier molecular orbital theory, but it is easily explained using the nucleophilic Fukui function-provided that one is willing to allow for the Fukui function to be negative. Negative Fukui functions emerge as key indicators of redox-induced electron rearrangements, where oxidation of an entire molecule (acetylene) leads to reduction of a specific region of the molecule (along the bond axis, between the carbon atoms). Remarkably, further oxidization of these substituted acetylenes (one can remove as many as four electrons!) causes the electron density along the C-C bond to increase even more. This work provides substantial evidence that the molecular Fukui function is sometimes negative and reveals that this is due to orbital relaxation.     

Correlated one-particle method: numerical results

In a previous paper a correlated one-particle method was formulated, where the effective Hamiltonian was composed of the Fock operator and a correlation potential. The objective was to define a correlated one-particle theory that would give all properties that can be obtained from a one-particle theory. The Fock-space coupled-cluster method was used to construct the infinite-order correlation potential, which yields correct ionization potentials (IP's) and electron affinities (EA's) as the negative of the eigenvalues. The model, however, was largely independent of orbital choice. To exploit the degree of freedom of improving the orbitals, the Brillouin-Brueckner condition is imposed, which leads to an effective Brueckner Hamiltonian. To assess its numerical properties, the effective Brueckner Hamiltonian is approximated through second order in perturbation. Its eigenvalues are the negative of IP's and EA's correct through second order, and its eigenfunctions are second-order Brueckner orbitals. We also give expressions for its energy and density matrix. Different partitioning schemes of the Hamiltonian are used and the intruder state problem is discussed. The results for ionization potentials, electron affinities, dipole moments, energies, and potential curves are given for some sample molecules.     

Comparison of DFT methods for molecular orbital eigenvalue calculations

We report how closely the Kohn-Sham highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) eigenvalues of 11 density functional theory (DFT) functionals, respectively, correspond to the negative ionization potentials (-IPs) and electron affinities (EAs) of a test set of molecules. We also report how accurately the HOMO-LUMO gaps of these methods predict the lowest excitation energies using both time-independent and time-dependent DFT (TD-DFT). The 11 DFT functionals include the local spin density approximation (LSDA), five generalized gradient approximation (GGA) functionals, three hybrid GGA functionals, one hybrid functional, and one hybrid meta GGA functional. We find that the HOMO eigenvalues predicted by KMLYP, BH&HLYP, B3LYP, PW91, PBE, and BLYP predict the -IPs with average absolute errors of 0.73, 1.48, 3.10, 4.27, 4.33, and 4.41 eV, respectively. The LUMOs of all functionals fail to accurately predict the EAs. Although the GGA functionals inaccurately predict both the HOMO and LUMO eigenvalues, they predict the HOMO-LUMO gap relatively accurately (approximately 0.73 eV). On the other hand, the LUMO eigenvalues of the hybrid functionals fail to predict the EA to the extent that they include HF exchange, although increasing HF exchange improves the correspondence between the HOMO eigenvalue and -IP so that the HOMO-LUMO gaps are inaccurately predicted by hybrid DFT functionals. We find that TD-DFT with all functionals accurately predicts the HOMO-LUMO gaps. A linear correlation between the calculated HOMO eigenvalue and the experimental -IP and calculated HOMO-LUMO gap and experimental lowest excitation energy enables us to derive a simple correction formula.     

Diagrammatic formulation of the second‐order many‐body multipartitioning perturbation theory

The second‐order multireference perturbation theory employing multiple partitioning of the many‐electron Hamiltonian into a zero‐order part and a perturbation is formulated in terms of many‐body diagrams. The essential difference from the standard diagrammatic technique of Hose and Kaldor concerns the rules of evaluation of energy denominators which take into account the dependence of the Hamiltonian partitioning on the bra and ket determinantal vectors of a given matrix element, as well as the presence of several two‐particle terms in zero‐order operators. The novel formulation naturally gives rise to a “sum‐over‐orbital” procedure of correlation calculations on molecular electronic states, particularly efficient in treating the problems with large number of correlated electrons and extensive one‐electron bases. ©1999 John Wiley & Sons, Inc. Int J Quant Chem 73: 395–401, 1999     

Modeling temporary anions in density functional theory: calculation of the Fukui function

Two approaches are investigated for modeling electron densities of temporary anions in density functional theory (DFT). Both rely on an artificial binding of the excess electron, in one case by a compact basis set and in the other by a potential wall. The key feature of the calculations is that the degree of binding is controlled in both cases by knowledge of the negative electron affinity of the corresponding neutral, approximated in terms of DFT local functional frontier orbital eigenvalues and vertical ionization potential, A=-(epsilon(LUMO)+epsilon(HOMO))-I. To illustrate the two approaches, Fukui functions for nucleophilic attack are determined in four molecules with increasingly negative electron affinities. They yield very similar results, which are notably different to those determined without artificial electron binding. The use of a potential wall has the attractive feature that large, diffuse basis sets can be used, avoiding the need for a compact basis, tailored to a particular molecule.     

Natural orbital Fukui function and application in understanding cycloaddition reaction mechanisms

A new condensed form of the Fukui function, the natural orbital Fukui function (NOFF), is proposed and derived from natural bond orbital occupancy. It is defined as the change in natural bond orbital occupancy upon electronic perturbation (electron addition to, or depletion from, a molecular system). Applying NOFF to a series of cycloaddition reactions (e.g., [4 + 2] and [2 + 1] cycloadditions) illustrates the effectiveness of the concept in interpreting bond breakage and formation mechanisms.     

Elimination, in electronic structure calculations, of redundant orbital products

We propose a direct method for reducing the dimension of the space of orbital products that occur, for example, in the calculation of time dependent density functional theory linear response and in Hedin's GW approximation to the electron propagator. We do this by defining, within the linear space of orbital products, a subspace of dominant directions that are associated with a certain eigenvalue problem. These directions span the entire linear space of products with an error that decreases approximately exponentially with their number. Our procedure works best for atomic orbitals of finite range and it avoids the use of extra sets of auxiliary fit functions.     

Bond fukui indices: Comparison of frozen molecular orbital and finite differences through mulliken populations

Bond Fukui functions and matrices are introduced for ab initio levels of theory using a Mulliken atoms in molecules model. It is shown how these indices may be obtained from first‐order density matrix derivatives without need for going to second‐order density matrices as in a previous work. The importance of taking into account the nonorthogonality of the basis in ab initio calculations is shown, contrasting the present results with previous work based on Hückel theory. It is shown how the extension of Fukui functions to Fukui matrices allows getting more insight into the nature of bond Fukui functions. All presently introduced indices respect the necessary normalization conditions and include the classical single atom condensed Fukui functions. © 2013 Wiley Periodicals, Inc.     

Theoretical study of the interaction between Pt(0) and MPH3 + fragments in complexes of the [Pt3 (��-CO)3(PH3)3]�CMPH3 + (M?=?Cu+, Au+, Ag+) type

Ab initio calculations suggest that a series of complexes of the [Pt₃(??-CO)₃(PH₃)₃]?CMPH₃ ⁺ type (M?=?Cu, Au, Ag) are stable. We have studied these complexes at the HF, MP2, B3LYP, and PBE levels of theory. The magnitude of the interaction energies and Pt₃?CM distances indicate a substantial covalent character of the bond, the latter being confirmed by orbital diagrams. The chemical bond is sensitive to electron correlation effects. In addition, the Fukui index of nucleophilic attack and electrophilicity index on the metal were used to explore possible sites where chemical reactivity may play a role.     

A direct,restricted-step,second-order MC SCF program for large scale ab initio calculations

A general purpose MC SCF program with a direct, fully second-order and step-restricted algorithm is presented. The direct character refers to the solution of an MC SCF eigenvalue equation by means of successive linear transformations where the norm-extended hessian matrix is multiplied onto a trial vector without explicitly constructing the hessian. This allows for applications to large wavefunctions. In the iterative solution of the eigenvalue equation a norm-extended optimization algorithm is utilized in which the number of negative eigenvalues of the hessian is monitored. The step control is based on the trust region concept and is accomplished by means of a simple modification of the Davidson—Liu simultaneous expansion method for iterative calculation of an eigenvector. Convergence to the lowest state of a symmetry is thereby guaranteed, and test calculations also show reliable convergence for excited states. We outline the theory and describe in detail an efficient implementation, illustrated with sample calculations.     

Activation Energy and Transition State Determination of the Olefin Insertion Process of Metallocene Catalysts Using a Semiempirical Molecular Orbital Calculation

The transition state of the olefin insertion process of metallocene catalysts can be determined by adopting the semiempirical PM3 model. In computational chemistry, the computational methods most employed are the ab initio method and density functional theory, which are very time consuming. The semiempirical molecular orbital method requires much less computational resources than the above methods. However, the accuracy and reliability of the semiempirical molecular orbital method remains to be determined. The PM3 model is the most recently developed the semiempirical molecular orbital method and can also be applied to transition metal calculations. This study is intended to investigate the reliability of computational results determined using semiempirical PM3 model on metallocene catalysts through comparison with published results on the density functional theory (DFT). The saddle point finding procedure is adopted to find the transition state of the ethylene insertion process of metallocene catalysts. Results on the geometry and energy trends of the ethylene insertion process of metallocene catalysts determined using the PM3 model are in good agreement with the DFT results. In addition, the saddle point of the potential energy surface of ethylene insertion is verified in accordance with the eigenvalue of the vibrational frequency spectrum. Correct eigenvalues indicate that the correct saddle point of the potential energy surface of ethylene insertion has been successfully located. Hence, the eigenvalue of the vibrational frequency spectrum is a valuable reference in terms of saddle point justification. Computational results and vibrational frequency spectrum analysis demonstrate that the PM3 model can be used to locate the correct saddle point of the potential energy surface. The results obtained using the PM3 model confirm that the eigenvalue of the transition state lies nearly on the vibrational frequency spectrum. The eigenvalues are also analyzed, providing a valuable reference for further studies of the transition state of olefin insertion of metallocene catalysts. The activation energies for the olefin insertion reaction are also studied for evaluation of the catalyst.     

Computing the Fukui function from ab initio quantum chemistry: approaches based on the extended Koopmans’ theorem

The extended Koopmans’ theorem is related to Fukui function, which measures the change in electron density that accompanies electron attachment and removal. Two approaches are used, one based on the extended Koopmans’ theorem differential equation and the other based directly on the expression of the ionized wave function from the extended Koopmans’ theorem. It is observed that the Fukui function for electron removal can be modeled as the square of the first Dyson orbital, plus corrections. The possibility of useful generalizations to the extended Koopmans’ theorem is considered; some of these extensions give approximations, or even exact expressions, for the Fukui function for electron attachment.     

Improved convergence of self-consistence procedures in the MC SCF theory

To obtain optimized orbitals within the MC SCF theory, the energy surface near a chosen point is approximated by a quadratic function of independent matrix elements of a small orthogonal orbital transformation. The method of a second-order one-electron Hamiltonian (OEH) is developed on the basis of this approximation. A procedure is proposed to define step coordinates, insuring a rapid descent along an average-energy surface also in the cases when the matrix of second energy derivatives has eigenvalues negative or close to zero. The results obtained in applying the OEH method for the calculation of ground and triplet states of uracile in the π-electron approximation are discussed. When a complete matrix of the second energy derivatives is used, the self-consistence procedure is quadratically convergent. An exponential, yet rapid enough convergence is provided by a simplified computation scheme neglecting cross derivatives.     

Electronegatlvity and the bonding character of Molecular Orbitals

The density-functional approach to Molecular Orbital theory shows that the chemical bonding potential is better described by orbital electronegativities than by ionization energies. This results from the fact that the electronic relaxation connected with ionization is not significant for the homolytic breaking of chemical bonds. Electronegativity, on the other hand, is an eigenvalue corresponding to the average potential seen by an electron as a molecular orbital changes into monocentric (atomic) orbitals.     

A Note on Nonlinear Parameters in Variational Methods

The dependence of the eigenvalues of matrices representing a quantum-chemical Hamiltonian in a model space on the nonlinear parameters of the trial functions is analyzed. Several theorems useful in determining the dependence of the matrix eigenvalues on the parameters are presented and their implications on the choice of the orbital basis sets are briefly discussed. A simple method of optimization of the parameter values is formulated.     

The relationship between the eigenvalues and eigenvectors of a similarity matrix and its associated Carbó index matrix

The eigenvalues and eigenvectors of a quantum similarity matrix are also generalized eigenvalues and eigenvectors of the associated matrix of Carbó indices. This establishes bounds on the spectrum of the Carbó index matrix; for example, a quantum similarity matrix is positive semidefinite if and only if the associated Carbó index matrix is also positive semidefinite. The generalized eigenvalue problem for the Carbó index matrix has a diagonal metric matrix on the right-hand-side. Every generalized eigenvalue problem can be written in this diagonal form (i.e., this form is not special to this application). This diagonally structure generalized eigenvalue problem is especially convenient because it can be converted to a conventional eigenvalue problem by a particularly simple partial Löwdin transformation.