Exciton Coupling and Conformational Changes Impacting the Excited State Properties of Metal Organic Frameworks

In recent years, the photophysical properties of crystalline metal-organic frameworks (MOFs) have become increasingly relevant for their potential application in light-emitting devices, photovoltaics, nonlinear optics and sensing. The availability of high-quality experimental data for such systems makes them ideally suited for a validation of quantum mechanical simulations, aiming at an in-depth atomistic understanding of photophysical phenomena. Here we present a computational DFT study of the absorption and emission characteristics of a Zn-based surface-anchored metal-organic framework (Zn-SURMOF-2) containing anthracenedibenzoic acid (ADB) as linker. Combining band-structure and cluster-based simulations on ADB chromophores in various conformations and aggregation states, we are able to provide a detailed explanation of the experimentally observed photophysical properties of Zn-ADB SURMOF-2: The unexpected (weak) red-shift of the absorption maxima upon incorporating ADB chromophores into SURMOF-2 can be explained by a combination of excitonic coupling effects with conformational changes of the chromophores already in their ground state. As far as the unusually large red-shift of the emission of Zn-ADB SURMOF-2 is concerned, based on our simulations, we attribute it to a modification of the exciton coupling compared to conventional H-aggregates, which results from a relative slip of the centers of neighboring chromophores upon incorporation in Zn-ADB SURMOF-2.


Introduction
Over the past two decades metal-organic frameworks (MOFs) have attracted significant attention due to their structural versatility and multitude of potential applications. They consist of metal/metal-oxo nodes connected by organic linker molecules and form highly regular, crystalline, and porous networks [1]. The huge variety of conceivable nodes and linkers has led to the synthesis of tens of thousands of different systems with specifically tuned properties [2]. Besides their established applications in areas like catalysis [3][4][5] gas storage [6][7][8] and gas separation [9,10] in recent years also their electronic and optical properties have attracted considerable interest [11][12][13][14][15]. Chromophores incorporated into MOFs have already been explored in sensing [16,17] and artificial light harvesting [18][19][20]. From a more Figure 1. Schematic structure of Zn-ADB SURMOF-2; grey squares represent the Zn-paddlewheel, which is bridged by four carboxylic acids of the ligand (see bottom left panel, Zn ions grey); blue rods represent the anthracene-dibenzoic acid linker, where the acidic H is shown for completeness; the carboxylic linkage groups in the paddle wheel and in the molecular skeletal formula are indicated with the same color code: black: carbon, red: oxygen; the tetragonal unit cell is marked with black lines, the unit cell axes are labelled as used in the text.
In the following, we will use quantum-mechanical simulations to provide an in-depth analysis of the excited state properties of ADB both as isolated chromophore and when π-stacked in a welldefined manner inside a coordinatively linked MOF. In particular, we will address the questions, (i) why there is a weak red-shift in the absorption spectrum despite the primarily H-aggregate type exciton coupling in the MOF and (ii) what triggers the red-shift of 0.24 eV of the MOF emission compared to the emission of individual chromophore in solution, which increases the energy difference between absorption and emission maxima from 0.45 eV in solution to 0.64 eV in the MOF. To answer these questions, we will analyze (a) the impact of the solvent, (b) the role played by the packing-motif in the MOF (H-vs. J-aggregates and the consequence of the relative arrangement of the chromophores), (c) the impact of changes in the molecular conformation, and (d) the nature of the excited states dominating absorption vs. emission processes in the MOF. In the following discussion, we will focus on the purely electronic properties of the excited states, acknowledging that for explaining certain details of the absorption and emission spectra of molecular aggregates one also has to include electron-phonon coupling effects [52], as, for example, reviewed for aggregates of paradistyrylbenzene by Gierschner et al. [53].

Molecular Properties and Solvent Effects
As a first step, we analyze the situation for an isolated ADB molecule. The transition to the lowest excited singlet state, S1, is optically allowed and the state is primarily described by an excitation from the highest occupied molecular orbital, HOMO, to the lowest unoccupied molecular orbital, LUMO. Both orbitals are largely localized on the central anthracene unit. This is also observed for the transition density associated with the S0→S1 transition (see Figure 2). It can be explained by a nearly complete breaking of conjugation between the anthracene core and the attached phenylenes, which is a consequence of an almost perpendicular arrangement of the π-planes in the respective units (with a twist angle of 84°). Consequently, when comparing the lowest excitation energies in anthracene and In the following, we will use quantum-mechanical simulations to provide an in-depth analysis of the excited state properties of ADB both as isolated chromophore and when π-stacked in a well-defined manner inside a coordinatively linked MOF. In particular, we will address the questions, (i) why there is a weak red-shift in the absorption spectrum despite the primarily H-aggregate type exciton coupling in the MOF and (ii) what triggers the red-shift of 0.24 eV of the MOF emission compared to the emission of individual chromophore in solution, which increases the energy difference between absorption and emission maxima from 0.45 eV in solution to 0.64 eV in the MOF. To answer these questions, we will analyze (a) the impact of the solvent, (b) the role played by the packing-motif in the MOF (H-vs. J-aggregates and the consequence of the relative arrangement of the chromophores), (c) the impact of changes in the molecular conformation, and (d) the nature of the excited states dominating absorption vs. emission processes in the MOF. In the following discussion, we will focus on the purely electronic properties of the excited states, acknowledging that for explaining certain details of the absorption and emission spectra of molecular aggregates one also has to include electron-phonon coupling effects [52], as, for example, reviewed for aggregates of para-distyrylbenzene by Gierschner et al. [53].

Molecular Properties and Solvent Effects
As a first step, we analyze the situation for an isolated ADB molecule. The transition to the lowest excited singlet state, S 1 , is optically allowed and the state is primarily described by an excitation from the highest occupied molecular orbital, HOMO, to the lowest unoccupied molecular orbital, LUMO. Both orbitals are largely localized on the central anthracene unit. This is also observed for the transition density associated with the S 0 →S 1 transition (see Figure 2). It can be explained by a nearly complete breaking of conjugation between the anthracene core and the attached phenylenes, which is a consequence of an almost perpendicular arrangement of the π-planes in the respective units (with a twist angle of 84 • ). Consequently, when comparing the lowest excitation energies in anthracene and in the isolated ADB molecule, one obtains an only rather small red-shift of 0.15 eV (see Table 1). The reason why there is any shift at all is a minor spread of the excited state onto the phenylenes, as can be inferred from the shapes of the frontier orbitals and the transition density (see Figure 2). Consistently, the electron withdrawing carboxylic acid substituents have an only very weak impact on the optical properties of ADB. This can be concluded from the data for diphenylanthracene (DPA) (see Table 1), which essentially coincide with those for ADB. Finally, it should be mentioned that the transition dipole for the lowest excited state in all molecules discussed in this paragraph is parallel to the short molecular axis of the anthracene moiety (cf. Figure 2).
Molecules 2020, 25, x FOR PEER REVIEW 4 of 20 in the isolated ADB molecule, one obtains an only rather small red-shift of 0.15 eV (see Table 1). The reason why there is any shift at all is a minor spread of the excited state onto the phenylenes, as can be inferred from the shapes of the frontier orbitals and the transition density (see Figure 2). Consistently, the electron withdrawing carboxylic acid substituents have an only very weak impact on the optical properties of ADB. This can be concluded from the data for diphenylanthracene (DPA) (see Table 1), which essentially coincide with those for ADB. Finally, it should be mentioned that the transition dipole for the lowest excited state in all molecules discussed in this paragraph is parallel to the short molecular axis of the anthracene moiety (cf. Figure 2).  Table 1. Experimental and calculated absorption and emission data for ADB in different solvents. For the sake of comparison, also gas-phase absorption data for diphenylanthracene (DPA) and anthracene-dibenzoic acid (ADB) are shown. The calculations have been performed at the PBE0/def2-TZVP//def2-SVP level and the reported values correspond to a purely electronic transition between the ground and the first excited state (S0 and S1, respectively). The oscillator strengths of the S0→S1 transitions are given in parentheses. Absorption (emission) data have been calculated for the S0 (S1) equilibrium geometries. The experimental values for the absorption and emission peaks in ethanol and for the MOF have been extracted from the experimental spectra shown in the Supplementary Material, which correspond to the spectra shown in [51] (including the respective Supporting Material). Here, in a first approximation, the calculated (vertical) excitation energies are associated with the positions of the maxima of the experimental spectra [54]. The position of the absorption peak in toluene was taken from [55]. For the sake of comparison with calculated data in later sections, the table also contains the experimental peak positions for Zn-ADB SURMOF-2 (i.e., chromophores incorporated into the MOF).   Table 1. Experimental and calculated absorption and emission data for ADB in different solvents. For the sake of comparison, also gas-phase absorption data for diphenylanthracene (DPA) and anthracene-dibenzoic acid (ADB) are shown. The calculations have been performed at the PBE0/def2-TZVP//def2-SVP level and the reported values correspond to a purely electronic transition between the ground and the first excited state (S 0 and S 1 , respectively). The oscillator strengths of the S 0 →S 1 transitions are given in parentheses. Absorption (emission) data have been calculated for the S 0 (S 1 ) equilibrium geometries. The experimental values for the absorption and emission peaks in ethanol and for the MOF have been extracted from the experimental spectra shown in the Supplementary Materials, which correspond to the spectra shown in [51] (including the respective Supporting Material). Here, in a first approximation, the calculated (vertical) excitation energies are associated with the positions of the maxima of the experimental spectra [54]. The position of the absorption peak in toluene was taken from [55]. For the sake of comparison with calculated data in later sections, the table also contains the experimental peak positions for Zn-ADB SURMOF-2 (i.e., chromophores incorporated into the MOF).

Absorption Experiment Calculation
Anthracene ( The strong localization of the S 1 state in the apolar center of the ADB molecule also explains, why including solvents of varying polarity in the calculations has hardly any impact on the lowest excited state (see Table 1). This theoretical finding is also consistent with experimental studies on DPA, which show that decreasing the polarity of the solvent from ε = 24 (ethanol) to ε = 2.4 (toluene) increases the excitation energy by at most 0.05 eV [55][56][57]. As a consequence, solvent effects are ignored for the present system and will not be accounted for in the following discussion.
As far as the emission properties of ADB are concerned, we calculate an energy of 2.78 eV for the purely electronic S 1 →S 0 transition (calculated for the S 1 equilibrium geometry). This corresponds to a rather large shift of~0.55 eV between the absorption and emission maxima. The magnitude of the shift can be explained by a reduced twist angle of the phenylene rings relative to the anthracene in the excited state equilibrium geometry (56 • in S 1 vs. 84 • in S 0 geometry), which results in an increase of conjugation. Overall, comparing the results in Table 1 shows that for the isolated molecule the simulations agree exceptionally well with the experiments. Thus, in the following we will focus on understanding the situation, when the ADB linkers are incorporated into the MOF.

Formation of H-and J-Aggregates in Zn-ADB SURMOF-2: Anthracene Dimers as Model Systems
For understanding the excited states of ADB incorporated into the MOF, as a first step, it is useful to analyze possible exciton couplings between the chromophores in the MOF structure. Depending on the relative alignment of the transition dipoles on adjacent chromophores, H-type or J-type aggregates are formed in a solid-state assembly. In H-type aggregates, the transition dipoles on adjacent chromophores are aligned in parallel, whereas J-type aggregates are formed with a head-to-tail alignment of the transition dipoles. In both aggregate types, the first excited state is shifted to lower energies compared to the isolated chromophore, but only in J-type aggregates this excited state is optically allowed, resulting in a red-shift of the absorption spectrum. In H-type aggregates, the first optically allowed state is typically found at energies higher than the parent state in the isolated chromophore, which results in a blue-shift [58].
For the following, discussion of the fundamental aspects of exciton coupling in Zn-ADB SURMOF-2 we will replace the ADB chromophores by anthracene units to simplify the interpretation. This does in no way affect the key conclusions regarding exciton coupling that will be discussed in this section and is justified by the strong localization of the excited state on the anthracene moiety in the isolated ADB molecule (vide supra).
Both aggregate types discussed above can be found in Zn-ADB SURMOF-2, as schematically shown in Figure 3: In a-direction, one observes J-type (green molecules) as well as H-type aggregates (blue molecules). The center-to center distance between the chromophores in both cases amounts to 19.80 Å. This is so large that the effect of exciton coupling becomes vanishingly small, resulting in a negligible splitting between the lowest excited states (see Table 2). For symmetry reasons, exactly the same situation as in a-direction is also found in b-direction (where now neighboring blue molecules form J-aggregates and neighboring green molecules form H-aggregates). The situation is fundamentally different in c-direction, where all nearest-neighbor molecules form H-aggregates. As in this direction the molecules are rather close, the impact of the coupling becomes noticeable. Thus, in the following we will be exclusively concerned with H-aggregates stacked in c-direction. Notably, here the center to center distance of 5.81 Å in the equilibrium structure of the MOF corresponds to a distance of 3.78 Å between neighboring π-planes of the anthracenes, as illustrated in Figure 3d. This is a consequence of a tilting of the anthracenes, which also results in a slip of the centers of neighboring chromophores parallel to the π-planes by 4.42 Å.  The distance between neighboring π-planes and the slip of the chromophores is also shown. Table 2. PBE0/def2-TZVP//def2-SVP calculated absorption energies and oscillator strengths (in parentheses) for the most relevant low-energy excitations for a single anthracene molecule, as well as for dimer and tetramer clusters along the three unit cell directions. For the H-aggregates, the S1 state and the lowest excited state with appreciable oscillator strength are listed. For the J-aggregate, besides the optically allowed S1 state, we also included the next state described by equivalent single-particle excitations (see text). The anthracene units in the model stacks are geometry-optimized at the singlemolecule level and then aligned in the same way (center-of-mass distance and tilt angles) as in the MOF optimized using periodic boundary conditions (cf., Figure 3d). For the c-stacked, H-aggregate-type dimer of anthracene, the splitting between the lowest excited state and the first state with appreciable oscillator strength amounts to 0.39 eV (see Table 2). Increasing the number of interacting chromophores somewhat increases the splitting (see last entry in Table 2), but considering the rather short distance between the neighboring π-planes, this splitting  Table 2. PBE0/def2-TZVP//def2-SVP calculated absorption energies and oscillator strengths (in parentheses) for the most relevant low-energy excitations for a single anthracene molecule, as well as for dimer and tetramer clusters along the three unit cell directions. For the H-aggregates, the S 1 state and the lowest excited state with appreciable oscillator strength are listed. For the J-aggregate, besides the optically allowed S 1 state, we also included the next state described by equivalent single-particle excitations (see text). The anthracene units in the model stacks are geometry-optimized at the single-molecule level and then aligned in the same way (center-of-mass distance and tilt angles) as in the MOF optimized using periodic boundary conditions (cf., Figure 3d). For the c-stacked, H-aggregate-type dimer of anthracene, the splitting between the lowest excited state and the first state with appreciable oscillator strength amounts to 0.39 eV (see Table 2). Increasing the number of interacting chromophores somewhat increases the splitting (see last entry in Table 2), but considering the rather short distance between the neighboring π-planes, this splitting still appears rather small. Moreover, in a typical H-aggregate one would expect the oscillator strength for the excitation into the S 1 state to be exactly zero for symmetry reasons.

Energy/eV (Oscillator Strength)
To rationalize these findings, as a first step the excited state structure of a cofacial (i.e., not slipped) anthracene model dimer shall be discussed: In the dimer, hybrid orbitals are formed from the molecular HOMOs and LUMOs of each of the molecules. In a single-particle picture, four excitations between these orbitals are possible (see Figure 4). For symmetry reasons, two of these excitations are optically allowed and two of them are forbidden (solid vs. dashed arrows in Figure 4). In the actual time dependent density functional theory (TD-DFT) calculations the single-particle excitations mix. This yields four excited states, which can be characterized by linear combinations of either the allowed or the forbidden single-particle excitations. This is shown in the first entry of Table 3 (slip 0.0 Å) for the two cofacial anthracene molecules. In the following, the four excited states will be denoted as S a , S b , S c , and S d . Here, S a refers to an excited state dominated by the positive linear combination of the forbidden single-particle transitions OS→UA and OA→US. In this context, O denotes to the highest occupied and U to the lowest unoccupied orbital of the dimer with a specific symmetry, where S and A specify, whether the dimer orbitals are symmetric (i.e., positive) or antisymmetric (i.e., negative) linear combinations of the orbitals of the individual molecules (see Figure 4). S b refers to the positive linear combination of the allowed transitions OA→UA and OS→US and S c and S d denote the negative linear combinations of the respective transitions (see last column of Table 3). In line with the involved single-particle excitations, transitions to states S a and S d are strictly optically forbidden (see oscillator strengths in Table 3), while excitations to states S b and S c are, in principle, optically allowed, although in our simulations the oscillator strength for excitations into S b are consistently much smaller than for excitations into S c . As a consequence, the position of S c determines the position of the first peak in the absorption spectrum.
Molecules 2020, 25, x FOR PEER REVIEW 7 of 20 still appears rather small. Moreover, in a typical H-aggregate one would expect the oscillator strength for the excitation into the S1 state to be exactly zero for symmetry reasons.
To rationalize these findings, as a first step the excited state structure of a cofacial (i.e., not slipped) anthracene model dimer shall be discussed: In the dimer, hybrid orbitals are formed from the molecular HOMOs and LUMOs of each of the molecules. In a single-particle picture, four excitations between these orbitals are possible (see Figure 4). For symmetry reasons, two of these excitations are optically allowed and two of them are forbidden (solid vs. dashed arrows in Figure  4). In the actual time dependent density functional theory (TD-DFT) calculations the single-particle excitations mix. This yields four excited states, which can be characterized by linear combinations of either the allowed or the forbidden single-particle excitations. This is shown in the first entry of Table  3 (slip 0.0 Å) for the two cofacial anthracene molecules. In the following, the four excited states will be denoted as Sa, Sb, Sc, and Sd. Here, Sa refers to an excited state dominated by the positive linear combination of the forbidden single-particle transitions OS→UA and OA→US. In this context, O denotes to the highest occupied and U to the lowest unoccupied orbital of the dimer with a specific symmetry, where S and A specify, whether the dimer orbitals are symmetric (i.e., positive) or antisymmetric (i.e., negative) linear combinations of the orbitals of the individual molecules (see Figure 4). Sb refers to the positive linear combination of the allowed transitions OA→UA and OS→US and Sc and Sd denote the negative linear combinations of the respective transitions (see last column of Table 3). In line with the involved single-particle excitations, transitions to states Sa and Sd are strictly optically forbidden (see oscillator strengths in Table 3), while excitations to states Sb and Sc are, in principle, optically allowed, although in our simulations the oscillator strength for excitations into Sb are consistently much smaller than for excitations into Sc. As a consequence, the position of Sc determines the position of the first peak in the absorption spectrum.    For the cofacial anthracene dimer (first entry in Table 3), the expected situation for a conventional H-aggregate is recovered: The optically forbidden S a state is lowest in energy and the first state with appreciable oscillator strength (S c ) lies 0.78 eV above the first excited state. To understand the different properties of the anthracene dimer in the H-aggregate configuration adopted in Zn-ADB SURMOF-2 discussed before, one has to consider the slip of the centers of neighboring chromophores in the MOF (Figure 3d): As discussed by Kazmaier and Hoffmann for simple model systems and for perylene [59] (and later found for a variety of organic semiconductors [60][61][62]), displacing the centers of coplanar molecules relative to each other results in a periodic variation of the splitting of the respective hybrid orbitals as a function of the displacement. This occurs due to the symmetry of the individual orbitals (see Figure 5a). Moreover, because of the decrease of the spatial overlap of the molecules with increasing displacement, the amplitude of the oscillations decreases. These changes of the orbital energies cause also variations in the energies of the above-mentioned four excited states, as shown in Figure 5b. These variations do not directly coincide with the variations for the orbital energies, which has two reasons: First, due to the different nodal patterns of the HOMO and LUMO (compare Figure 4), the slips at which the splitting between the HOMO and the HOMO − 1 vanishes differs from the slips at which the same occurs for the LUMO and LUMO + 1. Second, the dominant single-particle excitation describing specific excited states changes with the displacement, as shown by the filling of the symbols in Figure 5 (for details see figure caption).
Notably, the slip-induced changes in excitation energies are large enough for the order of the states to change as a function of the displacement. In fact, for a slip of 4.42 Å (the value obtained in Zn-ADB SURMOF-2, indicated by a vertical line in Figure 5) the optically weakly allowed S b state comes to lie lowest in energy instead of the strictly symmetry-forbidden S a state. This is reminiscent of the situation in crystals of dicyanodistyrylbenzene based molecules, where static symmetry-breaking renders the lowest excited state of H-aggregate coupled chromophores optically allowed [63]. Moreover, the splitting between the lowest excited state and the first state with appreciable oscillator strength (S c ) decreases by essentially a factor of two between the cofacial dimer (slip 0.0 Å) and the dimer in the Zn-ADB SURMOF-2 configuration (4.42 Å, see Table 3). This explains the somewhat unexpected excited state properties of the anthracene-dimer extracted from the Zn-ADB SURMOF-2 structure (vide supra). In passing we note that especially the change of the order of the states due to the slip will become relevant later, when discussing the emission properties of Zn-ADB SURMOF-2.
Molecules 2020, 25, x FOR PEER REVIEW 9 of 20 anthracene units. In particular, one has to study to what degree the conformations of the actual ADB units change upon incorporation into Zn-SURMOF2. Analogous dependence of the energies of the two lowest symmetryforbidden excitations (to states Sa and Sd) and symmetry-allowed excitations (to states Sb and Sc). Sa, Sb, Sc and Sd are described by superpositions of one-particle excitations according to Table 3. As detailed in the main text, Sa and Sd are dominated by the same single-particle excitations. The same applies to Sb and Sc. Therefore, these pairs are plotted using the same symbols. Filled symbols specify a larger weight for the OS→UA (in the case of Sa and Sd) and for the OA→UA single particle excitations (in the case of Sb and Sc), while open symbols specify that the respectively other single particle excitation dominates (OS→UA for Sa and Sd and OS→US for Sb and Sc). The vertical line indicates the slip of 4.42 Å that occurs in Zn-ADB SURMOF-2 (see Figure 3d).

Impact of the Chemical Linkages and The Solid-State Conformation on the Optical Properties of ADB in Zn-ADB SURMOF-2
As a first step towards answering that question, the impact of the bonding of the ADB chromophore to the metal nodes is assessed. For this purpose, we optimized a single ADB molecule suspended between two Zn-paddle wheels (pw) with the Zn and O atoms fixed to the positions they adopt in the periodic structure of the MOF discussed below. The Zn-nodes are additionally saturated by three acetate groups per paddlewheel (see Figure 6a). This structure in the following will be referred to as opt (pw-ADB-pw)1, where the subscript denotes that only a single ADB unit is considered and the superscript refers to a full geometry optimization (where only the Zn and O atoms are fixed). This geometry optimization yields a structure very similar to the isolated ADB molecule with bond angles within 0.5 degrees and bond lengths within 0.02 Å (cf. Supplementary Material, SI.2). The changes in bond-lengths are primarily triggered by fixing the Zn-Zn distance between the paddle wheels. Also, the nature and energy of the lowest excited state are very similar to those in the isolated molecule (cf. , Tables 1 and 4). This supports the finding from above that terminal substituents have essentially no impact on the lowest excited states of the chromophores. It also implies that for the present combination of chromophore and metal node, there is no "through-bond" electronic coupling (like in certain electrically conductive MOFs [14]) between adjacent chromophores in a-and bdirection. Analogous dependence of the energies of the two lowest symmetry-forbidden excitations (to states S a and S d ) and symmetry-allowed excitations (to states S b and S c ). S a , S b , S c and S d are described by superpositions of one-particle excitations according to Table 3. As detailed in the main text, S a and S d are dominated by the same single-particle excitations. The same applies to S b and S c . Therefore, these pairs are plotted using the same symbols. Filled symbols specify a larger weight for the OS→UA (in the case of S a and S d ) and for the OA→UA single particle excitations (in the case of S b and S c ), while open symbols specify that the respectively other single particle excitation dominates (OS→UA for S a and S d and OS→US for S b and S c ). The vertical line indicates the slip of 4.42 Å that occurs in Zn-ADB SURMOF-2 (see Figure 3d).
Still, there is one aspect of the calculations on the anthracene H-aggregates which is at variance with the experimental observation for Zn-ADB SURMOF-2: The simulations predict a blue shift of the absorption maximum (by 0.12 eV for the tetramer in Table 2 compared to the isolated molecule), while in the experiments a minor red shift of the absorption peak by 0.05 eV is observed (see Table 1). To understand that, one has to go beyond representing the chromophores in the simulations by anthracene units. In particular, one has to study to what degree the conformations of the actual ADB units change upon incorporation into Zn-SURMOF2.

Impact of the Chemical Linkages and The Solid-State Conformation on the Optical Properties of ADB in Zn-ADB SURMOF-2
As a first step towards answering that question, the impact of the bonding of the ADB chromophore to the metal nodes is assessed. For this purpose, we optimized a single ADB molecule suspended between two Zn-paddle wheels (pw) with the Zn and O atoms fixed to the positions they adopt in the periodic structure of the MOF discussed below. The Zn-nodes are additionally saturated by three acetate groups per paddlewheel (see Figure 6a). This structure in the following will be referred to as opt (pw-ADB-pw) 1 , where the subscript denotes that only a single ADB unit is considered and the superscript refers to a full geometry optimization (where only the Zn and O atoms are fixed). This geometry optimization yields a structure very similar to the isolated ADB molecule with bond angles within 0.5 degrees and bond lengths within 0.02 Å (cf. Supplementary Material, SI.2). The changes in bond-lengths are primarily triggered by fixing the Zn-Zn distance between the paddle wheels. Also, the nature and energy of the lowest excited state are very similar to those in the isolated molecule (cf., Tables 1 and 4). This supports the finding from above that terminal substituents have essentially no impact on the lowest excited states of the chromophores. It also implies that for the present combination of chromophore and metal node, there is no "through-bond" electronic coupling (like in certain electrically conductive MOFs [14]) between adjacent chromophores in aand b-direction.  Table 4. PBE0-D3/def2-TZVP//def2-SVP calculated absorption properties for the monomer opt (pw-ADB-pw)1, the dimer cut (pw-ADB-pw)2, and the tetramer cut (pw-ADB-pw)4 of node-linker-node units in MOF conformation stacked along c-direction. The superscripts in the system names denote, how the cluster has been generated. "opt" refers to a full geometry optimization, while "cut" refers to a system, generated by cutting out the central moiety (viz. tet2) of the fully optimized tetramer (see Figure 6c) and replicating it as a dimer or tetramer employing the periodicity of the bulk MOF. The description of the excited states and molecular orbitals in the dimer is analogous to that of the anthracene dimers in Table 3 to highlight the resemblance to the model system. Sa, Sb, Sc, and Sd are also the four lowest excited states of cut (pw-ADB-pw)2. For the tetramer, all transitions with oscillator strengths greater than 0.15 are shown; there, a nomenclature analogous to the dimer is no longer straightforward. To categorize the molecular orbitals nonetheless, the tetramer was grouped into three pairs of dimers and the symmetry of each pair is given. Again, S denotes a symmetric linear combination and A an antisymmetric one.  One of the crucial aspects not captured by the opt (pw-ADB-pw) 1 model system is, how the conformation of the chromophore is changed by the neighboring linkers in cdirection. To capture the influence of the neighbors, as a first step we optimized the structure of the 3D MOF employing periodic boundary conditions. For reasons explained in the Methods Section, this has been done using the PBE functional [64]. From the periodic structure we extracted a tetramer cluster repeated in c-direction (consisting of four ADB molecules bonded to two saturated Zn-paddle wheels). This cluster was then further optimized with the PBE0 functional [65] (like in the molecule-based simulations), fixing the positions of the Zn atoms. This is done to deal with geometries (in particular bond lengths) obtained at a level of theory consistent with the previously discussed simulations. The structure of that tetramer cluster, opt (pw-ADB-pw) 4 , is shown in Figure 6c. In this cluster, the geometries of the outermost (pw-ADB-pw) units are impacted by edge effects, but comparing the structures shown in Figure 6b,c, one sees that the two central units adopt an arrangement fully consistent with the periodic geometry optimizations. A more quantitative analysis of the results shows that all twist angles agree to within less than 1 • and that bond lengths are within 0.004 Å compared to the values from the periodic simulations (for more details see Supporting Information SI.4). Thus, we used one of these units (viz tet2 in Figure 6c) as the basic building block for the model systems used in the following. The structures constructed using this "cut-out" monomer will be denoted as cut (pw-ADB-pw) n , where the index n denotes the number of repeating units. Table 4. PBE0-D3/def2-TZVP//def2-SVP calculated absorption properties for the monomer opt (pw-ADB-pw) 1 , the dimer cut (pw-ADB-pw) 2 , and the tetramer cut (pw-ADB-pw) 4 of node-linker-node units in MOF conformation stacked along c-direction. The superscripts in the system names denote, how the cluster has been generated. "opt" refers to a full geometry optimization, while "cut" refers to a system, generated by cutting out the central moiety (viz. tet2) of the fully optimized tetramer (see Figure 6c) and replicating it as a dimer or tetramer employing the periodicity of the bulk MOF. The description of the excited states and molecular orbitals in the dimer is analogous to that of the anthracene dimers in Table 3 to highlight the resemblance to the model system. S a , S b , S c , and S d are also the four lowest excited states of cut (pw-ADB-pw) 2 . For the tetramer, all transitions with oscillator strengths greater than 0.15 are shown; there, a nomenclature analogous to the dimer is no longer straightforward. To categorize the molecular orbitals nonetheless, the tetramer was grouped into three pairs of dimers and the symmetry of each pair is given. Again, S denotes a symmetric linear combination and A an antisymmetric one. To better illustrate that, he OAAA and UAAA orbitals for cut (pw-ADB-pw) 4 are visualized in the Supplementary Material (SI.8). The geometric changes of all "cut" structures compared to optimized monomer structure, opt (pw-ADB-pw) 1 , arise from the impact of the neighboring ADB chromophores inside the MOF.

System
These changes primarily concern the twist of the anthracene and phenylene moieties of ADB relative to the plane in which the neighboring Zn atoms are arranged (compare Tables S3 and S4 in the Supplementary Materials). In the optimized monomer, opt (pw-ADB-pw) 1 , the phenylenes are essentially in the plane of the Zn atoms (see Figure 6a). Such a conformation is prevented in the 3D periodic structure by steric constraints, as can be inferred from the structure shown in Figure 6b. In fact, for phenylenes in the plane of the Zn atoms, the H atoms on neighboring rings would come much too close to each other. Consequently, the phenylenes are twisted by 24 • relative to the Zn plane in the periodic conformation. Also, the orientation of the anthracene units is significantly modified. While the anthracene plane is nearly perpendicular to the Zn-plane in opt (pw-ADB-pw) 1 (at 82 • ), the twist between the two planes is reduced to 42 • in the periodic structure. This reduced anthracene-Zn plane twist is primarily a result of van der Waals interactions between neighboring anthracenes trying to reduce the distance between the π-planes of the molecules (as shown in SI.3 contained in the Supplementary Materials). Most importantly, as a consequence of that also the angle between the phenylene and anthracene planes is reduced from 81 • in opt (pw-ADB-pw) 1 to 66 • in the periodic structure.
To study the impact of the changes in twist angles, we first discuss the properties of a single "cut" chromophore, cut (pw-ADB-pw) 1 , (see highlight in Figure 6c). The energy of the lowest excited state of cut (pw-ADB-pw) 1 is distinctly red-shifted (by 0.17 eV) compared to the fully optimized monomer, opt (pw-ADB-pw) 1 , and the oscillator strength is significantly increased (see Table 4). This can be attributed to the change of the conformation of the chromophore when incorporated into the MOF, where the main aspect is that the reduced twist between the phenylenes and the anthracene in the ADB unit results in an increased conjugation (see Supplementary Material SI.5).

The Final Absorption Spectrum of Zn-ADB SURMOF-2: Combining Conformational Changes and Aggregate Formation
In the previous sections we have shown that when assembling ADB chromophores into Zn-SURMOF-2, there are two competing effects regarding the change in the absorption spectrum: On the one hand, there is a red-shift of the lowest excited state in the MOF due to the conformational changes triggered by inter-linker interactions and a related increase of conjugation (see Section 2.3 and horizontal arrow in Figure 7). On the other hand, H-aggregate formation in c-direction causes a blue-shift of the first strongly allowed state, as discussed in Section 2.2 (see Table 2). This raises the question, how the combination plays out in the actual MOF, in which exciton coupling and conformational changes happen simultaneously. To address that, we constructed a dimer and a tetramer by assembling the "cut" monomer entities described above. This yielded cut (pw-ADB-pw) 2 and cut (pw-ADB-pw) 4 , where the former system has the advantage that its properties in terms of the nature of orbitals and excited states can be discussed in analogy to the situation of the anthracene dimers from Section 2.2. Indeed, it turns out that the nature and order of the excited states in cut (pw-ADB-pw) 2 are equivalent to those of the anthracene dimer with a slip of 4.42 Å (cf., Tables 3 and 4): In the TD-DFT simulations on cut (pw-ADB-pw) 2 , the lowest excited state displays S b character (the negative linear combination of OA→UA and OS→US transitions), while the state with the largest oscillator strength is the third excited state possessing S c character (the respective positive linear combination). Compared to the slipped anthracene dimer, the oscillator strengths of both states are increased, which is due to the spreading of the transition density onto the phenylene units (cf., Supplementary Material SI.5; see also comparison of the properties of anthracene and ADB in Table 1). This effect is a consequence of the reduced twist between the anthracene and the phenylene unit, when the ADB chromophores are incorporated into the MOF. As the oscillator strength associated with the S c state is by a factor of more than five higher than that of the S b state, the energy of S c determines the position of the first absorption peak. Due to the dimer formation, in cut (pw-ADB-pw) 2 this state is slightly blue-shifted by 0.08 eV compared to the corresponding monomer cut (pw-ADB-pw) 1 (see blue arrow Figure 7 and Table 4). Combining this blue shift by 0.08 eV with the red-shift by 0.17 eV between opt (pw-ADB-pw) 1 and cut (pw-ADB-pw) 1 due to conformational changes yields the overall red-shift of 0.09 eV between an isolated ADB molecule and cut (pw-ADB-pw) 2 . This is schematically shown by the black arrow in Figure 7.
A similar situation is obtained when calculating the excited states of the tetramer, cut (pw-ADB-pw) 4 . (see Table 4). Again, the lowest exited state is red-shifted compared to the monomer and states at higher energies dominate the absorption spectrum due to their larger oscillator strengths. The state with the highest oscillator strength amongst the first 20 excited states (S 14 ) is even somewhat further blue-shifted than in cut (pw-ADB-pw) 2 . Calculating a theoretical absorption spectrum from a superposition of Gaussian peaks with full widths at half maximum (FWHM) of 0.30 eV centered at the energies of the excited states of the tetramer and scaled by their oscillator strengths yields an absorption maximum at 3.29 eV (cf. Supplementary Material SI.7). This is only slightly higher than the experimental absorption maximum at 3.27 eV (see Table 1). In passing we note that an unambiguous signature of the weakly allowed S1 state in cut (pw-ADB-pw) 2 and cut (pw-ADB-pw) 4 cannot be identified in the experimental spectra, as discussed in more detail in the Supplementary Material (SI.9). this state is slightly blue-shifted by 0.08 eV compared to the corresponding monomer cut (pw-ADBpw)1 (see blue arrow Figure 7 and Table 4). Combining this blue shift by 0.08 eV with the red-shift by 0.17 eV between opt (pw-ADB-pw)1 and cut (pw-ADB-pw)1 due to conformational changes yields the overall red-shift of 0.09 eV between an isolated ADB molecule and cut (pw-ADB-pw)2. This is schematically shown by the black arrow in Figure 7. A similar situation is obtained when calculating the excited states of the tetramer, cut (pw-ADBpw)4. (see Table 4). Again, the lowest exited state is red-shifted compared to the monomer and states at higher energies dominate the absorption spectrum due to their larger oscillator strengths. The state with the highest oscillator strength amongst the first 20 excited states (S14) is even somewhat further blue-shifted than in cut (pw-ADB-pw)2. Calculating a theoretical absorption spectrum from a superposition of Gaussian peaks with full widths at half maximum (FWHM) of 0.30 eV centered at the energies of the excited states of the tetramer and scaled by their oscillator strengths yields an absorption maximum at 3.29 eV (cf. Supplementary Material SI.7). This is only slightly higher than the experimental absorption maximum at 3.27 eV (see Table 1). In passing we note that an unambiguous signature of the weakly allowed S1 state in cut (pw-ADB-pw)2 and cut (pw-ADB-pw)4 cannot be identified in the experimental spectra, as discussed in more detail in the Supplementary Material (SI.9). Overall, the above considerations explain, why the first absorption peak in Zn-ADB SURMOF-2 is not blue-shifted but rather red-shifted compared to the isolated ADB chromophore in solution in spite of the formation of H-aggregates. What remains to be explained is the significant shift of 0.64 eV between the absorption and emission maxima in Zn-ADB SURMOF-2.

Explaining the Red-Shifted Emission of ADB Molecules Incorporated into Zn-ADB SURMOF-2
For the isolated ADB molecule in solution, the rather large shift of 0.45 eV between absorption and emission maxima could be explained by a reduction of the twist angle between anthracene and phenylene units from 84 • to 56 • in the excited state equilibrium conformation (see Section 2.1). In the MOF that angle is already decreased to 66 • in the ground state, primarily due to the van der Waals interaction between neighboring chromophores (see Section 2.3). Moreover, a further planarization of the ADB linkers incorporated into Zn-ADB SURMOF-2 in the excited state is prevented by steric constraints due to the already tight packing of the ADB chromophores in c-direction in the ground state. Indeed, when optimizing the geometry of the (pw-ADB-pw) 4 tetramer in the S 1 electronic configuration, yielding S1 (pw-ADB-pw) 4 , the changes in tilt angles are only very minor (see Supporting Information SI.7). This applies in particular to the two central pw-ADB-pw units, where the S 1 state is primarily localized, as can be inferred from the excitation-induced changes in bond lengths and from the transition density shown in Figure 8.
As a consequence, one would expect a smaller shift between the absorption and emission maxima in the MOF. Indeed, when calculating the excited state properties of one of the two central pw-ADB-pw units cut from the optimized S 1 tetramer, S1,cut (pw-ADB-pw) 1 , one observes an only rather moderately shift of~0.3 eV (see monomer values in Tables 4 and 5). A similar shift is actually observed when comparing the energies of the lowest excited states for cut (pw-ADB-pw) 4 and S1 (pw-ADB-pw) 4 . The respective energies amount to 2.86 eV for ground-state conformation (see Table 4) and to 2.49 eV for the excited state conformation (see Table 5), where it should be stressed that both states have the same nature, being dominated by an excitation between equivalent orbitals (see Supplementary Material SI.8). This clearly shows that geometric relaxations in the excited state in Zn-ADB SURMOF-2 would result in a shift between absorption and emission maxima of only half the experimentally observed value of 0.64 eV.
state. Indeed, when optimizing the geometry of the (pw-ADB-pw)4 tetramer in the S1 electronic configuration, yielding S1 (pw-ADB-pw)4, the changes in tilt angles are only very minor (see Supporting Information SI.7). This applies in particular to the two central pw-ADB-pw units, where the S1 state is primarily localized, as can be inferred from the excitation-induced changes in bond lengths and from the transition density shown in Figure 8. Figure 8. Changes in bond lengths between the ground-and excited state configurations and shape of the transition density associated with the first excited state calculated to determine the localization of that excited state [66,67]. Panel (a) shows the difference of the bond lengths of bonds B1-B5 between ground and first excited state for (pw-ADB-pw)4 and for an isolated ADB molecule. The labelling of the bond lengths is explained in the inset. All structures were optimized at the PBE0/def2-SVP+D3 level. Panel (b) contains the transition density associated with the first excited state of S1 (pw-ADBpw)4 calculated with PBE0/def2-TZVP//def2-SVP+D3.
As a consequence, one would expect a smaller shift between the absorption and emission maxima in the MOF. Indeed, when calculating the excited state properties of one of the two central pw-ADB-pw units cut from the optimized S1 tetramer, S1,cut (pw-ADB-pw)1, one observes an only rather moderately shift of ~0.3 eV (see monomer values in Tables 4 and 5). A similar shift is actually observed when comparing the energies of the lowest excited states for cut (pw-ADB-pw)4 and S1 (pw- Figure 8. Changes in bond lengths between the ground-and excited state configurations and shape of the transition density associated with the first excited state calculated to determine the localization of that excited state [66,67]. Panel (a) shows the difference of the bond lengths of bonds B1-B5 between ground and first excited state for (pw-ADB-pw) 4 and for an isolated ADB molecule. The labelling of the bond lengths is explained in the inset. All structures were optimized at the PBE0/def2-SVP+D3 level. Panel (b) contains the transition density associated with the first excited state of S1 (pw-ADB-pw) 4 calculated with PBE0/def2-TZVP//def2-SVP+D3. Table 5. PBE0-D3/def2-TZVP//def2-SVP first excited-state properties for the monomer S1,cut (pw-ADB-pw) 1 and the tetramer S1 (pw-ADB-pw) 4 of node-linker-node units stacked along c-direction. The superscripts in the system definition denote, how the cluster has been generated. "S1" denotes a full geometry optimization in the first excited state, while "cut" refers to system generated by cutting out the central moiety, tet2, of the S 1 fully optimized tetramer analogous to Figure 6c. The description of the molecular orbital contributions of the S 1 optimized tetramer is analogous to that in Table 4. This suggests that the main features in the emission and absorption spectra are dominated by different electronically excited states: The key difference between them is that the absorption spectrum is influenced by all excited states and dominated by the state(s) with the highest oscillator strength (S 14 in the case of opt (pw-ADB-pw) 4 ). Conversely, what counts for the emission characteristics according to Kasha's rule are the properties of the lowest excited state, i.e., S 1 . This is in particular the case here, as due to the slip of the anthracene molecules, transitions between this state and the ground state are not optically forbidden in the H-aggregates of Zn-ADB SURMOF-2 (see Section 2.2) [63].

System
Therefore, to assess the shift between the maxima of the absorption and emission spectra, one has to compare the absorption maximum of cut (pw-ADB-pw) 4 (which we find at 3.29 eV as discussed in Section 2.4) and the S 1 energy of S1 (pw-ADB-pw) 4 (of 2.49 eV). This, indeed, yields a red-shift of 0.80 eV, which is consistent with the experimentally observed shift of 0.64 eV. The somewhat larger shift in the calculations occurs not only for Zn-ADB SURMOF-2, but also for ADB in solution (see Table 1). It is mostly due to a minor underestimation of the emission energy, as becomes evident, e.g., from the comparison between the calculated excited state properties and the experimental spectra in the Supplementary Material (SI.9).
In passing we note that even if the lowest excited state was forbidden in absorption, dynamic symmetry breaking due to geometry relaxations in the excited state could relax symmetry-selection rules [63,68]. This is, however, not the case in our TD-DFT simulations, as can be inferred from the reduced oscillator strength of the lowest excited state for the relaxed geometry of S1 (pw-ADB-pw) 4 compared to the ground-state conformation in cut (pw-ADB-pw) 4 . This is potentially a consequence of the delocalization of the exciton over two chromophores.
The above considerations show that the red-shifted emission of Zn-ADB SURMOF-2 is indeed a direct consequence of inter-chromophore interactions, where the subtleties of exciton coupling in slipped chromophores are crucial, while massive, excitation-induced conformational changes, as one would expect in classical excimers, do not play a role.
A final aspect that should be discussed is the possible consequence of the particularly low oscillator strength of 0.06 associated with the S 1 →S 0 emission transition of S1 (pw-ADB-pw) 4 . It implies that the radiative lifetime of Zn-ADB SURMOF-2 should be particularly long, in fact, much longer than the measured overall excited state lifetime of~4 ns [51] (which would then be determined by the non-radiative lifetime). It also suggests that in case there is some inhomogeneity in the sample including some less well-ordered regions, where the chromophores have a more isolated character, these regions will dominate the emission spectrum immediately after excitations due to the much larger oscillator strengths of isolated chromophores (see Table 1). Only, when the excited state populations in these regions have decayed or when the excitons have migrated to the fully crystalline parts of the samples, the red-shifted emission of the S 1 state in S1 (pw-ADB-pw) 4 will dominate. Such a red-shift of the emission with time has, indeed, been observed in the experiments on Zn-ADB SURMOF-2 [51].

Methods
For the investigation of the properties of Zn-ADB SURMOF-2, we performed Density Functional Theory, DFT, calculations on molecules, clusters, and crystalline, 3D periodic structures. All molecule and cluster-based calculations were performed with the ORCA 4.0.1 code [69] employing the PBE0 hybrid functional [65]. In the geometry optimizations of the ground and first excited state, we used the def2-SVP basis set [70] and Grimme's D3 dispersion correction [71]. Absorption and emission properties of the relaxed structures were calculated with the linear-response approach within time-dependent DFT (TD-DFT), employing the Tamm-Dancoff approximation (TDA) [72] and using the def2-TZVP basis set [70]. When studying absorption (emission) properties, the lowest twenty (six) excited states were calculated explicitly. To simulate the ADB molecule in solution, solvent effects for ethanol (ε = 24.3) and toluene (ε = 2.4) were included using the SMD continuum solvation model [73].
The periodic structure of the MOF was studied with the program FHI-Aims [74]. We employed a 2 × 2 × 8 k-point grid and Tight settings of the numeric atom-centered basis set. A detailed description of the used basis functions is given in the Supplementary Material (SI.1). The atomic ZORA correction was applied for treating relativistic effects [75]. Starting from the experimental XRD-data [51], the cell parameters of the tetragonal unit cell were fixed, while the atom positions within the cell were relaxed until the remaining forces were below 10 −3 eV/Å. The use of hybrid functionals in conjunction with periodic boundary conditions is prohibitively expensive for systems as large as the MOF studied here. Thus, in the FHI-Aims simulations we employed the Perdew-Burke-Ernzerhof (PBE) functional [64] with the Tkatchenko-Scheffler van der Waals correction [76], where the influence of the functional is addressed in detail in the Supplementary Materials.
Notably, excited state calculations employing DFT and periodic boundary conditions to date are impossible for systems as complex as the ones studied here. Moreover, to the best of our knowledge, no band structure code exist that would let us relax the relevant excited state geometries. Thus, to obtain consistent results, we used the periodic structures only as reference geometries for determining molecular conformations and for extracting starting structures for cluster optimizations.
Such clusters were built from (several) ADB linkers connected to two Zn-paddlewheels (extracted from the FHI-Aims calculations), with each paddlewheel coordinatively saturated with three acetate groups. During the cluster optimizations employing ORCA 4.0.1 and the PBE0 functional, the Zn and O atoms were kept fixed at the positions obtained in the periodic calculation. The geometries of the linkers were fully optimized. When studying clusters representing units of Zn-ADB SURMOF-2, no solvation model was employed, as experiments suggest a solvent free environment inside the MOF [51].

Conclusions
In this study, we provide a detailed explanation for the excited state properties of anthracene dibenzoic acid (ADB) both as isolated chromophore in solution as well as incorporated as linker into a metal-organic framework (Zn-ADB SURMOF-2). The latter provides a well-controlled arrangement of the chromophores relative to each other in the solid state. The comparison of the two situations (solution and porous solid) is facilitated by the observation that the solvent polarity has virtually no impact on the optical properties of ADB. Within the MOF, H-aggregate type coupling of the ADB entities is identified as the main type of excitonic interaction. Thus, one would expect a distinct blue-shift of the maximum of the absorption spectrum, which is, however, neither observed in the simulations nor in experiments [51]. One of the reasons for that is that in our simulations both the ordering of the excited states as well as the energetic splitting between excited states in the MOF are significantly modified compared to a conventional H-aggregate. This is due to a slip of the centers of neighboring chromophores relative to each other. A second, even more important aspect is that the geometric conformation of the ADB chromophores change considerably inside the MOF due to steric constraints resulting from the rigid network structure. This increases π-conjugation within the ADB linkers and triggers a red-shift of all excited states. The combination of the blue shift due to aggregate formation and the red shift due to conformational changes then leads to the a priori unexpected red-shift of the actual absorption spectrum.
Another striking feature of Zn-ADB SURMOF-2 is the huge shift of 0.64 eV between the maxima of the absorption and emission spectra. This is particularly surprising, as the comparably tight packing of the MOF-linkers in c-direction prevents a significant excitation-induced change of the conformation of the individual chromophores. This is in sharp contrast to the situation in solution, where the reduction of the twist between the anthracene and the phenylene units upon excitation provides a major contribution to the red-shift of the emission. There are also no significant changes in the relative arrangement of neighboring chromophores, which could modify the exciton coupling in the excited state. Instead, the highly red-shifted emission of Zn-ADB SURMOF-2 is a consequence of different excited states dominating the absorption and emission spectra. In absorption, primarily strongly allowed states count, which experience a very minor red-shift for the reasons discussed above. Conversely, for the emission, in line with Kasha's rule, he lowest excited state is the relevant one. Due to the exciton coupling between the strongly interacting linkers in the MOF, this state is strongly red-shifted and as a consequence of the relative slip of the π-planes of neighboring anthracenes, its calculated oscillator strength is small, but non-zero.
These considerations show that it is the subtle interplay of a variety of factors like inter-chromophore couplings, the reordering of states, and packing-induced conformational changes that determine the optical properties of MOFs and that disentangling these factors proofs difficult without performing suitable simulations.