Density functionals with asymptotic-potential corrections are required for the simulation of spectroscopic properties of materials

Five effects of correction of the asymptotic potential error in density functionals are identified that significantly improve calculated properties of molecular excited states involving charge-transfer character. Newly developed materials-science computational methods are used to demonstrate how these effects manifest in materials spectroscopy. Connection is made considering chlorophyll-a as a paradigm for molecular spectroscopy, 22 iconic materials as paradigms for 3D materials spectroscopy, and the VN− defect in hexagonal boron nitride as an example of the spectroscopy of defects in 2D materials pertaining to nanophotonics. Defects can equally be thought of as being “molecular” and “materials” in nature and hence bridge the relms of molecular and materials spectroscopies. It is concluded that the density functional HSE06, currently considered as the standard for accurate calculations of materials spectroscopy, should be replaced, in most instances, by the computationally similar but asymptotically corrected CAM-B3LYP functional, with some specific functionals for materials-use only providing further improvements.


Introduction
The asymptotic potential reects the energy required to take an electron from a system of interest and remove it to the surrounding vacuum. The most commonly used computational method for evaluating spectroscopic properties of molecules and materials, density-functional theory (DFT), can be implemented using functionals embodying many different levels of theory, from the local-density approximation (LDA), to generalised gradient approximations (GGA), to hybrid functionals that mix long-range Hartree-Fock exchange with local exchange, to range-corrected hybrid functionals designed to realistically represent the asymptotic potential, and beyond. This work focuses on the critical effects caused by the difference between basic LDA, GGA, meta-GGA, and hybrid functionals and functions with asymptotic correction: their ability to reliably calculate properties of the excited states of molecules and materials. If a spectroscopic transition involves electron transfer from one location to another, then the asymptotic potential error directly manifests in the calculated transition energies.
Two decades ago, the critical role of the asymptotic potential was recognised in regard to the evaluation of molecular spectroscopic properties. 1-23 A signicant benet from this was a gained ability to simulate excitation energies and exciton couplings in photosynthetic reaction centres. [24][25][26] Even for transitions not involving signicant electron transfer, the asymptotic potential error can induce large errors in properties such as calculated reorganization energies and their partitioning into the Huang-Rhys factors (electron-vibration coupling constants) that control spectral bandshapes as well as photochemical and photophysical reaction rates. [27][28][29] This occurs as such properties are controlled by the details of the entire excited-state manifold, not just the excited state of primary interest.
That such effects are also important for materials was recognised early on, 30 but functionals involving asymptotic correction were not implemented into appropriate soware packages. This situation changed recently with the implementation of the CAM-B3LYP asymptotically-corrected functional into both the Car-Parrinello Molecular Dynamics (CPMD) package 31 (gamma-point only) and the Vienna Ab initio Simulation Package (VASP) 32 (full implementation). In parallel, new density functionals appropriate only to materials are being developed that link the asymptotic-potential correction to the evaluation of the dielectric constant. [33][34][35] These methods have greatly improved the modelling of dielectric properties, as well as providing for improved predictions of band gaps. Such approaches lack the general applicability of established methods such as CAM-B3LYP, however, and involve extensive choices made between method options and possibly also their parameterisation. More broadly, many asymptotically corrected functionals have now been developed, including recent efforts, 23,36-44 reecting many and varied design choices. Whilst we consider mostly CAM-B3LYP herein, the elucidated basic principles are expected to be qualitatively representative of all asymptotically corrected functionals.
This work focuses on the basic understanding as to what asymptotic-potential correction does and why it is essential for studies on both molecules and materials. To do this, we focus on a system that can be regarded as being either a "molecule" or a "material": the defect spectroscopy of hexagon boron nitride (h-BN). This example is taken from the nanophotonics research eld, [45][46][47][48][49] an area that embeds into solid-state environments spatially localized properties, with applications ranging from quantum networks and quantum information to spin-photon interfaces. [49][50][51][52] Indeed, defects in h-BN have been associated with single-photon emission (SPE) [53][54][55][56] and optically detected magnetic resonance (ODMR). 57,58 Historically, only defects displaying magnetic properties have had their chemical natures determined, 57,59-65 but recently came the rst characterization based only on observed spectral properties. 66 This was made possible through extensive experimental characterization of composition, combined with computational spectroscopic predictions. In general, the role of computation has been important in all h-BN defect assignments. 61 For the calculation of the spectroscopic properties of defects in 2D materials, many difficult issues arise with both DFT and ab initio wavefunction approaches, with no method that is currently practical able to predict transition energies to within the desired accuracy of say AE 0.2 eV for all possible scenarios. 67 Of tested methods, the most generally reliable approach has been CAM-B3LYP, 7-9 which gives results within 0.5 eV of those from computationally demanding ab initio methods. 67 For both 2D and 3D materials, one of the most widely used density functionals is HSE06, 68,69 a method that also embodies range separation, 3 but uses this feature instead to enhance computational efficiency by removing the Hartree-Fock contributions at long range. 69 Hence a key deciency of the original PBE functional, 70 upon which HSE06 is based, is made more prominent. Currently, HSE06 is recognised as the standard for the prediction of band gaps of materials. 71 Nevertheless, for the case of defects in 2D materials, HSE06 has been established to make predictions that deviate signicantly from those of CAM-B3LYP and ab initio approaches 67,72 for problems that do not overtly involve charge transfer.
The example system considered herein is the V N À defect in a h-BN 2D layer. This defect consists of a nitrogen-atom vacancy that is negatively charged; it has been considered for a long time 73,74 as a possible contributor to observed h-BN spectroscopic properties, and of late also as a candidate for explaining some observed 58 ODMR, but always not all calculated and observed properties appear to match. Of interest herein, the orbitals associated with the defect are predicted to lie close to the h-BN conduction band, 64,74 with the result that chargetransfer transitions with low energies will occur. This defect therefore allows for the probing of the effect of the asymptoticpotential error on the properties of materials, in a way that directly links to known analogous molecular properties. First, Results Section (a) reviews the effects of the asymptotic potential on the perceived spectroscopy of an iconic molecular system, chlorophyll-a, 9,10 identifying 5 key effects. Section (b) then reviews how these effects manifest concerning predictions of the lowest-energy transition (bandgap) of simple materials, 32 materials historically perceived 71 as being well-treated without the need for asymptotic-potential correction. Comparison of predicted properties from modern materials-only functionals with asymptotic-potential correction to those of generally applicable functionals such as CAM-B3LYP are also reviewed therein, with the results contrasted to those of hybrid functionals including HSE06, as well as some other GGA and meta-GGA approaches.
This understanding is then used to interpret calculated properties of the spectroscopic manifold of V N À . We consider in Sections (c) and (d) molecular models of the defect (Fig. 1), and then 2D periodic models in Section (e) (Fig. 1). Calculations on model compounds for defects in 2D materials show very rapid convergence of calculated electronic properties with respect to increasing sample size for transitions localised within the defect. 65,66,72 Of signicance, model-compound calculations have also been shown to converge to the same results as obtained using analogous calculations on 2D periodic defect models. 72,75 Nevertheless, these generic results are not expected to apply to the charge-transfer transitions of V N À considered herein. In 3D materials, dielectric screening of long-range interactions can be critical to spectroscopic properties, but this effect does not apply here to the properties of 2D materials. 76 This is a critical factor leading to the found 65,66,72 rapid convergence of both molecular-model and 2D-model spectroscopic calculations with increasing sample size.

Results
(a) Overview of how corrections to the long-range potential can become essential for the qualitative interpretation of molecular spectroscopy Here, ways in which errors in the asymptotic potential manifest to introduce serious qualitative errors in molecular spectroscopy are reviewed, considering the iconic example of chlorophyll-a, arguably the world's most important chromophore. The major band systems of this and related molecules (porphyrins, hemes, etc.) were characterised by Gouterman in the 1960's, 77,78 in energy order, as: two $ red transitions, Q y then Q x , two $ blue transitions B x then B y , then a series of weak transitions starting at N x and N y . In these labels, x and y indicate absorption polarisation directions in the macrocycle plane. Ab initio spectral calculations verify this picture, no matter what approach is taken. 9,79,80 The (weak) Q x band proved difficult to assign as two x-polarised peaks are observed instead of the expected one, leading to two possible assignments that were debated continuously from 1982 (ref. 81 and 82) until 2013. 26 The resolution of this issue is that the Born-Oppenheimer approximation 83 fails for the Q bands, meaning that no classical spectroscopic model (e.g., using Huang-Rhys factors, 84 more general Franck-Condon factors, 85 and/or Herzberg-Teller interactions 86 ) can describe the observed spectroscopy: Q x absorption displays two maxima with a local minimum at the band centre. During the period in which the assignment of the spectrum of chlorophyll-a was intensely debated, DFT calculations depicted a scenario in which one of the two x-polarised bands was assigned to Q x and the other to N x , a state unconventionally predicted to occur at low energy. 87,88 The N bands are chargetransfer bands 9 that signicantly alter the molecular dipole, and hence are subject to the asymptotic-potential error in DFT. The error was established to be general and result in the underestimation of charge-transfer band energies by several eV. 10 The application of functionals such as CAM-B3LYP with asymptotic-potential correction resulted in predictions of N-band energies consistent with ab initio predictions and historical experimental interpretations. 9 Table 1 gives details of the CAM-B3LYP calculations for chlorophyll-a, tracing the effect of long-range correction by comparing these results to analogous B3LYP 89,90 ones and thence to HSE06 results. Highlighted are the differences between calculated and observed values, hampered somewhat as only bounds are available for the N bands and also for one other charge-transfer band involving peripheral-ring ketone excitation. Results for the popular uB97xD asymptotically corrected density functional 6 are also included and are similar to the CAM-B3LYP results. From the table, it is clear that functions without asymptotic correction like HSE06 and B3LYP give results that are signicantly different from those that do.
From Table 1, four key effects of asymptotic-potential correction are extracted and listed in Table 2, along with one additional effect. (1) Asymptotic correction signicant increases the bandgap between the HOMO and LUMO orbitals, with (2) a corresponding increase in the exciton binding energy of localised spectroscopic transitions. These two effects largely cancel so that (3) asymptotic correction does not greatly perturb localised transitions. Alternatively, for transitions involving charge transfer (4), the exciton binding energies are enhanced much less and hence charge-transfer transitions increase in energy. This creates order within the excited-state energy manifold. As a consequence  a Vertical excitation energies at the CAM-B3LYP/6-31G*/D3(BJ) groundstate geometry. 91 b Gas-phase, 92 with the Q x band centre (an absorption minimum) estimated using established methods: 26 (5), the vibronic couplings that occur between excited states become much better described, and hence the excited-state potential-energy surfaces become better described. This means that excited-state equilibrium geometries and spectroscopic reorganisation energies become better described, as well as vibrational line intensities arising from both Franck-Condon (Huang-Rhys) and Herzberg-Teller mechanisms. 27,28,91 (b) Overview of how corrections to the long-range potential can signicantly improve calculated properties of simple materials Concerning the spectroscopic properties of some simple materials, Table 3 provides a broad summary of the effects of asymptotic-potential correction, comparing known results for up to 22 materials for 6 functionals without asymptotic correction to those for 11 functionals with asymptotic correction. Some new calculations are performed for CAM-B3LYP, as detailed in Methods and ESI. † All results are listed in ESI Tables  S1 and S2, † with Table 3 providing a summary in terms of meanabsolute errors (MAD) and maximum errors (MAX). Firstly, Table 3 shows the MAD error in calculated lattice constants for the methods for which data is available. As previously noted, 32 the error for CAM-B3LYP is much less than that for HSE06, and indeed much less than that for B3LYP, indicating the importance of asymptotic correction. When predicting spectroscopic properties of materials, it may be necessary to optimise structures and hence the apparent reliability of CAM-B3LYP in this regard is of note.
In addition, Table 3 shows MAD and MAX errors for orbital bandgaps. This involves the comparison of calculated bandgaps with those deduced from experiment through estimations of the exciton binding energy and the assumption that the lowestobserved electronic transition corresponds to the HOMO to LUMO transition. The table is divided into methods with and without asymptotic correction: those without it, including HSE06, show large MAD and MAX errors, whilst those with it, including CAM-B3LYP, perform considerably better. This effect is demonstrated most clearly by comparison of the poor results for PBE0 (ref. 33) and B3LYP to the signicantly improved results for their asymptotically corrected counterparts, PBE0(a) 41 and CAM-B3LYP.
There is current interest in the development of meta-GGA functionals for the prediction of spectroscopic properties of materials. Available results for the 22 materials considered are Table 2 Known effects of asymptotic correction on calculated molecular spectroscopic properties Effect Description (1) Signicant increase in the HOMO-LUMO bandgap (2) Signicant increase in exciton binding energies for chargelocalised transitions, cancelling the increase in the HOMO-LUMO gap (3) Effects (1) and (2) combine to produce a signicant increase in the relative energy of charge-transfer transitions (4) Reduction in the differential error in state energies to provide an improved description of the excited-state manifold (5) Improved description of the shapes of potential energy surfaces, and hence reorganisation energies, Franck-Condon (Huang-Rhys) factors, Herzberg-Teller couplings, etc., as vibronic couplings with the excited-state manifold are much better represented  Table 3 were also designed considering the analysed data, whereas CAM-B3LYP is an old functional for which no materials properties were used in its construction. Methods like CAM-B3LYP are designed to reproduce directly observable quantities such as oxidation or reduction potentials and spectroscopic transition energies. Such quantities depend on orbital bandgaps, but include Coulomb and exchange interactions involving the associated charges as well. How this works for molecules was described in Table 2, with the results for materials from Table 3 and ESI Table S2 † clearly reecting effect (1): HOMO-LUMO bandgaps are signicantly increased in going from B3LYP to CAM-B3LYP or from PBE0 to BPE0(a).
Effect (2) depicts analogous changes to exciton binding energies, the prediction of which in materials CAM-B3LYP has been reported to be ve times more reliable than HSE06, 32 making effect (3) also important. Nevertheless, in materials such as silicon, no charge transfer is involved. Hence effect (2) has only a minor effect that is accurately predicted by CAM-B3LYP, 32 with the consequence that CAM-B3LYP overestimates the orbital bandgap and spectroscopic transition energy by 1.1 eV (Table 3 and ESI Table S2 †). Silicon is an unusual and (therefore) important material that has spawned modern generations of asymptotic corrected functionals ( Table 3) that focus on the correct description of its dielectric properties. 33,34,41,43,44 A contrasting example in ESI Table S2 † is LiF, for which both effects (2) and (3)  Prior to that, 122 the deduced exciton-binding energy was 1.00 AE 0.06 eV, leading to a perceived orbital bandgap of 13.60 AE 0.06 eV. Hence the CAM-B3LYP calculated exciton-binding energy is also in good agreement with the experimental data. In stark contrast, the exciton binding energy calculated by HSE06 is just 0.06 eV, too small by a factor of over 20. Lithium participates strongly in p back bonding and so charge-transfer is expected to be a signicant aspect of its properties, making essential its treatment by asymptotically corrected methods.
In a third example of note from ESI Table S2, † the orbital bandgap of a-Al 2 O 3 was initially estimated spectroscopically 123 to be at 8.8 eV, with improved analyses 106 leading to a revised value of 9.57 eV that is in agreement with conductivity measurements, 107 and an exciton binding energy of 0.13 eV. These later results are in good agreement with CAM-B3LYP predictions of an orbital bandgap of 9.42 eV and exciton binding energy of 0.20 eV. Alternatively, results for expensive state-of-the-art methods such as Bethe-Salpeter 124 (0.40 eV) and another GW-based approach 125 (0.37 eV) are signicantly too large, whilst HSE06 predicts just 0.010 eV.
Summarising these three examples, signicant differences between experiment and CAM-B3LYP calculated transition energies or orbital bandgaps can arise either (i) owing to shortcomings in the methodology, demanding use of improved modern long-range-corrected methods, 33,34,41,43,44 or (ii), when charge transfer is signicant, from shortcomings in experimental data interpretations. Owing to effect (2), asymptotically corrected methods, exampled here by CAM-B3LYP, are appropriate for the evaluation of exciton-binding energies. For LiF and a-Al 2 O 3 , CAM-B3LYP results agree to experiment to within ca. the experimental uncertainty, as has already been noted 32 for diamond, silicon, NaCl, MgO, 2D h-BN, and 3D h-BN. CAM-B3LYP therefore appears suitable for use in experimental data assignment, whereas asymptotically uncorrected methods, exampled here by HSE06, qualitatively fail to depict key spectroscopic quantities.

(c) The asymptotic-potential error in 2D defects perceived using TD-DFT and EOM-CCSD calculations on model defect clusters
The previous subsection concerned the properties of the lowestenergy transitions of 3D materials, but oen the questions of interest concern the assignment of many observed transitions within a manifold and hence the relative ordering of excited states becomes the most critical issue. Here, the relative ordering of transitions in model compounds depicting V N -À defect-associated transitions in h-BN are considered. Table 4 provides excited-state energies obtained using timedependent DFT (TD-DFT) 126 for the model compounds 1-4. TD-DFT is particularly well-suited to the study of defects as it only requires that its reference state be of mostly closed-shell nature, whereas most defect states are open shell, 67 which indeed is the case for V N À . These energies are shown diagrammatically in Fig. 2(a and b). Model 1 only supports excitations within the defect core, with larger models adding in additional excitations Table 4 Calculated vertical excitation energies for the V N À defect in h-BN, in eV, from the (1) 1 A 1 0 ground state of the model compounds 1-4 ( Fig. 1) that involve charge transfer to the h-BN conduction band. TD-CAM-B3LYP predicts that most transitions decrease in energy as the model size increases, with convergence quickly established. In contrast, TD-HSE06 predicts that most states continually decrease in energy as ring-size increases. The exceptions to this are the energies of (1) 1 E 0 and (1) 3 E 0 , which quickly converge. The electron densities of the key a 2 00 and e 0 orbitals involved in these transitions are shown in Fig. 2(d) and are all localized within the defect core. All other transitions, however, involve excitation to conduction-band orbitals, like the a 1 00 orbital shown in Fig. 2(d) that becomes occupied in the (1) 1 A 2 00 and (1) 3 A 2 00 states. Overall, these results parallel effects (2)-(4) that depict the impact of asymptotic-potential corrections on molecular spectroscopy, see Table 2.
The (1) 1 A 2 00 and (1) 3 A 2 00 states are highlighted as, for them, an occupied orbital is localized on the model boundary and hence could be considered as an artefact induced by using molecular cluster models. Indeed, for 2, all charge transfer bands are by necessity located on the boundary. Models like 3 and 4 support centrally located charge-transfer states, and, for 4, 17 chargetransfer transitions are predicted by TD-HSE06 at under 3 eV in energy, many of which involve more central destination orbitals. The charge-transfer effect depicted is therefore a robust prediction of the calculations; realistic predictions of its innite-sample limit are provided later using 2D periodicmodel simulations. Fig. 2(d) lists the orbital numbers of the illustrated orbitals, with #374 being the HOMO orbital and #375 the LUMO. The CAM-B3LYP ordering is indicated and identies (1) 1 A 2 00 and (1) 3 A 2 00 as resulting from the HOMO to LUMO transition; analogous HSE06 results are similar, see ESI. † Surprisingly, the low-energy states (1) 1 E 0 and (1) 3 E 0 arise from excitations deep into the unoccupied orbital space to orbitals #406 and #407. These are essentially defect orbitals, but their presence within the h-BN conduction band induces some mixing. As is oen found for defects, 67 orbital energy differences provide poor indications of state energy differences, particularly when charge transfer is involved. The Coulomb interactions between the electron and the hole involved in the excitation are critical to the energetics, interactions that scale with the asymptotic potential for charge-transfer transitions. Table 4 also shows transition energies for 1 evaluated using equation of motion coupled cluster (EOM-CCSD) theory, 127 another time-dependent approach considered to provide useful descriptions of the singlet and triplet manifolds of V N À . It does not suffer from anomalies concerning its description of longrange electrostatic or charge-transfer effects, and its predictions are similar to those of CAM-B3LYP, but sometimes different to those of HSE06.
(d) The asymptotic-potential error in 2D defects perceived using DFT and up to CCSD(T) calculations on model defect clusters As an alternative to TD-DFT, we now consider transition energies evaluated by calculating individual DFT energies, evoking the Kohn-Sham 128 and Gunnarsson-Lundqvist 129 theorems, for both the initial and nal states. These approaches can only be applied to only a limited selection of states, and is signicantly hampered for V N À as DFT and many ab initio approaches fail   Table 5 and illustrated in Fig. 2(c). Used are the CAM-B3LYP and HSE06 density functionals, as well as the ab initio approaches: Hartree-Fock theory (HF), 130 second-order Møller-Plesset (MP2) theory, 131 coupledcluster singles and doubles (CCSD) theory, [132][133][134] and this perturbatively corrected for triples excitations, CCSD(T). 135 For the local excitation (1) 3 E 0 , the ab initio methods appear to converge quickly as the treatment of electron correlation is systematically enhanced, giving results close to those of CAM-B3LYP, with HSE06 deviating slightly further. For the charge-transfer excitation (1) 3 A 2 00 , CAM-B3LYP deviates from CCSD by 0.4 eV, a large deviation, but one consistent with other worst-case predictions obtained using CAM-B3LYP. 67 On the other hand, HSE06 underestimates this energy by 0.9 eV. This is another example of the previously identied asymptotic-potential effects (2)-(4) from Table 2.
(e) The asymptotic-potential error in 2D defects perceived using DFT calculations on 2D-periodic models Transition energies can also be obtained for 2D-periodic models from energy differences, but it is difficult to get results for all but the lowest-energy state of each symmetry, and then only reasonable results can be expected if the two states of interest have minimal open-shell character. Results obtained for the (1) 3 B 1 component of the (1) 3 E 00 state are reported in Table 5, with the electron densities of the two key orbitals involved shown in Fig. 3 for P 75 (many more frontier orbitals, plus analogous results for P 53 and P 64 , are shown in ESI †). Owing to the use of boundary conditions in the periodic model, the defect, which should show both D 3h and its subgroup C 2v symmetry, cannot simultaneously display both 3-fold symmetry elements and the C 2v elements, with the utilized lattices retaining only C 2v symmetry. In this reduced symmetry, the (1) 1 A ' 1 ground state becomes (1) 1 A 1 , the HOMO orbital changes from a 2 00 to b 1 , and we focus on the a 1 component of the important e ' orbital from Fig. 2(d).
Comparing Fig. 3 and 2(d), we see that HSE06 and CAM-B3LYP both predict that the HOMO orbital of the singlet ground state has very similar form in both the molecular model 4 and the 2D model P 75 , being tightly localized inside the defect. For the 2D model, Fig. 3 shows that this orbital retains its form for spin-up electrons in the lowest triplet state. For the LUMO of the singlet ground state, Fig. 3 also shows similar results from HSE06 and CAM-B3LYP, depicting an orbital that is delocalized over the h-BN, but is most prominent in the box corners so as to maximize distance from the defect. The LUMO orbital becomes occupied in the dominant-spin component of the lowest-energy triplet state, and for it, Fig. 3 reveals very different characteristics predicted by HSE06 and CAM-B3LYP. CAM-B3LYP predicts a localized orbital similar to that predicted by the molecular model 4 (Fig. 2(d)), whereas HSE06 predicts an extremely delocalized orbital. Viewed from the perspective of the groundstate structure, what this indicates is that the electron-hole interaction energies perceived by CAM-B3LYP and by HSE06 are very different. This effect was noted previously for molecular spectroscopy as effect (2) in Table 2, and the HSE06 results are characteristic of signicant misrepresentation of the asymptotic potential. As Fig. 3 also indicates that effect (1) holds, with the HSE06 band gap being substantially less than that for CAM-B3LYP, effects (3) and (4) also hold and hence HSE06 signicantly underestimates the energies of charge-transfer transitions.
Expanding on the signicance of the electron-hole interaction energy, the extended results shown in ESI † indicate that both HSE06 and CAM-B3LYP predict that the (b 1 / a 1 ) excitation has a very large effect on the nature of most of the frontier orbitals. Hence simplistic ideas that the ground-state orbital band structure alone is sufficient to determine excited-state energetics clearly do not apply in this case. Just as the electron-hole interactions generated by the transition dominate the model-compound state properties revealed in Fig. 2, so it also controls the results of the 2D-periodic models. Errors in the asymptotic potential directly manifest in the locationdependence of the electron-hole interaction energy. As a result, HSE06 underestimates the energy of charge-transfer bands. For V N À1 , this effect is serious as the lowest-energy Table 5 Calculated vertical excitation energies for the V N À defect in h-BN, in eV, from the (1) 1 A 1 0 ground state a of the model compounds 1-4 (in D 3h symmetry), and the periodic layers P 53 , P 64 , and P 75 (in just C 2v symmetry as necessitated by the boundary conditions, utilizing either 1 Â 1 Â 1 or 1 Â 2 Â 2 k-points), obtained from state energy differences transition is incorrectly perceived as being a defect to conduction-band transition instead of an intra-defect transition. So, whereas 2D-periodic and molecular CAM-B3LYP simulations converge to very similar results (Table 5), HSE06 simulations do not.

Conclusions
The most important aspects of asymptotic correction to density functionals were identied through reviewing a key application in molecular spectroscopy and listed as ve effects in Table 2. Considering just bandgaps of simple materials, effects (1) and (2) were then shown to be signicant, despite these materials not having dominant charge-transfer character: asymptotic correction increases orbital bandgaps and exciton couplings to provide improved descriptions of the lowest-energy spectroscopic excitation. Effects (3) and (4), concerning the relative ordering of transitions of different charge-transfer types within the excited-state manifold, were then demonstrated to have profound effects on spectral assignment by considering the spectrum of the V N À defect in h-BN. When charge transfer is critical, methods without asymptotic correction fail to deliver qualitatively useful results; when charge-transfer character is not critical, asymptotic correction still leads to signicant quantitative improvements in calculated properties. Overall, correcting the asymptotic potential results in signicant increases to band gaps, differences that are compensated by analogous increases in exciton binding energies for charge-localised transitions, but le standing for charge-transfer bands. The net effects include: accurate prediction of band gaps, accurate predictions of exciton binding energies, accurate predictions of spectroscopic transition energies, and improved prediction of geometrical properties such as lattice constants. Related known effects include signicant improvements to calculated reorganisation energies, Huang-Rhys factors, and chemical and photochemical reaction rates, paving the way for the computational assignment of uninterpreted spectra.
HSE06 should be replaced as the default standard for the accurate and reliable prediction of materials spectroscopic properties by asymptotically corrected functionals. Approaches like CAM-B3LYP are now available in materials-science computational packages such as CPMD 31 and VASP 32 and require no more computational resources than does HSE06. 32 They offer signicant improvement for property prediction, in most cases. The exception to this concerns the prediction of the band gap of silicon, an unusual but extremely important Fig. 3 Shown for of the V N À1 defect in h-BN are the HOMO and LUMO relative electron densities, along a plane 0.8Å above the atoms, for orbitals #278 and #279 involved in the (1) 3 B 1 (b 1 / a 1 ) component of the (1) 3 E 00 lowest-energy state. These are based upon 2D model P 75 , using 1 Â 2 Â 2 k-points, and list each orbital energy, symmetry, and occupancy (see SI for many more frontier orbitals and analogous results for P 53 and P 64 ); for the triplet state, only densities for the spin-up component reflecting the electron majority are shown.
material. Density functionals with asymptotic correction, designed for materials only applications, are now also becoming widely available, [33][34][35] and offer solutions for silicon and related problems, whilst more generally applicable asymptotically correct methods are currently also being devised. 23,[36][37][38][39][40][41][42][43][44] Optimal methods that are independent of dimensionality will provide the opportunity for unied understanding of molecular and materials spectroscopies.

Methods
As shown in Fig. 1, calculations are performed for the model compounds 1, 2, 3, and 4 of V N À that comprise rings of B and N atoms surrounding the nitrogen vacancy, as well as for the model 2D periodic layers P 53 , P 64 , and P 75 . All molecular-model geometries are optimized using CAM-B3LYP/6-31G*, while the layers are optimized using HSE06; full details are provided in ESI. † Only vertical excitation energies are considered. The molecular calculations are performed using Gaussian-16 (ref. 136) B.01 with mostly the 6-31G* basis; 137 convergence of defect spectroscopic properties with respect to basis-set expansion is known to be very rapid. 72 The periodic-cell calculations are performed using VASP 5.4.4 (ref. 32, 138 and 139) (using "PREC ¼ HIGH", "PREFOCK ¼ NORMAL", PBE-PAW pseudopotentials 140 for B3LYP and HSE06 calculations and GW-PAW ones for CAM-B3LYP, with the implied basis sets (NGX, NGY, and NGZ, etc.) used for each calculation listed in the ESI †). 32 They are performed using either 1 Â 1 Â 1, 1 Â 2 Â 2, or 1 Â 3 Â 3 representations of the Brillion zone of (5 Â 3O3)R30 (P 53 ), (6 Â 4O3)R30 (P 64 ), and (7 Â 5O3)R30 (P 75 ) unit cells, with lattice vectors depicting an intrinsic h-BN BN bond length of 1.452Å. Fermi-level smearing is also applied, at an electronic temperature of 0.02 eV. Dipole corrections to mitigate the ctitious interactions between parallel h-BN planes are not used as they are only of order 0.01-0.03 eV. All singlet-state calculations are performed spin restricted, whereas all triplet-state calculations are spin unrestricted. Symmetry analyses 141 of the VASP excited states are listed in the ESI, † along with wavefunction projections that establish the relationship between the ground and excited states.
The new CAM-B3LYP results reported for lattice parameters (ESI Table S1 †), orbital bandgaps (ESI Table S2 †) and exciton binding energies (text) were also determined using VASP 5.4.4, 32,138,139 with optimised coordinates and calculation details for each material listed in the ESI. † PBE-PAW pseudopotentials 140 were again used for B3LYP and HSE06 calculations, with GW-PAW ones used for CAM-B3LYP. 32 The plane-wave basis sets were truncated at an energy of 850 eV, with many k-points sampled until a minimum was obtained in the orbital bandgap, see ESI Table S3. † Exciton-binding energies were determined for vertical excitation by direct evaluation of the energy of the triplet excited states corresponding to the observed singlet transitions, followed by application of (very minor) corrections for the singlet-triplet differences, with details given in ESI Table  S4. † The B3LYP geometry optimisations reported in Table 3 were performed using PREC ¼ HIGH and 12 Â 12 Â 12 k-points. These calculations differ from previously reported ones through the use of more modern pseudopotentials and the use of analytical functional derivatives. 32

Data availability
Signicant data is provided in ESI, † with other information availabe from the authors.

Author contributions
M. L., assisted by J. R. R., performed the calculations; J. R. R., assisted by M. L. and all others wrote the manuscript; R. K., R. D. A., M. J. F., and J. R. R. designed and supervised the work.

Conflicts of interest
There are no conicts to declare.