Density Functional Study of Tetraphenylporphyrin Long-Range Exciton Coupling

The performance of time-dependent density functional theory (TDDFT) for calculations of long-range exciton circular dichroism (CD) is investigated. Tetraphenylporphyrin (TPP) is used as a representative of a class of strongly absorbing chromophores for which exciton CD with chromophore separations of 50 Å and even beyond has been observed experimentally. A dimer model for TPP is set up to reproduce long-range exciton CD previously observed for a brevetoxin derivative. The calculated CD intensity is consistent with TPP separations of over 40 Å. It is found that a hybrid functional with fully long-range corrected range-separated exchange performs best for full TDDFT calculations of the dimer. The range-separation parameter is optimally tuned for TPP, resulting in a good quality TPP absorption spectrum and small DFT delocalization error (measured by the curvature of the energy calculated as a function of fractional electron numbers). Calculated TDDFT data for the absorption spectra of TPP are also used as input for a ‘matrix method’ (MM) model of the exciton CD. For long-range exciton CD, comparison of MM spectra with full TDDFT CD spectra for the dimer shows that the matrix method is capable of producing very accurate results. A MM spectrum obtained from TPP absorption data calculated with the nonhybrid Becke88–Perdew86 (BP) functional is shown to match the experimental brevetoxin spectrum ‘best’, but for the wrong reasons.


Introduction
The coupling of electronic excitations in molecular aggregates or in multi-chromophore molecules and metal complexes has long been an important research topic in chemistry and physics. [1,2] Of special interest are situations where the electric transition dipoles of the uncoupled excitations are geometrically arranged along a 'helical' path such that the coupled system exhibits electronic exciton circular dichroism (CD). [3][4][5][6][7] Exciton CD reveals a great amount of information about the relative orientation of individual chromophores with respect to each other and about the distance between them. Exciton CD can be very strong, even if the individual chromophores have no intrinsic chirality. In leading order, and at large separation of the chromophores, the rotatory strengths of the coupled transitions (the integrated intensities of the CD of individual excitations) are determined simply by the lengths, relative orientations, and relative distances of the electric transition dipole vectors as well as by the energies of the uncoupled transitions. [8][9][10] Equations for the 'matrix method' (MM) dipole coupling model [8] are provided in the Discussion section to illustrate the case.
Berova and collaborators have experimentally detected very long-range exciton CD with chromophore separations up to 50 for exciton coupling between two porphyrin based chromophores, for instance for tetraphenylporphyrin (TPP). [11] Tsubaki et al. have reported observations of TPP exciton CD of substituted chiral oligonaphthalenes at even larger distances of approximately 66 . [12,13] The strong coupling between TPP substituents in suitably derivatized biomolecules has allowed researchers to derive important information about their threedimensional structure. Examples where porphyrin-based exciton CD has been investigated or used in this context include the use of magnesium porphyrin to determine the absolute stereochemistry of chiral alcohols, [14] a theoretical study of Soret band coupling of bis-porphyrin derivatives using various exciton coupling models, [15] a determination of absolute stereochemistry of cyclic a-hydroxyketones with zinc-TPP tweezers, [16] and investigations of the conformational space of DNA. [17] A study of bis-porphyrin dimers of various derivatized biomolecules, including a brevetoxin bis-TPP derivative shown in Figure 1, has established the very long-range nature of the exciton CD. [11] In first-principles calculations of electronic absorption and CD spectra, the size of the system-dictating the number of basis functions, B-is one of the major limiting factors because of the scaling of the computational effort with B. Presently, the most frequently utilized method for electronic spectra calculations is (linear response) time-dependent density functional theory (TDDFT). [18][19][20] Reasonably efficient wave function-based methods for the treatment of exciton CD are available as The performance of time-dependent density functional theory (TDDFT) for calculations of long-range exciton circular dichroism (CD) is investigated. Tetraphenylporphyrin (TPP) is used as a representative of a class of strongly absorbing chromophores for which exciton CD with chromophore separations of 50 and even beyond has been observed experimentally. A dimer model for TPP is set up to reproduce long-range exciton CD previously observed for a brevetoxin derivative. The calculated CD intensity is consistent with TPP separations of over 40 . It is found that a hybrid functional with fully long-range corrected range-separated exchange performs best for full TDDFT calculations of the dimer. The range-separation parameter is opti-mally tuned for TPP, resulting in a good quality TPP absorption spectrum and small DFT delocalization error (measured by the curvature of the energy calculated as a function of fractional electron numbers). Calculated TDDFT data for the absorption spectra of TPP are also used as input for a 'matrix method' (MM) model of the exciton CD. For long-range exciton CD, comparison of MM spectra with full TDDFT CD spectra for the dimer shows that the matrix method is capable of producing very accurate results. A MM spectrum obtained from TPP absorption data calculated with the nonhybrid Becke88-Perdew86 (BP) functional is shown to match the experimental brevetoxin spectrum 'best', but for the wrong reasons.
well. [21] The attractiveness of TDDFT derives from the fact that it incorporates electron correlation at a computational cost that is comparable to Hartree-Fock (HF) theory, if hybrid functionals are used for the exchange. The scaling is formally of order B 4 but in practice often lower. Non-hybrid functionals afford better scaling yet. Still, the computational demand can be prohibitive if one attempts a reliable first-principles theoretical modeling of large systems that are of interest in biochemistry, not only because of the scaling as a function of B but also because of very large conformational spaces that may be encountered. The application of a dipole coupling model for well separated chromophores, or a more refined model including higher orders of the multipole expansion, offers a way out of this dilemma as far as exciton coupling is concerned. [15] For instance, exciton coupling models have been applied successfully in a recent study of the CD spectrum of bacteriorhodopsin. [22] Moreover, input data for such a model can be calculated from first principles for each chromophore separately, and the coupled excitations are then obtained from the lowerlevel coupling model at essentially negligible computational overhead. Herein, we explore this computational route for the long-range exciton CD of TPP pairs as encountered, for instance, in the aforementioned study of derivatized brevetoxin, [11] in conjunction with the 'matrix method' (MM). One aim of this work is to investigate the performance of such a two-level computational model by comparing the results of a dipolecoupling scheme using TPP monomer TDDFT data as input with full TDDFT dimer calculations.
TDDFT with popular standard hybrid and non-hybrid functionals may afford large errors in calculated excitation energies, if the excitation has an explicit or hidden charge-transfer (CT) character. [23][24][25][26][27] For multi-chromophore systems, the CT problem may also create large numbers of spurious low-energy CT excitations. [28,29] The CT problem of TDDFT can be corrected effectively by employing hybrid functionals with range-separated exchange, [30][31][32][33][34][35]26] in particular if the functional goes to pure HF exchange asymptotically (full long-range correction [LC]). However, the range-separation parameter in the exchange functional is strongly system dependent [36] and should therefore be determined system-specifically. Recently, there have been ways proposed of how to achieve a system-specific optimal 'tuning' of the range-separation parameter as well as other parameters in the exchange functional in an ab initio sense, based on cri-teria rooted in density functional theory (DFT). [27,[37][38][39][40][41][42][43] Another aim of this work is to investigate whether such a functional tuning is beneficial in the description of the TPP excitation spectrum and in the calculation of exciton coupling CD of TPP dimers. It is shown that the use of a system-specific optimally tuned range-separation parameter significantly improves the calculated TPP absorption spectrum compared with a universal parametrization. A LC functional is also shown to be the best choice for the study of long-range exciton coupling in full TDDFT calculations of the coupled system, as it suppresses spurious CT excitations. At large separations of the chromophores, it is shown that a simple dipole coupling model based on TDDFT monomer input data gives excellent agreement with long-range corrected TDDFT dimer spectra.

Computational Details
A geometry optimization of tetraphenylporphyrin (TPP) was performed with the Amsterdam Density Functional (ADF) program, [44][45][46] using the Becke88-Perdew86 (BP) functional, [47][48][49][50] dispersion corrections as devised by Grimme et al. [51] (see the ADF documentation regarding the DFT-D3 parametrization), and a combination of double-z polarized (C,N) and double-z (H) Slater-type basis sets from the ADF basis set library. Electronic absorption spectra of TPP were calculated with the NWChem program [52] using time-dependent Hartree-Fock (TDHF) [53,54] and TDDFT, [55][56][57] employing the 6-31G(d) Gaussiantype basis. The TDDFT calculations were carried out with the BP functional and with a fully long-range corrected hybrid functional with range-separated exchange based on the Perdew-Burke-Ernzerhof (PBE) exchange-correlation functional. [58,35,59] We use the acronym LC-PBE0 for the range-separated functional in this work. Details about the optimal tuning of this functional are provided in the Results and Discussion section. The tuning procedure was followed by a series of singlepoint energy calculations for TPP with integer and fractional electron numbers using a fractional occupations/fractional total electron number code implemented in NWChem by one of the authors. [41] For additional information on the functional tuning, please see Ref. [36] and our group's previous work on functional tuning. [41,40] The basis sets were chosen because they are computationally efficient, but at the same time they produce acceptable results for the TPP monomer benchmark (see below). The excitation spectrum is dominated by valence transitions that do not mandate diffuse functions in the basis set.
A TPP dimer configuration was derived from a structure of the brevetoxin derivative BTX-D(1,1,42)-1,42-diol bis(p-[10',15',20'-triphenyl-5'-porphyrinyl] benzoate) (referred to as BTX-D) as shown in Figure 1 and chosen to model exciton coupling of TPP dimers at large spatial separations. The geometry was inferred from Ref. [11], where a CD spectrum for BTX-D has been previously reported. A CD couplet in the spectral range of the TPP Soret band has been assigned to exciton cou- Figure 1. Brevetoxin substituted with TPP used in Ref. [11] to demonstrate long-range exciton CD.
pling of the respective TPP transitions at distances of up to 50 . [11] A structure of BTX-D was obtained from a conformational search using Spartan'08 and subsequent molecular mechanics optimizations. [60] The steroid linkage was then removed and each TPP unit was replaced with the BP optimized structures. The model affords an interchromophoric separation of approximately 42 . One of the TPP moieties was rotated along the centroid-centroid axis to match relative orientations as inferred from the experimental BTX-D CD spectrum. [11] The resulting geometry is referred to as the TPP dimer in this article. Our aim was not to explore the full conformational space of BTX-D but rather to obtain a reasonable separation between the TPP moieties for a low-energy conformer (because of the strong distance dependence of the exciton CD) in order to benchmark the matrix method (MM) results in comparison with full TDDFT calculations for the dimer. MM dipole-coupled CD spectra for the dimer were based on monomer spectral data calculated with TDHF, BP, and a global and tuned parametrization of LC-PBE0. Full TDDFT dimer calculations of CD spectra were performed with NWChem and Gaussian [61] with the 6-31G(d) basis. Dimer TDDFT calculations with the BP functional were performed with Gaussian due to convergence problems with the NWChem TDDFT code for large numbers of excitations (which are needed with BP in order to cover the experimentally accessible spectral range). Most of the CD spectra are based on the dipole-length representation of the rotatory strength. One of us has recently reported a TDDFT implementation in NWChem for CD spectra calculations with gauge-including atomic orbitals (GIAOs, or London orbitals) and with the dipole-velocity gauge. [62] This code was used to test that the dimer spectra reported herein are not contaminated by a gauge-origin dependence (see below and Ref. [62]). For more information on the MM, see work by Schellman et al. [8] and previous work by our group on exciton coupled CD spectra. [10] All calculated absorption and CD spectra were Gaussian broadened with a s value of 0.13 eV.

Results and Discussion
The geometry optimization of tetraphenylporphyrin (TPP) yielded a porphin backbone adopting a saddleback-shaped geometry and phenyl substituents with dihedral angles of approximately 66 8 between the porphin plane and the phenyl planes. The geometry is consistent with previous theoretical results [63][64][65] and with experimental data. [66] The remainder of this section is divided into three parts: 1) g-tuning of the LC-PBE0 functional for TPP; 2) an analysis of the TPP absorption spectrum calculated with various functionals, and comparisons with previously published calculated spectra and with experiment; 3) the analysis of the exciton CD for the TPP dimer model.

Tuning of the LC-PBE0 functional for TPP
A system-specific tuning of the LC-PBE0 functional for TPP was carried out by minimizing J' 2 , of Equation (1), as a function of the range-separation parameter g in the exchange functional. Here, N is the number of electrons for neutral TPP.
Specifically, we used the following range-separation of 1/r 12 in the exchange [Equation (2)]: [32] 1 The ionization potential (IP) is calculated as the difference of total energies, IP(N) = E(NÀ1)ÀE(N), and similarly for the corresponding (N + 1)-electron system. Further, e HOMO is the energy of the highest occupied molecular orbital (HOMO). For the elusive exact universal density functional, J' 2 = 0, the negative of the HOMO orbital energy would correspond exactly to the first ionization potential (IP). For an approximate functional, the idea is to determine the range-separation parameter such that an optimally small J' 2 is obtained, thereby giving physical meaning to the HOMO energy of both the N-electron and the (N + 1)-electron system. The process is thought to improve the fundamental gap. [42] The other two parameters in Equation (2), a = 0.25 and b = 0.75, were kept as in the originally proposed parametrization, [35] which affords g ¼ 0:3a À1 0 . We noted in previous work that tuning according to the IP criterion requires a + b = 1. [40,41] For a procedure to determine a to minimize the DFT delocalization error simultaneously with an optimal tuning of g per Equation (1), see Ref. [40]. The DFT delocalization error [67] can be quantified by the behavior of E(N) for fractional and integer electron numbers N. The exact energy should afford linear segments between integers, with slopes changing discontinuously [68] at integer values of N. Curvature in plots of E(N) is therefore indicative of delocalization error.
In Figure 2, J' 2 calculated for TPP is shown for varying g. Recent theoretical studies have shown [27,41] that as the size of a conjugated delocalized p system increases, the optimal range-separation parameter g tends to decrease. In turn, g is the reciprocal of a cut-off distance beyond which HF exchange in the functional starts to dominate. That is, a reduced value of g indicates that the delocalized p system benefits from local DFT exchange acting over a longer range of interelectronic distances. The optimal g * value for TPP is 0.105a À1 0 , that is, much below the originally proposed value of 0.3a À1 0 .
The behavior of E(N) as a function of a fractional electron number N is compared, in Figure 3, for the original parametrization of LC-PBE0 (a = 0.25, b = 0.75, g = 0.3) and the tuned version (g * = 0.105). We refer to the tuned functional as LC-PBE0g * from here on.
The stock parametrization of LC-PBE0 affords delocalization error for TPP, as demonstrated by the curvature in E(N). The negative curvature indicates that, for TPP, the value g = 0.3 gives too little delocalization. The behavior of E(N) for the tuned functional is much better. The curvature is nearly vanishing for both the electron-deficient and the electron-rich species. Because of the reasonably small curvatures obtained with g * at a = 0.25, we decided to forego a simultaneous optimization of a and g for TPP.

TPP excitation spectrum
Since the excitation spectrum of TPP in this work is used as the input for the matrix method (MM) coupling model, it is important to assess the quality of the spectrum calculated at various levels of theory. The currently accepted assignment of the experimental TPP spectrum, in reference to the Gouterman model, [69] is as follows: Weak absorption bands (Q-bands) at 1.86 eV (Q x ) and 2.26 eV (Q y ) are assigned to HOMO-to-LUMO and HOMO-to-LUMO + 1 transitions, respectively, in a molecular orbital (MO) representation. Isosurfaces of the relevant MOs are shown in Figure 4. A very intense absorption band around 3.06 eV, commonly referred to as the Soret band, is caused by a pair of transitions (B x , B y ) that are assigned to HOMOÀ1-to-LUMO and HOMOÀ1-to-LUMO + 1, respectively. The pairs of excitations responsible for the Q and B bands are calculated as 1 B 3u / 1 B 2u pairs. The symmetry labels for the excitations refer to the D 2h point group, even though the symmetry of the optimized TPP geometry is lower (C 2v ). However, the MO contributions to the excitations and the MO nodal patterns are similar to porphin as described by the Gouterman model. Note that the N transitions, assigned in previous and current works, are not considered in the Gouterman model, and therefore alternate criteria were taken from the literature. An excitation around 3.5 eV (N x ) has been assigned to a 1 B 3u excited state in Ref. [70] based on similarity-transformed equation-of-motion coupled-cluster singles and doubles (STEOM-CCSD) calculations on free base porphin. From Ref. [63], the N y transition was assigned to the corresponding porphin 1 B 2u excitation with dominant MO contributions from HOMOÀ3 to LUMO + 1 which is consistent with the calculations in this article. The spectral features of TPP above 4 eV are ignored in this work since the Q, B, and N transitions are the ones most important for the experimentally observed exciton coupling CD of TPP dimers.
In Table 1, the Q, B, and N transitions of TPP calculated by TDHF, and by TDDFT with various functionals, are compared with previously reported calculated [63] and with experimental  data. [71] The transition assignments were made based on the dominant contributions from occupied (occ) and unoccupied (unocc) MOs to the excitation transition densities, and in comparison with previous porphin [70,72] and TPP [63] calculations as well as the Gouterman model. The frontier orbitals obtained with Becke88-Perdew86 (BP) and LC-PBE0 are compared in Figure 4. Differences between the MOs obtained with different functionals are hardly visible in the isosurface plots. The nodal patterns of the frontier orbitals match those of the Gouterman model [69] and agree with previous DFT calculations. [63,73] Simulated spectra are shown in Figure 5. The THDF calculation yields Q bands which are red shifted with respect to the experiment, while the B and N bands are blue shifted. The assignment of the transitions, based on dominant contributions to the transition density matrix from occupied and unoccupied MOs, is in qualitative agreement with the Gouterman model. The N bands can be assigned as transitions from HOMOÀ4 to LUMO (N x ) and HOMOÀ4 to LUMO + 1 (N y ) where HOMOÀ4 of the HF calculation is equivalent to HOMOÀ3 for the DFT calculation of Ref. [63] and the present calculations for BP, LC-PBE0, and tuned LC-PBE0.
The BP functional yields blue shifted Q x and red shifted Q y transitions, with an assignment of these excitations that is in agreement with the Gouterman model. According to this model and the accepted assignment of the spectrum, the B x excitation should be lower in energy than B y . The BP functional incorrectly reverses the energies. This finding is consistent with a BP spectrum previously reported in Ref. [74]. An apparent N band may be attributed to a modestly intense pair of transitions at 3.7 eV, but the excitations do not have the expected MO contributions. With the BP functional, the excitations with strong HOMOÀ3-to-LUMO and HOMOÀ3-to-LUMO + 1 charac-ter expected for the N transitions appear at lower energy than the Soret peak, that is, below 3 eV. Moreover, the BP 'stick spectrum' in Figure 5 reveals many spurious excitations, in particular between 3 and 4 eV, with (mostly) low intensity. These transitions are not seen in the other spectra obtained with HF and asymptotically correct density functionals. The occurrence of these spurious transitions is likely a consequence of the incorrect long-range behavior of the BP exchange-correlation (XC) potential. In recent work, Baer, Kronik et al. noted that 'charge-transfer-like' excitations may occur in extended p chromophores, which render calculations with functionals that do not afford the correct long-range behavior suspicious. [39] The XC potentials of TDHF and LC-PBE0 (standard and tuned version) afford the correct long-range behavior, and therefore spurious transitions are not as prevalent. Clearly, the BP spectrum is seriously deficient despite the fact that the broadened absorption intensity resembles the experimental spectrum reasonably well.
The LC-PBE0 and LC-PBE0-g * functionals predict blue shifted Q, B, and N bands when compared to experiment (with the exception of Q y calculated with LC-PBE0). The tuning procedure is designed to improve the HOMO-LUMO energy difference toward the fundamental gap, and it has been shown that TDDFT response calculations based on the resulting functionals provide improved excitation energies as well [42] (where, similar to HF, the optical gap is calculated to be smaller than the HOMO-LUMO gap). From Table 1, the Q x transition is seen to be predicted worse for the tuned functional than for the stock parametrization. For all other transitions, however, the tuned LC-PBE0-g * functional gives excitation energies that are closer to the experimental band maxima. The g tuning also closes the gap between the B and N bands, which is favorable.
In Figure 6, the BP and tuned LC-PBE0 absorption spectra are compared to calculated spectra reported previously in the literature (GW + BSE [73] and TDDFT with the PBE0 functional [63] ) and to an experimental spectrum. [71] The TDHF spectrum is excluded because the excitation energies do not agree well with experiment. LC-PBE0 is also excluded because the tuned version yields an overall better spectrum. The GW + BSE and PBE0 spectra were generated by Gaussian line broadening (s = 0.13) based on excitation data taken from Ref. [73] and Ref. [63], respectively. For better comparison, the Soret maxima were scaled to the same value and the spectra are given in arbitrary units (a.u.).
The BP density functional is computationally rather inexpensive. However, as pointed out above, despite the apparent reasonable agreement with experiment, the deficiencies of the functional are clear when considering the assignments of the relevant transitions. For instance, the reader is reminded that the 'N' band in the BP spectrum is not the N-band proper. The computational savings from not including exact exchange (and a range separation of the exchange) are also not as advantageous as one might think because of the need to calculate a large number of excitations to cover the spectral range up to 4 eV (most of which are unphysical). This problem is exacerbated in the dimer calculations discussed below. The tuned LC-PBE0-g * functional performs well in comparison. We also note good agreement of the tuned hybrid with the GW + BSE calculation, apart from a modest blue shift.

Dimer calculations: TDHF, TDDFT, and the dipole coupling model
A TPP dimer was investigated in a relative arrangement approximating the geometry of BTX-D, as described in the computational details. The dimer model along with the relative orientation of two selected transition moments (B x ) is shown in Figure 7. This TPP dimer was used for all calculations discussed in this subsection. Given the large separation of the chromophores, it is important to investigate if the CD spectra calculated for the full dimer are affected by the gauge-origin problem of magnetic properties in finite-basis set calculations (see Ref. [62] and Ref. [75] for a detailed discussion of the problem in the context of CD spectra calculations, and for additional literature devoted to this topic). Figure 8 shows a comparison of TDHF calculations of the full dimer CD spectrum based on different rotatory-strength representations: dipole-length (standard basis), dipole-velocity, and dipole-length with use of gauge-including atomic orbital (GIAO) basis functions. The velocity and GIAO rotatory strengths are both origin invariant by design, but the velocitygauge may be affected more strongly by basis-set incomplete-ness. As the calculated CD spectra demonstrate, there is virtually no difference between the three CD spectra. Consequently, in order to save computational resources, all other dimer calculations have been performed with the dipole-length gauge without GIAOs.
In the matrix method (MM) coupling model, the individual excitations in individual chromophores are treated as quasiparticles (excitons) that interact electrostatically via their electric and magnetic multipole moments. For details and additional citations to original references not cited herein please see the Appendix of our previous work. [10] A coupled dimer with a single transition per monomer, for instance, is described by the Hamiltonian [Equation (3)], where the monomer excitations are labeled a, b, etc. and have energies E a , E b , etc. To lowest order in the multipole expansion, the coupling only considers the electric transition dipole moment vectors d a , d b , etc., which gives rise to the electrostatic interaction [Equation (4)]: Here, r ab = r b Àr a is the distance vector between the chromophore centers r a and r b (we use the center of nuclear charge), and r ab its length. Electrostatic units are used, where the square of an electric dipole divided by a volume yields an energy, and a magnetic dipole has the same unit as the electric dipole. In SI units, the equation for the potential carries an additional factor of (4pe 0 ) À1 . In the MM, the magnetic transition dipoles associated with the uncoupled excitations are obtained from the electric transition dipoles via Equation (5).
The excitation energies of the coupled system are the eigenvalues of H. To facilitate the next step, the eigenvectors of H are collected column-wise in a matrix C, and the uncoupled electric and magnetic transition dipoles are collected in matrices D 0 and M 0 , respectively. The coupled transition moments d i , m i for excitations i = 1, 2, etc. in the compound system are then obtained as the columns of the matrix products D 0 C and M 0 C. The rotatory strength of each coupled transition i is given as In the MM spectra reported in this subsection, the transition dipoles from a number M calculated excitations for the TPP monomer (using M = 25 for TDHF, LC-PBE0, and LC-PBE0-g * , and M = 100 for BP) have been used in the coupling model. That is, dipole interactions between pairs of transition dipoles centered on different TPP moieties were calculated and used in a 2 M by 2 M matrix model in order to obtain the coupled excitation energies and rotatory strengths for the dimer. Figure 9 displays the broadened CD spectra of the TPP dimer of Figure 7 obtained using the matrix method with calculated TPP monomer spectral data, and from full dimer calculations. The MM spectra generated from monomer data for a particular functional are indicated by dashed lines. For the theoretical methods that afford a correct long-range behavior of the XC potential (TDHF, and TDDFT with LC functionals), the agreement with the full dimer spectra is excellent. Some deviations occur near the high-energy cutoff of the spectra. We consider this an artifact of the upper-energy cutoff for the TDDFT dimer spectra, made necessary by the memory limitations of Figure 6. A) Calculated absorption spectra compared to previously reported theoretical absorption spectra: [63,73] BP (c), tuned LC-PBE0 (a), GW + BSE (c), PBE0 (a). B) Selected calculated spectra compared to experiment: [71] BP (c), tuned LC-PBE0 (a), experiment (c). The intensities were scaled such that the Soret peaks have the same intensities.  Figure 11. Phenyl groups were replaced with hydrogen for clarity, and interchromphoric distance was reduced for the plot. the code used. That is, the upper energy range of the broadened dimer spectra may change in intensity if additional, possibly intense, excitations were available to generate the broadened spectra. For the exciton CD in the energy range of the calculated TPP Soret bands, the excellent agreement of the MM spectra with the full TDDFT spectra (HF and LC functionals) demonstrates that this couplet is indeed attributable to longrange dipole interactions of the intense TPP B transitions, along with some dipole-coupling contributions from the Q and N transitions.
The BP dimer spectrum is clearly affected strongly by the TDDFT CT problem. The dimer calculation afforded a large number of spurious transitions; over 150 excitations were needed to cover the rather small energy range shown in Figure 9. For comparison, 100 excitations in the BP monomer spectrum, many of which are already spurious, reach up to 5.1 eV. Some of the deficiencies noted for the TPP monomer spectrum are not hidden anymore in the dimer. A MM CD spectrum constructed from the Q, B, and higher energy ('N') transitions of the BP calculation (Figure 9 B) leaves out spurious transitions and shows that such a coupling would at least yield a qualitatively correct spectrum. In contrast, at 42 separation of the TPP moieties, the LC-PBE0 calculations (standard and tuned versions) afford the expected physically correct picture: The exciton CD of the dimer just above 3 eV is caused by simple electrostatic coupling, mainly, of the TPP B transition dipoles.
The deficiencies in the BP dimer TDDFT CD spectrum therefore have several origins. For instance, the MM(BP) spectrum displays a low-energy tail of the first (negative) CD band that is not visible in the dimer calculations with the other functionals. As pointed out above, the excitations that would be expected to form the TPP N band are obtained as pre-Soret features in the BP monomer spectrum (see also Ref. [74]). It is the dipole coupling of these excitations that creates the intense lowenergy tail in the MM(BP) CD spectrum. The full dimer TDDFT spectrum, however, is not even in agreement with the dipole coupling model based on deficient BP monomer data. The charge-transfer problem/wrong asymptotic behavior of the XC potential creates many additional spurious excitations, along with exciton CD from coupling among these and more physical transitions. These problems render the BP dimer CD spectrum altogether unreliable. For truly long-range coupling of excitations, it is possible that any 'pure' functional component with a wrong asymptotic behavior causes problems akin to those found for BP. This would also affect popular functionals such as B3LYP or PBE0 which afford only 20 and 25 % HF exchange, globally, in the exchange part of the functional. In comparison, the LC hybrid functionals appear unproblematic, as expected.  Figure 10 displays MM spectra obtained from BP and LC-PBE0-g * next to an experimental CD spectrum for BTX-D (see Computational Details). Based on the intensity of the simulated CD spectra, considering that computed intensities may typically deviate from experiment by factors of two, the magnitude of the BTX-D circular dichroism is consistent with exciton coupling of the TPP Soret (B) transitions at distances of 40-50 . Additional features seen in the experimental BTX-D spectrum but not in the MM(LC-PBE0-g * ) spectrum may be potentially attributed to excitations in the steroid backbone, vibronic fine structure of the excitations, the conformational flexibility of BTX-D, and possibly interactions of BTX-D with the solvent. The calculated spectra have been broadened empirically, with a uniform Gaussian width independent of the excitation energy. The full TDDFT dimer spectra in Figure 9 as well as the MM spectra indicate that the exciton couplet from the Soret bands can be nonconservative, for instance due to coupling of additional transitions with the B transitions. Moreover, if the broadening of the spectrum increases significantly at higher energies (shorter wavelengths), the observed trough and peak heights of the exciton CD couplet would differ. It is possible that such mechanisms lead to the much weaker observed CD intensity above 3 eV (below about 420 nm) in the experimental spectrum, compared to the negative CD band between 2.9 and 3 eV. However, other factors cannot be ruled out at this point. A full modeling of the spectrum would require exploration of the conformational space of BTX-D along with solvent effects, which is beyond the scope of this work. We note in passing that force-field based dynamics and conformer searches may bias calculations toward certain configurations. [76] Careful benchmarks, possibly in comparison to ab initio molecular dynamics, appear to be in order.
An interesting point arising from the comparison in Figure 10 is the apparently 'best' agreement of the MM(BP) spectrum with the experimental CD of BTX-D. Based on the analysis presented above, the low-energy tail in the MM(BP) spectrum is generated by exciton coupling of spurious pre-Soret N-type transitions in the monomer. Therefore, this calculated MM(BP) spectrum looks more similar to experiment than the other dimer spectra presented in this section for entirely wrong reasons. More accurate calculations with LC functionals expose this problem.
In the remainder of this section we consider additional aspects of the dimer exciton CD that can be extracted from the dipole coupling model. To facilitate the discussion, consider first a simplified MM setup with two degenerate uncoupled transitions, E a = E b = E, at a distance r, with identical electric transition dipole magnitudes d a = d b = d. For a situation where the two dipole vectors form angles of a a and a b , respectively, with r ab as shown in Figure 11, and adopt a dihedral angle of b along r ab , the potential V ab = V ba = V works out to be as given by Equation (6).
The excitation energies for the coupled system are E 1/2 = E AE V. The corresponding rotatory and dipole strengths are given by Equations (7 a) and (7 b).
As expected, the rotatory strengths change sign as the dihedral angle changes sign because of the sin b term. The sin a a/b terms also indicate that the rotatory strengths vanish if one of the dipoles points along the inter-chromophore vector.  The CD couplets seen experimentally for TPP dimers have been assigned to coupling between the intense Soret transitions, which is consistent with our dimer calculations. The transition dipole vectors for the B transitions lie in the porphin plane. Assume a parallel TPP stack, that is a a = a b = 90 8. If only two transitions were coupled, one B x or B y from each monomer, the dihedral angle adopted in the dimer model of Figure 7 would produce a couplet with the opposite sign of the full TDDFT or MM spectra, and opposite to experiment. A MM coupling model for both monomer transitions, B x and B y , results in a 4 4 Hamiltonian matrix. As additional data demonstrate (see the Supporting Information), for a 'stackeddisk' arrangement, the exciton CD vanishes if the monomer transitions are degenerate or nearly degenerate. As long as the excitations are degenerate, the CD couplet is negligible in a one-side-tilted stacked-disk arrangement (a a = 90 8 and a b ¼ 6 90 8). A substantially stronger exciton CD is created for relative monomer orientations with both a angles deviating substantially from 90 8 and, optionally, a substantial energetic splitting between the B x and B y transitions of the monomer. These conditions are fulfilled for the dimer model shown in Figure 7. Therefore, the single TPP dimer conformer used to represent an approximate average BTX-D structure affords a physically reasonable setup, geometrically and electronically as far as the TPP moieties are concerned, for a study of the long-range exciton CD of porphyrin-based chromophore dimers.

Conclusions
Tetraphenylporphyrin (TPP) has been used in this work as a representative for a strongly absorbing chromophore for which measurements of long-range exciton circular dichroism (CD; at interchromophore distances of up to 50 ) have been previously reported. Two main questions were addressed: 1) How accurate is an exciton coupling model such as the matrix method (MM) with input from time-dependent density functional theory (TDDFT) for the description of exciton CD at such spatial chromophore separations? 2) What types of density functionals are best suited for full TDDFT calculations of longrange exciton CD, and what functionals are best suited to generate reliable input data for a lower-level model such as the MM? In the context of these questions, an 'optimal tuning' of hybrid functionals with range-separated exchange has been addressed as well.
The results were found to be unambiguous: For long-range exciton coupling of large p chromophores, the use of a fully long-range corrected hybrid functional with range-separated exchange is beneficial for both 1) and 2). When reliable input data are generated from first principles methods, the matrix method based on a simple electric-dipole coupling is seen to produce good-quality exciton CD spectra for TPP dimers at spatial separations inferred from experimental data. For chromophores with extended p systems, it appears to be highly beneficial to optimally tune the range-separation parameter, as it is strongly system-specific. For the TPP test case, the tuning produced a much smaller range-separation parameter than typically used in global parametrizations of range-separated hybrids. The tuning resulted in a more accurate absorption spectrum, and significantly reduced DFT delocalization error as evidenced by small curvatures of E(N) when plotted for electron numbers (N) around the neutral system. These findings are in line with recent work by Baer, Kronik, et al. on the excitations of extended p systems [39,42] and recent findings by our group on the performance of range-separated hybrids in calculations of the optical rotation of helicenes. [41] The disastrous dimer TDDFT results obtained with the 'pure' functional BP demonstrate that systems with spatially separated chromophores potentially pose problems for all functionals that retain a fraction of approximate local DFT exchange asymptotically. Even a MM spectrum based on BP data for the TPP monomer was found to be strongly deficient, due to exciton coupling in the dimer of spurious TPP N bands appearing energetically below the B transitions. Curiously, this spectrum agrees 'best' with an experimental TPP exciton coupling CD spectrum reported in Ref. [11] for a brevetoxin derivative.