Bond Orders and Their Relationships with Cumulant and Unpaired Electron Densities

Natural orbitals and two-particle correlators as tools for the analysis of effective exchange couplings in solids

Using generalizations of spin-averaged natural orbitals and two-particle charge correlators for solids, we investigate the electronic structure of antiferromagnetic transition-metal oxides with a fully self-consistent, imaginary-time GW method. Our findings disagree with the Goodenough-Kanamori (GK) rules that are commonly used for the qualitative interpretation of such solids. First, we found a strong dependence of the natural orbital occupancies on momenta, contradicting GK assumptions. Second, along the momentum path, the character of natural orbitals changes. In particular, the contributions of oxygen 2s orbitals are important, which has not been considered in the GK rules. To analyze the influence of the electronic correlation on the values of effective exchange coupling constants, we use both natural orbitals and two-particle correlators and show that electronic screening modulates the degree of superexchange by stabilizing the charge-transfer contributions, which greatly affects these coupling constants. Finally, we give a set of predictions and recommendations regarding the use of density functional, Green's function, and wave-function methods for evaluating effective magnetic couplings in molecules and solids.

Cumulants as the variables of density cumulant theory: A path to Hermitian triples

We study the combination of orbital-optimized density cumulant theory and a new parameterization of reduced density matrices in which the variables are the particle-hole cumulant elements. We call this combination OλDCT. We find that this new Ansatz solves problems identified in the previous unitary coupled cluster Ansatz for density cumulant theory: the theory is now free of near-zero denominators between occupied and virtual blocks, can correctly describe the dissociation of H₂, and is rigorously size-extensive. In addition, the new Ansatz has fewer terms than the previous unitary Ansatz, and the optimal orbitals delivered by the exact theory are the natural orbitals. Numerical studies on systems amenable to full configuration interaction show that the amplitudes from the previous ODC-12 method approximate the exact amplitudes predicted by this Ansatz. Studies on equilibrium properties of diatomic molecules show that even with the new Ansatz, it is necessary to include triples to improve the accuracy of the method compared to orbital-optimized linearized coupled cluster doubles. With a simple iterative triples correction, OλDCT outperforms other orbital-optimized methods truncated at comparable levels in the amplitudes, as well as coupled cluster single and doubles with perturbative triples [CCSD(T)]. By adding four more terms to the cumulant parameterization, OλDCT outperforms CCSDT while having the same O(V⁵O³) scaling.

Assessing the orbital-optimized unitary Ansatz for density cumulant theory

The previously proposed Ansatz for density cumulant theory that combines orbital-optimization and a parameterization of the 2-electron reduced density matrix cumulant in terms of unitary coupled cluster amplitudes (OUDCT) is carefully examined. Formally, we elucidate the relationship between OUDCT and orbital-optimized unitary coupled cluster theory and show the existence of near-zero denominators in the stationarity conditions for both the exact and some approximate OUDCT methods. We implement methods of the OUDCT Ansatz restricted to double excitations for numerical study, up to the fifth commutator in the Baker-Campbell-Hausdorff expansion. We find that methods derived from the Ansatz beyond the previously known ODC-12 method tend to be less accurate for equilibrium properties and less reliable when attempting to describe H₂ dissociation. New developments are needed to formulate more accurate density cumulant theory variants.

Topological population analysis and pairing/unpairing electron distribution evolution: Atomic B3+ cluster bending mode, a case study

Local and non-local topological treatment of electronic distributions are applied to a simple out of equilibrium case of an electron-deficient three-atom cluster, B3⁺. The bending movement is described in detail through the onset and disappearance of critical points defining two kinds of molecular structures, characterizing a transition state (TS) and predicting two stable equilibrium geometries. All points in this rich evolution and the structural change in the out of equilibrium conformations has been featured and distinguished by the behavior of the population magnitudes and of the paired and unpaired electron densities within the non-local and local points of view of the topological formalism. The unpaired or electron hole density appears as relevant in both versions, the non-local or integrated one, in which it is sometimes called free-valence and also for its complementary counterpart, the local one, to describe and to quantify the interatomic interactions. The stability of the cluster B₃⁺ is characterized in terms of a topologically defined ring structure and the highest total two- and three-center populations, thus showing the role of the geometry, the covalence, and the complex patterns. Consideration of the electron correlation effects constitutes the basement of the results gathered, thus displaying their influence in the formation and breaking of boron bonding interactions.

Photochemical oxidation of di-n-butyl phthalate in atmospheric hydrometeors by hydroxyl radicals from nitrous acid

The photochemical oxidation of di-n-butyl phthalate (DBP) by ^•OH radicals from nitrous acid (HONO) in atmospheric hydrometeors was explored by two techniques, steady-state irradiation, and laser flash photolysis (LFP). The effects of atmospheric liquid parameters on DBP transformation were systematically evaluated, showing that DBP does not react with HONO directly and ^•OH-initiated reactions are crucial steps for consumption and transformation of DBP. Two reaction channels are operative: ^•OH addition and hydrogen atom abstraction. The overall rate constant for the reaction of DBP with ^•OH is 5.7 × 10⁹ M^-1 s^-1, and its specific rate constant for addition is 3.7 × 10⁹ M^-1 s^-1 determined by using laser flash photolysis technique. Comparing the individual reaction rate constant for aromatic ring addition with the total rate constant, the majority of the ^•OH radicals (about 65%) attack the aromatic ring. The major transformation products were identified by GC-MS, and the trends of their yields derived from both ring addition and H-abstraction with time are discussed. These results provide important insights into the photochemical transformation of DBP in atmospheric hydrometeors and contribute to atmospheric aerosol chemistry.

Do Organometallic CH4-Me(+p) Adducts and X4H(+) (X = P, As) Clusters Undergo Two-Electron Three-Center Interactions? Some Aspects of Discussion

Most of the systems possessing true two-electron three-center interactions are electron deficient compounds like boron hydrids, closo-boranes, and some organic ions such as butonium cations. In this work, we perform a detailed study of the electron distribution for two different types of systems to which likewise interactions has been adjudicated: organometallic CH4-Me(+p) (p = 1, 2) adducts with Me, alkaline and earth alkaline metallic ions of Li, Na, K, Be, Mg, Ca in their stable gaseous phase and X4H(+) (X = P, As) simple clusters. For this purpose, topological analysis of the electron density decomposed into its effectively paired and unpaired contributions has been carried out looking for complex interactions.

Pairing and unpairing electron densities in organic systems: two-electron three center through space and through bonds interactions

Two-electron three-center bonding interactions in organic ions like methonium (CH5(+)), ethonium (C2H7(+)), and protonated alkanes n - C4H11(+) isomers (butonium cations) are described and characterized within the theoretical framework of the topological analysis of the electron density decomposition into its effectively paired and unpaired contributions. These interactions manifest in some of this type of systems as a concentration of unpaired electron cloud around the bond paths, in contrast to the well known paradigmatic boron hydrids in which it is not only concentrated close to the atomic nucleus and the bond paths but out of them and over the region defined by the involved atoms as a whole. This result permits to propose an attempt of classification for these interactions based in such manifestations. In the first type, it is called as interactions through bonds and in the second type as interactions through space type.

An efficient generalized polyelectron population analysis in orbital spaces: the hole-expansion methodology

We present relations leading to an efficient generalized population analysis in orbital spaces of usual delocalized molecular orbital wave functions. Besides the calculation of the diagonal elements of the reduced density matrices of any order, one can also calculate efficiently the probabilities (or, in general, the weights) of various occupation schemes of local electronic structures, by using generalized density operators referring to both electrons and electron holes. Within this population analysis, correlated molecular orbital wave functions can be used, and there are no restrictions to the number of the analyzed electrons and electron holes. It is based on the hole-expansion methodology, according to which a given electronic population is expanded in terms involving only electron holes, which as shown, can be calculated very efficiently; usual difficulties arising from the necessity to handle extremely large local determinantal basis sets are avoided, without introducing approximations. Although an emphasis is given for populations in the basis of orthogonal orbital spaces (providing probabilities), the case of nonorthogonal ones is also considered in order to show the connection of the generalized populations and the traditional weights obtained from valence-bond wave functions. Physically meaningful populations can be obtained by using natural orbitals, such as the natural atomic orbitals (NAOs) (orthogonal orbitals) or the pre-NAO's (nonorthogonal orbitals); numerical applications for pyrrole molecule are presented in the basis of these natural orbitals.

Energy decompositions according to physical space partitioning schemes: Treatments of the density cumulant

This article is a continuation of our previous paper on schemes of energy decompositions of molecular systems in the real space [D. R. Alcoba et al., J. Chem. Phys. 122, 074102 (2005)] now using correlated state functions. We study, according to physical arguments, the appropriate management of the density cumulant arising from the second-order reduced density matrix at correlated level, whose contributions can be assigned to one-center or to two-center terms in the energy partitioning. Our treatments are applied within two physical space partitioning schemes: the Bader partitioning into atomic basins and the fuzzy atom procedure. The results obtained in selected molecules are analyzed and discussed in detail.

Control of delocalization and structural changes by means of an electric field

The strength and, mainly, the direction of a static electric field can be used to control delocalization effects occurring in a non-polar pi-system. The delocalization energy, the weights, and the probabilities of some local electronic structures, the behavior of electron pairs, and the electronic fluctuations are considered and examined in cis-butadiene, used as model system. The effects of the electric field are detected and evaluated in the basis of natural orbital spaces appropriate to investigate the behavior of one- and poly-electron distributions. The consequences of modifying the delocalization effects on structural changes are also investigated. Full geometry optimizations in both Hartree-Fock and MP2 levels show that the changes in bond lengths, guided by the changes of the behavior of the electronic assembly, can be controlled by means of the electric field.

Laplacian Field of the Effectively Unpaired Electron Density: Determination of Many-Body Effects on Electron Distributions

This work carries out the study of the Laplacian field function of the electron density L(r) = -nabla2rho(r) splitted in two contributions rho(r) = rho(p)(r) + rho(u)(r), which correspond to the effectively paired and effectively unpaired electron densities, respectively. The visualization of the concentration and depletion of these fields and their spatial localization show no contribution of the effectively unpaired electrons to the conventional bonding among two centers, but the field -nabla2rho(u)(r) provides an interesting structure. We also study the reliability of the information contained in the partitioning of this electron density field function for describing nonclassical bondings as the three-center two-electron ones.

A study of the partitioning of the first-order reduced density matrix according to the theory of atoms in molecules

This work describes a simple spatial decomposition of the first-order reduced density matrix corresponding to an N-electron system into first-order density matrices, each of them associated to an atomic domain defined in the theory of atoms in molecules. A study of the representability of the density matrices arisen from this decomposition is reported and analyzed. An appropriate treatment of the eigenvectors of the matrices defined over atomic domains or over unions of these domains allows one to describe satisfactorily molecular properties and chemical bondings within a determined molecule and among its fragments. Numerical determinations, performed in selected molecules, confirm the reliability of our proposal.

Determination of Three-Center Bond Indices from Population Analyses: A Fuzzy Atom Treatment

This work proposes the use of the treatment referred to as fuzzy atoms to describe three-center bond indices within studies of electron population analysis. A simple manipulation of our algorithms reported previously for describing multicenter bondings enables us to introduce this methodology in our mathematical framework, providing suitable numerical determinations of three-center bond indices, two-center bond ones, and electron atomic populations. The results, obtained in selected systems, are discussed and compared to those arising from other procedures of population analysis.

Correlated holes and their relationships with reduced density matrices and cumulants

This paper describes a matrix formulation for the correlated hole theory within the framework of the domain-averaged model in many electron systems (atoms, molecules, condensed matter, etc.). General relationships between this quantity and one-particle reduced density matrices for any independent particle or correlated state functions are presented. This formulation turns out to be suitable for computational purposes due to the straightforward introduction of cumulants of two-particle reduced density matrices within the quantum field structure. Numerical calculations in selected simple molecular systems have been performed in order to determine preliminary correlated values for such a quantity.

Energy decompositions according to physical space partitioning schemes

This work describes simple decompositions of the energy of molecular systems according to schemes that partition the three-dimensional space. The components of those decompositions depend on one and two atomic domains thus providing a meaningful chemical information about the nature of different bondings among the atoms which compose the system. Our algorithms can be applied at any level of theory (correlated or uncorrelated wave functions). The results reported here, obtained at the Hartree-Fock level in selected molecules, show a good agreement with the chemical picture of molecules and require a low computational cost in comparison with other previously reported decompositions.

Determination of atomic valence indices from population analyses at correlated level

This work carries out a study of atomic valences within molecular systems based on Mulliken and topological population analyses at correlated level. The use of the unpaired electron densities leads to suitable relationships between valences, free valence indices, and bond indices, which turn out to be quite useful for computational purposes. The results arising from both methods at correlated and uncorrelated levels are compared in a large series of chemical compounds. Several interesting conclusions are drawn out and analyzed in detail.