Does the radical GPRI strongly depend on the population scheme? A comparative study to predict radical attack on unsaturated molecules with the radical general-purpose reactivity indicator.

The reactivity of 22 unsaturated molecules undergoing attack by a methyl radical (⋅CH) have been elucidated using the condensed radical general-purpose reactivity indicator (condensed radical GPRI) appropriate for relatively nucleophilic or electrophilic molecules. Using the appropriate radical GPRI equation for electrophilic attack or nucleophilic radical attack, seven different population schemes were used to assign the most reactive atoms in each of the 22 molecules. The results show that the condensed radical GPRI is sensitive to the population scheme chosen, but less sensitive than the radical Fukui function.

Bioaccumulative chemicals are either too hard or too soft: Conceptual density functional theory as a screening tool for emerging pollutants.

Conceptual density functional theory (CDFT) descriptors were computed to predict the environmental fate of approximately 6,000 widely used industrial chemicals. CDFT descriptors aligned with a molecule's possible bioaccumulation mechanism, i.e.

Warning! The negative divergence of the stress-tensor does not always yield the Ehrenfest force.

It has been assumed that the negative divergence of all stress tensors in common use yields the same force. This work finds that this is untrue, and, in fact, can vary wildly. We demonstrate this for the hydrogen atom, the one-particle isotropic harmonic oscillator, and a particle in an infinite spherical well where the exact density, pair-density, and the first order reduced density matrix are known for ground and excited states without any approximation.

An information-theoretic approach to basis-set fitting of electron densities and other non-negative functions.

The numerical ill-conditioning associated with approximating an electron density with a convex sum of Gaussian or Slater-type functions is overcome by using the (extended) Kullback-Leibler divergence to measure the deviation between the target and approximate density. The optimized densities are non-negative and normalized, and they are accurate enough to be used in applications related to molecular similarity, the topology of the electron density, and numerical molecular integration. This robust, efficient, and general approach can be used to fit any non-negative normalized functions (e.

Predicting the Mass Spectra of Environmental Pollutants Using Computational Chemistry: A Case Study and Critical Evaluation.

Organic pollutants can be identified by comparing their electron ionization (EI) mass spectra with those in libraries or obtained from authentic standards. Nevertheless, libraries are incomplete; standards may be unavailable or too costly, or their synthesis may be too time-consuming. This study evaluates the performance of quantum chemical electron ionization mass spectrometry (QCEIMS) vis-à-vis competitive fragmentation modeling (CFM) for suspect screening and unknown identification.

Molecular QTAIM Topology Is Sensitive to Relativistic Corrections.

The topology of the molecular electron density of benzene dithiol gold cluster complex Au -S-C H -S'-Au' changed when relativistic corrections were made and the structure was close to a minimum of the Born-Oppenheimer energy surface. Specifically, new bond paths between hydrogen atoms on the benzene ring and gold atoms appeared, indicating that there is a favorable interaction between these atoms at the relativistic level. This is consistent with the observation that gold becomes a better electron acceptor when relativistic corrections are applied.

The general setting for the zero-flux condition: The lagrangian and zero-flux conditions that give the heisenberg equation of motion.

Generalizing our recent work on relativistic generalizations of the quantum theory of atoms in molecules, we present the general setting under which the principle of stationary action for a region leads to open quantum subsystems. The approach presented here is general and works for any Hamiltonian, and when a reasonable Lagrangian is selected, it often leads to the integral of the Laplacian of the electron density on the region vanishing as a necessary condition for the zero-flux surface. Alternatively, with this method, one can design a Lagrangian that leads to a surface of interest (though this Lagrangian may not be, and indeed probably will not be, "reasonable").

A reference-free stockholder partitioning method based on the force on electrons.

We argue that when one divides a molecular property into atom-in-a-molecule contributions, one should perform the division based on the property density of the quantity being partitioned. This is opposition to the normal approach, where the electron density is given a privileged role in defining the properties of atoms-in-a-molecule. Because partitioning each molecular property based on its own property density is inconvenient, we design a reference-free approach that does not (directly) refer atomic property densities.

Relativistic (SR-ZORA) quantum theory of atoms in molecules properties.

The Quantum Theory of Atoms in Molecules (QTAIM) is used to elucidate the effects of relativity on chemical systems. To do this, molecules are studied using density-functional theory at both the nonrelativistic level and using the scalar relativistic zeroth-order regular approximation. Relativistic effects on the QTAIM properties and topology of the electron density can be significant for chemical systems with heavy atoms.

Using the general-purpose reactivity indicator: challenging examples.

We elucidate the regioselectivity of nucleophilic attack on substituted benzenesulfonates, quinolines, and pyridines using a general-purpose reactivity indicator (GPRI) for electrophiles. We observe that the GPRI is most accurate when the incoming nucleophile resembles a point charge. We further observe that the GPRI often chooses reactive "dead ends" as the most reactive sites as well as sterically hindered reactive sites.

Kinetic and electron-electron energies for convex sums of ground state densities with degeneracies and fractional electron number.

Properties of exact density functionals provide useful constraints for the development of new approximate functionals. This paper focuses on convex sums of ground-level densities. It is observed that the electronic kinetic energy of a convex sum of degenerate ground-level densities is equal to the convex sum of the kinetic energies of the individual degenerate densities.

In pursuit of negative Fukui functions: molecules with very small band gaps.

A justification for the likely presence of negative Fukui functions in molecules with small band gaps is given, and a computational study performed to check whether molecules with small band gaps have negative Fukui functions to a chemically significant extent is reported. While regions with negative Fukui functions were observed, significantly negative values for the atom-condensed Fukui functions were not observed.

The Ehrenfest force topology: a physically intuitive approach for analyzing chemical interactions.

Modified ANO-RCC basis sets are used to determine twelve molecular graphs of the Ehrenfest force for H2, CH4, CH2O, CH3NO, C2H2, C2H4, C3H3NO, N4H4, H2O, (H2O)2, (H2O)4 and (H2O)6. The molecular graphs include all types of topological critical points and a mix of bonding types is chosen to include sigma-, π- and hydrogen-bonding. We then compare a wide range of point properties: charge density, trace of the Hessian, eigenvalues, ellipticity, stiffness, total local energy and the eigenvectors are calculated at the bond critical points (BCPs) and compared for the Ehrenfest, QTAIM and stress tensor schemes.

In pursuit of negative Fukui functions: examples where the highest occupied molecular orbital fails to dominate the chemical reactivity.

In our quest to explore molecules with chemically significant regions where the Fukui function is negative, we explored reactions where the frontier orbital that indicates the sites for electrophilic attack is not the highest occupied molecular orbital. The highest occupied molecular orbital (HOMO) controls the location of the regions where the Fukui function is negative, supporting the postulate that negative values of the Fukui function are associated with orbital relaxation effects and nodal surfaces of the frontier orbitals. Significant negative values for the condensed Fukui function, however, were not observed.

Quantum theory of atoms in molecules: results for the SR-ZORA Hamiltonian.

The quantum theory of atoms in molecules (QTAIM) is generalized to include relativistic effects using the popular scalar-relativistic zeroth-order regular approximation (SR-ZORA). It is usually assumed that the definition of the atom as a volume bounded by a zero-flux surface of the electron density is closely linked to the form of the kinetic energy, so it is somewhat surprising that the atoms corresponding to the relativistic kinetic-energy operator in the SR-ZORA Hamiltonian are also bounded by zero-flux surfaces. The SR-ZORA Hamiltonian should be sufficient for qualitative descriptions of molecular electronic structure across the periodic table, which suggests that QTAIM-based analysis can be useful for molecules and solids containing heavy atoms.

How ambiguous is the local kinetic energy?

The local kinetic energy and the closely related local electronic stress tensor are commonly used to elucidate chemical bonding patterns, especially for covalent bonds. We use three different approaches-transformation properties of the stress tensor, quasiprobability distributions, and the virial theorem from density-functional theory-to clarify the inherent ambiguity in these quantities, discussing the implications for analyses based on the local kinetic energy and stress tensor. An expansive-but not universal-family of local kinetic energy forms that includes the most common choices and is suitable for both chemical-bonding and atoms-in-molecule analysis is derived.

Relationships between the third-order reactivity indicators in chemical density-functional theory.

Relationships between third-order reactivity indicators in the closed system [N, v(r)], open system [mu, v(r)], and density [rho(r)] pictures are derived. Our method of derivation unifies and extends known results. Among the relationships is a link between the third-order response of the energy to changes in the density and the quadratic response of the density to changes in external potential.

Predicting the reactivity of ambidentate nucleophiles and electrophiles using a single, general-purpose, reactivity indicator.

We recently proposed a new reactivity indicator, termed the "general-purpose reactivity indicator", Xi, which describes not only the classical reactivity paradigms, but also describes reactions that are neither frontier-orbital nor electrostatically controlled. This indicator was proposed to be especially useful for reactants with multiple reactive sites, especially if the nature of the reactivity at those sites was different. This suggests that this reactivity indicator is especially appropriate for ambidentate molecules; this paper confirms this hypothesis.

Conceptual Density-Functional Theory for General Chemical Reactions, Including Those That Are Neither Charge- nor Frontier-Orbital-Controlled. 2. Application to Molecules Where Frontier Molecular Orbital Theory Fails.

This paper examines cases where frontier molecular orbital theory is known to fail, specifically electrophilic aromatic substitution reactions on isoquinoline and borazarophenanthrenes. While we are able to explain the experimental regioselectivity preferences for isoquinoline without too much difficulty, describing the regioselectivity of the borazarophenanthrenes is much more challenging. This is attributed to the fact that these molecules lie between the electrostatic (or charge) control and electron-transfer (or frontier-orbital) control paradigms.

Conceptual Density-Functional Theory for General Chemical Reactions, Including Those That Are Neither Charge- nor Frontier-Orbital-Controlled. 1. Theory and Derivation of a General-Purpose Reactivity Indicator.

A new general-purpose reactivity indicator is derived. Unlike existing indicators, this indicator can describe the reactivity of molecules that lie between the electrostatic (or charge) control and electron-transfer (or frontier-orbital) control paradigms. Depending on the parameters in the indicator, it describes electrostatic control (where the electrostatic potential is the appropriate indicator), electron-transfer control (where the Fukui function's potential is the appropriate indicator), and intermediate cases (where linear combinations of the electrostatic potential and the Fukui function's potential are appropriate indicators).

Indices for predicting the quality of leaving groups.

The inherent quality of leaving groups in chemical reactions is related to their ionization potential and electron affinity using a quadratic model for the dependence of the energy on the number of electrons. A good leaving group for nucleophilic substitution/elimination reactions is one where the difference in energy between the system with the "optimum" number of electrons and the anion is small. Similarly, a good leaving group for electrophilic substitution/elimination reactions is one where the difference in energy between the system with the optimum number of electrons and the cation is small.