Next Article in Journal / Special Issue
First Principles Molecular Dynamics Study of Catalysis for Polyolefins: the Ziegler-Natta Heterogeneous System.
Previous Article in Journal / Special Issue
Chemical Reactivity Dynamics and Quantum Chaos in Highly Excited Hydrogen Atoms in an External Field: A Quantum Potential Approach
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Can Unrestricted Density-Functional Theory Describe Open Shell Singlet Biradicals?

Department of Theoretical Chemistry, G¨oteborg University, Reutersgatan 2, S-413 20 G¨oteborg, Sweden
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2002, 3(4), 360-394; https://0-doi-org.brum.beds.ac.uk/10.3390/i3040360
Submission received: 29 August 2001 / Accepted: 17 December 2001 / Published: 25 April 2002
(This article belongs to the Special Issue Application of Density Functional Theory)

Abstract

:
Unrestricted density functional theory (UDFT) can be used for the description of open-shell singlet (OSS) biradicals provided a number of precautions are considered. Biradicals that require a two-determinantal wave function (e.g. OSS state of carbenes) cannot be described by UDFT for principal reasons. However, if the overlap between the open-shell orbitals is small (the single electrons are located at different atomic centers) errors become small and, then, the principal failure of UDFT in these cases is not apparent and may even be disguised by the fact that UDFT has the advantage of describing spin polarization better than any restricted open shell DFT method. In the case of OSS biradicals with two- or multiconfigurational character (but a one-determinantal form of the leading configuration), reasonable results can be obtained by broken-symmetry (BS)-UDFT, however in each case this has to be checked. In no case is it reasonable to lower the symmetry of a molecule to get a suitable UDFT description. Hybrid functionals such as B3LYP perform better than pure DFT functionals in BS-UDFT calculations because the former reduce the self-interaction error of DFT exchange functionals, which mimics unspecified static electron correlation effects, so that the inclusion of specific static electron correlation effects via the form of the wavefunction becomes more effective.

1 Introduction

Biradicals such as p-didehydrobenzene [1] or α,3-didehydrotoluene [2, 3] have attracted considerable attention in the last ten years because they are the intermediates of the Bergman reaction [4] of enediynes and determine the biological activity of the latter compounds, in particular their potential as antitumor drugs. Their ground state is a singlet rather than a triplet [1, 2, 3, 6, 7], which seems to tune their reactivity in such a way that DNA is attacked by these biradicals in a regioselective fashion [8]. Since singlet open-shell biradicals are difficult to detect and to investigate by experimental means [1, 2, 3, 6, 7, 8, 9], most of their properties have been determined by quantum chemical calculations, which have become an indispensable tool in this connection [1, 3, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23]. Beside high-quality calculations carried out with wave function theory (WFT) using Coupled Cluster methods such as CCSD(T) [24] or BD(T) [25], work for the larger enediyne systems is increasingly carried out with density functional theory (DFT) [26, 27]. DFT has become an attractive tool for the investigation of enediyne biradicals because DFT and in particular unrestricted DFT (UDFT) lead to reasonable descriptions of biradicals [1, 11, 12, 13, 17, 18, 19, 20, 21, 22, 23] and in this way provides a convenient way of circumventing the use of much more expensive ab initio methods in these cases.
Actually, singlet open-shell biradicals can possess multideterminant (multireference) character, which implies that single determinant methods such as Hartree-Fock (HF) theory or correlation-corrected methods based on HF will fail to describe these molecules in an appropriate way. Reliable descriptions can only be obtained if a multireference approach such as MRCI or MRCC [28] is used or alternatively a single-determinant method covering a high amount of dynamic electron correlation is applied.  Methods such as CCSD(T) or BD(T) will partially compensate for the deficiencies of a single-determinant starting wave function and in this way provide a reasonable account of the properties of biradicals [10, 12, 16, 18, 29].
In view of these considerations it is astonishing that a single-determinant method such as Kohn-Sham DFT [26] has provided reasonable descriptions for some of the biradicals investigated. Cremer, Kraka, and co-workers obtained reasonable energies, geometries, and one-electron properties for m- and p-didehydrobenzene biradicals using UDFT [1, 3, 11, 12, 13, 29]. This observation is in line with the experience made by other authors when using UDFT for the description of open shell singlet (OSS) systems [18, 19, 20, 21, 22, 30, 31]. Gräfenstein and co-workers [11] analyzed the performance of UDFT in the case of p-didehydrobenzene in detail and showed that UDFT actually covers some static electron correlation in spite of its single-determinant nature. The success of DFT in describing singlet biradicals has led other authors to extend the application of this method also to systems where DFT (either unrestricted or restricted, UDFT or RDFT) should fail for principal reasons. Such attempts have have been the basis for confusing statements about the applicability of DFT. In this connection, those DFT investigations are particularly problematic that seem to lead to correct results. They are used to declare UDFT the method of choice for all chemical problems involving biradicals.
In this work, we will critically discuss the question to which extent DFT and in particular UDFT can be used for a reliable description of biradicals with multireference character. For this purpose, we will consider simple biradicals such as methylene in its OSS biradical state (1B1, 1), α, 3-didehydrotoluene (1A″, 2) or 1, 4-didehydrobenzene (1Ag, 3, Scheme 1). These electronic systems are typical representatives for the large classes of carbenes (silylenes, etc.) and didehydroannulenes, which are currently under intense investigations either to synthesize carbenes with increased life-times and special properties for synthesis or to design new antitumor drugs. The size of the biradicals experimentally interesting (e.g. diarylcarbenes or the p-didehydrobenzenes generated from naturally occurring enediynes such as dynemicin, calicheamicin, etc) exclude quantum chemical descriptions with high level ab initio methods and, therefore, DFT is amply used in these cases mostly without the chance of verifying DFT results by other more accurate methods. We have chosen molecules 1, 2, and 3 as suitable examples for the experimentally relevant biradicals because DFT results for 1, 2, and 3 can be compared with available experimental data and the results of high level ab initio calculations. In this work, we seek some basic answers on the performance of DFT in the case of biradicals with multireference character rather than presenting DFT descriptions for a large number of interesting systems.
According to a general scheme of classifying multireference systems suggested by Gräfenstein and Cremer [32] one can distinguish between type 0, type I, type II, and type III systems (Figure 1). A closed shell molecule without any multireference character such as benzene (type 0 systems) can be represented by a single determinant Φ0 (in Figure 1, the 0, 1, etc. denote Φ0, Φ1, etc.) with weight factor w = 1 if just the Hartree potential for the electron interactions is used within Kohn-Sham DFT. By switching on the real electron-electron interaction with the help of the parameter λ (λ increasing from 0 to 1, adiabatic connection scheme [33], Figure 1), configuration state functions Ψ1, Ψ2, etc. corresponding to excited states are mixed into the wave function of the type 0 system thus describing dynamic electron correlation. For most type 0 systems it is sufficient to express excited configuration by just one determinant (Ψ1 = Φ1, Ψ2 = Φ2, etc.), in which an occupied orbital is replaced by an unoccupied orbital with additional nodal planes so that left-right, in-out, angular, and other correlation effects can be described. Since these higher configuration state functions possess weight factors of small but finite magnitude, the weight factor of Ψ0 = Φ0 becomes smaller than 1, but Ψ0 still dominates the total wave function. For λ = 0 → 1, the energy decreases in the way that electron correlation is covered by an increasingly admixture of excited configurations Ψ1 = Φ1, Ψ2 = Φ2, etc. to the ground state wave function Ψ0 = Φ0 (Figure 1). The difference E(λ = 1) − E(λ = 0) is the correlation energy EC, i.e. the energy change due to electron correlation effects (Figure 1).
Type II systems are multireference systems with strong static (non-dynamic) electron correlation [32]. This is reflected by the fact that for λ = λx an (or several) excited configuration state function(s) start(s) to make (a) large contribution(s) to the true wave function, which is reflected by a large increase in the corresponding weight factor(s) (Figure 1). The function Ψ0 does no longer dominate the wave function. The strong mixing in of a higher configuration is reflected by a strong decrease of the energy at λx. The difference EC covers now both the dynamic and the static electron correlation energy.
Ozone is a typical example for a type II system, because the description of its ground state wave functions requires several configuration state functions Ψi. RDFT as a single determinant method is no longer able to describe such a molecule correctly, which is reflected by the fact that the RDFT energy is considerably above the true energy as indicated in Figure 1. The UDFT energy, however, provides often a reasonable approximation for the true energy (Figure 1) because it can cover some multi-configurational effects. Since there is no guarantee that this is really the case for any type II system, it is the goal of this work to clarify for which type II systems UDFT represents a reasonable approximation.
There is a continuous transition from type 0 to type II systems, which implies that one can encounter molecules with weak, medium or strong static electron correlation effects. Type I systems do not take an intermediate position because their wave function is of different nature than that of either type 0 or type II systems [32]. The ground state configuration Ψ0 still dominates the true wave function for λ = 0 → 1 as in the case of type 0 systems, however Ψ0 is already for λ = 0 a two-determinantal wave function where the two determinants Φ0 and Φ′0 contribute with equal weight factors w = w′. Typical representatives of type I systems are 1 in its 1B1 OSS or 2 in its 1A″ OSS state. Hence, if UDFT is able to describe a type II system such as 3 in a reasonable way, this does not automatically imply that a reasonable UDFT description is also obtained for type I systems.
The necessity of a two-determinantal wave function in the case of type I systems does not imply any long-range (static) electron correlation; on the contrary type I systems such as the 1B1 OSS state of methylene are molecules with just dynamic electron correlation, i.e. excited configurations Ψ1 = Φ1, Ψ2 = Φ2, etc. contribute to the true wave function just weakly similar as in the case of type 0 systems. Nevertheless, RDFT fails in these cases because as a single-determinantal method it cannot describe a single-configurational but two-determinantal wave function. Also, UDFT does not present a useful alternative for type I systems as is indicated in Figure 1. The fact that UDFT can describe a type II system such as 3 in a reasonable way does not automatically imply that a reasonable UDFT description is obtained for type I systems.
Similarly as in the case of type 0 systems, type I systems can adopt multiconfigurational character and one has to consider that Ψ0 and excited configuration state functions Ψ1, Ψ2, etc. substantially contribute to the true wave function where each configuration state function is built up by two or more determinants with equal weight. Electron systems that require such a multi-configurational multi-determinantal wave function are denoted as type III systems (Figure 1). Biradical 2 is a representative of a type III system and as such the most difficult biradical to be described. Again, there is a continuous transition from type I to type III systems. In view of the particular nature of the type I and type III wave functions, we have to clarify whether DFT can describe type I and type III systems at all.
It is noteworthy that for the classification of atoms and molecules into type 0, I, II, ands III systems elements and techniques of WFT are used and connected with the DFT description. In this connection, however it is important to realize that electron correlation is differently defined in WFT and DFT. In the former case, it is defined as the difference between the exact energy of the interacting many-electron system and the energy of a suitable reference system, which is normally given by the HF description. Both WFT and Kohn-Sham DFT use the energy of a single-determinantal wave function as reference. This reference function is optimized with respect to the total energy for HF while in KS theory the reference function is optimized to provide the best ρ(r), which means that the reference energy is more positive for KS than for HF; hence the exact KS correlation energy is larger (by its absolute value) than the exact correlation energy in WFT. Also it has to be considered that in DFT the kinetic energy is just calculated for the non-interacting electrons. Hence the total correlation energy covers a kinetic energy correction due to the correlated movement of the electrons.

2 Ways of Describing Biradicals by Unrestricted Hartree-Fock Theory

Since the procedures used to describe biradicals by UDFT originate in HF theory [34], we will briefly review how UHF is used in this connection. The initial guess for an OSS state is normally chosen to be of closed-shell nature. The spin symmetry must be broken explicitly if the singlet biradical is to be described by UHF. This is done by manipulating orbitals ϕR = R and ϕS = S that describe the biradical character of the system. As a rule, R and S will be HOMO and LUMO of the closed-shell guess. They are shown schematically for 1 and 2 in Figure 2a.
The unrestricted wave function takes the following form:
Ijms 03 00360 i001
where A is the antisymmetrizer, Ψ core U denotes the (closed-shell) core part, and Ψ open U describes the open-shell part. We use the symbol U rather than UHF to indicate that Eq. (1) and the following equations are equally valid for UHF and UDFT. Actually, Ψ core U will no longer represent an exact closed-shell function in the final unrestricted wave function due to the independent optimization of the α and β spin orbitals. There are two ways to build Ψ open U in the initial guess. If the HOMO and LUMO are of σ and π character or, more generally, belong to different irreducible representations of the molecular symmetry, these two ways will lead to qualitatively different results.
Permuted-Orbitals (PO) UDFT. The ordering of the orbitals is changed for one of the spin orientations; typically, the HOMO and LUMO are exchanged for β spin orientation (see Figure 2a and Figure 2b):
Ijms 03 00360 i002
where non-overlined and overlined symbols denote α and β spin orbitals, respectively. In the special case of molecules 1 and 2, Eq. (2) takes the form
Ijms 03 00360 i003
As shown in Figure 2b for 1, the resulting initial state is a mixture of a σπ singlet (S) and a σπ triplet (T) component with equal weights,
Ijms 03 00360 i004
Ijms 03 00360 i005a
Ijms 03 00360 i005b
i.e. the spin symmetry of the closed-shell initial state is broken. Regarding spatial symmetry, the PO-UHF (PO-UDFT) reference state will belong to a one-dimensional irreducible representation, which is antisymmetric with respect to the mirror plane (1,3B1 for 1, 1,3A″ for 2), in distinction to the closed-shell initial guess, which belongs to the totally symmetric irreducible representation (1A1 for 1, 1A′ for 2).
Broken-Symmetry (BS) UDFT. In BS-UDFT, the initial guess is constructed as follows:
Ijms 03 00360 i006a
Ijms 03 00360 i006b
Ijms 03 00360 i006c
with a rotation angle θ that is optimized during the SCF iterations. After the optimization, ϕr and ϕs are the localized counterparts of ϕR and ϕS, respectively (see 3 in Figure 2a). They resemble the Generalized-Valence-Bond (GVB) orbitals as is shown in Figure 2b. The reference state Ψ open BS−U can be rewritten as
Ijms 03 00360 i007
The resulting state is a mixture of the | σ σ ¯ 〉 and | π π ¯ 〉 singlet states and the| σ π ¯ T (MS = 0) triplet state, i.e. the spin symmetry is broken as in the case of PO-UHF (PO-UDFT). With respect to spatial symmetry, the BS-UHF (BS-UDFT) reference state does not belong to any irreducible representation but is part of a reducible representation that consists of a 1A1 and a 3B1 contribution for 1 and a 1A′ and a 3A″ contribution for 2. However, it is important to note that the BS-UHF (BS-UDFT) reference function is not completely asymmetric but represents an irreducible representation of a mixed spin-space symmetry group where all reflections at the mirror plane are combined with a simultaneous flip of all spins in the system. The PO-UHF (PO-UDFT) wave function is not symmetric with respect to this combined symmetry group.
The BS-UHF (BS-UDFT) wave function (7) covers in principle the GVB wave function and, therefore, is able to describe two-configurational problems such as biradicals where the rotational angle θ determines the amount of biradical character. Since for a given problem it is sometimes difficult to determine the optimized θ value and the exact form of (7), we use the natural orbital occupation numbers (NOON) of the mixing orbitals to calculate the biradical character covered by the BS-UHF (BS-UDFT) wave function.
Both for PO-UHF (PO-UDFT) and for BS-UHF (BS-UDFT), the wave function belongs to an irreducible representation of an appropriately chosen symmetry group, and the breaking of the spin symmetry does not imply a reduction of the molecular symmetry, i.e. the molecule remains mirror-symmetric. If the biradical does not possess any symmetry as in the OSS state of 90°-rotated alkene 4, the spatial symmetry of the PO-UHF (PO-UDFT) wave function vanishes while the spin symmetry remains. The BS-UHF (BS-UDFT) wave function obtains, beside the three components of Eq. (7), a contribution of the state | R S ¯ S. The symmetry of the BS-UDFT wave function with respect to the combined space-spin symmetry group vanishes. Eventually, the PO-UHF (PO-UDFT) and BS-UHF (BS-UDFT) wave functions are no longer qualitatively different.

3 Description of Open Shell Singlet Biradicals by DFT: Principal Shortcomings of UDFT

Kohn-Sham DFT relates the real electronic system (atom or molecule) to a reference system of non-interacting particles that has the same electron density distribution ρ(r) as the real electronic system The total energy E of the real system can then be represented (compare with Figure 1) as
Ijms 03 00360 i008a
Ijms 03 00360 i008b
where EX is the exchange energy, EC the correlation energy, Ts the kinetic energy, Vext the external interaction energy, and VH the Coulomb interaction energy. In distinction to E(λ = 0) and EX, the correlation energy EC cannot be calculated exactly. Thus, practical DFT calculations have to use approximate expressions for EC that rest upon the description of the correlation hole in a homogeneous or weakly inhomogeneous electron gas A similar approximation is made for EX. Even though EX can be calculated exactly, the approximate expressions for EX within the exchange-correlation energy are superior in accuracy due to a cancellation of errors between EX and EC.
Standard Kohn-Sham DFT calculations rely on assumption (1) and assumption (2) (see e.g [26, 27]):
  • The reference state is a single-determinant state. This assumption is essential for the calculation of EX.
  • The correlation hole in the real system is described reasonably well by the model correlation hole from the homogeneous or weakly inhomogeneous electron gas. This assumption is essential for the calculation of EC.
In singlet biradicals of type I, II, or III, one or both of these assumptions are not fulfilled. Therefore, we have to investigate whether a DFT description of these systems is still possible and in particular whether these systems can be described with UDFT.
The reference wave function (i.e. the wave function for λ = 0 in the adiabatic connection scheme [33]) of a type-I singlet biradical takes the following form:
Ijms 03 00360 i009
where the open-shell part ΨTD is now represented by two Kohn-Sham determinants Φ0 and Φ′0 possessing identical weight factors (Figure 1) according to Eq. (5a):
Ijms 03 00360 i010a
or in the special case of molecules 1 and 2 (Figure 2b)
Ijms 03 00360 i010b
Assumption (1) of standard Kohn-Sham DFT is not fulfilled for type I systems (OSS states of 1, 2 and similar biradicals) and, therefore, neither RDFT nor UDFT is suitable for describing these systems. A proper DFT description of OSS biradicals implies that the functional EX is modified. This does not simply concern the question of what approximate functional has to be chosen but the question of the appropriate form of the EX functional. The exact exchange energy of ΨTD can be written in the form
Ijms 03 00360 i011a
Ijms 03 00360 i011b
The exchange energy of the triplet component | R S ¯ T (MS = 0) is identical to that of the triplet components |RS〉 (MS = 1) and |RS〉 (MS = −1), which is readily expressed in UDFT as is EX (| R S ¯ 〉) for function (2) in the spin-asymmetric state. KRS is the exchange integral between orbitals ϕR and ϕS. Hence, Eq. (11) can be used as a starting point for constructing KS schemes for type I systems, and a number of such schemes have been suggested in the literature [35, 36, 37]. The restricted open-shell singlet (ROSS) scheme by Gräfenstein, Kraka, and Cremer [35] starts from Eq. (11a). Frank and co-workers [36] developed a scheme based on Eq. (11b). The Restricted Open-shell KS (ROKS) formalism by Filatov and Shaik [37] is a generalization of the approach of Frank and co-workers [36] and makes it possible to treat type I systems as well as other low-spin open-shell states. The three methods differ slightly in the way the correlation energy is determined. However, these differences are of minor importance for the computational results and we refer to the original publications for details [35, 36, 37].
Eq. (11b) reminds of the sum rule that was used by several authors [38] to calculate the OSS energy of type I systems. However, the formalisms of Ref.s [36] and [37] do not employ independent DFT calculations for the states |RS〉 (MS = 1 to represent | R S ¯ T) and | R S ¯ 〉 of Eq. (11b). Instead, Eq. (11b) is used to set up one KS scheme and the system is described based on a single Kohn-Sham DFT calculation.
Table 1 summarizes the results of ROSS, ROKS, and UDFT calculations for the three lowest states (3B1, 1A1, and 1B1) of 1. Dunning’s correlation-consistent polarized valence triple-zeta basis set cc-pVTZ [39] was used in the calculations. The GGA (generalized gradient approximation) functional BLYP [40] and the hybrid functional B3LYP [41] were employed to obtain representative values. The ROSS calculations were carried out with the program package COLOGNE 99 [42], the ROKS calculations with a local version of CADPAC 5 [43], and the standard-DFT calculations with the Gaussian 98 program package [44]. Results from CAS-DFT calculations [32, 45] are also included in Table 1 to have an appropriate reference at the DFT level that includes multi- and not just two-configurational effects in a systematic manner. For each state, the calculated NOONs for the σ and π frontier orbitals of 1 are given where an integer number indicates that the NOON value is fixed by the method used.
Since there is no experimental value available for the 1B13B1 energy splitting of 1, Bauschlicher’s [46] configuration interaction (CI) value of 33.4 kcal/mol, which is derived from a series of calculations with varying basis set, was used as a suitable reference value. Here and in the following, positive values for singlet-triplet splittings indicate a triplet ground state, negative ones a singlet ground state. For the 1A13B1 splitting, the experimental value of McKellar and co-workers [47] corrected for vibrational effects [48] was used. Experimental [49, 50] and calculated geometries are also compared in Table 1.
DFT provides reasonable descriptions for the two lowest states of 1, which can compete with high level WFT results. For example, the BLYP and B3LYP results for the 1A13B1 splitting obtained with restricted DFT methods are 9.4 (REKS: 8.8 kcal/mol, Table 1), and 10.5 kcal/mol, respectively, which have to be compared with an experimental value of 9.1 kcal/mol. UDFT exaggerates the splitting (Table 1, numbers in parentheses) somewhat because the triplet state is better described than the singlet state. CAS-DFT underestimates the splitting, however the value of 6.9 kcal/mol is still better than that of CASPT2 [51] and other high level calculations [48]. The performance of DFT is a direct reflection of the fact that the biradical character of the 1A1 state of 1 is rather small (8.7 % at GVB-DFT; 8.8 % at CAS-DFT, see Table 1). Hence, an RDFT description of this state (0 % biradical character) does not lead to a large error provided it covers sufficient dynamic electron correlation.
The situation changes when describing the OSS state 1B1 of 1 and the 1B13B1 splitting. ROSS and ROKS lead to reasonable values (30.4 and 30.8 kcal/mol, Table 1), which are in line with both the CI (33.4 kcal/mol) and the CAS-DFT result (35.2 kcal/mol, Table 1). PO-UDFT, in contrast, gives a 1B13B1 splitting of 10.1 (BLYP) and 11.4 kcal/mol (B3LYP, Table 1), respectively, i.e. the splitting is more than 20 kcal/mol off the reference value. Also, the 1B1 rather than the 1A1 state is predicted to be the lowest singlet state of 1. For BS-UDFT, the situation deteriorates even more in the way that the 1B13B1 splitting is predicted to be 4-5 kcal/mol (Table 1). The geometry of the 1B1 state is now characterized by unreasonably long CH bonds (1.09 and 1.10 Å, Table 1) and a HCH angle (114°) more than 30° smaller than predicted by ROSS or ROKS. The NOON values (1.53 and 0.47, Table 1) reveal that BS-UDFT underestimates the biradical character by 53 % while PO-UDFT actually provides the correct value of 100 %.
By inspection of Eq.s (4) and (11) as well as the diagrams for type I systems in Figure 1, the dramatic failure of UDFT in the case of 1 becomes understandable. The PO-UDFT energy of the 1B1 state of 1 has to be too low by at least the magnitude of the exchange integral KRS, which is 15.6 kcal/mol in the case of 1. We conclude that even a qualitatively correct PO-UDFT-description of 1B1-1 is not possible. A correct description of the exchange interaction between the unpaired electrons can only be provided by ROSS and ROKS (REKS) or the more general CAS-DFT method. In the case of 1, already a two-configurational WFT description (covering the two-determinantal character of the wave function for a type I system) obtained for example with ROHF for open-shell singlet state [53c] and combined with the DFT correlation energy (ROHF-DFT [53]) already leads to reasonable excitation energies (25.8 kcal/mol, Table 1; with the 6-31G(d, p) basis 29.6 kcal/mol are obtained). At this level of theory, the 1A1 state is described by GVB-DFT [53], which includes nondynamic electron correlation effects thus providing for this state a better description than the ROHF-DFT description can provide for the triplet state. Consequently, the singlet(1A1)-triplet(3B1) splitting is strongly underestimated at this level of theory (1.2 kcal/mol, Table 1). CAS-DFT provides a more balanced account of the excitation energies of 1 although the underestimation of the 1A13B1 splitting seems to be typical of WFT-DFT combination methods and has been attributed to the double-counting of correlation effects [45]. The failure of BS-UDFT in the description of 1 will be discussed in detail in chapter 5.

4 Why does UDFT Give Reasonable Results for Certain Open Shell Singlet Biradicals?

Table 2 lists results of DFT calculations for the lowest 1A″ and 3A″ states of 2. The corresponding geometries are shown in Figure 3a. Experiments [3, 7] indicate that the 1A″ state is the ground state of 2 and that the singlet-triplet splitting should be smaller than 5 kcal/mol. High-level CI calculations [21b, 54] suggest a singlet-triplet splitting of about −2 kcal/mol, i.e. the 1A″ below the 3A″ state in line with the experimental observations.
¿From the discussion in chapter 3, one should expect that UDFT fails to describe the relative energies of the singlet and triplet states of 2. However, PO-UDFT predicts a singlet-triplet splitting of -0.5 and -1.0 kcal/mol, respectively, while both ROSS and ROKS calculations give erroneous singlet-triplet ordering as reflected by positive singlet-triplet splittings (0.8 and 0.9 kcal/mol, Table 2). Although these results seem to contradict the discussion in chapter 2, they become understandable in view of considerations (a) and (b):
(a) According to Eq. (11a), the basic failure of PO-UDFT to correctly describe the wave function of the OSS state of a type I system becomes more apparent the larger the exchange integral KRS is. If the two unpaired electrons are well separated as in the case of 2, KRS will become small and the basic misconception of UDFT is disguised. For 2, KRS is equal to 0.42 kcal/mol at B3LYP/cc-pVTZ and the artificial lowering of the UDFT energy of the OSS state is negligible.
(b) It is known [32, 45, 55] that the 1A″ and 3A″ states of 2 show spin-polarization effects in the π-electron system, which require a multi-configuration description, i.e. strictly speaking 2 is a type III rather than a type I system. In the 1A″ state, this leads to stabilizing αα and ββ spin interactions at the single electron centers while for the 3A″ state one of these interactions is replaced by a destabilizing single-center αβ interaction. Contrary to the one-center biradical 1, biradical 2 can only be correctly described if the method used covers these spin polarization effects. ROSS and ROKS as restricted methods have to use identical α and β spin orbitals and, therefore, they are not able to cover spin polarization effects whereas UDFT can describe these effects due to the flexibility of the spin-unrestricted orbitals.
The advantages of (spin symmetry broken) PO-UDFT prove more important than the inappropriate description of EX for the description of the singlet-triplet splitting of 2 or similar biradicals. This will no longer be true if magnetic properties or in general properties that depend on the spin density distribution are considered. A description that both accounts for static and dynamic correlations such as CAS-DFT [32, 45] is in any case preferable in such a situation because it is based on the correct wave function, describes exchange more completely, and covers also spin polarization effects. In the case of 2, CAS-DFT gives the singlet-triplet splitting of 2 as -3 kcal/mol [45] in line with difference-dedicated configuration-interaction (DDCI2) calculations [54]. GVB-DFT [53] overestimates the stability of the singlet state similar in the case of 1 (Table 2), which underlines the necessity of a balanced and, therefore, multi- rather than two-configurational approach.

5 UDFT Description of Type II Biradicals

Biradical 3 represents a type II biradical, the wave function of which can be set up with the help of Eq. (12).
Ijms 03 00360 i012
The configurations | R R ¯ 〉 and | S S ¯ 〉 in Eq. (12) are non-equivalent but nearly degenerate, i.e neither |C2| nor C1 is small. The biradical character of wave function (12) depends on the ratio |C2/C1|, which becomes clear when expressing (12) in terms of the (localized) GVB orbitals ϕr and ϕs.
Ijms 03 00360 i013
The first term in Eq. (13) gives the ionic terms, which become less important in the way the importance of the doubly excited state | S S ¯ 〉 increases. The second term describes the singlet biradical state, which gains importance with decreasing ionic terms. Hence, the ratio |C2/C1| is a direct indicator for the biradical character and, by this, determines the NOON values of frontier orbitals ϕR and ϕS.
The prototype of a type II system is the stretched H2 molecule, where ϕR and ϕS are the and * orbitals. Type II systems possess long-range non-dynamic electron correlations, which contradicts assumption (2) of standard Kohn-Sham DFT. An RDFT treatment of these systems, although formally a correct approach, will in practice seriously underestimate electron correlation and yield too high an energy for the singlet state [11, 30] as is indicated in Figure 1. One approach to treat type II systems with DFT is to start from a two-configuration reference function that covers the static correlation effects and to use DFT for the description of dynamic correlation effects only as is done in the GVB-DFT approach [53, 56], or the two-configuration (TC)-DFT methods suggested by Borowski et al. Filatov and Shaik [58] developed the restricted-ensemble KS (REKS) method, which represents the total energy of a multi-configuration system as a weighted average over an ensemble of single-determinant states. In the case of type I systems, the REKS method becomes equivalent to the ROKS method.
BS-UDFT provides an alternative way of describing type II systems (compare, Eq.s 12 and 7). The electrons in localized orbitals ϕr and ϕs (Figure 2a) of the BS-UDFT wave function (6a) avoid each other in a similar way as in the type II reference state (12). As was discussed in detail by Gräfenstein and co-workers [11], the static correlation energy is in this way expressed as a part of the exchange energy. Hence, UDFT can lead to a reasonable description of type II biradicals provided a number of points are considered.
(a) Hybrid functional such as B3LYP in combination with BS-UDFT cover static electron correlation effects in a better way than GGA functionals such as BLYP. This becomes obvious when calculating the NOON values of the BS-UDFT description of 3: BLYP predicts a biradical character of just 51 % and by this 20 % below the corresponding B3LYP value (Table 3) or the value of high level WFT calculations [10]. The low biradical character is related to the fact that BLYP leads to a more stable RDFT description of 3 than B3LYP [11]. Hence, BLYP emphasizes the closed-shell singlet structure more than B3LYP, which is also indicated by the calculated geometry of the 1Ag state of 3 (Figure 3b): CC bond length alternation is enhanced, which is indicative of enhanced through-bond coupling of the single electrons and a reduction of the biradical character [12, 13b].
Recently, we have demonstrated with the help of electron density studies using ab initio densities as suitable reference densities that the self-interaction error of commonly used exchange functionals leads to the inclusion of nondynamic electron correlation effects in an unspecified way thus suppressing an efficient introduction of static electron correlation via a two-(multi-) configurational wave function. Since hybrid functionals such as B3LYP reduce the self-interaction error, static electron correlation effects mimicked by the exchange functional play a smaller role and, accordingly, can be more effectively introduced via the two-(multi-) configurational form of the wavefunction [59].
(b) The use of hybrid functionals leads to another advantage. Eq. (7) shows that the UDFT wave function for the OSS state is a superposition of singlet and triplet functions. For spin-unrestricted WFT, it was suggested that the influence of the triplet contribution is compensated by calculating the energy ES of the biradical state as an average using the sum formula (14) [38]
Ijms 03 00360 i014
Ijms 03 00360 i015
where ET is the energy of the corresponding triplet state and S ^ 2 is the operator of the total spin. Application of this formula for pure DFT makes little sense, it becomes however useful as soon as hybrid functionals are used. As discussed in Ref. [11], Eq. (14) can be exploited to improve BS-UDFT results for type II systems.
(c) The magnitude of the singlet-triplet splitting provides information about the reliability of the BS-UDFT results, which are reasonable in the case of small splittings (< 3 kcal/mol), but become unreliable for large splittings (> 5 kcal/mol). Appropriate test calculations have to clarify this point in each case.
Table 3 shows RDFT, UDFT, REKS, GVB-DFT, and CAS-DFT results for the lowest singlet (1Ag) and triplet (3B1u) states of 3. The corresponding geometries are shown in Figure 3b. It is known from experiment [6] that 3 possesses a singlet ground state, the singlet-triplet splitting being −3.5± 0.5 kcal/mol. RDFT gives a singlet-triplet splitting of 13.1 kcal/mol for B3LYP and -0.1 kcal/mol for BLYP, respectively. Hence, the 1Ag state is above or at equal level as the 3B1u state contrary to experimental observations. Inspection of Figure 1 reveals that this failure of DFT results from the fact that within the single determinant approach the influence of a second (or additional) configuration(s) is not considered (NOON values of 2 and 0 for the b1u-HOMO and the ag-LUMO, Table 2) and, therefore, the stabilization of the 1Ag state resulting from static electron correlation is missed. BS-UDFT partially compensates this error by adding 72 % biradical character to the DFT wave function (NOON values: 1.28 and 0.72, Table 3).
The B3LYP functional is more sensitive to indicate this basic failure of RDFT and, therefore, RDFT/B3LYP deviates more from the true description and becomes more unstable than RDFT/BLYP. This has lead to the conclusion [17, 18, 19, 20] that hybrid functionals are inappropriate for a description of 3 and related systems. Actually, if the basic error in the RDFT description is corrected by using BS-UDFT (Figure 1) it will turn out that hybrid functionals such as B3LYP combined with a sufficiently large basis set lead in general to a higher accuracy than GGA functionals such as BLYP [11]. BS-UDFT/B3LYP gives a singlet-triplet splitting of -2.6 kcal/mol, which is close to the CAS-DFT result of -2.5 kcal/mol and also close to the experimental value (Table 3) [6]. BS-UDFT/BLYP similar as REKS/BLYP exaggerate the stability of the singlet state and by this give too negative splittings (-4.5 and -5.7 kcal/mol, Table 3). REKS/BLYP itself underestimates the biradical character of the 1Ag-state (28 %, NOON values 1.72 and 0.28, Table 3) by 30-40 % relative to high level ab initio methods [10]. Consequently, the REKS/BLYP geometry of 3 (Figure 3b) resembles more a HF than a CCSD(T) geometry [10].
Use of the sum formula can lead to even more negative singlet-triplet splittings, which suggests that application of the empirical correction has to be decided from case to case. We note in this connection that the sum formula is rather useful if BS-UDFT energies are compared with RDFT values within a reaction system [11]. The largest singlet-triplet splitting is obtained by GVB-DFT, which again is a consequence of the the inadequate description of the triplet state by ROHF-DFT.
A PO-UDFT calculation for the singlet state of 3 gives at B3LYP/cc-pVTZ an energy 38.9 kcal/mol above the triplet state, i.e. about 40 kcal/mol off the correct value. The reason is that PO-UDFT is basically inappropriate to describe type II systems: the PO-UDFT KS wave function contains no contribution of the sought type-II singlet state but describes the type I 1B1u (OSS) state of 3. This state has a relatively high energy because it suffers from ionic terms, which lead to large Coulomb repulsion (see Eq. 16).
Ijms 03 00360 i016
Clearly, PO-UDFT must not be used for the description of type II systems. Also, the use of RDFT is erroneous in the case of a type II system with large biradical character. Of course, one could imagine a situation in which the biradical character of a given system is not known (unknown |C2|) and therefore RDFT is applied. By checking the stability of the RDFT solution [60] it becomes obvious that UDFT has to be applied in the case of 3 and an erroneous description can be avoided in this way.

6 Use and Misuse of UDFT

The motivation for using BS-UDFT results from the necessity of describing type II systems in the most cost-efficient way. This can lead to problems and the usefulness of this approach has to be checked from case to case. By no means is it justified to derive from a successful BS-UDFT description of type II systems the validity of this description in the case of type I biradicals. Nevertheless this has been done in the literature [23] and, therefore, these attempts deserve special consideration.
Inspection of Eq. (7) reveals that Ψ open BS−UDFT does not contain the state of interest, i.e. | R S ¯ S. Instead, states are contained in the BS-UDFT wave function, for which the orbitals ϕR or ϕS are doubly occupied. It depends on the relative energies of the | R R ¯ 〉, | S S ¯ 〉, and | R S ¯ T states and the overlap between the | R R ¯ 〉 and | S S ¯ 〉 states (which relates to the matrix element KRS) whether the BS-UDFT state is higher or lower in energy than the | R S ¯ 〉 UDFT state.
For 1, the | R R ¯ 〉 = | σ σ ¯ 〉 state is the 1A1 state, which is lower in energy than the 1B1 state to be described. This explains why BS-UDFT results in a lower energy for the 1B1 state, which is 5 kcal/mol below the PO-UDFT state thus increasing the UDFT error to 25 kcal/mol (Table 1). For 2, in contrast, the | R R ¯ 〉 and | S S ¯ 〉 states (Figure 2a) are higher in energy than both the | R S ¯ T and | R S ¯ S states, which is due to the higher Coulomb repulsion energy in the | R R ¯ 〉 and | S S ¯ 〉 states. As a consequence, the BS-UDFT state has a higher energy than the PO-UDFT state. One gets with the B3LYP functional a BS-UDFT singlet-triplet splitting of 22.5 kcal/mol, i.e. the energy of the BS-UDFT state is by 23.5 kcal/mol higher than that of the PO-UDFT state (Table 2). The biradical character (58 %, Table 2) is strongly underestimated, which is in line with the fact that the BS-UDFT wave function (7) gives ionic terms now too much weight.
A stability test [60] reveals that the BS-UDFT solution is internally unstable with respect to a rotation of the orbitals ϕr and ϕs (lowest eigenvalue λ of the stability matrix: -0.0636, see also Figure 4). This orbital rotation describes a breaking of the symmetry with respect to the combined spin-space symmetry group, which was discussed for the BS-UDFT state in chapter 2. This explains observations made by Schreiner and Prall [23]: If 2 and 3 are described with BS-UDFT employing an exactly planar (i.e. mirror-symmetric) starting geometry, one will obtain the BS-UDFT solution characterized by a high relative energy and a geometry that is unstable with respect to a breaking of the mirror symmetry. Upon distorting the starting geometry of 2 and 3 slightly toward non-planarity, the symmetry of the BS-UDFT initial guess is broken (Figure 4). The sought state | R S ¯ 〉 enters into the initial guess, and the final BS-UDFT solution takes the form of | R S ¯ 〉. The energy is lowered and the (correct) Cs symmetry of 2 is reduced to C1 symmetry.
Schreiner and Prall [23] concluded that 2 may be slightly non-planar. Furthermore, they argued that pure DFT functionals are more appropriate for the description of biradicals than hybrid functionals because some of the discrepancies observed for the planar starting geometry occur only for the B3LYP functional but not for the BLYP and BPW91 functionals. The discussion of BS-UDFT shows that both of these claims are erroneous and have to be rejected. The correct choice of the UDFT initial guess (PO-UDFT rather than BS-UDFT) confirms that the geometry of 2 is planar [61]. As in the case of the RDFT description of 3, the results for B3LYP change more strongly with an incorrect use of DFT than those for pure DFT functionals, whereas DFT, when applied correctly, gives reasonable or even better results for hybrid functionals.
In summary, we state that UDFT, if correctly applied, is in many cases an appropriate method for the description of singlet biradicals while it will fail in some cases for principal reasons. It is possible to decide in advance whether a UDFT description is appropriate and, if so, which kind of UDFT approach has to be used:
  • For a type I system, UDFT fails for principal reasons and should not be used unless special precautions are taken. For example, PO-UDFT can be used in some cases if the overlap between the open-shell orbitals is small, i.e. the two unpaired electrons sit at different, preferably non-neighboring atoms. In this situation the exchange integral KRS will be negligible and the major deficiency of UDFT is disguised so that a reasonable description seems to be obtained. The PO-UDFT description of rotated alkene 4 represents such a case The use of BS-UDFT in these cases is always erroneous.
  • For a type I system with considerable overlap between the open-shell orbitals, UDFT is always inappropriate. This applies e.g. to the OSS state of carbenes, where the two unpaired electrons are located both at the same atom. Correct results can only be obtained with ROSS [35] or ROKS [37].
  • For type II systems, BS-UDFT will give a reasonable description provided the corresponding singlet-triplet splitting is small (< 5 kcal/mol) while PO-UDFT is inappropriate because it describes an excited OSS state of the type II system.
  • Most difficult is always the description of type III systems because they possess a multiconfigurational-multideterminantal wave function. For a symmetric biradical such as 2, the PO-UDFT provides qualitatively correct energies and geometries because of a small value of KRS (related to a small singlet-triplet splitting) and the fact that UDFT describes spin polarization better than any restricted open-shell DFT method. BS-UDFT will always be erroneous as reflected by an unreasonably high OSS energy. Also, ROSS and REKS fail in these cases because they cannot describe spin-polarization correctly.
  • Neither type I, type II or type III systems can be described by RDFT. This trivial fact was overlooked for some time (see, e.g. Ref. [20]). In the case that the biradical character of a given molecule is not known one can avoid an erroneous use of RDFT by making an appropriate stability test [11, 59] of the RDFT energy.
  • The calculation of NOON values helps to determine the biradical character of a molecule and to find the appropriate method for description. This procedure is regularly used when analyzing CASSCF wave functions, however its usefulness in connection with BS-UHF or BS-UDFT descriptions is hardly exploited.
  • Hybrid functionals such as B3LYP perform better than pure DFT functionals when describing OSS biradicals simply because they reduce the self-interaction error of commonly used DFT exchange functionals. Hence, the inclusion of static electron correlation via the wave function is more effective than for GGA functionals [59].
  • Lowering the molecular symmetry for the initial guess (e.g., breaking the mirror symmetry) may remedy an incorrect choice of the UDFT method (i.e. PO-UDFT or BS-UDFT). However, this is neither necessary nor reasonable. Firstly, both the erroneous breaking of the initial symmetry and the transition from a BS-UDFT to a PO-UDFT wave function (or vice versa) require unnecessary computational expenses. Secondly, the final geometry and electronic structure in such an optimization will still be slightly asymmetric, which may lead to qualitatively incorrect results for symmetry-dependent quantities (e.g., infrared intensities).
  • The differences between PO-UDFT and BS-UDFT vanish in the case of C1 symmetry of the biradical to be investigated. The question which of these approaches is more reasonable can only be answered by inspection of the resulting UDFT orbitals. However, it is always dangerous to consider the quality of calculated molecular properties as an appropriate criterion for the reliability of the UDFT approach chosen because this may disguise basic deficiencies of PO- or BS-UDFT.
  • In any case it is advisable to use CAS-DFT as a more complete method to check the reliability of a given UDFT description. GVB-DFT is not reliable enough to fulfill this task as was demonstrated for type II and type III systems.
It is still a widespread belief that the reliability of UDFT results should be tested with the help of the expectation value 〈 S ^ 2 〉. Various authors [11, 62, 63] have pointed out that 〈 S ^ 2 〉 is of little use in this connection: a large value of 〈 S ^ 2 〉 may be obtained in a case in which the singlet-triplet splitting is small and by this the BS-UDFT description actually rather reliable (see above). Apart from this one has to emphasize that 〈 S ^ 2 〉 in DFT has no longer the same meaning as in WFT. a) The KS wave function corresponds to the situation of noninteracting electrons, which differs considerably from the many-particle state represented by the KS density. Hence, one cannot expect that 〈 S ^ 2 〉 calculated from the KS Slater determinant has any physical significance. -b) The KS wave function can only describe one-particle quantities correctly while the expectation value 〈 S ^ 2 〉 is a two-particle quantity. Hence, the value of 〈 S ^ 2 〉 calculated with the KS Slater determinant is not correct and an expectation value for S ^ 2 that deviates from S(S + 1) does not necessarily indicate spin contamination.

Acknowledgment

This work was supported by the Swedish Natural Science Research Council (NFR). Calculations were done on the supercomputers of the Nationellt Superdatorcentrum (NSC), Linköping, Sweden. The authors thank the NSC for a generous allotment of computer time.

References

  1. Marquardt, R.; Balster, A.; Sander, W.; Kraka, E.; Cremer, D.; Radziszewiski, J. G. p-benzyne. Angew. Chem. Int. Ed. 1998, 37, 955–958. [Google Scholar] [CrossRef]
  2. (a) Myers, A. G. Proposed structure of the neocarzinostatin chromophore-Methyl thioglycolate adduct - a mechanism for the nucleophilic activation of neocarzinostatin. Tetrahedron Lett. 1987, 28, 4493–4496. [Google Scholar] (b) Myers, A. G.; Proteau, P. J.; Handel, T. M. Stereochemical assignment of neocarzinostatin chromophore - structures of neocarzinostatin chromophore methyl thioglycolate adducts. J. Am. Chem. Soc. 1988, 110, 7212–7214. [Google Scholar] (c) Myers, A. G.; Dragovich, P. S.; Kuo, E. Y. Studies on the thermal generation and reactivity of a class of (σ, π)-1,4-biradicals. J. Am. Chem. Soc. 1992, 114, 9369–9386. [Google Scholar]
  3. Sander, W.; Wandel, H.; Bucher, G.; Gräfenstein, J.; Kraka, E.; Cremer, D. α, 3-didehydro-5-methyl-6-hydroxytoluene: Matrix isolation of a diradical related to the neocarzinostatin chromophore. J. Am. Chem. Soc. 1998, 120, 8480–8485. [Google Scholar] [CrossRef]
  4. (a) Jones, R. R.; Bergman, R. G. p-Benzyne. Generation as an intermediate in a thermal isomerization reaction and trapping evidence for the 1,4-benzenediyl structure. J. Am. Chem. Soc. 1972, 94, 660–661. [Google Scholar] (b) Bergman, R. G. Reactive 1,4-dehydroaromatics. Acc. Chem. Res. 1973, 6, 25–31. [Google Scholar] (c) Lockhart, T. P.; Comita, P. B.; Bergman, R. G. Kinetic evidence for the formation of discrete 1,4-dehydrobenzene intermediates. Trapping by inter- and intramolecular hydrogen atom transfer and observation of high-temperature CIDNP. J. Am. Chem. Soc. 1981, 103, 4082–4090. [Google Scholar] (d) Lockhart, T. P.; Bergman, R. G. Evidence for the reactive spin state of 1,4-dehydrobenzenes. J. Am. Chem. Soc. 1981, 103, 4091–4096. [Google Scholar] (e) Gleiter, R.; Kratz, D. Conjugated enediynes - An old topic in a different light. Angew. Chem. Int. Ed. 1993, 32, 842–845. [Google Scholar] (f) Turro, N. J.; Evenzahav, A.; Nicolaou, K. C. Photochemical analog of the bergman cycloaromatization reaction. Tetrahedron Lett. 1994, 35, 8089–8092. [Google Scholar] (g) Kuwatani, Y.; Ueda, I. Synthesis of new Sworski-type dehydroannulenes - 3,4-benz-1,2,5,6,8,9,12,13-octadehydro[14]annulenes. Angew. Chem. Int. Ed. 1995, 34, 1892–1894. [Google Scholar]
  5. For reviews see: (a) Enediyne Antibiotics as Antitumor Agents; Borders, D. B.; Doyle, T. W. (Eds.) Marcel Dekker: New York, 1995. (b) Neocarzinostatin: The Past, present, and Future of an Anticancer Drug; Maeda, H.; Edo, K.; Ishida, N. (Eds.) Springer: New York, 1997. (c) Nicolaou, K. C.; Smith, A. L. Molecular design, chemical synthesis, and biological action of enediynes. Acc. Chem. Res. 1992, 25, 497–503. [Google Scholar] (d) Nicolaou, K. C.; Dai, W.-M. Chemistry and biology of the enediyne anticancer antibiotics. Angew. Chem. Int. Ed. 1991, 30, 1387–1416. [Google Scholar] (e) Pogozelski, W. K.; Tullius, T. D. Oxidative strand scission of nucleic acids: Routes initiated by hydrogen abstraction from the sugar moiety. Chem. Rev. 1998, 98, 1089–1107. [Google Scholar] (f) Maier, M. E.; Boße, F.; Niestroj, A. J. Design and synthesis of dynemicin analogs. Eur. J. Org. Chem. 1999, 1, 1–13. [Google Scholar] (g) Thorson, J. S.; Shen, B.; Whitwam, R. E.; Liu, W.; Li, Y.; Ahlert, J. Enediyne biosynthesis and self-resistance: A progress report. Bioorganic Chemistry 1999, 27, 172–188. [Google Scholar] (h) Wisniewski Grimssom, J.; Gunawardena, G. U.; Klingberg, D.; Huang, D. The chemistry of enediynes, enyne allenes and related compounds. Tetrahedron 1996, 19, 6453–6518. [Google Scholar] (i) Fallis, A. G. 1998 Alfred Bader Award Lecture - Tangents and targets: The synthetic highway from natural products to medicine. Can. J. Chem. 1999, 77, 159–177. [Google Scholar] (j) Caddick, S.; Delisser, V. M.; Doyle, V. E.; Khan, S.; Avent, A. G.; Vile, S. Studies toward the synthesis of natural and unnatural dienediynes 1. Approaches to a functionalised bicyclic ring system. Tetrahedron 1999, 55, 2737–2754. [Google Scholar]
  6. Wenthold, P. G.; Squires, R. R.; Lineberger, W. C. Ultraviolet photoelectron spectroscopy of the o-, p-, and p-benzyne negative ions. Electron affinities and singlet-triplet splittings for o-, m-, and p-benzyne. J. Am. Chem. Soc. 1998, 120, 5279–5290. [Google Scholar] [CrossRef]
  7. (a) Wenthold, P. G.; Wierschke, S. G.; Nash, J. J.; Squires, R. R. α, 3-dehydrotoluene - Experimental and theoretical evidence for a singlet ground-state. J. Am. Chem. Soc. 1993, 115, 12611–12612. [Google Scholar] (b) Logan, C. F.; Ma, J. C.; Chen, P. The photoelectron-spectrum of the α, 3-dehydrotoluene biradical. J. Am. Chem. Soc. 1994, 116, 2137–2138. [Google Scholar]
  8. (a) Hoffner, J. H.; Schottelius, M. J.; Feichtinger, D.; Chen, P. Chemistry of the 2,5-didehydropyridine biradical: Computational, kinetic, and trapping studies toward drug design. J. Am. Chem. Soc. 1998, 120, 376–385. [Google Scholar] (b) Logan, C. F.; Chen, P. Ab initio calculation of hydrogen abstraction reactions of phenyl radical and p-benzyne. J. Am. Chem. Soc. 1996, 118, 2113–2114. [Google Scholar] (c) Schottelius, M. J.; Chen, P. 9,10-dehydroanthracene: p-benzyne-type biradicals abstract hydrogen unusually slowly. J. Am. Chem. Soc. 1996, 118, 4896–4903. [Google Scholar] (d) Chen, P. Design of diradical-based hydrogen abstraction agents. Angew. Chem. Int. Ed. 1996, 35, 1478–1480. [Google Scholar]
  9. Roth, W. R.; Hopf, H.; Horn, C. The energy well of diradicals. 5. 1,3,5-Cyclohexatriene-1,4-diyl and 2,4-cyclohexadiene-1,4-diyl. Chem. Ber. 1994, 127, 1765–1779. [Google Scholar] [CrossRef]
  10. (a) Kraka, E.; Cremer, D. CCSD(T) investigation of the bergman cyclization of enediyne - Relative stability of o-didehydrobenzene, m-didehydrobenzene, and p-didehydrobenzene. J. Am. Chem. Soc. 1994, 116, 4929–4936. [Google Scholar] (b) Kraka, E.; Cremer, D. Ortho-benzyne, meta-benzyne, and para-benzyne - A comparative CCSD(T) investigation. Chem. Phys. Lett. 1993, 216, 333–340. [Google Scholar]
  11. Gräfenstein, J.; Hjerpe, A. M.; Kraka, E.; Cremer, D. An accurate description of the Bergman reaction using restricted and unrestricted DFT: Stability test, spin density, and on-top pair density. J. Phys. Chem. 2000, 104, 1748–1761. [Google Scholar] [CrossRef]
  12. Kraka, E.; Cremer, D. Computer design of anticancer drugs. A new enediyne warhead. J. Am. Chem. Soc. 2000, 122, 8245–8264. [Google Scholar] [CrossRef]
  13. (a) Kraka, E.; Cremer, D. Structure and stability of enediynes containing heteroatoms - A quantum chemical investigation. Theochem 2000, 506, 191–211. [Google Scholar] (b) Kraka, E.; Cremer, D. The para-didehydropyridine, para-didehydropyridinium, and related biradicals - A contribution to the chemistry of enediyne antitumor drugs. J. Comp. Chem. 2001, 22, 216–229. [Google Scholar]
  14. Koga, N.; Morokuma, K. Comparison of biradical formation between enediyne and enyne allene - abinitio CASSCF and MRSDCI study. J. Am. Chem. Soc. 1991, 113, 1907–1911. [Google Scholar] [CrossRef]
  15. (a) Wenthold, P. G.; Paulino, J. A.; Squires, R. R. The absolute heats of formation of ortho-benzyne, meta-benzyne, and para-benzyne. J. Am. Chem. Soc. 1991, 113, 7414–7415. [Google Scholar] (b) Wenthold, P. G.; Squires, R. R. Biradical thermochemistry from collision-induced dissociation threshold energy measurements - Absolute heats of formation of ortho-benzyne, meta-benzyne, and para-benzyne. J. Am. Chem. Soc. 1994, 116, 6401–6412. [Google Scholar] (c) Wierschke, S. G.; Nash, J. J; Squires, R. R. A multiconfigurational SCF and correlation-consistent CI study of the structures, stabilities, and singlet-triplet splittings of o-benzyne, m-benzyne, and p-benzyne. J. Am. Chem. Soc. 1993, 115, 11958–11967. [Google Scholar]
  16. (a) Lindh, R.; Persson, B. J. Ab-initio study of the Bergman reaction - the autoaromatization of hex-3-ene-1,5-diyne. J. Am. Chem. Soc. 1994, 116, 4963–4969. [Google Scholar] (b) Lindh, R.; Lee, T. J.; Bernhardsson, A.; Persson, B. J.; Karlström, G. Extended ab-initio and theoretical thermodynamics studies of the Bergman reaction and the energy splitting of the singlet o-benzynes, m-benzynes, and p-benzynes. J. Am. Chem. Soc. 1995, 117, 7186–7194. [Google Scholar] (c) Lindh, R.; Ryde, U.; Schütz, M. On the significance of the trigger reaction in the action of the calicheamicin γ I 1 anti-cancer drug. Theor. Chem. Acta 1997, 97, 203–210. [Google Scholar]
  17. (a) Cramer, C. J.; Nash, J. J.; Squires, R. R. A reinvestigation of singlet benzyne thermochemistry predicted by CASPT2, coupled-cluster and density functional calculations. Chem. Phys. Lett. 1997, 277, 311–320. [Google Scholar] (b) Cramer, C. J.; Debbert, S. Heteroatomic substitution in aromatic sigma biradicals: the six pyridynes. Chem. Phys. Lett. 1998, 287, 320–326. [Google Scholar]
  18. Cramer, C. J. Bergman, aza-Bergman, and protonated aza-Bergman cyclizations and intermediate 2,5-arynes: Chemistry and challenge to computation. J. Am. Chem. Soc. 1998, 120, 6261–6269. [Google Scholar] [CrossRef]
  19. (a) Schreiner, P. R. Monocyclic enediynes: Relationships between ring sizes, alkyne carbon distances, cyclization barriers, and hydrogen abstraction reactions. Singlet-triplet separations of methyl-substituted p-benzynes. J. Am. Chem. Soc. 1998, 120, 4184–4190. [Google Scholar] (b) Schreiner, P. R. Cyclic enediynes: relationship between ring size, alkyne carbon distance, and cyclization barrier. Chem. Commun. 1998, 4, 483–484. [Google Scholar]
  20. Chen, W.-C.; Chang, N.-Y.; Yu, C.-H. Density functional study of Bergman cyclization of enediynes. J. Phys. Chem. A 1998, 102, 2584–2593. [Google Scholar] [CrossRef]
  21. (a) Engels, B.; Lennartz, C.; Hanrath, M.; Schmittel, M.; Strittmatter, M. Regioselectivity of biradical cyclizations of enyne-allenes: Influence of substituents on the switch from the Myers-Saito to the novel C2–C6 cyclization. Angew. Chem. Int. Ed. 1998, 37, 1960–1963. [Google Scholar] (b) Engels, B.; Hanrath, M. A theoretical comparison of two competing diradical cyclizations in enyne-allenes: The Myers-Saito and the novel C2–C6 cyclization. J. Am. Chem. Soc. 1998, 120, 6356–6361. [Google Scholar]
  22. Cramer, C. J.; Squires, R. R. Quantum chemical characterization of the cyclization of the neocarzinostatin chromophore to the 1,5-didehydroindene biradical. Organic Lett. 1999, 1, 215–218. [Google Scholar] [CrossRef]
  23. Schreiner, P. R.; Prall, M. Myers-Saito versus C2–C6 (“Schmittel”) cyclizations of parent and monocyclic enyne-allenes: Challenges to chemistry and computation. J. Am. Chem. Soc. 1999, 121, 8615–8627. [Google Scholar] See also (b) Prall, M.; Wittkopp, A.; Schreiner, P. R. Can fulvenes form from enediynes? A systematic high-level computational study on parent and benzannelated enediyne and enyne-allene cyclizations. J. Phys. Chem. A 2001, 105, 9265–9274. [Google Scholar] (c) Prall, M.; Wittkopp, A.; Fokin, A. A.; Schreiner, P. R. Substituent effects on the Bergman cyclization of (Z)-1,5-hexadiyne-3-enes: A systematic computational study. J. Comp. Chem. 2001, 22, 1605–1614. [Google Scholar]
  24. Raghavachari, K.; Trucks, G. W.; Pople, J. A.; Head-Gordon, M. A 5th-order perturbation comparison of electron correlation theories. Chem. Phys. Lett. 1989, 157, 479–483. [Google Scholar] [CrossRef]
  25. Handy, N. C.; Pople, J. A.; Head-Gordon, M.; Raghavachari, K.; Trucks, G. W. Sizeconsistent Brueckner theory limited to double substitutions. Chem. Phys. Lett. 1989, 164, 185–192. [Google Scholar] [CrossRef]
  26. Kohn, W.; Sham, L. J. Self-Consistent Equations Including Exchange and Correlation Effects. Phys. Rev. 1965, 140, A 1133–1138. [Google Scholar] [CrossRef]
  27. For reviews on DFT methods, see for example (a) Parr, R. G.; Yang, W. International Series of Monographs on Chemistry 16: Density-Functional Theory of Atoms and Molecules; Oxford University Press: New York, 1989. [Google Scholar] (b) Density Functional Methods in Chemistry; Labanowski, J. K.; Andzelm, J. W. (Eds.) Springer: Heidelberg, 1990. (c) Theoretical and Computational Chemistry, Vol. 2, Modern Density Functional Theory - A Tool For Chemistry; Seminario, J. M.; Politzer, P. (Eds.) Elsevier: Amsterdam, 1995. (d) Chemical Applications of Density Functional Theory; ACS Symposium Series 629; Laird, B. B.; Ross, R. B.; Ziegler, T. (Eds.) American Chemical Society: Washington, DC, 1996. (e) Lecture Notes in Physics, Density Functionals: Theory and Applications; Joubert, D. (Ed.) Springer: Hei-delberg, 1997. (f) Recent Advances in Computational Chemistry, Vol. 1, Recent Advances in Density Functional Methods, Part II; Chong, D. P. (Ed.) World Scientific: Singapore, 1997. (g) Electronic Density Functional Theory, Recent Progress and New Directions; Dobson, J. F.; Vignale, G.; Das, M. P. (Eds.) Plenum Press: New York, 1998. (h) Gill, P. Encyclopedia of Computational Chemistry; Schleyer, P. v. R., Allinger, N. L., Clark, T., Gasteiger, J., Kollman, P. A., Schaefer, H. F., III, Schreiner, P. R., Eds.; Wiley: Chichester, UK, 1998; Vol. 1, p. 678. [Google Scholar]
  28. Szalay, P. G. Modern Ideas in Coupled-Cluster Methods; Bartlett, R. J., Ed.; World Scientific: Singapore, 1997; p. 81. [Google Scholar]
  29. Kraka, E.; Cremer, D.; Bucher, G.; Wandel, H.; Sander, W. A CCSD(T) and DFT investigation of m-benzyne and 4-hydroxy-m-benzyne. Chem. Phys. Lett. 1997, 268, 313–320. [Google Scholar] [CrossRef]
  30. For a recent review, see Bally, T.; Borden, W. T. Reviews Comp. Chem.; Lipkowitz, K. B., Boyd, D. B., Eds.; Wiley: New York, 1999; Vol. 13, p. 1. [Google Scholar]
  31. For related work, see (a) Ziegler, T.; Rauk, A.; Baerends, E. J. On the calculation of multiplet energies by the Hartree-Fock-Slater method. Theor. Chim. Acta 1977, 43, 261–271. [Google Scholar] (b) Noodleman, L.; Case, D. A. Density-functional theory of spin polarization and spin coupling in iron-sulfur clusters. Adv. Inorg. Chem. 1992, 38, 423–470. [Google Scholar] (c) Lovell, T.; McGrady, J. E.; Stranger, R.; MacGregor, S. Optimized structures of bimetallic systems: A comparison of full- and broken-symmetry density functional calculations. Inorg. Chem. 1996, 35, 3079–3080. [Google Scholar] (d) Noodleman, L.; Post, D.; Baerends, E. J. Symmetry breaking and ionization from symmetry equivalent inner shells and lone pairs in Xα theory. Chem. Phys. 1982, 64, 159–166. [Google Scholar] (e) Cramer, C. J.; Dulles, F. J.; Giesen, D. J.; Almlöf, J. Density-functional theory - Excited-states and spin annihilation. Chem. Phys. Lett. 1995, 245, 165–170. [Google Scholar]
  32. Gräfenstein, J.; Cremer, D. Can density functional theory describe multi-reference systems? Investigation of carbenes and organic biradicals. Phys. Chem. Chem. Phys. 2000, 2, 2091–2103. [Google Scholar] [CrossRef]
  33. (a) Langreth, D. C.; Perdew, J. P. The exchange correlation energy of a metallic surface. Solid State Commun. 1975, 17, 1425–1429. [Google Scholar] (b) Gunnarsson, O.; Lundqvist, B. I. Exchange and correlation in atoms, molecules, and solids by the spin-density-functional formalism. Phys. Rev. B 1976, 13, 4274–98. [Google Scholar]
  34. (a) Fukutome, H. Unrestricted Hartree-Fock theory and its applications to molecules and chemical reactions. Int. J. Quantum Chem. 1981, 20, 955–1065. [Google Scholar] (b) Yamaguchi, K. Symmetry and broken symmetry in molecular orbital (MO) descriptions of unstable molecules. Generalized MO theoretical studies on 1,3-dipolar species. J. Mol. Struct. (THEOCHEM) 1983, 12, 101–120. [Google Scholar]
  35. Gräfenstein, J.; Cremer, D.; Kraka, E. Density functional theory for open-shell singlet biradicals. Chem. Phys. Lett. 1998, 288, 593–602. [Google Scholar] [CrossRef]
  36. Frank, I.; Hutter, J.; Marx, D.; Parrinello, M. Molecular dynamics in low-spin excited states. J. Chem. Phys. 1998, 108, 4060–4069. [Google Scholar] [CrossRef]
  37. Filatov, M.; Shaik, S. Spin-restricted density functional approach to the open-shell problem. Chem. Phys. Lett. 1998, 288, 689–697. [Google Scholar] Application of spin-restricted open-shell Kohn-Sham method to atomic and molecular multiplet states. J. Chem. Phys. 1999, 110, 116–125.
  38. Lim, M. H.; Worthington, S. E.; Dulles, F. J.; Cramer, C .J. Chemical Applications of Density Functional Theory; ACS Symposium Series 629; Laird, B. B., Ross, R. B., Ziegler, T., Eds.; American Chemical Society: Washington, DC, 1996; p. 402. [Google Scholar]
  39. Dunning, T. H., Jr. Gaussian-basis sets for use in correlated molecular calculations. 1. The atoms boron through neon and hydrogen. J. Chem. Phys. 1989, 90, 1007–1023. [Google Scholar] [CrossRef]
  40. (a) Becke, A. D. Density-functional exchange-energy approximation with correct asymptotic-behavior. Phys. Rev. A 1988, 38, 3098–3100. [Google Scholar] (b) Lee, C.; Yang, W.; Parr, R. P. Development of the Colle-Salvetti correlation-energy formula into a functional of the electron-density. Phys. Rev. B 1988, 37, 785–789. [Google Scholar]
  41. (a) Becke, A. D. Density-functional thermochemistry. 3. The role of exact exchange. J. Chem. Phys. 1993, 98, 5648–5652. [Google Scholar] See also (b) Stephens, P. J.; Devlin, F. J.; Chabalowski, C. F.; Frisch, M. J. Ab-initio calculation of vibrational absorption and circular-dichroism spectra using density-functional force-fields. J. Phys. Chem. 1994, 98, 11623–11627. [Google Scholar]
  42. Kraka, E.; Gräfenstein, J.; Gauss, J.; Reichel, F.; Olsson, L.; Konkoli, Z; He, Z; Cremer, D. COLOGNE 99; Göteborg University: Göteborg, 1999. [Google Scholar]
  43. Amos, R. D. CADPAC 5: The Cambridge Analytic Derivatives Package, issue 5, University of Cambridge, 1992. with contributions from Alberts, I. L.; Andrews, J. S.; Colwell, S. M.; Handy, N. C.; Jayatilaka, D.; Knowles, P. J.; Kobayashi, R.; Koga, N.; Laidig, K. E.; Maslen, P. E.; Murray, C. W.; Rice, J. E.; Sanz, J.; Simandiras, E. D.; Stone, A. J.; Su, M. D.
  44. Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Zakrzewski, V. G.; Montgomery, J. A., Jr.; Stratmann, R. E.; Burant, J. C.; Dapprich, S.; Millam, J. M.; Daniels, A. D.; Kudin, K. N.; Strain, M. C.; Farkas, O.; Tomasi, J.; Barone, V.; Cossi, M.; Cammi, R.; Mennucci, B.; Pomelli, C.; Adamo, C.; Clifford, S.; Ochterski, J.; Petersson, G. A.; Ayala, P. Y.; Cui, Q.; Morokuma, K.; Malick, D. K.; Rabuck, A. D.; Raghavachari, K.; Foresman, J. B.; Cioslowski, J.; Ortiz, J. V.; Stefanov, B. B.; Liu, G.; Liashenko, A.; Piskorz, P.; Komaromi, I.; Gomperts, R.; Martin, R. L.; Fox, D. J.; Keith, T.; Al-Laham, M. A.; Peng, C. Y.; Nanayakkara, A.; Gonzalez, C.; Challacombe, M.; Gill, P. M. W.; Johnson, B.; Chen, W.; Wong, M. W.; Andres, J. L.; Gonzalez, C.; Head-Gordon, M.; Replogle, E. S.; Pople, J. A. Gaussian 98, Revision A.5, Gaussian, Inc.: Pittsburgh PA, 1998.
  45. Gräfenstein, J; Cremer, D. The combination of density functional theory with multiconfiguration methods - CAS-DFT. Chem. Phys. Lett. 2000, 316, 569–577. [Google Scholar]
  46. (a) Bauschlicher, C. W. Accurate ab initio calculations for the 1B11A1 separation in methylene. Chem. Phys. Lett. 1980, 74, 273–276. [Google Scholar] (b) Bauschlicher, C. W.; Langhoff, S. R. On the 1A13B1 separation in CH2 and SiH2. J. Chem. Phys. 1987, 87, 387–391. [Google Scholar]
  47. McKellar, A. R. W.; Bunker, P. R.; Sears, T. J.; Evenson, K. M.; Saykally, R. J.; Langhoff, S. R. Far infrared laser magnetic resonance of singlet methylene: singlet-triplet perturbations, singlet-triplet transitions, and the singlet-triplet splitting. J. Chem. Phys. 1983, 79, 5251–64. [Google Scholar] [CrossRef]
  48. Bettinger, H. F.; Schleyer, P. v. R.; Schreiner, P. R.; Schaefer, H. F., III. The Encyclopedia of Computational Chemistry; Schleyer, P. v. R., Allinger, N., Clark, T., Gasteiger, J., Kollman, P. A., Schaefer, H. F., III, Schreiner, P. R., Eds.; Wiley-Interscience: Chichester, 1998; p. 183. [Google Scholar]
  49. Petek, H.; Nesbitt, D. J.; Darwin, D. C.; Ogilby, P. R.; Moore, C. B.; Ramsay, D. A. Analysis of CH2 A1A1 (1,0,0) and (0,0,1) Coriolis-coupled states, A1A1X 3B1 spin orbit coupling, and the equilibrium structure of CH2 X 1A1 state. J. Chem. Phys. 1989, 91, 6566–6578. [Google Scholar] [CrossRef]
  50. Jensen, P.; Bunker, P. R. The potential surface and stretching frequencies of X3B1 methylene (CH2) determined from experiment using the Morse oscillator-rigid bender internal dynamics Hamiltonian. J. Chem. Phys. 1988, 89, 1327–1332. [Google Scholar] [CrossRef]
  51. Andersson, K.; Roos, B. Multiconfigurational 2nd-order perturbation-theory - A test of geometries and binding-energies. Int. J. Quant. Chem. 1993, 45, 591–607. [Google Scholar] [CrossRef]
  52. Bobrowicz, F. W.; Goddard, W. A., III. Methods of Electronic Structure Theory, Modern Theoretical Chemistry; Vol. 3, Schaefer, III, Ed.; Plenum Press: New York, 1977; p. 79. [Google Scholar]
  53. (a) Kraka, E. Homolytic dissociation-energies from GVB-LSDC calculations. Chem. Phys. 1992, 161, 149–153. [Google Scholar] (b) Kraka, E.; Cremer, D.; Nordholm, S. Molecules in Natural Science and Biomedicine; Maksic, Z., Eckert-Maksic, M., Eds.; Ellis Horwood: Chichester, 1971; p. 351. [Google Scholar] (c) In the case of the OSS state of 1 or 2, the GVB-DFT description simplifies to a ROHF-DFT description. See, Davidson, E. R. Spin-restricted open-shell self-consistent-field theory. Chem. Phys. Lett. 1973, 21, 565–567. [Google Scholar]
  54. Cabrero, J.; Ben-Amor, N.; Caballol, R. Singlet-triplet gap in α-n-dehydrotoluene and related biradicals: An ab initio configuration interaction study. J. Phys. Chem. A 1999, 103, 6220–6224. [Google Scholar] [CrossRef]
  55. Wenthold, P. G.; Wierschke, S. G.; Nash, J. J.; Squires, R. R. Biradical thermochemistry from collision-induced dissociation threshold energy measurements 2. Experimental and theoretical studies of the mechanism and thermochemistry of formation of α, n-dehydrotoluene biradicals from gas-phase halide elimination reactions. J. Am. Chem. Soc. 1994, 116, 7378–7392. [Google Scholar]
  56. Wu, W.; Shaik, S. VB-DFT: A nonempirical hybrid method combining valence bond theory and density functional energies. Chem. Phys. Lett. 1999, 301, 37–42. [Google Scholar] [CrossRef]
  57. Borowski, B.; Jordan, K. D.; Nichols, J.; Nachtigall, P. Investigation of a hybrid TCSCF-DFT procedure. Theor. Chem. Acc. 1998, 99, 135–140. [Google Scholar] [CrossRef]
  58. Filatov, M.; Shaik, S. A spin-restricted ensemble-referenced Kohn-Sham method and its application to diradicaloid situations. Chem. Phys. Lett. 1999, 304, 429–437. [Google Scholar] [CrossRef]
  59. Cremer, D. Mol. Phys. in press.
  60. (a) Bauernschmitt, R.; Ahlrichs, R. Stability analysis for solutions of the closed shell Kohn-Sham equation. J. Chem. Phys. 1996, 104, 9047–9052. [Google Scholar] (b) Seeger, R.; Pople, J. A. Self-consistent molecular orbital methods. XVII. Constraints and stability in Hartree-Fock theory. J. Chem. Phys. 1977, 66, 3045–3050. [Google Scholar]
  61. Wenthold, P. G.; Lipton, M. A. A density functional molecular orbital study of the C2–C7 and C2–C6 cyclization pathways of 1,2,4-heptatrien-6-ynes. The role of benzannulation. J. Am. Chem. Soc. 2000, 122, 9265–9270. [Google Scholar] [CrossRef]
  62. Pople, J.; Gill, P.M.W.; Handy, N.C. Spin-unrestricted character of Kohn-Sham orbitals for open-shell systems. Int. J. Quant. Chem. 1995, 56, 303–305. [Google Scholar] [CrossRef]
  63. Gräfenstein, J.; Cremer, D. On the diagnostic value of S ^ 2 in Kohn-Sham density functional theory. Mol. Phys. 2001, 99, 981–989. [Google Scholar] [CrossRef]
Scheme 1. Schematic representation of the molecules under discussion.
Scheme 1. Schematic representation of the molecules under discussion.
Ijms 03 00360 g005
Figure 1. Schematic representation of wave function and energy of type 0, type I, type II, and type III electronic systems as described by RDFT or UDFT with approximate exchange-correlation functionals (dotted lines) or exact Kohn-Sham DFT (solid line). Numbers 0, 1, 2, etc. denote configuration state functions Ψ0, Ψ1, Ψ2, etc. or determinants Φ0, Φ1, etc. The weights w of the configuration state functions in the true wave function are schematically shown (first row of diagrams) in dependence of the parameter λ that stepwise increases the electron correlation energy Ec from zero (λ = 0) to its true value (λ = 1). The corresponding changes in the Kohn-Sham energy EKS and in the exact energy E(exact) is schematically given in the second row of diagrams.
Figure 1. Schematic representation of wave function and energy of type 0, type I, type II, and type III electronic systems as described by RDFT or UDFT with approximate exchange-correlation functionals (dotted lines) or exact Kohn-Sham DFT (solid line). Numbers 0, 1, 2, etc. denote configuration state functions Ψ0, Ψ1, Ψ2, etc. or determinants Φ0, Φ1, etc. The weights w of the configuration state functions in the true wave function are schematically shown (first row of diagrams) in dependence of the parameter λ that stepwise increases the electron correlation energy Ec from zero (λ = 0) to its true value (λ = 1). The corresponding changes in the Kohn-Sham energy EKS and in the exact energy E(exact) is schematically given in the second row of diagrams.
Ijms 03 00360 g001
Figure 2a. Schematic representation of HOMO and LUMO for biradicals 1 and 2. For 3, HOMO and LUMO are given in canonical (top) and localized (bottom) representation.
Figure 2a. Schematic representation of HOMO and LUMO for biradicals 1 and 2. For 3, HOMO and LUMO are given in canonical (top) and localized (bottom) representation.
Ijms 03 00360 g002a
Figure 2b. Schematic representation of the construction of PO-UDFT and BS-UDFT wave function in the case of 1. The changes in the orbital form are indicated at the top (PO-UDFT) and at the bottom (BS-UDFT).
Figure 2b. Schematic representation of the construction of PO-UDFT and BS-UDFT wave function in the case of 1. The changes in the orbital form are indicated at the top (PO-UDFT) and at the bottom (BS-UDFT).
Ijms 03 00360 g002b
Figure 3. Calculated CC bond lengths (in Å) for (a) α, 3-didehydrotoluene (2) and (b) 1, 4-didehydrobenzene (3).
Figure 3. Calculated CC bond lengths (in Å) for (a) α, 3-didehydrotoluene (2) and (b) 1, 4-didehydrobenzene (3).
Ijms 03 00360 g003
Figure 4. Schematic representation of open-shell orbitals ϕR = R and ϕS = S for 2 after mixing within the BS-UDFT approach (right side, middle), after repeated mixing because of internal instability (right side, bottom), and after mixing caused by geometry distortions (right side, top). Atoms Cα, C1, C2, and C3 are represented by dots. The geometry distortion is strongly exaggerated to represent orbital mixing effects.
Figure 4. Schematic representation of open-shell orbitals ϕR = R and ϕS = S for 2 after mixing within the BS-UDFT approach (right side, middle), after repeated mixing because of internal instability (right side, bottom), and after mixing caused by geometry distortions (right side, top). Atoms Cα, C1, C2, and C3 are represented by dots. The geometry distortion is strongly exaggerated to represent orbital mixing effects.
Ijms 03 00360 g004
Table 1. Energies, geometries and NOON values for the three lowest states of methylene (1). a)
Table 1. Energies, geometries and NOON values for the three lowest states of methylene (1). a)
MethodFunctionalE or ∆ ErCHHCHNOON(a1)NOON(b1)
3B1 State
RODFTB3LYP−39.166681.077133.911
BLYP−39.141391.084134.411
ROKSBLYP−39.142141.083134.611
PO-UDFTB3LYP−39.168741.077135.01.0001.000
BLYP−39.143961.083135.51.0001.000
ROHF-DFT bLYP−39.146081.07133.911
CAS(6,6)-DFT cCS−39.120671.102131.01.0001.000
Expt. d1.070134.011
1A1 State
RDFTB3LYP e10.5 (11.8)1.110101.620
BLYP e9.4 (10.4)1.120100.920
REKSBLYP8.81.119101.11.9800.020
GVB-DFT bLYP6.31.110101.61.9860.014
CAS(6,6)-DFT cCS6.91.138100.01.9090.088
Expt. d9.11.111101.9N/AN/A
1B1 State
ROSS-DFTB3LYP30.41.072147.511
BLYP29.81.078148.611
ROKSBLYP30.81.077148.611
PO-UDFTB3LYP11.41.077138.11.0001.000
BLYP10.11.083138.31.0001.000
BS-UDFTB3LYP5.31.101114.01.5280.472
BLYP4.41.083114.31.5250.475
ROHF-DFT bLYP25.81.072147.511
CAS(6,6)-DFT cCS35.21.131134.11.0001.000
CI f35.41.077141.5N/AN/A
a) Absolute energies E in hartree, relative energies ∆E in kcal mol−1, bond lengths in Å, bond angles in degrees. The natural orbital occupation numbers (NOON) are given for the two highest occupied MOs. Integer values denote NOON values fixed by the method used, floating point numbers are given when other than integer values are possible. All calculations were done with Dunning’s cc-pVTZ basis set [39] unless otherwise noted. The CAS-DFT calculations were carried out with the Colle-Salvetti (CS) functional.
b) Geometries taken from DFT. cc-pVTZ basis at DFT geometries.
c) From Ref. [45].
d) Experimental geometries from Ref. [49] (1A1) and Ref. [50] (1A1: 1.075 Å, 133.9 deg; and 3B1). E(1A1) − E(3B1) from Ref. [47], converted from T0 to Te according to Ref. [48].
e) Energy differences in parentheses are relative to the PO-UDFT energy of the 3B1 state.
f) E(1B1) − E(3B1) estimated by extrapolations to large basis sets, Ref. [46a]; E(1A1) − E(3B1) from SOCI+Q/ANO calculations, Ref. [46b].
Table 2. Singlet and triplet energies for α-3-didehydrotoluene (2). a)
Table 2. Singlet and triplet energies for α-3-didehydrotoluene (2). a)
NOON values
MethodFunctionalE(3A″)E(1A′)3A″, a’3A″, a”1A′, a’1A′, a”
ROSS B3LYP -270.31971 0.8 1 1 1 1
BLYP−270.199310.91111
REKSBLYP−270.206130.91111
PO-UDFTB3LYP−270.32322−1.01.0001.0001.0001.000
BLYP−270.20152−0.51.0001.0001.0001.000
BS-UDFTB3LYP−270.3232222.51.0001.0001.4230.577
GVB-DFT bLYP−270.11665−7.41111
CAS(8,8)-DFT cCS−270.03909−2.61111
DDCI2 d−2.0
Expt. e−5... 01111
a) Absolute energies E in Hartree, relative energies ∆E in kcal mol−1. All calculations were done with Dunning’s cc-pVTZ basis set [39] unless otherwise noted. The method is given for the description of the OSS state while for the triplet state either RODFT or UDFT was used. The natural orbital occupation numbers (NOON) are given for the two highest occupied MOs. Integer values denote NOON values fixed by the method used, floating point numbers are given when other than integer values are possible.
b) Triplet calculated with ROHF-DFT. ROSS Geometries used. See Ref. [53].
c) Ref. [32]. The Colle-Salvetti (CS) functional and the 6-31G(d,p) basis set was used.
d) Ref. [54].
e) Ref. [7].
Table 3. Singlet and triplet energies for 1,4-didehydrobenzene (3). a)
Table 3. Singlet and triplet energies for 1,4-didehydrobenzene (3). a)
NOON values
MethodFunctionalE(3Blu)E(1Ag)3Blu, blu3Blu, ag1Ag, blu1Ag, ag
RDFTB3LYP−230.9525413.11120
BLYP−230.85665−0.11120
REKSBLYP−230.86028−5.7111.7230.277
BS-UDFTB3LYP−230.95254−2.61.0001.0001.2810.719
BLYP−230.85606−4.51.0001.0001.4920.508
BS-UDFT + SFbB3LYP−230.95254−4.9
BLYP−230.85606−7.3
PO-UDFTB3LYP−230.9525438.91.0001.0001.0001.000
GVB-DFTLYP−230.77758−10.4111.1790.821
CAS(8,8)-DFTcCS−230.70904−2.51.0001.0001.2420.758
Expt. d−3.5 ± 0.511not known
a) Absolute energies E in Hartree, relative energies ∆E in kcal mol−1. All calculations were done with Dunning’s cc-pVTZ basis set [39] unless otherwise noted. The method is given for the description of the singlet state while for the triplet state either RODFT or UDFT was used. The natural orbital occupation numbers (NOON) are given for the two highest occupied MOs. Integer values denote NOON values fixed by the method used, floating point numbers are given when other than integer values are possible. UDFT results are from Ref. [11].
b) BS-UDFT + sum formula.
c) Ref. [32]. The Colle-Salvetti (CS) functional and the 6-31G(d,p) basis set was used.
d) Ref. [6], converted to Te according to Ref. [11].

Share and Cite

MDPI and ACS Style

Gräfenstein, J.; Kraka, E.; Filatov, M.; Cremer, D. Can Unrestricted Density-Functional Theory Describe Open Shell Singlet Biradicals? Int. J. Mol. Sci. 2002, 3, 360-394. https://0-doi-org.brum.beds.ac.uk/10.3390/i3040360

AMA Style

Gräfenstein J, Kraka E, Filatov M, Cremer D. Can Unrestricted Density-Functional Theory Describe Open Shell Singlet Biradicals? International Journal of Molecular Sciences. 2002; 3(4):360-394. https://0-doi-org.brum.beds.ac.uk/10.3390/i3040360

Chicago/Turabian Style

Gräfenstein, Jürgen, Elfi Kraka, Michael Filatov, and Dieter Cremer. 2002. "Can Unrestricted Density-Functional Theory Describe Open Shell Singlet Biradicals?" International Journal of Molecular Sciences 3, no. 4: 360-394. https://0-doi-org.brum.beds.ac.uk/10.3390/i3040360

Article Metrics

Back to TopTop