Generalized-active-space pair-density functional theory: an efficient method to study large, strongly correlated, conjugated systems

Predicting ground- and excited-state properties of open-shell organic molecules by electronic structure theory can be challenging because an accurate treatment has to correctly describe both static and dynamic electron correlation.


Introduction
Open-shell organic systems are important in several areas of chemical science including magnetism, 1-3 photovoltaic materials, 4-6 photochemistry, 7,8 light-powered devices such as switching and sensing devices, 2,9 and biological, 10,11 biomimetic, 12,13 and organometallic 14,15 catalysis. However, studying their electronic structure is challenging because they oen exhibit strong electron-correlation, i.e., their electronic structure is strongly multicongurational due to near-degeneracy correlation, oen called static correlation 16 (systems of this type will be called multireference systems, and those without strong static correlation will be called single-reference systems).
Quantitative treatments of open-shell systems require one to describe not only static correlation but also dynamic correlation, which arises from electrons avoiding one another at short range to minimize repulsion or correlating their motion at large distances leading to dispersion and dispersion-like forces.
Organic molecules with open-shell singlet character are prime examples of strongly correlated systems. Oligoacenes (compounds consisting of several linearly fused benzene rings; see Fig. 1) have singlet ground states that develop increasing open-shell character with increasing length and are considered as prime test cases to study the performance of electronic structure methods for open-shell systems. [17][18][19][20][21][22][23][24][25][26][27][28][29][30] These open-shell singlets may be considered to be diradicals for oligoacenes of moderate length, and polyradicals for still longer ones, with about two unpaired (or signicantly partially unpaired) spins (one a, one b) for each ve or so rings. 17,27 In the last few decades, oligoacenes and their derivatives have also become of great interest for applications due to their charge transport properties, [31][32][33][34] complex excited state dynamics, [35][36][37][38] and electronic structure. 39 A potentially transformative process for improving the efficiency of photovoltaic devices is singlet ssion, in which a singlet exciton decays into two triplet excitons. [35][36][37] Long oligoacenes are especially promising candidates for singlet ssion due to their favourable crystal packing and electronic structure. 40,41 The efficiency of singlet ssion is strongly correlated with polyradical or diradical character of the ground state 42,43 and with the singlet-triplet (ST) gap of the molecule. 44 While short oligoacenes mostly have closed-shell ground states, longer acenes (beyond pentacene) develop open-shell character 17 and the ground state can in principle be either a singlet or a triplet, with both having multireference character. Predicting the ground spin state and the singlet-triplet energy splitting is essential to understanding the singlet ssion mechanism and to designing new materials.
An excellent overview of the electronic structure of oligoacenes and a review of the various theoretical treatments up to 2010 has been provided by Bettinger. 45 Angliker et al. extrapolated available experimental UV-vis data and predicted all oligoacenes beyond nonacene to have a triplet ground state. 39 This was supported by initial Kohn-Sham density functional theory (KS-DFT) calculations that considered only closed-shell singlet states, 46 but when open-shell singlet states were considered employing a broken-symmetry formalism, KS-DFT calculations predicted singlet diradical states as ground states for oligoacenes beyond pentacene. 47 Several electronic structure calculations including density matrix renormalization group (DMRG) calculations with a Pariser-Parr-Pople (PPP) Hamiltonian 48 and with an ab initio Hamiltonian, 17 two electron reduced density matrix (vRDM) calculations, 19 and coupled cluster calculations 20,21 also predicted all acenes up to dodecacene to have a singlet ground state, i.e., the diradicals and polyradicals are antiferromagnetically coupled. More recently, Ibeji and Ghosh employed spin-ip methods to compute the singlet-triplet gaps of oligoacenes and through extrapolation of their data, they have shown that there is no singlet-triplet crossover for innite chain acenes. 26 Diradicals and polyradicals are strongly correlated, and one of the reasons behind these differences in the prediction of the spin ground state of long oligoacenes is that there have been no electronic structure methods that can reliably treat systems that are both strongly correlated and large. Two candidate methods are KS-DFT and coupled cluster theory. The KS-DFT method is affordable for large systems, and it represents openshell singlets using broken-symmetry Slater determinants with electron spin component M S equal to 0 and having two (or more) unpaired electrons with opposite spins in different molecular orbitals. These broken-symmetry states are a mixture of a singlet and a triplet state; if one were able to use an exact exchange-correlation functional, KS-DFT could give the exact results for the singlet using a broken-symmetry solution, but available approximate functionals tend to give an energy between that of the singlet and triplet. Although several methods have been advanced to extract the singlet state energy using the broken-symmetry solution, [49][50][51][52] and although they are sometimes successful, they are not reliable. Coupled cluster theory can also in principle yield the correct singlet energy, but conventional coupled cluster theory uses a single-conguration reference wave function, and to make a treatment with such a reference state reliable for multireference systems oen requires including triple and higher excitations, which is impractical for medium-sized and larger systems. Another way to treat biradicals and polyradicals with lower excitation levels is to use spin-ip (SF) approaches 26,53,54 that take single-determinantal high spin states as the reference state.
An alternative with greater prospects for success in treating multireference systems is a method that adds dynamical correlation to a multiconguration self-consistent eld (MCSCF) wave function 55 that has the correct spin symmetry and is variationally optimized for both the orbitals used in the conguration state functions (CSFs) and the coefficients of the CSFs in a conguration interaction expansion of the wave function. There are several possible ways to choose the congurations in the CI expansion to try to balance accuracy (which can demand a large number of congurations) and affordability (which can demand a smaller number). Many of these are best explained by rst considering the complete active space self-consistent eld (CASSCF) method, 56 which is a special case of MCSCF in which a full conguration interaction expansion (i.e., one including all possible ways to assign the electrons to orbitals) of the wave function is constructed over a specied active space of n electrons and m orbitals with all the other orbitals either doubly occupied or vacant. However, the number of CSFs in this full-CI-within-awindow wave function increases exponentially with the activespace size and this approach already reaches its practical limit for closed-shell singlets when n ¼ 18 and m ¼ 18 and for openshell states with similar n and m (depending on the state). This means that CASSCF calculations including all oligoacene valence-p orbitals in the active space are possible for at most up to tetracene. Therefore methods like the generalized active space SCF (GASSCF), 57,58 the restricted active space SCF (RASSCF), 59-62 the occupation restricted multiple active space (ORMAS), 63 and split-GAS 64 have been developed to remove many of the less important CSFs from the full CI. For example, in the GASSCF model, the active orbitals are placed into distinct subspaces, and accumulated minimum and maximum electron occupation numbers are applied to each subspace. Within a subspace, one includes all possible spin-and symmetry-adapted CSFs that can be constructed with this occupation number constraint plus, optionally, a restricted set of intersubspace excitations. This allows larger active spaces (larger n and m) than one can treat with CASSCF, while keeping the numbers of CFSs affordable.
Since an MCSCF calculation includes only a small fraction of the dynamic correlation, one must include the remaining dynamic correlation, which is necessary for chemical accuracy, by a post-SCF step. Until recently, the most affordable method that uses an MCSCF wave function as a reference for adding dynamic correlation in a post-SCF step has been complete active space second order perturbation theory (CASPT2) 65,66 or restricted active space second order perturbation theory (RASPT2). 61 Recently the capability has been developed to also carry out GASPT2 (ref. 67) calculations, although with the approximation that excitations into the orthogonal complement of the GAS space in the CAS space are omitted from the perturbed wave function. For this reason, the GASPT2 results presented in the following have only a subset of the second-order terms and hence are incomplete. However, it is still useful to report the GASPT2 results, because they are the only affordable calculations that one can perform at present using a perturbation treatment on top of a GASSCF wave function. Although these PT2 methods are capable of providing high accuracy, 68 they are not suitable for large systems due to the rapid increase of computational cost and memory requirements.
As an alternative, we have recently developed a method called multiconguration pair-density functional theory (MC-PDFT). 69 It is a new kind of density functional theory in which the electronic energy is calculated from the kinetic energy, density, and on-top pair density of a multicongurational wave function. The MC-PDFT method requires a reference calculation to generate an on-top pair density, which is the probability of nding two electrons at a given point in space. Whereas the density functional in KS-DFT depends on the spin-up electron density and the spin-down electron density and is called the exchangecorrelation functional, the density functional in MC-PDFT depends on the total electron density and on-top pair density and is called the on-top density functional. MC-PDFT has been successfully used to compute excitation energies, 70-72 barrier heights, 73 and transition metal energetics. 74 In these prior calculations, we have employed several different on-top density functionals that have usually given similar results. In the present article, we use only one of these, namely tPBE, which was dened in our rst MC-PDFT paper. 69 In general, any method that generates a two-body density matrix can be used as an MC-PDFT reference, although so far applications have been limited to using CASSCF and GASSCF wave functions, for which one may view MC-PDFT as a post-MCSCF method. For calculations based on GASSCF, a special case of GASSCF called the separated pair (SP) approximation has been dened and proven to be successful. 75 SP is the special case in which no more than two orbitals are included in any GAS subspace and in which intersubspace excitations are excluded. The SP approximation typically leads to far fewer congurations than CASSCF for the same n and m, and it has been successful for ground-state calculations. For excitation energies in the present paper, we explore alternative GASSCF partitions based on frontier orbitals.
MC-PDFT using CASSCF or GASSCF reference wave functions has been tested only for small systems. In the present study, we employed MC-PDFT to compute the singlet-triplet energy difference for oligoacenes ranging in length from naphthalene to dodecacene. A reasonable choice of active space for a p-conjugated system is an unabridged valence-p active space, with each occupied orbital having a correlating orbital, but a CASSCF with this active space reaches its limit at tetracene (18 electrons in 18 orbitals). On the other hand, within the GASSCF formalism, active spaces as large as 50 electrons in 50 orbitals can be employed. We here test the performance of MC-PDFT for various GASSCF wave functions with the aim of providing a general prescription on how to perform MC-PDFT calculations on large p-conjugated systems.

Computational methods
The geometries of the oligoacenes for both the singlet and triplet states were optimized using the B3LYP 76 exchangecorrelation functional and the 6-31G(d,p) basis set. For comparison with the present MC-PDFT calculations, we also carried out KS-DFT calculations with the 6-31+G(d,p) basis set with the PBE0 (ref. 77) exchange-correlation functional as well as with the PBE, 78 one on which both PBE0 and tPBE are based. All KS-DFT calculations were performed with the unrestricted broken-symmetry approach (labeled "variational" in previous work 52 ) using Gaussian 09. 79 All CASPT2 and GASPT2 calculations were performed using the standard empirical IPEA 80 shi value of 6.80 eV (0.25 a.u.) and an imaginary shi 81 of 5.44 eV (0.2 a.u.). All MC-PDFT calculations were performed with the tPBE 69 on-top functional. We have also performed PBE, 74 tBLYP 69 and BLYP 74 calculations for hexacene to see the functional dependence of the MC-PDFT results. Results with these other functionals are reported in the ESI. † The 6-31+G(d,p) basis set was used for the CASSCF, GASSCF, CASPT2, GASPT2, and tPBE calculations. For all these calculations, Cholesky decomposition via resolution of the identity, 82 was used to facilitate the computation of the twoelectron integrals. All CASSCF, GASSCF, CASPT2, GASPT2, and tPBE calculations were performed in a locally modied version of Molcas 8.1. 83 Oligoacenes ranging from naphthalene (n-acene with n ¼ 2) to dodecacene (n-acene with n ¼ 12) were investigated. All calculations were performed by imposing D 2h symmetry. In this point group, all singlet states considered belong to the 1 A g irreducible representation, and the triplet states belong to the 3 B 3u irreducible representation. Adiabatic singlet-triplet energy gaps were determined as electronic energy differences between the energies calculated at the respective optimized geometries of the singlet ground state and the triplet state, whereas the vertical singlet-triplet gaps are calculated at the optimized singlet geometries. To calculate singlet-triplet gaps for innite chain acenes for tPBE, we have tted our data to an exponential decay of the form a + b exp(Àcx) as done in ref. 26.
Several CASSCF and GASSCF active space choices were tested, as discussed next.

CASSCF active space choices
For the CASSCF calculations, we use the notation CAS(n,m) where n is the number of active electrons, and m is the number of orbitals in the active space. In the paper, we report the performance of MC-PDFT and CASPT2 for a minimal active space, i.e., CAS(2,2). Larger active spaces, in particular CAS (4,4) and CAS (8,8) were also considered; all results obtained for the larger complete active spaces are reported in the ESI. †

GASSCF active spaces
All GASSCF active spaces include all valence p orbitals in the active space, but four GASSCF partitions are examined in which different choices are made for the accumulated minimum and maximum electron occupation numbers (Fig. 2).
Frontier partition with one set of interspace excitations (FP-1). In this partition, GAS1 includes the HOMO and LUMO orbitals, GAS2 includes all other occupied p orbitals, and GAS3 includes all unoccupied p* orbitals. Single and double excitations are allowed within GAS1. Moreover, single and double excitations are allowed from GAS2 to GAS3, and no excitations are allowed in or out of GAS1.
Double-frontier partition with one set of interspace excitations (DFP-1). In this partition, GAS1 includes the HOMO and LUMO orbitals, GAS2 includes the HOMOÀ1 and LUMO+1 orbitals, GAS3 includes all the other occupied p orbitals, and GAS4 includes all the remaining unoccupied p* orbitals. In this case single and double excitations are allowed within GAS1 and within GAS2; no intersubspace excitations are allowed for GAS1 and GAS2; single and double excitations are allowed from GAS3 to GAS4.
Wide-frontier partition with one set of interspace excitations (WFP-1). In this partition, GAS1 includes the HOMO, LUMO, HOMOÀ1, and LUMO+1 orbitals, GAS2 includes all the other occupied p orbitals and GAS3 include all the other unoccupied p* orbitals. Thus GAS1 has four orbitals, and GAS2 and GAS3 each have (n À 4)/2 orbitals, where n is the number of valence p orbitals (n equals the number of carbon atoms in the oligoacene). All excitations are allowed within GAS1; single and double excitations are allowed from GAS2 to GAS3 and no intersubspace excitations are allowed for GAS1.
Wide-frontier partition with three sets of interspace excitations (WFP-3). In this partition, labelling of GAS1 and GAS2 are switched with respect to the WFP-1 case, and GAS3 remains the same, but now single and double excitations are allowed from GAS1 to GAS2, GAS1 to GAS3, and GAS2 to GAS3. The minimum and maximum number of electrons in GAS1 are n À 6 and n À 4 respectively; and the minimum and maximum numbers of electrons in GAS1 + GAS2 are n À 2 and n respectively (so the maximum excitations into GAS3 is 2). All excitations are allowed within GAS2.

Highest-level available literature estimates and asymptotic limit
Literature values for vertical and adiabatic singlet-triplet gaps in oligoacenes from previous electronic structure calculations are reported in Table 1. The CCSD(T) values are obtained starting from restricted Hartree-Fock reference wave functions. 21 Such CCSD(T) calculations become less reliable as oligoacene size increases because of the increasingly strongly correlated character. The DMRG values 17 are based on complete active space conguration interaction without orbital optimization and with small basis sets, and for hexacene and smaller acenes they include three p* orbitals for each p orbital (for octacene and larger acenes they use the same active space as used here, i.e., one p* orbital for each p orbital, but with the very small STO-3G basis set, and these results are not used for comparison here). The vRDM results 19 are based on an active space containing all valence p electrons (as in some of the present calculationssee below). But neither DMRG nor vRDM attempt to recover all the dynamic correlation energy, which makes the accuracy of these results questionable. The multireference Moller-Plesset (MRMP) results 84 for anthracene and tetracene are based on active spaces with 12 electrons in 12 orbitals where six p and six p* orbitals are included. Even though MRMP is capable of recovering dynamic correlation energy, such calculations with balanced active spaces are not affordable for larger acenes. Multireference conguration interaction singles and doubles with correction for quadruple and higher excitations (MRCISD + Q) 85 results are obtained using reference states obtained from complete active space calculations using an active space of 8 electrons in 8 orbitals. The +Q correction approximately corrects for the underestimate of dynamic correlation energy by MRCISD as it only includes single and double excitations. We take the results in the vertical gap column and the average adiabatic gap column as highestlevel available literature estimates, and in the following we report average absolute deviations from these values as mean unsigned deviations (MUDs).
In ref. 26, the authors t the best estimates of the adiabatic singlet-triplet gaps to an exponential decay of the form a + b exp(Àcx), from which the extrapolated singlet-triplet gaps for innitely long polyacenes are 5.06 kcal mol À1 for the best estimates and 5.37 kcal mol À1 for their best spin-ip results. We will compare our long chain singlet-triplet gap values to the value of 5.06 kcal mol À1 .

Results and discussion
We rst discuss the CASPT2 and tPBE results using the CAS(2,2) active space (Table 2; for additional CAS(4,4) and CAS (8,8) results, see Tables S1 and S2, † respectively). Both CASPT2 and tPBE predict a singlet ground state for all oligoacenes. From nonacene to dodecacene, the vertical singlet-triplet gap is lower than the adiabatic one, which cannot occur physically but it must be recalled that these energy differences are all computed using xed B3LYP geometries (to facilitate comparison with prior predictions), i.e., these are not vertical excitation energies for each different theory's ground-state geometry. This inversion of the energy ordering is observed for other active spaces as well, and it would be interesting in future work to optimize all geometries at the CASSCF, CASPT2, and tPBE levels; at present, however, this is not technically feasible. In any case, the reported MUD is meaningful since all vertical excitation values in the table (both our new calculations and the literature values) are computed for the same geometries. Table 2 shows that both CASPT2 and tPBE with the (2,2) CASSCF active space agree remarkably well with the literature values up to decacene; however, tPBE predicts increasingly small gaps for undecacene and dodecacene, which is not in agreement with the best available literature values. Extrapolation of the CASPT2 and tPBE results give limiting singlet-triplet gap of 7.1 kcal mol À1 and 1.3 kcal mol À1 respectively for an innite acene compared to the value of 5.1 kcal mol À1 noted above as a previous best estimate. We note that prior studies have indicated that it is important to add other p orbitals into the active space, 22 especially for longer oligoacenes, to fully account for diradical and ultimately polyradical character. A systematic way to do this is to include all occupied and unoccupied valence p orbitals in the active space. However, such a large active space is not affordable beyond tetracene in the CASSCF formalism; instead, we do this in the present work with the GASSCF frontier partitions dened above. Table 3 reports the number of CSFs for the GASSCF and corresponding CASSCF wave functions, including all valence p orbitals; this table shows a tremendous  Complete results for the other CASSCF active spaces, for the other two partitions of the valence p GASSCF spaces, and for the KS-DFT calculations are given in the electronic (ESI †), and mean unsigned deviations for the CASSCF(2,2) active space, for all four frontier partitions of the valence p GASSCF spaces and for KS-DFT will be compared below. Table 4 shows that tPBE with the FP-1 partition agrees extremely well with literature values; in particular, the MUDs are respectively 1.5 kcal mol À1 and 1.6 kcal mol À1 for vertical and adiabatic singlet-triplet gaps. It is also encouraging to notice that tPBE predicts an asymptotic value of 7.6 kcal mol À1 compared to the reference value of 5.1 kcal mol À1 . The incomplete GASPT2 calculations and the tPBE calculations predict very similar results for naphthalene, but they differ more for anthracene and tetracene. Table 5 shows that for WFP-3 calculations, which include more congurations in the CI expansion because single and double excitations between GAS spaces are allowed, tPBE again agrees reasonably well with the literature values, with MUDs of 5.7 kcal mol À1 and 4.1 kcal mol À1 for vertical and adiabatic singlet-triplet gaps, respectively. These results are a little worse than with the smaller CSF list of FP-1 when compared to the literature values. WFP-3 predicts the asymptotic gap of 1.9 kcal mol À1 which is narrower than FP-1. Note that GASPT2 calculations could only be performed up to tetracene with the WFP-3 partition because of unaffordable resource demand (maximum available runtime of 240 hours) for the larger systems. Nevertheless, it is interesting that tPBE and GASPT2 agree quite well with each other for the rst three systems. Note that GASSCF    With wave function methods, one would expectall other things being equalmore accurate results with a larger CI expansion in the active space, but it is well known that this is not always the case because it is equally important or more important that the active space includes the qualitatively important near-degeneracy effects in a well-balanced way as that it is large. The same balance is important in MC-PDFT, and the results with the larger conguration lists are not always more accurate. With that caveat, we note that for both widefrontier partitions (WFP-1 and WFP-3), tPBE systematically predicts smaller vertical and adiabatic singlet-triplet gaps for longer acenes with respect to the literature values. This is not surprising, though, since the DMRG and v2RDM literature values were all obtained without taking into account the total dynamic correlation. In our calculations, dynamic correlation is recovered using the on-top functional of pair-density functional theory. It has been predicted previously that the inclusion of dynamic correlation effects would decrease the singlet-triplet gaps by a few kcal mol À1 . 17 Thus, while the reported MUD for the WFP-3 partition is larger than that for the FP-1 partition, the predictions may be more accurate for the former and may indeed constitute the set of reference values against which other models should be compared. From a chemical standpoint, it is particularly signicant that neither the tPBE calculations with the (2,2) CASSCF space nor the tPBE calculations with any of the frontier-orbital partitions of active spaces including all p orbitals show a singlet-triplet crossover up to dodecacene. Table S8 in the ESI † shows results with other on-top functionals, and it is encouraging that Table S8 † shows that the results also do not depend strongly on the choice of on-top functional. Table 6 compares MC-PDFT mean unsigned deviations for CASSCF(2,2) reference wave functions and all four GASSCF partitions of the valence p active space. It is encouraging that the results do not depend strongly on the MCSCF wave function. It is interesting that the results are even slightly better on average for the frontier partitions involving smaller CSF lists, especially when we recall the caveats associated with the dependence of the results on the size of the conguration list and with the unknown accuracy of the literature values.
The percentage contribution of the most dominant conguration (the nominal HF conguration) for the singlet state decreases along the series (Fig. 3) suggesting the presence of signicant multicongurational character for longer oligoacenes. The oscillations in Fig. 3 are reminiscent of that for the fundamental excitation gap of oligoacenes, that Korytár et al. 86 have attributed to a Dirac-like cone in the band structure. However, the smooth decrease in the HF conguration weight observed with the more complete WFP-3 suggests that the observed oscillations with fewer CSFs may alternatively be an artefact of a small active space.
Several studies have discussed the diradical and polyradical character of long oligoacenes. 17,19,22,27 This character can be analyzed in terms of the occupation numbers of the GASSCF natural orbitals (NOs). Fig. 4 shows the occupation number of the HONOÀ1, HONO, LUNO, and LUNO+1 for different oligoacenes for the FP-1 and WFP-3 partitions of the active space. Hereif n denotes the number of valence p electrons -HONO denotes the NO with occupation number n/2 when the orbitals are ordered by decreasing occupation number, and LUNO denotes orbital (n/2) + 1; HONOÀ1 denotes orbital (n/2) À 1, and LUNO+1 denotes orbital (n/2) + 2. The LUNO occupation number usually increases along the series, but again the trend is not monotonic. No symmetry switch occurs between the HONO and LUNO up to dodecacene, contrary to the observations in ref. 19 and 27. Interestingly for WFP-3, the occupation number of LUNO+1 also tends to increase, although at a slower rate than the LUNO occupation. An increasing occupation of the LUNO+1 is a signature of growing polyradical character. 17 Again, signicantly reduced oscillatory behaviour with WFP-3 is possibly attributable to the more complete nature of this active space. The adiabatic singlet-triplet gaps obtained for tPBE using CAS(2,2), FP-1, and WFP-3 active spaces, and the literature values are compared in Fig. 5. With both CAS(2,2) and FP-1 reference wave functions, tPBE results match reasonably well with the literature values, while with WFP-3 reference wave functions, tPBE predicts systematically smaller adiabatic singlet-triplet gaps than the literature values, which may reect better accounting for the differential effects of correlation on the two state energies.

Concluding remarks
In the present paper, we used the recently developed MC-PDFT method to compute the singlet-triplet gaps in oligoacenes ranging from naphthalene to dodecacene. We showed that the tPBE ontop density functional performs well with both minimal and valence-p active spaces. Using generalized active space reference wave functions, we were able to perform calculations with orbitally optimized active spaces up to 50 electrons and 50 orbitals and with dynamic correlation included by an on-top density functional. Even though it is not easy to generalize prescriptions for partitioning the active space, because the best partitions are expected to be system-and problem-dependent, the various frontier partitions used here may be useful for other similar problems as well as for oligoacenes. Our rst conclusion is that the results are moderately sensitive to active space choice. Our second conclusion is that even small active spaces can yield good results and active spaces with sparse GAS partitioning can yield useful accuracy. However, it is not possible to draw further conclusions at this stage because, in the absence of experimental measurements, we have no way to assess the accuracy of the various high-level available literature values, and thus it is not possible to determine which active space provides the most accurate results. To amplify on that point, it may well be that one of the computational protocols employed herein is the most accurate of all methods to date, but a validation of that possible hypothesis will require measurement or demonstrably more complete calculations.
The GASSCF method has great potential to generate physically meaningful wave functions for large systems, and MC-PDFT based on these wave functions can provide chemical accuracy. Overall this study provides a framework in terms of both wave function theory and density functional theory that can be extended to study ground-and excited-state properties of organic electronic materials and other problems involving large, strongly correlated, conjugated systems.