Use of a unitary wavefunction in the calculation of static electronic properties

A unitary coupled cluster method is advocated in this paper for the calculation of static properties. Corresponding to the perturbed Hamiltonian H() including the relevant static property, a suitable unitary wavefunction is envisaged. It is shown that a specific nonvariational model of calculating various order static properties utilising this unitary ansatz results in simplifications compared to the previous Coupled Cluster Theories using only hole-particle excitation parameters formulated for this purpose.NCL Communications No. 3533     

On the redistribution of electrons for chemical reaction systems

A regional density-functional theory is formulated and applied to the study of ground-state electron redistributions during the course of a chemical reaction. If for a given increment of the reaction process, accumulation of electrons occurs in a certain region of space, then it is called the dynamic acceptor region, denoted by P. The complement is called the dynamic donor region, denoted by Q. The regional energy itself is determined as a unique functional of the electron density of the total system. The regional transfer potentials are defined in such a way that they add to give the total chemical potential, and their values along the reaction coordinate are found to be different between P and Q. The difference between the regional transfer potentials is shown to provide the driving force for electron transfer from Q to P. A characteristic coordinate for following electron transfer and an associated excitation potential are introduced. The excitation potential is a measure of regional virtual excitation due to regional interactions. The regional transfer potential gives the local character of electron transferability, while the excitation potential gives the global character. The theory encompasses the concepts of regional hardness and softness and sheds light on the HSAB principle.     

ORBKIT: A modular python toolbox for cross‐platform postprocessing of quantum chemical wavefunction data

ORBKIT is a toolbox for postprocessing electronic structure calculations based on a highly modular and portable Python architecture. The program allows computing a multitude of electronic properties of molecular systems on arbitrary spatial grids from the basis set representation of its electronic wavefunction, as well as several grid‐independent properties. The required data can be extracted directly from the standard output of a large number of quantum chemistry programs. ORBKIT can be used as a standalone program to determine standard quantities, for example, the electron density, molecular orbitals, and derivatives thereof. The cornerstone of ORBKIT is its modular structure. The existing basic functions can be arranged in an individual way and can be easily extended by user‐written modules to determine any other derived quantity. ORBKIT offers multiple output formats that can be processed by common visualization tools (VMD, Molden, etc.). Additionally, ORBKIT possesses routines to order molecular orbitals computed at different nuclear configurations according to their electronic character and to interpolate the wavefunction between these configurations. The program is open‐source under GNU‐LGPLv3 license and freely available at https://github.com/orbkit/orbkit/ . This article provides an overview of ORBKIT with particular focus on its capabilities and applicability, and includes several example calculations. © 2016 Wiley Periodicals, Inc.     

Applications of a local grid method for modeling chemical dynamics at a mean-field level

A local grid method proposed earlier is used to model chemical dynamical events in more than one dimension. Two different mean-field routes are applied to model problems representing dynamics of isomerization, H⁺-ion transfer, energy transfer, etc. The methods are seen to work with equal facility for both time-dependent and time-independent potentials. © 1996 John Wiley & Sons, Inc.     

The localized electrons detector as an ab initio representation of molecular structures

The local single particle momentum is proposed as a localized‐electrons detector (LED) that provides a direct three‐dimensional representation of bonding interactions in molecules. It is given exclusively in terms of the electron density and its gradient. We show that the graphical representation of bonding interactions given by LED is consistent with the local curvatures of the electron density as given by the eigenvalues of the Hessian matrix, according to a local symmetry classification of the critical points here introduced. LED consistently complements the topological analysis of the electron density given by the quantum theory of atoms in molecules, by providing a graphical representation of the symmetry of the bonding interactions in molecular systems. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem 110:2418–2425, 2010     

Unpaired electrons at the second‐order reduced density matrix level: Covalent bonding,and coulomb and fermi correlations in closed shell systems

The usual one‐electron populations in atomic orbitals of closed shell systems are split into unpaired and paired at the (spin‐dependent) second‐order reduced density matrix level. The unpaired electron in an orbital is defined as the “simultaneous occurrence of an electron and an electron hole of opposite spins in the same spatial orbital,” which for simplicity is called “electropon.” The electropon population in a given orbital reveals whether and to what degree the Coulomb correlations, and hence, the chemical bonding between this orbital and the remaining orbitals of the system are globally favorable or unfavorable. The interaction of two electropons in two target orbitals reveals the quality (favorable or unfavorable) and the strength of the covalent bonding between these orbitals; this establish a bridge between the notion of “unpaired electrons” and the traditional covalent structure of valence‐bond (VB) theory. Favorable/unfavorable bonding between two orbitals is characterized by the positive/negative (Coulomb) correlation of two electropons of opposite spins, or alternatively, by the negative/positive (Fermi) correlation of two parallel spin electropons. A spin‐free index is defined, and the relationship between the electropon viewpoint for chemical bonding and the well‐known two‐electron Coulomb and Fermi correlations is established. Benchmark calculations are achieved for ethylene, hexatriene, benzene, pyrrole, methylamine, and ammonia molecules on the basis of physically meaningful natural orbitals. The results, obtained in the framework of both orthogonal and nonorthogonal population analysis methods, provide the same conceptual pictures, which are in very good agreement with elementary chemical knowledge and VB theory. © 2013 Wiley Periodicals, Inc.     

Linear scaling local correlation approach for solving the coupled cluster equations of large systems

A linear scaling local correlation approach is proposed for approximately solving the coupled cluster doubles (CCD) equations of large systems in a basis of orthogonal localized molecular orbitals (LMOs). By restricting double excitations from spatially close occupied LMOs into their associated virtual LMOs, the number of significant excitation amplitudes scales only linearly with molecular size in large molecules. Significant amplitudes are obtained to a very good approximation by solving the CCD equations of various subsystems, each of which is made up of a cluster associated with the orbital indices of a subset of significant amplitudes and the local environmental domain of the cluster. The combined effect of these two approximations leads to a linear scaling algorithm for large systems. By using typical thresholds, which are designed to target an energy accuracy, our numerical calculations for a wide range of molecules using the 6-31G or 6-31G* basis set demonstrate that the present local correlation approach recovers more than 98.5% of the conventional CCD correlation energy.     

Electron correlation studies by means of local-scaling transformations and electron-pair density functions

Electron correlation is one long standing problem of computational electronic structure theory. Even more, with the advent of the density functional theory and, in particular, with its Kohn–Sham implementation, the separation of the non-dynamical and dynamical components of the electron correlation has became an unavoidable requirement towards construction of reliable exchange-correlation functionals. In this paper, we address the analysis of the separation of the non-dynamical and dynamical electron correlation effects from two complementary viewpoints, namely, analysis of the correlation energy components and the analysis of the electron-pair density. The former approach will make use of the local-scaling transformations and the latter will be based on the study of intracule and extracule densities.Work supported by grant 9/UPV-00203.215-13527/2001of the Office of Universities and Research of the The Goverment of the Basque Country and, by grant BQU2001-0208 of the Spanish Ministry of Education and Science.Partially supported by grant G-97000741 of CONICIT of Venezuela.     

Superlinear scaling in master-slave quantum chemical calculations using in-core storage of two-electron integrals

We describe the implementation of a parallel, in-core, integral-direct Hartree-Fock and density functional theory code for the efficient calculation of Hartree-Fock wave functions and density functional theory. The algorithm is based on a parallel master-slave algorithm, and the two-electron integrals calculated by a slave are stored in available local memory. To ensure the greatest computational savings, the master node keeps track of all integral batches stored on the different slaves. The code can reuse undifferentiated two-electron integrals both in the wave function optimization and in the evaluation of second-, third-, and fourth-order molecular properties. Superlinear scaling is achieved in a series of test examples, with speedups of up to 55 achieved for calculations run on medium-sized molecules on 16 processors with respect to the time used on a single processor.     

Effects of electron correlations in a surface-molecule model for the adiabatic electrochemical electron transfer reactions with allowance for the electrostatic repulsion of electrons on an effective orbital of metal

Within the framework of a surface-molecule model for the adiabatic electrochemical electron transfer reactions, exact expressions for the adiabatic free energy surfaces are obtained and the diagrams of kinetic modes are constructed with allowance made for the electrostatic repulsion between electrons with the opposite spin projection both on the valence orbital of the reactant and on the effective electron orbital of the metal. It is shown that taking into account the electrostatic repulsion on the effective orbital of the metal and the correlation effects connected with it is very substantial for a number of electrochemical electron-transfer reactions and leads not only to an alteration of the activation free energies but also to qualitatively different forms of adiabatic free energy surfaces in some regions of values of the model’s parameters.     

Dynamic behavior of the IRC in chemical laser systems

The ab initio potential energy surface of the HF elimination reaction of CH₂FOH is surveyed. In order to elucidate the mechanism of the vibrational excitation of the product HF, the vibrational frequency correlation diagram is analyzed along the IRC.     

Exact relations between the electron density and external potential for systems of interacting and noninteracting electrons

The differential virial theorem (DVT) is an explicit relation between the electron density ρ( r ), the external potential, kinetic energy density tensor, and (for interacting electrons) the pair function. The time‐dependent generalization of this relation also involves the paramagnetic current density. We present a detailed unified derivation of all known variants of the DVT starting from a modified equation of motion for the current density. To emphasize the practical significance of the theorem for noninteracting electrons, we cast it in a form best suited for recovering the Kohn–Sham effective potential v_s( r ) from a given electron density. The resulting expression contains only ρ( r ), v_s( r ), kinetic energy density, and a new orbital‐dependent ingredient containing only occupied Kohn–Sham orbitals. Other possible applications of the theorem are also briefly discussed. © 2012 Wiley Periodicals, Inc.     

Organic mixed-valence compounds: a playground for electrons and holes

Mixed-valence (MV) compounds are excellent model systems for the investigation of basic electron-transfer (ET) or charge-transfer (CT) phenomena. These issues are important in complex biophysical processes such as photosynthesis as well as in artificial electronic devices that are based on organic conjugated materials. Organic MV compounds are effective hole-transporting materials in organic light emitting diodes (OLEDs), solar cells, and photochromic windows. However, the importance of organic mixed-valence chemistry should not be seen in terms of the direct applicability of these species but the wealth of knowledge about ET phenomena that has been gained through their study. The great variety of organic redox centers and spacer moieties that may be combined in MV systems as well as the ongoing refinement of ET theories and methods of investigation prompted enormous interest in organic MV compounds in the last decades and show the huge potential of this class of compounds. The goal of this Review is to give an overview of the last decade in organic mixed valence chemistry and to elucidate its impact on modern functional materials chemistry.     

Electron‐pair density decomposition for core–valence separable systems

The electron pair density of a core‐valence separable system can be decomposed into three parts: core‐core, core‐valence, and valence‐valence. The core‐core part has a Hartree‐Fock like structure. The core‐valence part can be written as Γ_cv (1,2) = γ_c (1,1)γ_v (2,2) ? γ_c (1,2)γ_v (2,1) + γ_c (2,2)γ_v (1,1) ? γ_c (2,1)γ_v (1,2), where only the 1‐matrices from the core and valence orbitals contribute. The valence‐valence part is left to be determined from the reduced frozen‐core type wave function, which often contains the essential information on the electron correlation and the chemical bond. We demonstrate the analysis to the ground state of negative ion Li^? and 2¹Σ_u⁺ excited state of the Li₂ molecule. © 2012 Wiley Periodicals, Inc.     

Evaluation of the electron momentum density of crystalline systems from ab initio linear combination of atomic orbitals calculations

Alternative techniques are presented for the evaluation of the electron momentum density (EMD) of crystalline systems from ab initio linear combination of atomic‐orbitals calculations performed in the frame of one‐electron self‐consistent‐field Hamiltonians. Their respective merits and drawbacks are analyzed with reference to two periodic systems with very different electronic features: the fully covalent crystalline silicon and the ionic lithium fluoride. Beyond one‐electron Hamiltonians, a post‐Hartree–Fock correction to the EMD of crystalline materials is also illustrated in the case of lithium fluoride. © 2012 Wiley Periodicals, Inc.