Paul W. Ayers

The physical basis of the hard/soft acid/base principle

The dependence of molecular properties on the chemical hardness is explored. As the chemical hardness of a molecule increases, its size and polarizability typically decrease and its charge and electronegativity typically increase. On the basis of these properties, the interaction energy between hard and soft acids and bases is quantified, and the physical basis of the global and local hard/soft acid/base (HSAB) principles is elucidated.

Variational principles for describing chemical reactions: Condensed reactivity indices

Two recent papers [P. W. Ayers and R. G. Parr, J. Am. Chem. Soc. 122, 2010 (2000); 123, 2007 (2001)] have shown how variational principles for the energy may be used to derive and elucidate the significance of the chemical reactivity indices of density-functional theory. Here, similar ideas are applied, yielding a systematic, mathematically rigorous, and physically sound approach to condensed reactivity indices. First, we use the variational principle for the energy to derive an expression for the condensed Fukui function index in terms of the condensed hardness kernel. Next, we address an important open problem pertaining to condensed reactivity indices: when (if ever) is the condensed Fukui function for an atom in a molecule negative? In particular, our analysis confirms the observation, hitherto based only on computational evidence, that the Hirshfeld partitioning is optimal for obtaining non-negative Fukui functions. We also hypothesize that the strong diagonal dominance of the condensed hardness kern...

Chemical Reactivity Descriptors for Ambiphilic Reagents: Dual Descriptor, Local Hypersoftness, and Electrostatic Potential

The second-order response of the electron density with respect to changes in electron number, known as the dual descriptor, has been established as a key reactivity indicator for reactions like pericyclic reactions, where reagents accept and donate electrons concurrently. Here we establish that the dual descriptor is also the key reactivity indicator for ambiphilic reagents: reagents that can act either as electrophiles or as nucleophiles, depending on the reaction partner. Specifically, we study dual atoms (which are proposed to act, simultaneously, as an electron acceptor and an electron donor), dual molecules (which react with both electrophiles and nucleophiles, generally at different sites), and dual ion-molecule complexes (which react with both cations and anions). On the basis of our analysis, the dual atom (an Al(I) that has been purported to be dual in the literature) is actually pseudodual in the sense that it does not truly accept electrons from a nucleophiles; rather, it serves as a conduit through which an electrophile can donate electrons to the attached aromatic ring. For understanding dual ion-molecule complexes, it helps to understand that the dual descriptor makes a key contribution to the long-range portion of the quadratic hyperpolarization. In all cases, a complete description of the reactivity of the ambiphilic reagent requires considering both an orbital-based descriptor of electron transfer (the dual descriptor or the local hypersoftness) and the electrostatic potential. The local hypersoftness strongly resembles the dual descriptor.

Critical thoughts on computing atom condensed Fukui functions.

Different procedures to obtain atom condensed Fukui functions are described. It is shown how the resulting values may differ depending on the exact approach to atom condensed Fukui functions. The condensed Fukui function can be computed using either the fragment of molecular response approach or the response of molecular fragment approach. The two approaches are nonequivalent; only the latter approach corresponds in general with a population difference expression. The Mulliken approach does not depend on the approach taken but has some computational drawbacks. The different resulting expressions are tested for a wide set of molecules. In practice one must make seemingly arbitrary choices about how to compute condensed Fukui functions, which suggests questioning the role of these indicators in conceptual density-functional theory.

An elementary derivation of the hard/soft-acid/base principle

The hard/soft-acid/base (HSAB) principle indicates that hard acids prefer binding to hard bases (often forming bonds with substantial ionic character) while soft acids prefer binding to soft bases (often forming bonds with substantial covalent character). Though the HSAB principle is a foundational concept of the modern theory of acids and bases, the theoretical underpinnings of the HSAB principle remain murky. This paper examines the exchange reaction, wherein two molecules, one the product of reacting a hard acid and a soft base and the other the product of reacting a soft acid with a hard base, exchange substituents to form the preferred hard-hard and soft-soft product. A simple derivation shows that this reaction is exothermic, proving the validity of the HSAB principle. The analysis leads to the simple and conceptually appealing conclusion that the HSAB principle is a driven by simple electron transfer effects.

Density-based energy decomposition analysis for intermolecular interactions with variationally determined intermediate state energies.

The first purely density-based energy decomposition analysis (EDA) for intermolecular binding is developed within the density functional theory. The most important feature of this scheme is to variationally determine the frozen density energy, based on a constrained search formalism and implemented with the Wu-Yang algorithm [Q. Wu and W. Yang, J. Chem. Phys. 118, 2498 (2003)]. This variational process dispenses with the Heitler-London antisymmetrization of wave functions used in most previous methods and calculates the electrostatic and Pauli repulsion energies together without any distortion of the frozen density, an important fact that enables a clean separation of these two terms from the relaxation (i.e., polarization and charge transfer) terms. The new EDA also employs the constrained density functional theory approach [Q. Wu and T. Van Voorhis, Phys. Rev. A 72, 24502 (2005)] to separate out charge transfer effects. Because the charge transfer energy is based on the density flow in real space, it has a small basis set dependence. Applications of this decomposition to hydrogen bonding in the water dimer and the formamide dimer show that the frozen density energy dominates the binding in these systems, consistent with the noncovalent nature of the interactions. A more detailed examination reveals how the interplay of electrostatics and the Pauli repulsion determines the distance and angular dependence of these hydrogen bonds.

Atoms in molecules, an axiomatic approach. I. Maximum transferability

Central to chemistry is the concept of transferability: the idea that atoms and functional groups retain certain characteristic properties in a wide variety of environments. Providing a completely satisfactory mathematical basis for the concept of atoms in molecules, however, has proved difficult. The present article pursues an axiomatic basis for the concept of an atom within a molecule, with particular emphasis devoted to the definition of transferability and the atomic description of Hirshfeld.

A New Mean-Field Method Suitable for Strongly Correlated Electrons: Computationally Facile Antisymmetric Products of Nonorthogonal Geminals.

We propose an approach to the electronic structure problem based on noninteracting electron pairs that has similar computational cost to conventional methods based on noninteracting electrons. In stark contrast to other approaches, the wave function is an antisymmetric product of nonorthogonal geminals, but the geminals are structured so the projected Schrödinger equation can be solved very efficiently. We focus on an approach where, in each geminal, only one of the orbitals in a reference Slater determinant is occupied. The resulting method gives good results for atoms and small molecules. It also performs well for a prototypical example of strongly correlated electronic systems, the hydrogen atom chain.

Strategies for computing chemical reactivity indices

Abstract. Two recent articles [(2000) J Am Chem Soc 122: 2010, (2001) J Am Chem Soc 123: 2007] have explored electron-density-based and external-potential-based chemical reactivity indices. In this article, methods are presented for computing these indices from the output of a Kohn–Sham density functional theory calculation.

Density-functional theory calculations with correct long-range potentials

A variational method for forcing the exchange-correlation potential in density-functional theory (DFT) to have the correct asymptotic decay is developed. The resulting exchange-correlation potentials are much improved while the total energies remain essentially the same, compared with conventional density-functional theory calculations. The highest occupied orbital energies from the asymptotically corrected exchange-correlation potentials are found to provide significantly more accurate approximations to the ionization potential (for a neutral molecule) and the electron affinity (for an anion) than those from conventional calculations, although the results are usually inferior to direct methods by computing energy differences. Extending recent results from exchange-only DFT, we show that exact exchange-correlation potential is nonuniform asymptotically. Correcting the asymptotic decay of approximate exchange-correlation potentials towards the exact functional form binds the highest occupied orbitals for atomic and molecular anions, which supports the use of DFT calculations for negatively charged molecular species. With this technique, even hybrid functionals have local exchange-correlation potentials, effectively removing the largest objection to including these functionals in the panoply of Kohn-Sham DFT methods.