{'fraction_non_alphanumeric': 0.026874373690443654, 'fraction_numerical': 0.017946615650906442, 'mean_word_length': 4.624594363791631, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 6, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Electronically Excited States in Poly(p-phenylenevinylene): Vertical Excitations and Torsional Potentials from High-Level Ab Initio Calculations

Ab initio second-order algebraic diagrammatic construction (ADC(2)) calculations using the resolution of the identity (RI) method have been performed on poly-(p-phenylenevinylene) (PPV) oligomers with chain lengths up to eight phenyl rings. Vertical excitation energies for the four lowest π–π* excitations and geometry relaxation effects for the lowest excited state (S1) are reported. Extrapolation to infinite chain length shows good agreement with analogous data derived from experiment. Analysis of the bond length alternation (BLA) based on the optimized S1 geometry provides conclusive evidence for the localization of the defect in the center of the oligomer chain. Torsional potentials have been computed for the four excited states investigated and the transition densities divided into fragment contributions have been used to identify excitonic interactions. The present investigation provides benchmark results, which can be used (i) as reference for lower level methods and (ii) give the possibility to parametrize an effective Frenkel exciton Hamiltonian for quantum dynamical simulations of ultrafast exciton transfer dynamics in PPV type systems.


INTRODUCTION
Poly-(p-phenylenevinylene) (PPV, Scheme 1a) plays a paradigmatic role in understanding the electronically excited states and the excitation energy transfer (EET) in conjugated organic molecules to be used in electronic devices aiming at applications in photovoltaics and electroluminescence. 1−5 Understanding the EET at a molecular level is of fundamental importance for a successful design of efficient photovoltaic devices. Recent experiments have shown fascinating features of coherent PPV dynamics 5 and a rationalization in terms of structural dynamics along the PPV chain has been suggested. 6 It is well documented that a structural relaxation in connection with the evolution of the electronically excited states plays an important role for the EET dynamics. 7 These processes are coupled in a complex way because at least several electronic states are involved and nonadiabatic effects have to be considered due to crossings of these states depending on their structural evolution.
Electronic structure calculations have the potential to clarify the above-addressed questions but are facing severe problems because of the large molecular sizes to be handled in addition to the already extremely difficult task of calculating reliably the required electronic states. In view of this situation, several semiempirical calculations have been performed using the collective electronic oscillator (CEO) method combined with the Austin Model 1 (AM1) 8 applied to PPV and to polyphenylene ethynylene units 9 and the Pariser−Parr−Pople (PPP) π electron Hamiltonian 10 for surface hopping dynamics. The vibronic structure of the lowest optical transitions in PPV has been studied using a monoexcited configuration interaction for the electronic transition and an empirical description of the electron phonon coupling. 11 The torsional dependence of the potential energy surface for the first excited state around the vinylene single bonds was investigated by means of Zerner's intermediate neglect of differential overlap (ZINDO) method 12 and comparison was made with time-dependent density functional theory (TDDFT) results. The TDDFT method has been used also to study the singlet−triplet splitting in oligomers of α-thiophenes, p-phenylenes, PPV, and ladder-type oligophenylenes. 13 Furthermore, the localization of the electronic excitation along the polymeric PPV chains in dependence of different functionals was investigated. 14 In a subsequent publication by the same authors 15 environmental polarization has been included noting that the localization of charged polarons was sensitive to environmental polarization whereas neutral states were less affected. Following a different approach, a diabatic Hamiltonian for application to wavepacket dynamics in PPV driven by torsional modes along the vinylene single bonds and the bond length alternation has been developed. 16 Intrachain vs interchain energy transfer dynamics has been studied in conjugated polymers on the basis of Forster-type approaches 17  Ab initio investigations are much more costly than semiempirical or TDDFT calculations and, therefore, the focus concentrated on the smaller oligomers of PPV starting with complete active space perturbation theory to second-order (CASPT2) calculations on the electronic spectrum of stilbene 19 to symmetry-adapted cluster-configuration interaction (SAC-CI) investigations up to four PV units. 20 Torsional motions around the vinylene single and double bonds and the location of conical intersections have been investigated at complete active space self-consistent field (CASSCF) and CASPT level for stilbene. 21 Whereas these methods have been shown to provide reliable data for electronically excited states of πconjugated oligomers, their scope in terms of oligomer sizes is quite limited because of the drastically growing computational cost with increasing chain length. To satisfy the urgent demand for ab initio investigations for significantly larger oligomer length, the approximate coupled cluster method to secondorder (CC2) 22 has been shown to be very useful in calculations on methylene-bridges oligofluorenes 23 and oligo-p-phenylenes. 24 The combination with the resolution of the identity (RI) method 25 allowed computationally efficient calculations on several excited states and the availability of analytic energy gradients for excited states 26 proved to be an especially interesting feature. A systematic benchmark investigation on the capabilities of the CC2 method can be found in ref 27. The related second-order algebraic construction method (ADC(2)) method 28,29 gives similar results in comparison to CC2; 29,30 it has, however, the advantage over CC2 that the excited states are obtained as eigenvalues of a hermitian matrix whereas in coupled-cluster response the excitation energies are obtained as eigenvalues of a non-Hermitian Jacobi matrix. 29 It is, however, to be expected that the second-order character of the method can lead to artifacts. Recent benchmark equation of motion excitation energy coupled-cluster (EOMEE-CC) calculations on DNA nucleobases 31 and comparison with CC2 led to the conclusion that CC2 reproduced π−π* excitations remarkably well in comparison to the higher-level EOMEE-CC approach with singles and doubles and noniterative triples (CCSD(T)) method. CC2 showed noteworthy deficiencies primarily in describing n−π* states, a type of excitation not relevant for the UV spectrum of PPV. Because of the aforementioned similarity of ADC(2) to CC2 it can be expected that the former method will show similar good performance for π−π* states.
One major goal of this work is to provide reliable and consistent information on the excited state properties of PPV oligomers such as vertical excitation energies, torsional potentials, and defect localization in the S 1 state for extended oligomer sizes using the aforementioned RI-ADC(2) method. Even though the TDDFT method would be clearly preferred because of its computational efficiency, we decided for the ADC(2) method because it is free of any ambiguities of choosing an appropriate density functional. Such choices between different functionals will be critical because it has been shown that the defect localization (self-trapping of the exciton) strongly depends on the amount of Hartree−Fock exchange included in the selected functional. 14 PPV consists of alternating phenylene and vinylene units, which are connected through π-conjugation. The structure of the excited states is determined by the interplay between states deriving from phenylene and vinylene units as well as by defects in the molecular structure. The nature of the excited states in PPV is still a matter of debate. Thus, our second goal is to get more insight into the electronic mechanism of interaction between different subunits and to characterize the electronic nature of the excitons using previously developed methods based on the analysis of transition density matrices. 32 For that purpose two types of oligomers representing PPV have been chosen: (i) the phenyl end-capped version denoted (PV) n P and (ii) the vinylene end-capped oligomer denoted (PV) n . Here P and V denote phenyl and vinyl units, respectively. The purpose of these two choices was to investigate the effect of the different chain terminations on geometry and UV spectra. Aside from the electronic nature of the excitons, geometric distortions are important. The phenomenon of exciton−phonon coupling and its relation to exciton localization has been studied by several groups. 10,33−35 Moreover, it has been pointed out that defects breaking the electronic coupling are of highest interest. 5,36,37 Torsional coordinates around the junctions between the vinylene and phenylene subunits have been studied in this work as they play a prominent role as has been shown in simulations of the vibrational broadening of the UV spectrum 12 and in the photodynamics following the absorption process. 10,38 Additionally, we have investigated the bond length alternation (BLA) 10 computed from bond length differences for the interring junction and within the phenylene units relevant for the trapping of defects. Finally, this work is also aiming at the preparation of the basis to fit parameters to our ab initio data to be used in diabatic Hamiltonians 16,39 and subsequent wavepacket dynamics simulations.

COMPUTATIONAL DETAILS
All excited state calculations were performed at the ADC(2) level. 28,29 Ground state geometry optimizations were carried out using the Møller−Plesset perturbation theory to secondorder (MP2). 40 In both cases the RI approximation 29,41−43 and the SV(P) 44 basis was used. All vertical excitations were computed using the SV(P) and TZVP 45 basis sets; for the smaller oligomers the TZVPP basis 41 was used as well. For (PV) 7 P and (PV) 8 torsional potential curves and the geometry optimization in the S 1 state were performed with the split valence SV basis for reasons of computational economy by removing the polarization functions from the SV(P) basis set.
For a compact characterization of the molecular structure the following two BLA parameter definitions were used following ref 10 Here the d i R refer to bond distances within a ring R whereas the d i J refer to bond lengths in the vinyl junction J. The labeling of the bonds within one unit is shown in Scheme 2.
The symmetry of the planar (PV) n P structures is C 2h and that of (PV) n is C s . Rigid torsional potential curves (i.e., in the absence of geometry optimization for the remaining coordinates) in C 2 symmetry were computed for (PV) n P by considering the four torsions labeled T 1 , T 2 , T 3 , and T 4 (Scheme 3) individually. The torsions were defined as linear combinations of individual torsional angles τ as given, e.g., for T 1 as τ 23,24,25 + τ 23,24,25′ . The indices characterizing the angles refer to the bond numbers shown in Scheme 3. For simplicity, only τ 23,24,25 is given in the flowing text to characterize T 1 . Analogous definitions were adopted for the torsions T 2 to T 4 . To preserve C 2 symmetry, the symmetry equivalent bond on the other side of the molecule was twisted as well. The same definitions (without C 2 symmetry equivalence) were used also for (PV) 8 .
An extrapolation of the S 1 excitation energy to infinite chain lengths was performed according to the formula with n defined in Scheme 1. E(N) is the excitation energy at this chain length. The two parameters E 0 and α represent the excitation energy of an isolated oscillator and the coupling strength between neighboring oscillators, respectively. After taking the square of eq 3, the resulting formula may be readily used for a linear regression analysis against cos(π/(N + 1)). The excited states were analyzed according to a recently developed scheme 32 based on previous work by Tretiak and Mukamel 48 and Luzanov and Zhikol. 49 For that purpose the transition density matrix D 0α[AO] between the ground state and excited state α expressed in the atomic orbital (AO) basis is considered. Furthermore, the system is partitioned into fragments A, B, ... Then the charge transfer number from fragment A to B for the transition to state α is computed by summation over basis functions a and b located on the respective fragments A, B, ... where the overlap matrix S [AO] is used to account for nonorthogonality of the AOs in the sense of a Mulliken population analysis. For more details see ref 32. The charge transfer numbers take a more concrete meaning when the excited state is viewed as an electron-hole pair with respect to the ground state. Then Ω AB α can be understood as the probability of simultaneously finding the hole on fragment A and the electron on fragment B.
In the present study the analysis was carried out on the basis of PV fragments, where the formal cuts were performed through the vinyl double bonds. This choice was taken to have analogous fragments representing the overall symmetry of the molecule. Following ref 32 the ADC(2) singly excited cluster amplitudes were chosen to represent the transition density matrix because they dominate the expansion of the ADC(2) wave function.

RESULTS AND DISCUSSION
3.1. Vertical Excitations. The molecular structure of the (PV) n P (n = 1−7) oligomers is presented in Scheme 1a. The planar ground state structures were optimized at the MP2/ SV(P) level using C 2h symmetry followed by the calculation of the four lowest vertical singlet excitation energies (1 1 B u , 2 1 A g , 2 1 B u , 3 1 A g ). The excitation energies for the oligomer series computed at ADC(2)/TZVP level are plotted in Figure 1a as a function of the inverse number of double bonds N (eq 4). In the range from n = 4 to n = 7 an almost linear decrease between ΔE and 1/N is found for all states. The S 1 (1 1 B u ) energy was fitted and extrapolated according to eq 5 using the data for (PV) n P (n = 2, ..., 7) with parameters E 0 = 9.775 eV and α = 0.456. This fit works remarkably well with a squared correlation shows that the three lowest excited states S 1 −S 3 lead to practically the same asymptotic values for N → ∞, which suggests that these states derive from the same exciton band. In contrast, the fourth state (3 1 A g ) stays distinct at higher energies even at longer chain lengths with a polymer limit of 3.72 eV. Comparison of the fit linear in 1/N with the Kuhn fit (eq 5) for the S 1 state shows the well-known underestimation of the asymptotic limit by the former approach. 47 For the first members n = 2 to n = 4 in the (PV) n P series the energy curve for the 2 1 B u state is almost horizontal and starts with a linear decrease only for n ≥ 5.
The oscillator strengths are plotted in Figure 1b. This figure shows that the lowest excited state (1 1 B u ) carries almost all the oscillator strength and that there is a strong increase with growing chain length. This increase in the oscillator strength is linear at the beginning of the series (n = 3−5) and continues with a significantly enhanced rise for larger n. For longer oligomers (starting with n = 5) also the 2 1 B u state possesses a nonnegligible, but still comparatively small, oscillator strength. The A g states are dark for symmetry reasons. The direction of the transition dipole moments is shown in Figure 1S of the Supporting Information. They are in parallel to the molecular plane, pointing along the long molecular axis.
To explore the basis set effect of the vertical excitation energies in more detail, calculations using larger basis sets have been performed. The results are collected in Table 1. For n ≤ 4, SV, SV(P), TZVP, and TZVPP basis sets have been used whereas for the larger oligomers (n = 5, 6, and 7) calculations have been carried out only with the SV, SV(P), and TZVP basis sets. Increasing the size of the basis set quite generally decreases the vertical excitation energies of all four excited states in relation to the ground state. The energy decreases by about 0.3 eV for all oligomers and all four states computed when going from SV to TZVP. It should be noted at this point that the SV basis does not contain any d-functions on C and that the SV(P) and SVP sets include one d-set on C. Extension to TZVPP (inclusion of an additional d-set and on one f-function) decreases the excitation energies further by about 0.1 eV for the n = 1, 2, 3, and 4 cases computed. Even though the basis set dependence of vertical excitation energies with respect to the ground state amounts to about 0.4−0.5 eV, it is also noted that the difference in the energies between the excited states remains almost constant within ∼0.1 eV for all basis sets investigated (e.g., see the results for n = 7 in Table 1). This is an important finding justifying the use of the relatively small SV basis for comparing the different excited states.
The computed vertical excitation energies show quite good agreement with the experimental results available for n = 1 to n = 4. 50 Except for n = 1 the TZVPP results are lower than the experimental values by 0.1−0.25 eV. The increase of the basis set will certainly increase the difference somewhat more. The value for the gasphase vertical excitation obtained from extrapolation of experimental oligomer spectra in solution is 3.25 eV (as computed from ref 47, Table 1, as ΔE vert = E 00 + ΔE eq (abs) − ΔE solv ). Our value extrapolated to infinite chain length using the TZVP basis is 2.91 eV. The increase of the basis set to TZVPP for n = 1−4 reduced the vertical excitation  energy by another 0.1 eV and additional, smaller basis set effects could reduce the vertical excitation energy further, giving in summary an estimated value of 2.6−2.7 eV at ADC(2) level in the complete basis set limit. The discrepancies to the experimental value of 3.25 eV can be partly ascribed to inadequacies present in the ADC(2) approach, but possibly also to the reduction of the conjugation length in the PPV polymer due to kinks in the chains resulting from the flat S 0 potential (see below) around the planar structures and a corresponding increase in the observed excitation energy.

Geometry Relaxation in Ground and Excited
States. To characterize the geometries optimized for S 0 and S 1 , the BLA parameters as defined in eqs 1 and 2 were used. The results for the junctions and rings are shown in Figure 2a,b, respectively, and the underlying bond distances are displayed in Figure 2S (Supporting Information). In the ground state (Figure 2a), all the vinylene junctions show a significant bond length alternation (e.g., d 4 ∼ 0.10 Å), indicating the pronounced difference between the single and double bonds. In the excited state d J decreases significantly toward the center of the chain with a minimal value for d 4 ≈ 0.025 Å. The d R BLA values for the phenylene units ( Figure 2b) show a different trend. In the ground state the bond length alternation is very small (0.02 Å), representing undistorted phenyl rings. In the S 1 state a pronounced quinoid distortion, shortening of the bonds 5 and 6 and an increasing the bonds 1−4 in the phenylene units, is found. Again, this distortion is pronounced at the center, which can be understood as a trapping of the exciton. A similar picture of defect localization in the S 1 state was given by Sterpone and Rossky 10 on the basis of combined PPP/force field calculations. The TDDFT study of Nayyar et al. 14 reflects the strong dependence of the defect localization on the density functional used. The general gradient approximation (GGA) functional PBE does not result in any significant defect confinement; increasing the amount of Hartree−Fock exchange leads to enhanced localization. The present ADC(2) results confirm the strongly localized picture with an extension of the trapping over about four repeat units.
The minimum to minimum S 0 /S 1 transition is 3.246 eV, as compared to the vertical excitation energy of 3.445 eV ( Table 1, SV(P) basis). The resulting reduction of 0.199 eV is somewhat smaller than the experimental estimate ΔE eq of 0.28 eV. 47 3.3. Torsional Potentials. Besides the BLA, the torsional interring modes around the vinylene single bonds are being considered as most important for the excitonic coupling of different PPV units. 16 Rigid torsional potential curves have been computed for the ground and the first four excited states of (PV) 7 P along the four torsional angles, T 1 , T 2 , T 3 and T 4 (Scheme 3). The results were obtained by varying the torsional angles between the adjacent rings in 30°steps. Ground state torsions are shown in Figure 3. All four torsional curves are very similar. For T 1 , T 2 , and T 3 the planar geometry displayed in Scheme 3 (torsional angle 180°) is the most stable one, showing a very flat potential surface until 150°. For T 4 , the geometry with a torsional angle of 150°is most stable by a small margin of 0.011 eV. In all cases the maximum lies at 90°w ith a barrier of 0.30 eV for T 1 and the other barrier heights within 0.03 eV.
Torsional potential curves computed for the first four excited states are shown in Figure 4 for the angles T 1 −T 4 . For all torsions and for all four excited states the planar geometry with a dihedral angle of 180°is the most stable one. The barrier to the torsional rotation is located at an angle of 90°in all cases. At the planar geometries the four excited states investigated are well separated. But as the torsional angle changes, for T 1 and T 2 the energies of the S 1 and S 2 states and S 3 and S 4 states approach each other. At 90°they are pairwise almost degenerate. For T 3 , the S 2 and S 3 states become almost degenerate at 90°and for T 4 the same situation happens for the S 3 and S 4 states while at the same time the other two states remain separate from each other. A rationalization of this situation will be given on the basis of the characteristics of the orbitals investigated in the following paragraph and the discussion in section 3.4.  The two highest occupied and the two lowest unoccupied orbitals of the planar (PV) 7 P oligomer are depicted in Figure 5.
Inspection of this figure shows that the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) have their maximal amplitudes in the middle of the chain whereas the HOMO−1 and LUMO+1 orbitals have two regions of maximal amplitudes in each of the halves of the chain. From Table 1S of the Supporting Information it is observed that for the planar geometry the lowest electronic excitation is dominated by a HOMO → LUMO transition followed by a smaller contribution from HOMO−1 to LUMO1. On twisting the planar geometry around any of the torsional angles, the participation of the HOMO−1 to LUMO+1 transition increases and the contribution of the HOMO → LUMO transition simultaneously decreases. Finally, for 90°, the contributions from HOMO−LUMO and HOMO−1 → LUMO+1 transitions are almost the same for all of the four torsional angles. Thereafter, the importance of the HOMO → LUMO transition slowly increases again. Plots of the frontier orbitals for T 1 (Figure 6) and T 2 ( Figure 3S, Supporting Information) at 90°show for the four orbitals HOMO−1 to LUMO+1 two regions separated by an extended nodal area in the center. For T 3 ( Figure 4S, Supporting Information), the HOMO and LUMO are again concentrated in the center whereas the (HOMO−1 and LUMO+1) pair is concentrated at the chain ends. In case of T 4 ( Figure 5S, Supporting Information), the HOMO and LUMO are again localized in the middle of the chain and there is negligible participation of the phenyl rings at the ends of the chain.
To examine the effect of the addition of an extra terminal vinyl group, rigid T 1 potential energy curves have been calculated for (PV) 8 . The resulting curves (Figure 7) are very similar to the ones obtained for torsion T 1 in (PV) 7 P ( Figure  4). Again, the first excited state has its global minimum at a planar geometry and the states S 1 and S 2 approach each other along the torsion and become closest at 90°. Thus, removing    the equivalence of the two chain halves in (PV) 7 P due to the C 2 symmetry by introducing the terminal vinyl group has only a minor influence on the torsional potentials.
The effect of geometry relaxation has been tested for the T 1 rotation in S 0 and S 1 . These curves were computed by fixing the torsional angle T 1 and optimizing the remaining geometry parameters using a C 2 symmetry restraint. Results are displayed in Figure 8. For the ground state, the geometry with a dihedral angle of 30°is the most stable configuration. The stabilization with respect to the planar structure is, however, quite small and amounts to 0.03 eV. The barrier height with respect to the 90°b arrier is 0.25 eV. For comparison, the barrier computed from the rigid rotation is 0.30 eV. The structure with the dihedral angle of 150°is only 0.01 eV higher in energy than the most stable structure at 30°. For the first excited state, the planar geometries are the most stable ones and the overall shape of the torsional curve is similar to the one of the rigid T 1 rotation ( Figure 4). The entire relaxed S 1 torsional curve is stabilized due to the geometry optimization in comparison to the one computed for the rigid rotation. This stabilization amounts to 0.105 eV for the planar structure at 180°and to 0.07 eV at 90°. As a result, the rotational T 1 barrier is increased by the geometry optimization by a0.035−0.335 eV.
3.4. Character of the Excited States. To get a more detailed insight into the excited states, an analysis of the charge transfer numbers (eq 3) was carried out. The results are presented in Figure 9. In each of these diagrams the position of the hole and electron are plotted along the vertical and horizontal axes, respectively, with eight boxes in each direction representing the eight PV units. In these boxes the value of Ω AB α , which amounts to the probability of simultaneously finding the hole on fragment A and the electron on fragment B, is coded in grayscale. The length of the exciton is seen along the diagonal of the plot, whereas the electron hole separation (charge transfer character) is represented by the off-diagonal width. At the planar equilibrium geometry the three lowest excited states can be seen as rather tightly bound excitons with small off-diagonal width. The Ω AB α plot for S 2 differs from that for S 1 by a region of low amplitude in the center. This difference originates from a nodal plane in the excitonic wave function. Similarly, the Ω AB α plot for S 3 indicates two nodal planes. An interpretation of this nodal progression has been given by Wu et al. 51 in terms of a particle-in-a-box model. Thus, these three states can be seen as belonging to the same exciton band that would arise at infinite chain length. This fact is also consistent with the observation made above that the excitation energies of these three states converge to the same value at infinite chain length ( Figure 1). The torsion around T 1 at an angle of 90°creates a break in the middle of the chain (see the MOs in Figure 6) dividing the molecule into two equivalent weakly coupled fragments. S 1 /S 2 and S 3 /S 4 form two pairs of states of very similar excited state structures. These pairs could each be seen as the positive and negative linear combinations of two localized fragment states. The small gap between each of these pairs of states (Figure 4, T 1 ) can be directly identified with the splitting due to excitonic coupling. For the next two torsions (T 2 , T 3 ) at 90°, the oligomer is effectively split into three parts, two border units and the center unit. For the 90°s tructure around T 2 , S 1 and S 2 can be identified as a pair of excitons on the outer fragments, showing a respective small energy splitting in Figure 4. The Ω AB α plots for S 3 and S 4 exhibit a completely different excited state structure. This suggests that their near-degeneracy in Figure 4 is only accidental. For the case of T 3 the near-degenerate S 2 and S 3 states can be approximately identified as a pair of excitons. In the case of T 4 , only the two outer PV units are twisted and thereby decoupled. This leads to a central conjugated planar region of the remaining chain. Therefore, the first four excited states have similar appearances to the planar equilibrium structure. The T 4 structure at 90°can be seen as (PV) 5 P with two additional Figure 8. Relaxed potential energy curves for the ground state and the first excited states of (PV) 7 P for the torsional angle T 1 . The ground state energy minimum at τ = 30°is taken as reference. Figure 9. Plots of the charge transfer numbers Ω AB α (eq 6) for the first four singlet excited states of (PV) 7 P at its planar geometry and after 90°torsion around the T 1 , T 2 , T 3 , and T 4 angle, respectively. The axis of each square numbers the eight phenylenevinylene units (going from left to right and bottom to top, respectively). The value of Ω AB α is coded in grayscale tone. decoupled end groups and the low splitting between S 3 and S 4 can be understood to derive just from the same type of accidental degeneracy that was present for these two states in planar (PV) 5 P (Figure 1).

CONCLUSION
The lowest four π−π* transitions in the series of PPV oligomers have been investigated by means of the ab initio RI-ADC(2) method. It has been shown that this method is well suited to describe vertical excitations within the π system of PPV and that also excited state geometry optimizations can be performed reliably. In contrast to the TDDFT approach, 14 ADC(2) does not suffer from an intrinsic bias toward defect localization or delocalization in the excited state of PPV. Moreover, due to the resolution of the identity approach, efficient calculations of sufficiently large oligomers sizes could be performed, which are of direct use for describing polymer properties for the excited state. Two important types of structural parameters characterizing the excited state relaxation of PPV have been investigated in detail. The BLA displays the characteristic quinoid distorsions found by means of geometry relaxation in the lowest excited singlet state of conjugated polymers, which could be quantified at high computational level by our calculations. 8,10,24 The torsional modes give insight into the excitonic coupling for several excited states, which has, to our knowledge, not been investigated in detail so far. Combined with the analysis of transition density matrices, an overall picture of the relation of these excited states has been given that should provide a useful benchmark for calculations using lower level, but computationally more efficient, methods.
The present results show the most likely local relaxation pathway after Franck−Condon excitation, the self-trapping 8 of the exciton in the S 1 state by BLA relaxation as shown in Figure  3. However, vibrational energy relaxation is not immediate, and coherent BLA oscillations may be observed on a time scale of tens to hundreds of femtoseconds. 52 The EET dynamics, which will compete with the trapping process, can be anticipated from the splitting of the exciton on torsions along the vinylene C−C bonds ( Figure 9) and the concomitant excitonic couplings. A strong coupling between electron and phonon motions can be expected both for the trapping and for the EET process. The analysis of transition densities indicates for the planar (PV) 7 P molecule the emergence of an excitonic band structure. A significant change in the excited states followed from torsions around the vinylene junctions. Consequently, the qualitative picture of the excitonic interactions depends significantly on the location of the torsion within the oligomer. In the extreme case of a 90°torsion around the central vinylene bond, electronic decoupling and the formation of pairs of equivalent excitonic states could be observed. For torsions around the other vinylene C−C bonds a more complex pattern deriving from the presence of three decoupled fragments was found. In a subsequent work, 39 the present electronic structure results will be employed to parametrize an effective Frenkel exciton Hamiltonian. 39 This ab initio-based diabatic model Hamiltonian will be used as the basis for quantum dynamical simulations of ultrafast exciton transfer dynamics in PPV type systems.  Figure 1S, computed bond lengths in Figure 2S, and MO plots for 90°t orsions around T 2 −T 4 are depicted in Figures 3S−5S. This material is available free of charge via the Internet at http:// pubs.acs.org.