Usefulness of the Colle–Salvetti model for the treatment of the nondynamic correlation

Organic Emitters Showing Excited-States Energy Inversion: An Assessment of MC-PDFT and Correlation Energy Functionals Beyond TD-DFT

The lowest-energy singlet (S1) and triplet (T1) excited states of organic conjugated chromophores are known to be accurately calculated by modern wavefunction and Time-Dependent Density Functional Theory (TD-DFT) methods, with the accuracy of the latter heavily relying on the exchange-correlation functional employed. However, there are challenging cases for which this cannot be the case, due to the fact that those excited states are not exclusively formed by single excitations and/or are affected by marked correlation effects, and thus TD-DFT might fall short. We will tackle here a set of molecules belonging to the azaphenalene family, for which research recently documented an inversion of the relative energy of S1 and T1 excited states giving rise to a negative energy difference (ΔEST) between them and, thereby, contrary to most of the systems thus far treated by TD-DFT. Since methods going beyond standard TD-DFT are not extensively applied to excited-state calculations and considering how challenging this case is for the molecules investigated, we will prospectively employ here a set of non-standard methods (Multi-Configurational Pair Density Functional Theory or MC-PDFT) and correlation functionals (i.e., Lie–Clementi and Colle–Salvetti) relying not only on the electronic density but also on some modifications considering the intricate electronic structure of these systems.

Multiconfiguration Pair-Density Functional Theory

Kohn-Sham density functional theory with the available exchange-correlation functionals is less accurate for strongly correlated systems, which require a multiconfigurational description as a zero-order function, than for weakly correlated systems, and available functionals of the spin densities do not accurately predict energies for many strongly correlated systems when one uses multiconfigurational wave functions with spin symmetry. Furthermore, adding a correlation functional to a multiconfigurational reference energy can lead to double counting of electron correlation. Multiconfiguration pair-density functional theory (MC-PDFT) overcomes both obstacles, the second by calculating the quantum mechanical part of the electronic energy entirely by a functional, and the first by using a functional of the total density and the on-top pair density rather than the spin densities. This allows one to calculate the energy of strongly correlated systems efficiently with a pair-density functional and a suitable multiconfigurational reference function. This article reviews MC-PDFT and related background information.

Hartree-Fock implementation using a Laguerre-based wave function for the ground state and correlation energies of two-electron atoms

An implementation of the Hartree-Fock (HF) method using a Laguerre-based wave function is described and used to accurately study the ground state of two-electron atoms in the fixed nucleus approximation, and by comparison with fully correlated (FC) energies, used to determine accurate electron correlation energies. A variational parameter A is included in the wave function and is shown to rapidly increase the convergence of the energy. The one-electron integrals are solved by series solution and an analytical form is found for the two-electron integrals. This methodology is used to produce accurate wave functions, energies and expectation values for the helium isoelectronic sequence, including at low nuclear charge just prior to electron detachment. Additionally, the critical nuclear charge for binding two electrons within the HF approach is calculated and determined to be Z^HF_C=1.031 177 528.This article is part of the theme issue 'Modern theoretical chemistry'.

The role of topology in organic molecules: origin and comparison of the radical character in linear and cyclic oligoacenes and related oligomers

We discuss the nature of electron-correlation effects in carbon nanorings and nanobelts by a combined approach based on FT-DFT and RAS-SF methods.

A hierarchy of static correlation models

It is commonly accepted in the scientific literature that the static correlation energy, Estat, of a system can be defined as the exact correlation energy of its valence electrons in a minimal basis. Unfortunately, the computational cost of calculating the exact correlation energy within a fully optimized minimal basis grows exponentially with system size, making such calculations intractable for all but the smallest systems. However, analogous to single-reference methods, it is possible to systematically approximate both the treatment of electron correlation and flexibility of the minimal basis to reduce computational cost. This yields a hierarchy of methods for calculating Estat, ranging from coupled cluster methods in a minimal atomic basis up to full valence complete active space methods with a minimal molecular orbital basis constructed from a near-complete atomic orbital basis. By examining a variety of dissociating diatomics, along with equilibrium and transition structures for polyatomic systems, we show that standard coupled cluster models with minimal atomic basis sets (e.g., STO-3G) offer a convenient and cost-effective hierarchy of black box estimates for Estat in small- to medium-sized systems near their equilibrium geometries. To properly describe homolytic bond dissociation, it is better to use a more flexible basis set expansion so that each atomic orbital can effectively adapt to its molecular environment.

Combining two-body density correlation functionals with multiconfigurational wave functions using natural orbitals and occupation numbers

We propose a procedure that combines multiconfigurational (MC) wave functions with two-body density correlation functionals by transforming the latter into functionals of the MC natural orbitals and occupation numbers. The method is tested with the spectroscopic constants of a set of 11 diatomics, the diradical-involved automerization barrier of cyclobutadiene, the energy difference between triplet and open-shell singlet states in He and the methylene molecule, and the magnetic coupling constants of several systems, such as NiO, KNiF(3), K(2)NiF(4), La(2)CuO(4), alpha-4-dehydrotoluene, 1,1('),5,5(')-tetramethyl-6,6(')-dioxo-3,3(')-biverdazyl, [Cu(2)Cl(6)](-2), copper(II) acetate monohidrate and H-He-H. The procedure is applied to the Colle-Salvetti [Theor. Chim. Acta 37, 329 (1975); 53, 55 (1979)], functional and to a size-consistent functional depending on the on-top pair density (F1-5-N(eff)). On average, the best results are provided by the transformed F1-5-N(eff) [J. Chem. Phys. 114, 2022 (2001)] functional.

Torsional potential of 4,4'-bipyridine: ab initio analysis of dispersion and vibrational effects

Ab initio calculations using restricted Hartree-Fock, second-order Møller-Plesset perturbation theory (MP2), density-functional theory (DFT), and coupled-cluster methods have been done to obtain the torsional potential-energy profile of the aza-aromatic molecule 4,4'-bipyridine. The torsional potential is evaluated adiabatically by fixing the normalized sum of the dihedral angles through the C-C inter-ring bond at several values along the torsional path and relaxing the remaining degrees of freedom. Previous discrepancies between MP2 and DFT internal rotation barrier heights are removed, and seen to be mostly due to the underestimation of the dispersion energy in the coplanar conformer by MP2 when using relatively small basis sets. The calculations indicate that the barrier height between the twisted global minimum and the 0 degrees conformer is around 1.5-1.8 kcal mol-1 while that corresponding to the 90 degrees one is about 2.0-2.2 kcal mol-1. This same relative energy ordering of the coplanar and perpendicular conformers was experimentally derived from nuclear magnetic resonance (NMR) measurements of 1H dipolar couplings on 4,4'-bipyridine solutions in a nematic liquid crystal, although the barrier heights are much lower than those estimated from NMR experiments in the gas phase. The DFT infrared spectrum and zero-point vibrational energy corrections to the torsional energy profile have also been calculated, the latter having a small influence on the torsional potential-energy profiles.

Valence bond theory for chemical dynamics

This essay provides a perspective on several issues in valence bond theory: the physical significance of semilocal bonding orbitals, the capability of valence bond concepts to explain systems with multireferences character, the use of valence bond theory to provide analytical representations of potential energy surfaces for chemical dynamics by the method of semiempirical valence bond potential energy surfaces (an early example of specific reaction parameters), by multiconfiguration molecular mechanics, by the combined valence bond-molecular mechanics method, and by the use of valence bond states as coupled diabatic states for describing electronically nonadiabatic processes (photochemistry). The essay includes both ab initio and semiempirical approaches.

Analytic gradients for the spin-conserving and spin-flipping equation-of-motion coupled-cluster models with single and double substitutions

Analytic gradient expressions for the spin-conserving and spin-flipping equation-of-motion coupled-cluster models with single and double substitutions are derived using a Lagrangian approach for the restricted and unrestricted Hartree-Fock references, both for the case of all orbitals being active in correlated calculations and for the frozen core and/or virtual orbitals. Details of the implementation within the Q-CHEM electronic structure package are discussed. The capabilities of the new code are demonstrated by application to cyclobutadiene.