Unrestricted Natural Orbitals Research Articles

Automatic Construction of the Initial Orbitals for Efficient Generalized Valence Bond Calculations of Large Systems.

We propose an efficient general strategy for generating initial orbitals for generalized valence bond (GVB) calculations which makes routine black-box GVB calculations on large systems feasible. Two schemes are proposed, depending on whether the restricted Hartree-Fock (RHF) wave function is stable (scheme I) or not (scheme II). In both schemes, the first step is the construction of active occupied orbitals and active virtual orbitals. In scheme I, active occupied orbitals are composed of the valence orbitals (the inner core orbitals are excluded), and the active virtual orbitals are obtained from the original virtual space by requiring its maximum overlap with the virtual orbital space of the same system at a minimal basis set. In scheme II, active occupied orbitals and active virtual orbitals are obtained from the set of unrestricted natural orbitals (UNOs), which are transformed from two sets of unrestricted HF spatial orbitals. In the next step, the active occupied orbitals and active virtual ones are separately transformed to localized orbitals. Localized occupied and virtual orbital pairs are formed using the Kuhn-Munkres (KM) algorithm and are used as the initial guess for the GVB orbitals. The optimized GVB wave function is obtained using the second-order self-consistent-field algorithm in the GAMESS program. With this procedure, GVB energies have been obtained for the lowest singlet and triplet states of polyacenes (up to decacene with 96 pairs) and the singlet ground state of two di-copper-oxygen-ammonia complexes. We have also calculated the singlet-triplet gaps for some polyacenes and the relative energy between two di-copper-oxygen-ammonia complexes with the block-correlated second-order perturbation theory based on the GVB reference.

Full-valence density matrix renormalisation group calculations on meta-benzyne based on unrestricted natural orbitals. Revisit of seamless continuation from broken-symmetry to symmetry-adapted models for diradicals

ABSTRACTIn this work, we show that the natural orbitals of unrestricted hybrid density functional theory (UHDFT) can be used as the active space orbitals to perform multireference (MR) calculations, for example, the density matrix renormalisation group (DMRG) and Mukherjee-type (Mk) MR coupled-cluster (CC) method. By including a sufficiently large number of these natural orbitals, full-valence (FV) active space can be identified without recourse of the expensive self-consistent procedures for DMRG-SCF. Several useful chemical indices are derived based on the occupation numbers of the natural orbitals for seamless continuation from broken-symmetry (BS) to symmetry-adapted (SA) methods. These procedures are used on 1,3-didehydrobenzene (meta-benzyne) to calculate its singlet (S)-triplet (T) gap. We compare our results to available experiments and computational results obtained by several other groups. We see our procedures as a seamless bridge between single-reference BS methods, such as UHDFT, and the SA MR methods, such as FV DMRG and MkMRCC.

Symmetry Breaking and Broken Ergodicity in Full Configuration Interaction Quantum Monte Carlo.

The initiator full configuration interaction quantum Monte Carlo method (i-FCIQMC) is applied to the binding curve of N2 in Slater-determinant Hilbert spaces formed of both canonical restricted Hartree-Fock (RHF) and symmetry-broken unrestricted Hartree-Fock (UHF) orbitals. By explicit calculation, we demonstrate that the technique yields the same total energy for both types of orbital but that as the bond is stretched, FCI expansions expressed in unrestricted orbitals are substantially more compact than their restricted counterparts and more compact than those expressed in split-localized orbitals. These unrestricted Hilbert spaces, however, become nonergodic toward the dissociation limit, and the total wave function may be thought of as the sum of two weakly coupled, spin-impure, functions whose energies are nonetheless very close to the exact energy. In this limit, it is a challenge for i-FCIQMC to resolve a spin-pure wave function. The use of unrestricted natural orbitals is a promising remedy for this problem, as their expansions are more strongly weighted toward lower excitations of the reference, and they provide stronger coupling to higher excitations than do UHF orbitals.

Semiempirical UNO–CAS and UNO–CI: Method and Applications in Nanoelectronics

Unrestricted Natural Orbital-Complete Active Space Configuration Interaction, abbreviated as UNO-CAS, has been implemented for NDDO-based semiempirical molecular-orbital (MO) theory. A computationally more economic technique, UNO-CIS, in which we use a configuration interaction (CI) calculation with only single excitations (CIS) to calculate excited states, has also been implemented and tested. The class of techniques in which unrestricted natural orbitals (UNOs) are used as the reference for CI calculations is denoted UNO-CI. Semiempirical UNO-CI gives good results for the optical band gaps of organic semiconductors such as polyynes and polyacenes, which are promising materials for nanoelectronics. The results of these semiempirical UNO-CI techniques are generally in better agreement with experiment than those obtained with the corresponding conventional semiempirical CI methods and comparable to or better than those obtained with far more computationally expensive methods such as time-dependent density-functional theory. We also show that symmetry breaking in semiempirical UHF calculations is very useful for predicting the diradical character of organic compounds in the singlet spin state.

The unrestricted natural orbital-restricted active space method: methodology and implementation

The full configuration interaction method in the space of fractionally occupied unrestricted natural orbitals (UNO-CAS method) is extended to excited states as well as to strongly correlated and reactive systems with large active spaces. This is accomplished by␣using restricted active space (RAS) wave functions introduced by Olsen et al. [(1988) J Chem Phys 89: 2185] and using the UNOs without the expensive orbital optimization step. In RAS, the space of active orbitals is subdivided into three groups: a group with essentially doubly occupied orbitals (RAS1), the usual CAS space (RAS2), and a space with weakly occupied active orbitals (RAS3). We select these spaces on the basis of the occupation numbers of the UNOs. All possible electron distributions are allowed in the usual CAS space, but the number of vacancies is limited in RAS1 and the number of electrons is limited in RAS3. We discuss an efficient algorithm for generating a RAS wave function. This is based on the Handy-Knowles determinantal expansion with an addressing scheme adopted for the restricted expansion. Results for both ground and excited states of azulene and free base porphyrin are presented.