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

A variance-based optimization for determining ground and excited N-electron wave functions within the doubly occupied configuration interaction scheme

This work describes optimizations of N-electron system wave functions by means of the simulated annealing technique within the doubly occupied configuration interaction framework. Using that technique, we minimize the energy variance of a Hamiltonian, providing determinations of wave functions corresponding to ground or excited states in an identical manner. The procedure that allows us to determine electronic spectra can be performed using treatments of restricted or unrestricted types. The results found in selected systems, described in terms of energy, spin, and wave function, are analyzed, showing the performance of each method. We also compare these results with those arising from more traditional approaches that minimize the energy, in both restricted and unrestricted versions, and with those obtained from the full configuration interaction treatment.

Extraction of local spin-coupled states by second quantized operators

We present a methodology for analyzing chemical bonds embedded in the electronic wave function of molecules, especially in terms of spin correlations or so-called "local spin." In this paper, based on biorthogonal second quantization, the spin correlation functions of molecules are naturally introduced, which enables us to extract local singlet and local triplet elements from the wave function. We also clarify the relationship between these spin correlations and traditional chemical concepts, i.e., resonance structures. Several chemical reactions, including the intramolecular radical cyclization and the formation of preoxetane, are demonstrated to verify the analysis method numerically.

Variational determination of the two-electron reduced density matrix within the doubly occupied configuration interaction framework: Treatments of triplet N-electron systems

This work performs variational determinations of two-electron reduced density matrices corresponding to open-shell N-electron systems within the framework of the doubly occupied configuration interaction treatment, traditionally limited to studies of closed-shell systems. The procedure has allowed us to describe satisfactorily molecular systems in triplet states following two methods. One of them adds hydrogen atoms at an infinite distance of the triplet system studied, constituting a singlet supersystem. Energies and reduced density matrices of the triplet system are obtained by removing the contributions of the added atoms from the singlet supersystem results. The second procedure determines variationally the two-electron reduced density matrices corresponding to the triplet systems by means of adequate couplings of basis-set functions. Both models have been managed by imposing N-representability conditions on the reduced density matrix calculations. Results obtained from these methods for molecular systems in triplet ground states are reported and compared with those provided by benchmark methods.

An analysis of valence electronic structure from a viewpoint of resonance theory: Tautomerization of formamide and diazadiboretidine

The resonance theory is still very useful in understanding the valence electron structure. However, such a viewpoint is not usually obtained by general-purpose quantum chemical calculations, instead requires rather special treatment such as valence bond methods. In this study, we propose a method based on second quantization to analyze the results obtained by general-purpose quantum chemical calculations from the local point of view of electronic structure and analyze diazadiboretidine and the tautomerization of formamide. This method requires only the "PS"-matrix, consisting of the density matrix (P-matrix) and overlap matrix, and can be computed with a comparable load to that of Mulliken population analysis. A key feature of the method is that, unlike other methods proposed so far, it makes direct use of the results of general-purpose quantum chemical calculations.

Block deformation analysis: Density matrix blocks as intramolecular deformation density

Block deformation analysis as deformation density of atomic orbitals is introduced to analyze intramolecular interactions. In this respect, density matrix blocks in terms of natural atomic orbitals are employed to find interacting and noninteracting multicenter subsystem and extract the corresponding deformation density. Eigenanalysis of this deformation density is performed to result eigenvalues and eigenorbitals as displaced charge due to the intramolecular interaction and orbital space responsible for charge reorganization, respectively, that possesses advantages of other methods, simultaneously. It is applied to several small molecules, different types of carbon allotropes including zero-, one-, and two-dimensional nanostructures, and challenging systems such as ortho-hydrogen atoms in planar biphenyl. Results highly correlate with delocalization and Wiberg bond indices and show that eigenvalues of block deformation analysis deserved to be considered as bonding index.

Unrestricted treatment for the direct variational determination of the two-electron reduced density matrix for doubly occupied-configuration-interaction wave functions

This work extends to the unrestricted orbital approach the procedure described in our previous report [Alcoba et al., J. Chem. Phys. 148, 024105 (2018)] for determining variationally the two-electron reduced density matrix arising from doubly occupied-configuration-interaction wave functions by imposing two- and three-index N-representability conditions. An analysis of the numerical results obtained in selected systems, from both restricted and unrestricted treatments, allows one to assess the performance of these methodologies as well as to show the influence of the P, Q, G, T1, and T2 positivity conditions. We highlight the satisfactory results obtained within the unrestricted scheme.

Direct variational determination of the two-electron reduced density matrix for doubly occupied-configuration-interaction wave functions: The influence of three-index N-representability conditions

This work proposes the variational determination of two-electron reduced density matrices corresponding to the ground state of N-electron systems within the doubly occupied-configuration-interaction methodology. The P, Q, and G two-index N-representability conditions have been extended to the T1 and T2 (T2') three-index ones and the resulting optimization problem has been addressed using a standard semidefinite program. We report results obtained from the doubly occupied-configuration-interaction method, from the two-index constraint variational procedure and from the two- and three-index constraint variational treatment. The discussion of these results along with a study of the computational cost demanded shows the usefulness of our proposal.

Spin contamination-free N-electron wave functions in the excitation-based configuration interaction treatment

This work deals with the spin contamination in N-electron wave functions provided by the excitation-based configuration interaction methods. We propose a procedure to ensure a suitable selection of excited N-electron Slater determinants with respect to a given reference determinant, required in these schemes. The procedure guarantees the construction of N-electron wave functions which are eigenfunctions of the spin-squared operator Sˆ(2), avoiding any spin contamination. Our treatment is based on the evaluation of the excitation level of the determinants by means of the expectation value of an excitation operator formulated in terms of spin-free replacement operators. We report numerical determinations of energies and 〈Sˆ(2)〉 expectation values, arising from our proposal as well as from traditional configuration interaction methods, in selected open-shell systems, in order to compare the behavior of these procedures and their computational costs.

Study of Electron Delocalization in 1,2-, 1,3-, and 1,4-Azaborines Based on the Canonical Molecular Orbital Contributions to the Induced Magnetic Field and Polyelectron Population Analysis

The electron delocalization in 1,2-azaborine, 1,3-azaborine, and 1,4-azaborine is studied using canonical molecular orbital contributions to the induced magnetic field (CMO-IMF) method and polyelectron population analysis (PEPA). Contour maps of the out-of-plane component of the induced magnetic field (Bz(ind)) of the π system show that the three azaborines, in contrast with borazine, sustain much of benzene's π-aromatic character. Among them, 1,3-azaborine exhibits the strongest π delocalization, while 1,4-azaborine is the weakest. Contour maps of Bz(ind) for individual π orbitals reveal that the differentiation of the magnetic response among the three isomers originates from the π-HOMO orbitals, whose magnetic response is governed by rotational allowed transitions to unoccupied orbitals. The low symmetry of azaborines enables a paratropic response from HOMO to unoccupied orbitals excitations, with their magnitude depending on the shape of interacting orbitals. 1,3-Azaborine presents negligible paratropic contributions to Bz(ind) from HOMO to unoccupied orbitals transitions, where 1,2- and 1,4-azaborine present substantial paratropic contributions, which lead to reduced diatropic response. Natural bond orbital (NBO) analysis employing PEPA shows that only the 1,3-azaborine contains π-electron fully delocalized resonance structures.