Interpolation of property-values between electron numbers is inconsistent with ensemble averaging

Quantum states of physical domains in molecular systems: A three-state model approach

The physical regions (domains or basins) within the molecular structure are open systems that exchange charge between them and, consequently, house a fractional number of electrons (net charge). The natural framework describing the quantum states for these domains is the density matrix (DM) in its grand-canonical version, which corresponds to a convex expansion into a set of basis states of an integer number of electrons. In this report, it is shown that the solution for these quantities is supported by the DM expansion into three states of different numbers of particles: the neutral and two (edge) ionic states. The states and the average number of particles in the domains (fractional occupation population) are determined by the coefficients of the expansion in terms of the fundamental transference magnitudes, revealing the donor/acceptor character of the domains by which the quantum accessible states are discussed.

1,3 Dipolar Cycloaddition of Münchnones: Factors behind the Regioselectivity

The 1,3 dipolar cycloaddition reactions of münchnones and alkenes provide an expedite synthetic way to substituted pyrroles, an exceedingly important structural motif in the pharmaceutical and material science fields of research. The factors governing their regioselectivity rationalization are not well understood. Using several approaches, we investigate a set of 14 reactions (featuring two münchnones, 12 different alkenes, and two alkynes). The Natural Bond Theory and the Non-Covalent Interaction Index analyses of the noncovalent interaction energies fail to predict the experimental major regioisomer. Employing global cDFT descriptors or local ones such as the Fukui function and dual descriptor yields similarly inaccurate predictions. Only the local softness pairing, within Pearson's Hard and Soft Acids and Bases principle, constitutes a reliable predictor for the major reaction product. By taking into account an estimator for the steric effects, the correct regioisomer is predicted. Steric effects play a major role in driving the regioselectivity, as was corroborated by energy decomposition analysis of the transition states.

Electronic properties of amino acids and nucleobases: similarity classes and pairing principles from chemical reactivity indices

We present a new classification scheme for amino acids and nucleobases based on the electronic properties of the individual molecules. Using chemical reactivity indices such as electronegativity, electrophilicity, and chemical hardness, we can identify similarities and differences between each class of amino acids and nucleobases. Notable differences emerge in particular with regard to high, neutral or low electronegativity as well as different combinations of chemical hardness. Our approach allows us to relate these insights to the properties of the side groups in terms of a unique reference scheme. We further show that hydrophobic differences between amino acids are rather negligible in the context of electronic properties. Our classification scheme also rationalizes the occurrence of distinct stable nucleobase pairs and clearly emphasizes certain differences between individual molecules. The stability and abundant occurrence of Watson-Crick nucleobase pairs is further discussed in the context of the minimum electrophilicity principle.

Molecular interactions from the density functional theory for chemical reactivity: Interaction chemical potential, hardness, and reactivity principles

In the first paper of this series, the authors derived an expression for the interaction energy between two reagents in terms of the chemical reactivity indicators that can be derived from density functional perturbation theory. While negative interaction energies can explain reactivity, reactivity is often more simply explained using the “|dμ| big is good” rule or the maximum hardness principle. Expressions for the change in chemical potential (μ) and hardness when two reagents interact are derived. A partial justification for the maximum hardness principle is that the terms that appear in the interaction energy expression often reappear in the expression for the interaction hardness, but with opposite sign.

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

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

From adaptive resolution to molecular dynamics of open systems

ABSTRACT

We provide an overview of the Adaptive Resolution Simulation method (AdResS) based on discussing its basic principles and presenting its current numerical and theoretical developments. Examples of applications to systems of interest to soft matter, chemical physics, and condensed matter illustrate the method's advantages and limitations in its practical use and thus settle the challenge for further future numerical and theoretical developments.

Exploring chemical space with alchemical derivatives: alchemical transformations of H through Ar and their ions as a proof of concept

Alchemical derivatives have been used previously to obtain information about transformations in which the number of electrons is unchanged. Here an approach for combining changes in both the number of electrons and the nuclear charge is presented.

Electronegativity-a perspective

Electronegativity is a very useful concept but it is not a physical observable; it cannot be determined experimentally. Most practicing chemists view it as the electron-attracting power of an atom in a molecule. Various formulations of electronegativity have been proposed on this basis, and predictions made using different formulations generally agree reasonably well with each other and with chemical experience. A quite different approach, loosely linked to density functional theory, is based on a ground-state free atom or molecule, and equates electronegativity to the negative of an electronic chemical potential. A problem that is encountered with this approach is the differentiation of a noncontinuous function. We show that this approach leads to some results that are not chemically valid. A formulation of atomic electronegativity that does prove to be effective is to express it as the average local ionization energy on an outer contour of the atom's electronic density.

How predictive could alchemical derivatives be?

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

Benchmark values of chemical potential and chemical hardness for atoms and atomic ions (including unstable anions) from the energies of isoelectronic series

We present benchmark values for the electronic chemical potential and chemical hardness from reference data for ionization potentials and electron affinities. In cases where the energies needed to compute these quantities are not available, we estimate the ionization potential of the metastable (di)anions by extrapolation along the isoelectronic series, taking care to ensure that the extrapolated data satisfy reasonable intuitive rules to the maximum possible extent. We also propose suitable values for the chemical potential and chemical hardness of zero-electron species. Because the values we report are faithful to the trends in accurate data on atomic energies, we believe that our proposed values for the chemical potential and chemical hardness are ideally suited to conceptual studies of chemical trends across the periodic table. The critical nuclear charge (Z critical) of the isoelectronic series with 2 < N < 96 has also been reported for the first time.

New Fukui, dual and hyper-dual kernels as bond reactivity descriptors

We define three new linear response indices with promising applications for bond reactivity using the mathematical framework of τ-CRT (finite temperature chemical reactivity theory). The τ-Fukui kernel is defined as the ratio between the fluctuations of the average electron density at two different points in the space and the fluctuations in the average electron number and is designed to integrate to the finite-temperature definition of the electronic Fukui function. When this kernel is condensed, it can be interpreted as a site-reactivity descriptor of the boundary region between two atoms. The τ-dual kernel corresponds to the first order response of the Fukui kernel and is designed to integrate to the finite temperature definition of the dual descriptor; it indicates the ambiphilic reactivity of a specific bond and enriches the traditional dual descriptor by allowing one to distinguish between the electron-accepting and electron-donating processes. Finally, the τ-hyper dual kernel is defined as the second-order derivative of the Fukui kernel and is proposed as a measure of the strength of ambiphilic bonding interactions. Although these quantities have never been proposed, our results for the τ-Fukui kernel and for τ-dual kernel can be derived in zero-temperature formulation of the chemical reactivity theory with, among other things, the widely-used parabolic interpolation model.

Going beyond the three-state ensemble model: the electronic chemical potential and Fukui function for the general case

The analytical working equations for the chemical potential and the Fukui function for the case of any number of ground and excited states is presented.