Search

The computational approach to the Hirshfeld [Theor. Chim. Acta 44, 129 (1977)] atom in a molecule is critically investigated, and several difficulties are highlighted. It is shown that these difficulties are mitigated by an alternative, iterative version, of the Hirshfeld partitioning procedure. The iterative scheme ensures that the Hirshfeld definition represents a mathematically proper information entropy, allows the Hirshfeld approach to be used for charged molecules, eliminates arbitrariness in the choice of the promolecule, and increases the magnitudes of the charges. The resulting “Hirshfeld-I charges” correlate well with electrostatic potential derived atomic charges. [ABSTRACT FROM AUTHOR]

The relaxation of the Pauli principle associated with density scaling is examined. Scaling the density has been investigated in the development of density functional computational methods with higher accuracy. Scaling the density by ρ(r)ζ=ρ(r)/ζ reduces the number of electrons to M=N/ζ when ζ>1. The minimum kinetic energy of the scaled density, Ts[ρ/ζ], can be scaled back to the N-electron system by multiplying the M-electron Kohn-Sham-type occupation numbers by ζ to produce Tζ[ρ]. This relaxes the Pauli principle when the orbital occupation numbers are greater than 1 in the N-electron system. The effects of antisymmetry on solutions to the Kohn-Sham equations are examined for Ne and the Be isoelectronic series. The changes in Tζ[ρ] and the exchange energy Exζ[ρ] when ζ is varied show that these two quantities are inextricably linked. [ABSTRACT FROM AUTHOR]

The Woodward-Hoffmann rules for pericyclic reactions, a fundamental set of reactivity rules in organic chemistry, are formulated in the language of conceptual density functional theory (DFT). DFT provides an elegant framework to introduce chemical concepts and principles in a quantitative manner, partly because it is formulated without explicit reference to a wave function, on whose symmetry properties the Woodward-Hoffmann [J. Am. Chem. Soc. 87, 395 (1965)] rules are based. We have studied the initial chemical hardness response using a model reaction profile for two prototypical pericyclic reactions, the Diels-Alder cycloaddition of 1,3-butadiene to ethylene and the addition of ethylene to ethylene, both in the singlet ground state and in the first triplet excited state. For the reaction that is thermally allowed but photochemically forbidden, the initial hardness response is positive along the singlet reaction profile. (By contrast, for the triplet reaction profile, a negative hardness response is observed.) For the photochemically allowed, thermally forbidden reaction, the behavior of the chemical hardness along the initial stages of the singlet and triplet reaction profiles is reversed. This constitutes a first step in showing that chemical concepts from DFT can be invoked to explain results that would otherwise require invoking the phase of the wave function. [ABSTRACT FROM AUTHOR]

We provide a rigorous proof that the Hohenberg-Kohn theorem holds for spin densities by extending Lieb’s Legendre-transform formulation to spin densities. The resulting spin-density-functional theory resolves several troublesome issues. Most importantly, the present paper provides an explicit construction for the spin potentials at any point along the adiabatic connection curve, thus providing a formal basis for the use of exchange-correlation functionals of the spin density in the Kohn-Sham density-functional theory (DFT). The practical implications of this result for unrestricted Kohn-Sham DFT calculations is considered, and the existence of holes below the Fermi level is discussed. We argue that an orbital’s energy tends to increase as its occupation number increases, which provides the basis for a computational algorithm for determining the occupation numbers in Kohn-Sham DFT and helps explain the origin of Hund’s rules and holes below the Fermi level. [ABSTRACT FROM AUTHOR]

A comprehensive analysis is presented for the acid-base double-exchange reaction as well as the associated acid-displacement and base-displacement “half-reactions” with the goal of elucidating the meaning of the hard/soft acid/base (HSAB) principle and the conditions for its validity. When electron-transfer effects are important and other effects are negligible, the HSAB principle is driven by the surpassing stability of the soft acid/soft base product. When electrostatic effects dominate the reactivity, the HSAB principle is driven by the surpassing stability of the hard acid/hard base product. Because electron-transfer effects favor soft/soft interactions, while electrostatic effects favor hard/hard interactions, acid-base exchange reactions may be used to determine whether a reagent’s reactivity is dominated by electron-transfer or by electrostatic effects. Because electron-transfer and electrostatic considerations separately favor the HSAB principle whenever the electronic chemical potentials of the acids and bases involved in the reaction are similar, our analysis provides strong support for the HSAB principle. The electronic chemical potential measures the intrinsic strength of acids and bases. [ABSTRACT FROM AUTHOR]

Given a general, N-particle Hamiltonian operator, analogs of the Hohenberg-Kohn theorem are derived for functions that are more general than the particle density, including density matrices and the diagonal elements thereof. The generalization of Lieb’s Legendre transform ansatz to the generalized Hohenberg-Kohn functional not only solves the υ-representability problem for these entities, but, more importantly, also solves the N-representability problem. Restricting the range of operators explored by the Legendre transform leads to a lower bound on the true functional. If all the operators of interest are incorporated in the restricted maximization, however, the variational principle dictates that exact results are obtained for the systems of interest. This might have important implications for practical work not only for density matrices but also for density functionals. A follow-up paper will present a useful alternative approach to the v- and N-representability problems based on the constrained search formalism. [ABSTRACT FROM AUTHOR]

Classical dynamics can be described with Newton’s equation of motion or, totally equivalently, using the Hamilton-Jacobi equation. Here, the possibility of using the Hamilton-Jacobi equation to describe chemical reaction dynamics is explored. This requires an efficient computational approach for constructing the physically and chemically relevant solutions to the Hamilton-Jacobi equation; here we solve Hamilton-Jacobi equations on a Cartesian grid using Sethian’s fast marching method [J. A. Sethian, Proc. Natl. Acad. Sci. USA 93, 1591 (1996)]. Using this method, we can—starting from an arbitrary initial conformation—find reaction paths that minimize the action or the time. The method is demonstrated by computing the mechanism for two different systems: a model system with four different stationary configurations and the H+H2→H2+H reaction. Least-time paths (termed brachistochrones in classical mechanics) seem to be a suitable chioce for the reaction coordinate, allowing one to determine the key intermediates and final product of a chemical reaction. For conservative systems the Hamilton-Jacobi equation does not depend on the time, so this approach may be useful for simulating systems where important motions occur on a variety of different time scales. © 2004 American Institute of Physics. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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 kernel is sufficient for the non-negativity of the Fukui function. Errors in the partitioning of molecules into atoms and inadequate treatment of correlation are pinpointed as the most likely causes of negative condensed Fukui functions. We conclude by noting that the condensed Fukui functions are, in some respects, more appropriate indicators of a molecular site’s reactivity than the Fukui function itself. © 2002 American Institute of Physics. [ABSTRACT FROM AUTHOR]

Examines the axiomatic basis for the concept of an atom within a molecule. Difficulty in providing a satisfactory mathematical basis; Definition of transferability; Description of Hirshfeld atom in a molecule.

Discusses the alternative definition of exchange-correlation charge in density functional theory (DFT). Numerical evaluation of the exchange-correlation charge by using the Zhao-Morrison-Parr method; Radial distribution of exchange-correlation charge and atomic density for atom neon; Comparison of accurate exchange-correlation charge for atom argon.

Approximate generalized overlap amplitudes (GOAs), also called Feynman–Dyson amplitudes, have been calculated from a full configuration interaction (CI) wave function for the ground state of beryllium using the extended Koopmans’ theorem (EKT). The GOAs were also calculated from the ground state CI wave function for Be and the ground state and excited state CI wave functions for Be+. The EKT GOAs are nearly identical to the corresponding CI GOAs for the lower 2S states of Be+ and for the Be+ 1s2s2 state which has a large GOA occupation number. There are many CI GOAs for which there is not a corresponding EKT GOA. This may be due in part to the limited size of the basis set and to the inability to include natural spin orbitals with small occupation numbers in the EKT calculations. © 1995 American Institute of Physics. [ABSTRACT FROM AUTHOR]

The atoms in molecules associated with the Hirshfeld partitioning minimize the generalized Hellinger-Bhattacharya distance to the reference pro-atom densities. Moreover, the reference proatoms can be chosen by minimizing the distance between the pro-molecule density and the true molecular density. This provides an alternative to both the heuristic "stockholder" and the mathematical information-theoretic interpretations of the Hirshfeld partitioning. These results extend to any member of the family of f-divergences. [ABSTRACT FROM AUTHOR]

We show that the hard/soft acid/base principle holds when electron-transfer effects are dominant and the weaker acid and stronger base are harder than the other acidic and basic reagents. In this case the preference of strong acids for strong bases and weak acids for weak bases reinforces the preference of hard acids for hard bases and soft acids for soft bases. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

We provide new arguments supporting the validity of the maximum hardness and the minimum electrophilicity principles, considering the overall change of these descriptors in a charge-transfer reaction. We analyze two cases: (a) how the reactivity is affected when we perturb one reagent, keeping the other constant; (b) how the hardness and electrophilicity change when we treat the interaction between the reagents as a perturbation. [ABSTRACT FROM AUTHOR]

In methods like geminal-based approaches or coupled cluster that are solved using the projected Schrödinger equation, direct computation of the 2-electron reduced density matrix (2-RDM) is impractical and one falls back to a 2-RDM based on response theory. However, the 2-RDMs from response theory are not N-representable. That is, the response 2-RDM does not correspond to an actual physical N-electron wave function. We present a new algorithm for making these non-N-representable 2-RDMs approximately N-representable, i.e., it has the right symmetry and normalization and it fulfills the P-, Q-, and G-conditions. Next to an algorithm which can be applied to any 2-RDM, we have also developed a 2-RDM optimization procedure specifically for seniority-zero 2-RDMs. We aim to find the 2-RDM with the right properties which is the closest (in the sense of the Frobenius norm) to the non-N-representable 2-RDM by minimizing the square norm of the difference between this initial response 2-RDM and the targeted 2-RDM under the constraint that the trace is normalized and the 2-RDM, Q-matrix, and G-matrix are positive semidefinite, i.e., their eigenvalues are non-negative. Our method is suitable for fixing non-N-representable 2-RDMs which are close to being N-representable. Through the N-representability optimization algorithm we add a small correction to the initial 2-RDM such that it fulfills the most important N-representability conditions. [ABSTRACT FROM AUTHOR]

A recent paper [P. W. Ayers, J. Chem. Phys122, 141102 (2005)] considered the hard/soft acid/base exchange reaction, showing that the products associated with the hard/soft acid/base rule (in which the hard acid and hard base are bound, as are the soft acid and soft base) have lower energy than the alternative (in which the hard acid and soft base would have been bound and similarly the soft acid and hard base). Here we show that the maximum hardness principle also predicts this result. Unlike the previous derivation, we do not need to make any assumptions about the relative strength of the acids and bases. [ABSTRACT FROM AUTHOR]

We have implemented the sweep algorithm for the variational optimization of SU(2) x U(1) (spin and particle number) invariant matrix product states (MPS) for general spin and particle number invariant fermionic Hamiltonians. This class includes non-relativistic quantum chemical systems within the Born-Oppenheimer approximation. High-accuracy ab initio finite field results of the longitudinal static polarizabilities and second hyperpolarizabilities of one-dimensional hydrogen chains are presented. This allows to assess the performance of other quantum chemical methods. For small basis sets, MPS calculations in the saturation regime of the optical response properties can be performed. These results are extrapolated to the thermodynamic limit. [ABSTRACT FROM AUTHOR]

Two recent papers [S. Liu, P. W. Ayers, and R. G. Parr, J. Chem. Phys. 111, 6927 (1999); A. Go¨rling, Phys. Rev. Lett. 83, 5459 (1999)] have stated that ∫∇[sup 2]ν[sub xc](r)dr=4π, where ν[sub xc](r) is the exchange-correlation potential of density functional theory. Here, we derive this sum rule and related rules such as ∫∇[sup 2]ν[sub x](r)dr=4π and ∫∇[sup 2]ν[sub c](r)dr=0, where ν[sub x](r) and ν[sub c](r) are the exchange and correlation components of ν[sub xc](r). Using similar methods, we derive the sum rule for the “screening” portion of the exchange-correlation potential and also “generalized” sum rules for ν[sub c](r) and the “response” portion of the exchange-correlation potential, v[sub xc][sup response](r). From the sum rule for v[sub xc][sup response](r), we deduce the asymptotic decay of the density response of the hole-correlation function. We conclude by discussing the probable utility of these results for the development of new exchange-correlation functionals. © 2001 American Institute of Physics. [ABSTRACT FROM AUTHOR]

The energy of a many-particle system is not convex with respect to particle number for r-k interparticle repulsion potentials if k > log34 ≈ 1.262. With such potentials, some finite electronic systems have ionization potentials that are less than the electron affinity: they have negative band gap (chemical hardness). Although the energy may be a convex function of the number of electrons (for which k = 1), it suggests that finding an analytic proof of convexity will be very difficult. The bound on k is postulated to be tight. An apparent signature of non-convex behavior is that the Dyson orbital corresponding to the lowest-energy mode of electron attachment has a vanishingly small amplitude., (© 2024 Author(s). Published under an exclusive license by AIP Publishing.)

Condensed-to-atom Fukui functions which reflect the atomic reactivity like the tendency susceptible to either nucleophilic or electrophilic attack demonstrate the bonding trend of an atom in a molecule. Accordingly, Fukui functions based concepts, that is, bonding reactivity descriptors which reveal the bonding properties of molecules in the reaction were put forward and then applied to pericyclic and cluster reactions to confirm their effectiveness and reliability. In terms of the results from the bonding descriptors, a covalent bond can readily be predicted between two atoms with large Fukui functions (i.e., one governs nucleophilic attack while the other one governs electrophilic attack, or both of them govern radical attacks) for pericyclic reactions. For Si n O m clusters' reactions, the clusters with a low O atom ratio readily form a bond between two Si atoms with big values of their Fukui functions in which they respectively govern nucleophilic and electrophilic attacks or both govern radical attacks. Also, our results from bonding descriptors show that Si-Si bonds can be formed via the radical mechanism between two Si atoms, and formations of Si-O and O-O bonds are possible when the O content is high. These results conform with experimental findings and can help experimentalists design appropriate clusters to synthesize Si nanowires with high yields. The approach established in this work could be generalized and applied to study reactivity properties for other systems.

Using results from atomic spectroscopy, we show that there are two types of flat-planes conditions. The first type of flat-planes condition occurs when the energy as a function of the number of electrons of each spin, Nα and Nβ, has a derivative discontinuity on a line segment where the number of electrons, Nα + Nβ, is an integer. The second type of flat-planes condition occurs when the energy has a derivative discontinuity on a line segment where the spin polarization, Nα - Nβ, is an integer, but does not have a discontinuity associated with an integer number of electrons. Type 2 flat planes are rare-we observed just 15 type 2 flat-planes conditions out of the 4884 cases we tested-but their mere existence has implications for the design of exchange-correlation energy density functionals. To facilitate the development of functionals that have the correct behavior with respect to both fractional number of electrons and fractional spin polarization, we present a dataset for the chromium atom and its ions that can be used to test new functionals.

Variational second order density matrix theory under "two-positivity" constraints tends to dissociate molecules into unphysical fractionally charged products with too low energies. We aim to construct a qualitatively correct potential energy surface for F(3)(-) by applying subspace energy constraints on mono- and diatomic subspaces of the molecular basis space. Monoatomic subspace constraints do not guarantee correct dissociation: the constraints are thus geometry dependent. Furthermore, the number of subspace constraints needed for correct dissociation does not grow linearly with the number of atoms. The subspace constraints do impose correct chemical properties in the dissociation limit and size-consistency, but the structure of the resulting second order density matrix method does not exactly correspond to a system of noninteracting units.

A variational optimization of the second-order density matrix under the P-, Q-, and G-conditions was carried out for a set of diatomic 14-electron molecules, including N(2), O(2) (2+), NO(+), CO, and CN(-). The dissociation of these molecules is studied by analyzing several chemical properties (dipole moments, population analysis, and bond indices) up to the dissociation limit (10 and 20 A). Serious chemical flaws are observed for the heteronuclear diatomics in the dissociation limit. A careful examination of the chemical properties reveals that the origin of the dissociation problem lies in the flawed description of fractionally occupied species under the P-, Q-, and G-conditions. A novel constraint is introduced that imposes the correct dissociation and enforces size consistency. The effect of this constraint is illustrated with calculations on NO(+), CO, CN(-), N(2), and O(2)(2+).

In analogy to the electrophilicity, we define potentialphilicity indicators that represent energetically favorable ways to change the external potential of a molecule at fixed electron number. Similarly, we define a potentialphobicity to represent the least favorable way to change the external potential of a molecule. The resulting indicators should be useful for describing how molecular geometries change and predicting favorable and unfavorable ways for a reagent to approach a molecule. The linear response function enters plays a very important role in this approach, analogous to the role of the hardness for the electrophilicity or the hardness kernel for the Fukui function. The mathematical properties of the response function and its implications for these reactivity indicators are discussed in depth.