Introduction

Top Abstract Introduction Results and Discussion Conclusion References

The science and technology of electronic and optoelectronic devices based on organic -conjugated polymers is evolving rapidly. Numerous high-performance photonic devices fabricated from organic polymers and molecular crystals have been made (1, 2), including light-emitting diodes (LEDs) (2-8) and electrochemical cells (9), display panels (10-12), photovoltaic cells (13-15), photodetectors (16-18), transistors (19-21), light-emitting field-effect transistors (22), biosensors (23), imaging devices (24, 25), and solid-state lasers (26-29). The rational design of novel plastic materials with enhanced functionalities requires insight into their electronic structure, charge and energy transport, and photoexcitation dynamics. This study has become a research frontier that involves many challenges for theory, experiment, and synthesis; however, a sufficient fundamental understanding is yet to be achieved.

The theoretical description of -conjugated molecular systems is challenging because of electronic correlation effects and strong electron-phonon interactions. It is well understood from the early days of polyacetylene studies (1, 30) that coupling of the electrons to the nuclear degrees of freedom results in a very rich photophysics of solitons, polarons, and excitons and makes an important distinction of "soft" condensed matter from rigid solids based on semiconductor or crystalline metal materials. On the other hand, electron correlation effects have a dramatic impact on the electronic structure, leading, e.g., to the inverse ordering of 2A_g and 1B_u excited states in finite polyenes (31). The early studies conducted for finite systems in real space (32, 33) or for periodic lattices using one-dimensional band theories (k space) (34) were based on simple empirical Hamiltonian models such as Su-Schrieffer-Heeger (SSH) (30, 35) and made an invaluable contribution to our qualitative understanding of such phenomena.

However, an accurate description of molecular chemical structures including defects, disorder, and chemical substitutions is beyond an empirical treatment but possible with quantum-chemical methods. The development of powerful computers and versatile program packages such as the Gaussian suite (36) leads to the growing importance of quantum chemistry for the understanding underlying fundamental electronic phenomena and for providing accurate predictions of relevant observables. For ground-state adiabatic potential surfaces (e.g., geometry optimization, transition-state searches, and molecular and reaction dynamics) these approaches are currently well developed. In most cases, ground-state wavefunctions can be well represented by a single Slater determinant (37), and this approximation allows the treatment of fairly large molecular systems. However, the computation of excited states requires accounting for electronic correlation effects by using, for example, a multiconfigurational representation (38), which is a numerically expensive procedure. Nevertheless, accurate computations of excited-state structures, UV-visible spectra, and nonlinear polarizabilities at the fixed optimal ground-state geometry are currently routine even for large molecules with sizes up to several hundred atoms (38, 39). Configuration interaction combined with semiempirical Hamiltonian models such as the intermediate neglect of differential overlap/spectroscopy (40) or the time-dependent density functional theory (41) are the most popular approaches for excited-states calculations, the former method being less expensive than the latter.

Computation of the excited-state dynamical processes that involve nuclear degrees of freedom still remains a computational challenge and is available only for small molecules. The recently developed excited-state molecular dynamics (ESMD) approach (42-44), schematically illustrated in Fig. 1, represents a reasonable compromise between accuracy and numerical effort and allows adiabatic photoexcitation dynamics to be followed on picosecond time scales in molecules with hundreds of atoms. This approach is based on the collective electronic-oscillator method (39), which combines the random-phase approximation for a many-electron problem (which includes essential electronic correlations) with a semiempirical Hamiltonian model [such as intermediate neglect of differential overlap/spectroscopy (45) or Austin model 1 (46)]. The collective electronic oscillator allows accurate prediction of spectroscopic observables by using fast Krylov-space diagonalization techniques to calculate excited states with minimal numerical effort (39). The ESMD utilizes the collective electronic oscillator approach to calculate gradients of the excited-state potential energy surface (forces) on the fly. These quantum-mechanical forces are used further to propagate classical Newtonian equations of motion for the nuclei to model the time-dependent vibronic evolution (42). This approach therefore is applicable within the Born-Oppenheimer approximation and allows treatment of adiabatic photoexcited dynamics.

View larger version (17K): [in this window] [in a new window]   Fig. 1.   Geometry of cis-polyacetylene and schematic representation of molecular dynamics propagation. The excited-state potential energy Ee(q) as a function of nuclear coordinates q, displacements , vertical absorption A, and fluorescence F frequencies.

In this article we investigate the excited-state electron-nuclear dynamics of a polyacetylene chain (which is a simple -conjugated polymer and serves as a prototype of a one-dimensional electronic system with strong coupling to nuclear degrees of freedom). Our focus is on the formation of photoexcited "breathers," i.e., multiquanta vibronic states. A breather is a dynamic state representing a spatially localized, time-persistent nonlinear excitation. Extensive studies based on SSH-type models predicted the formation of breather-like states in polyenes after either photoexcitation (32, 33, 47-49) or electronic ground-state perturbations of sufficient energy (50). A clear physical picture emerges, for example, for photoexcitation dynamics in trans-polyacetylene (which has two degenerate ground-state conformations): An electron-hole pair (exciton) created by an absorbed photon decays to a charged soliton-antisoliton pair. These quasiparticles rapidly become spatially separated on the time scale of a few tens of femtoseconds, whereas a part of the initial photoexcitation energy, which is left behind, forms a localized breather (33, 47). In polymers with a nondegenerate ground state such as cis-polyacetylene, this picture is more complicated: The exciton cannot decay into a fully separated soliton-antisoliton pair because of confinement effects. Instead, the composite dynamics of an electron, a hole and a breather, develops (48). In these models, the size of spatial localization strongly depends on the empirical Hamiltonian parameters (30, 35). These model investigations also attributed peaks near the   * edge region in the photoinduced absorption as a spectroscopic signature of breathers (47, 48, 51). However, it took almost two decades to develop ultrafast sub-5-fs experimental spectroscopic techniques to monitor real-time vibronic dynamics of conjugated polymers that may support the existence of breathers (52).

Although it is straightforward to treat practically any organic molecular structure (within certain size limits) with the ESMD method (39, 42, 43), we focus here on polyacetylene. This choice is for two reasons: (i) the ability to compare these calculations with early studies of breathers calculated within the empirical Hamiltonian models (32, 33, 47-49) and (ii) recent experimental evidence for breather excitations in polyacetylene (52). Our results reproduce all the essential features of the early studies and predict characteristic properties that do not depend on empirical Hamiltonian parameters. In addition, our results provide all the information necessary for modeling a variety of spectroscopic data, which allows us to connect these simulations with recent experiments (52) and to predict additional spectroscopic signatures of breathers. Finally, we study the dynamics in the presence of a structural defect (an acetylene bond) and chemical substitutions (halogenization) to understand the effect of chain imperfections on the photoexcitation dynamics.

	Results and Discussion

Top Abstract Introduction Results and Discussion Conclusion References

We used the Austin model 1 Hamiltonian for all calculations presented in this article. Austin model 1 has been parameterized specifically to reproduce ground-state properties (such as chemical energy, geometry, dipole moment, etc.) in a variety of molecules (46). In addition, this model provides a reasonable accuracy for excited-state properties in combination with the random-phase approximation (53, 54). The calculations start from the optimal geometry of a linear cis-polyacetylene oligomer obtained with the standard GAUSSIAN 98 computational package (36), i.e., the ground-state potential energy E_g(q) is minimal in the space of nuclear coordinates q, which span the 3N  6-dimensional space, with N being the total number of atoms in the molecule. We further use the ESMD computational package to follow photoexcitation adiabatic dynamics up to a 0.5-ps time scale (42). This approach calculates the excited-state potential energy as E_e(q) = E_g(q) + (q). Here (q) is the excitation frequency to the lowest 1B_u (band-gap) state of polyacetylene. The program uses numerical derivatives of E_e(q) with respect to each nuclear coordinate q_i (i = 1, ... , 3N  6) to calculate forces and subsequently to step along the excited-state hypersurface using these gradients. A standard Verlet molecular dynamics algorithm (55) has been used for propagation of the Newtonian equations of motion. No dissipative processes are included here, so the total energy (initial ground-state energy plus the energy of the absorbed photon) is conserved in the simulations. In a separate simulation we introduced an artificial dissipative force into the equations of motion to obtain the minimum of the excited-state potential energy E_e(q) corresponding to the relaxed geometry.

To analyze excited-state dynamics of cis-polyacetylene we calculate oligomers (with free boundary conditions) long enough (n = 120 carbon atoms) to mimic the infinite chain limit [typically spectroscopic observables in polyenes saturate to the bulk limit at 60-80 carbon atom chain lengths (38, 39)]. To monitor the vibrational dynamics we focus on the bond-length alternation parameter, which reflects an uneven distribution of the  electrons over the bonds (Peierls distortion) and is therefore strongly coupled to the electronic system. The bond-length alternation parameter r_n is defined as

[ 1 ]

where l_n is the bond length between the nth and n + 1th atoms along the chain. In addition, we define a local compression parameter s_n as

[ 2 ]

where x_n is the deviation of nth atom from its equilibrium position. The s_n parameter reflects the changes of the lattice constant. r_n and s_n have been defined in a similar fashion for the SSH Hamiltonian (47, 48).

To follow the electronic dynamics we further use a real-space two-dimensional representation of the transition density matrices (39). These matrix elements reflect the changes in electronic density and coherences and therefore characterize the properties of an exciton or bound electron-hole pair created after photoexcitation (39, 56). Transition densities are defined as

[ 3 ]

where (c_j) are creation (annihilation) operators of an electron at the jth atomic orbital, and |g (|) is the ground (excited)-state many-electron wavefunction (39, 56). The diagonal elements ()_jj represent the net charge induced in the jth atomic orbital by the external field. The off-diagonal elements ()_ij with i  j represent the joint probability amplitude of finding an electron and a hole located at the ith and jth atomic orbitals, respectively. Thus, a transition density plot depicts probabilities of an electron moving from one molecular position (horizontal axis) to another (vertical axis) after electronic excitation.

Fig. 2 shows the variation of r_n and s_n, and Fig. 3 displays the correspondent transition densities at specific times up to 100 fs. At the moment of photoexcitation (t = 0) the molecule in the ground state has a constant dimerization (r_n = 0.095 Å) and s_n is vanishing along the polymer chain (i.e., there is no variation in lattice constant). The effects from the molecular ends are small and local. Therefore, after the photoexcitation, the exciton created is initially delocalized along the entire chain (Fig. 3, t = 0). The diagonal direction of the plot shows the coordinate of the exciton's center of mass measured by the delocalization size L_d, which is distributed evenly among all carbon atoms (except the molecular ends) at t = 0. The size of the exciton (maximum distance between an electron and a hole), which is characterized by the off-diagonal extent of the plot (measured by the coherence size L_c), is 20-25 carbon atoms and is much smaller than the oligomer length.

View larger version (28K): [in this window] [in a new window]   Fig. 2.   Variation of bond-length alternation [rn (Å)] and compression [sn (Å)] parameters in (polyacetylene-cis)60 oligomer along the polymer's backbone at t = 0, 16, 32, 44, 60, 70, 84, 100, and  fs after photoexcitation. t = 0 () corresponds to absorption (fluorescence) events for the molecule being at the ground (excited)-state equilibrium geometry.

View larger version (41K): [in this window] [in a new window]   Fig. 3.   (Upper) Contour plots of the transition densities between the ground state (|0) and the first excited 1Bu state (|1) in (polyacetylene-cis)60 oligomer at t = 0, 16, 32, 44, 60, 70, 84, 100, and  fs after photoexcitation. The axes are labeled by the carbon atoms along the chain. The color map is given on the bottom of the color plots. (Lower) Variation of the diagonal (Ld) and off-diagonal (Lc) coherence sizes in carbon atoms with time. Equilibrium values of these quantities at excited-state relaxed geometry (t = ) are also shown. Lc and Ld define, respectively, the sizes of the exciton and localization of the exciton on the polymer backbone.

Because of strong coupling to vibrational degrees of freedom (electron-phonon coupling) an exciton rapidly distorts the lattice in the middle of the chain (Fig. 3, t = 16 fs), locally reducing the bond-length alternation (Fig. 2, t = 16 fs), and localizes itself in this region on the time scale of 20 fs. This exciton self-trapping process is typical for one-dimensional conjugated polymers (30, 42, 43). When an exciton distorts the lattice, dynamical vibrational excitations (phonons) are created, appearing as waves in r_n and s_n on the edges of the exciton potential well (Fig. 2, t = 16 fs). The subsequent dynamics can be qualitatively described as: the phonon "waves" propagate in opposite directions, reflect from the chain ends, and finally meet in the middle of the chain (Fig. 2, t = 32 fs), pulling an exciton out of its well and delocalizing the excitation (Fig. 3, t = 32 fs). The exciton, in turn, attempts to localize again, creating more phonons, i.e., the energy is exchanging between electronic and vibrational degrees of freedom. Subsequently, the magnitudes of r_n and s_n increase, and the nuclear dynamics becomes more complex because more vibrational modes are involved (Fig. 2, t = 44, 60, 70, 84, and 100 fs). In agreement with refs. 47 and 48, we also notice that s_n roughly behaves as a derivative of r_n, even though these are formally unrelated parameters corresponding to slow (acoustic) and fast (optical) phonons, respectively. The electronic dynamics in turn mirrors structural deformations by exciton delocalization in the regions with the lower bond-length alternation (Fig. 3, t = 44, 60, 70, 84, and 100 fs). This variation of diagonal delocalization is a characteristic "breathing" pattern. The overall electron-nuclear dynamics is a nonlinear and nonperiodic (on relevant femtosecond-picosecond time scales) process, because all vibrational degrees of freedom are coupled to electronic structure and to each other. The most significant structural changes occur in the middle of the chain, which can be envisioned as a localized and persistent vibrational excitation or a breather, i.e., we observe coupled nonlinear dynamics of collective electronic (electron and hole) and lattice (breather) excitations. A similar process has been observed in simulations with a simple SSH model (48). However, there is an important distinction from these early studies: The nuclear dynamics of a real polymer is far richer and more complex than the dynamics of one-dimensional lattice. Finally, if dissipation is introduced, the exciton is localized in the middle of the chain in a region with reduced bond-length alternation over 30-carbon-atoms length, which formally corresponds to the long-time limit (Figs. 2 and 3, t =  fs).

To quantitatively characterize the time scales of the exciton dynamics, we calculated the diagonal (L_d) and off-diagonal (L_c) exciton coherence sizes as respective inverse participation ratios (following ref. 39). Variation of L_d and L_c with time is shown in Fig. 3 Bottom. These quantities vary quasiperiodically with typical periods of 34 fs (L_d) and 18 fs (L_c) determined from their power spectra (not shown). We also notice that the magnitude of L_d changes considerably (by 50%), whereas L_c does not vary significantly, i.e., photoexcitation dynamics affects the spatial localization of an electron-hole pair but not the exciton size.

To relate these time scales to molecular vibrations, we calculated vibrational normal modes for the ground-state equilibrium geometry and expanded the excited-state structural changes in the basis set of these coordinates to obtain dimensionless displacements _i (i = 1, ... , 3N  6). We find that many fast and slow nuclear modes have significant displacements and therefore are dynamically coupled to the electronic system. The variation in the displacements for the most pronounced molecular vibrations are shown in Fig. 4. Vibrational modes with frequencies ₁ = 1,825.8 and '₁ = 1,826.8 cm¹ are related to "in-phase" and "out-of-phase" C==C stretches and overall have positive and negative displacements, respectively, corresponding to an elongated double bond. These modes most resemble phonons with k = 0 and 2 momenta in the long-chain limit, respectively, with typical periods of 18.3 fs. Displacements corresponding to these nuclear motions show quasiperiodic variations around their relaxed values. There are several other quasidegenerate C==C stretching modes related to the optical phonon band in the long-chain limit. We observe that C==C stretching modes corresponding to phonons with higher momenta have consequently smaller displacements (note that phonons with odd momenta do not contribute to the dynamics because the inversion symmetry is preserved during the simulations). Another fast nuclear motion corresponding to the CC stretch (₂ = 1,079.8 cm¹) has an overall negative displacement (shorter single bond) and a period of 30.9 fs. In addition, several vibrational modes corresponding to the slow motion with typical periods of a few hundreds femtoseconds are also strongly coupled to the electronic system. Variation of the displacements of two slow modes with frequencies 61 and 101 cm¹ are shown in Fig. 4. The same types of vibrational modes with similar couplings to the electronic system have been identified in a study of excited-state potentials of short polyenes (n = 8-20 carbon atoms) by using a configuration interaction/semiempirical approach (57).

View larger version (22K): [in this window] [in a new window]   Fig. 4.   Variation of dimensionless displacements along selected vibrational normal modes with time. These nuclear motions strongly couple to the electronic degrees of freedom and can be interpreted as C==C (1 = 1,825.8 cm1 and '1 = 1,826.8 cm1) and CC (2 = 1,079.8 cm1) stretches (fast modes), polymer backbone stretch (3 = 101.2 cm1), and wagging motion (4 = 61.2 cm1) (slow modes). Equilibrium values of these quantities at the excited-state relaxed geometry (t = ) are also shown.

As expected, the L_c variation is related to C==C stretches, which reflects the dependence of the exciton size on the bond-length alternation measuring the degree of Peierls dimerization. However, we notice that the L_d period of 34 fs does not correspond to any normal mode with a significant displacement, i.e., the L_d variation is related to a collective localized vibrational excitation (breather) that cannot be represented by a linear vibrational phonon mode.

We next analyze the time dependence of several observables (Fig. 5) to study possible spectroscopic signatures of the vibronic dynamics described above. Fig. 6 shows power spectra that identify the relevant frequency components. We expect that the signatures of C==C stretching motions, which show the strongest coupling to the electronic system (Fig. 4), will be the most pronounced. The excited-state potential energy E_e varies with a period of 9 fs, which is related to a doubled frequency of the C==C stretching motion. This could be interpreted even within a simple model that assumes that the excited-state energy has an elastic harmonic term Ku²/2, with K being the spring constant and u  sin(₁t) the atomic displacement. It then is clear that E_e will vary with a frequency of 2₁.

View larger version (26K): [in this window] [in a new window]   Fig. 5.   Variation of electronic potential energy Ee, transition frequencies of the first five electronic excited states, 1  5, ground 0(0) and excited 1(0) state static polarizabilities, and transition dipole µ01 between the ground state (|0) and the first excited 1Bu state (|1) with time. Equilibrium values of these quantities at the excited-state relaxed geometry (t = ) are also shown.

View larger version (17K): [in this window] [in a new window]   Fig. 6.   Power spectra of electronic potential energy Ee, transition frequency 1, excited-state static polarizability 1(0), and transition dipole µ01.

Variations of transition frequencies are also coupled to C==C stretching motions that change the bond-length alternation with an 18-fs period (Fig. 5). It is instructive to draw a connection with solid-state band theories, which show that the band-gap energy (which corresponds to the ₁ transition frequency) has a component u from the Peierls distortion, with  being the electron-phonon coupling constant and u  sin(₁t). We notice that the magnitude of frequency variations decreases for the higher energy electronic states, which are more localized as compared with the band-gap state and therefore less sensitive to the bond-length alternation variations. In addition, the power spectrum of ₁ has weaker frequency components at 34 and 70 fs (Fig. 6), which cannot be associated with any specific vibrational mode. Similarly, the ground [(0)] and excited [(1)] state polarizabilities (calculated as the second derivatives of the ground- and excited-state energies, respectively) have an 18-fs period varying with phase opposite to that of the transition frequencies. This can be rationalized by recalling that in a single-state approximation, the polarizability is inversely proportional to the square of the transition energy  = f/², with f being the oscillator strength. The power spectra of polarizabilities also contain 34- and 70-fs components (Fig. 6).

We next calculate the transition dipole moment between ground and excited states as the expectation value of the dipole operator µ: µ₀₁ = Tr(µ₁). This quantity, and subsequently the correspondent oscillator strength (not shown), are dominated by an 34-fs period (see Fig. 5), which corresponds to the breather excitation, and weaker 18- and 70-fs components. For the delocalized band-gap excitations each segment of the polymer chain along the L_d length contributes almost equally to the oscillator strength, which leads to overall scalings µ₀₁ and f₀₁  n. Thus, changes in the diagonal delocalization size L_d cause variation of the transition dipole and the oscillator strength. Vibronic spectroscopy of small molecules is typically described by using the Condon approximation, which implies that the transition dipole does not change in the course of photoexcitation dynamics. However, the significant changes in the transition dipole magnitude we observe in the present case show that the Condon approximation is not always applicable for extended molecular systems. Variation of the transition dipole is an important signature of the breather vibrational excitations, which may be detected spectroscopically.

In summary, the power spectra of all spectroscopic observables (transition frequency, transition dipole, and polarizability) contain clear frequency components with periods of 18, 34, and 70 fs. Only the 18-fs harmonic is related to C==C stretching vibrational normal mode (₁  1,826 cm¹), whereas the other two harmonics are not associated with the vibrational spectrum, with 34 fs being the breather excitation.

To explore the effect of defects and chain imperfections on photoexcitation dynamics we simulated a cis-polyacetylene oligomer that contains an alkynyl (triple-bond) chemical defect. The resulting variations of r_n, s_n, and transition densities are shown in Fig. 7 Upper. We first notice that the triple bond is a very weak defect that does not break the conjugation but locally increases the bond-length alternation. Therefore, initially (Fig. 7, t = 0) an exciton is still delocalized along the entire molecule, which is similar to pristine polyacetylene (compare with Fig. 3, t = 0). However, the dynamics is very different: Immediately an exciton moves to the most conjugated part of the chain, and all subsequent vibronic dynamics, which includes breather formation, occurs on this segment with weak penetration of an excitation to the neighboring segment (Fig. 7). This is the natural relaxation of our photoexcitation to the lowest   * electronic state, which lies in the largest conjugation segment of the chain. The defect provides a strong impedance mismatch, localizing the excited energy. Of course, photoexcitation to higher * electronic orbital is possible with different dynamics. Examination of r_n and s_n confirms this conclusion. We also calculated chains with stronger defects that break the conjugation (such as sp³ defects). The results (not shown) indicate that these defects greatly enhance the localization, effectively breaking the chain into well separated pieces and limiting the dynamics to the longest conjugated region.

View larger version (45K): [in this window] [in a new window]   Fig. 7.   (Upper) The effect of alkynyl (triple-bond) defect in (polyacetylene-cis)50 oligomer on photoexcitation dynamics. (Left) Contour plots of the transition densities between the ground state (|0) and the first excited 1Bu state (|1) at t = 0, 44, 70, and  fs. (Right) Variation of bond-length alternation [rn (Å)] and compression [sn (Å)] parameters along the polymer backbone at t = 0, 44, 70, and  fs. (Lower) The effect of halogenization (chlorine substitutions) in (polyacetylene-cis)50 oligomer on photoexcitation dynamics. (Left) Contour plots of the transition densities between the ground state (|0) and the first excited 1Bu state (|1) at t = 0, 55, 83, and  fs. (Right) Variation of bond-length alternation [rn (Å)] and compression [sn (Å)] parameters along the polymer backbone at t = 0, 55, 83, and  fs. On the color plots the axes are labeled by the carbon atoms along the chain, and the color map is as described in the Fig. 3 legend.

In our final numerical experiment, we modeled the effects of chemical substitutions by chlorinating one third of the chain, i.e., calculating a H(CH)₃₃  (CCl)₁₇Cl oligomer. The resulting variations of r_n, s_n, and transition densities are shown in Fig. 7 Lower. Even though the bond-length alternation is increased along the chlorinated segment (Fig. 7, r_n at t = 0), the electronegativity of Cl effectively attracts an exciton. At t = 0 an excitation is already essentially localized on the substituted end. Subsequent photoexcitation dynamics occurs on the chlorinated segment. We observe strong electron coupling to C==C stretching modes at 1,620 cm¹ (21 fs) on the chlorinated segment; whereas C==C stretching modes at 1,820 cm¹ spatially localized along the unsubstituted part of the chain do not couple to an exciton. Subsequently, the breather period in this system is 40 fs.

	Conclusion

Top Abstract Introduction Results and Discussion Conclusion References

To conclude, we modeled here photoexcitation dynamics of conjugated cis-polyacetylene oligomers using an ESMD quantum-chemical approach. The results show several characteristic features of this dynamics. (i) An excitation moves to the largest conjugated segment of the polymer chain (where it has the lowest electronic energy) by means of short-range coherent localization (42) or long-range Förster transfer mechanism (58). (ii) An exciton creates phonon excitations and significant local distortions of the lattice. The C==C stretching modes exhibit the strongest coupling to the electronic system among all nuclear motions. These phonons form a complex nonlinear localized vibrational excitation (breather) that cannot be associated with any specific single linear vibrational normal mode. The subsequent dynamics is spatially localized and appears as an exchange of the energy between an exciton and a breather. (iii) The breather excitation has distinct spectroscopic signatures in the power spectra of all spectroscopic observables (such as transition frequencies and transition dipole moments, oscillator strengths), namely, a characteristic component with a period of 34 fs, which is roughly a double period of the C==C stretching motion. We observe that the breather component has the strongest signature in observables related to the intensities (transition dipoles and oscillator strengths), whereas energy variables (such as transition energies) are dominated by vibrational mode frequency components (in particular the C==C stretching mode with a period of 18.3 fs).

The magnitudes of the transient absorption spectra are related to the respective transition dipole moments and therefore contain the breather vibrational component in their power spectra, which should allow identification of breather excitations experimentally (52). Our calculated breather period of 34 fs agrees well with previous estimates for polyacetylene of 38 fs (47, 48) and 30 fs (33) obtained with the SSH models and with a recent experimental value of 44 fs (52). The breather observed in the experiment has a short lifetime of 50 fs. We can understand this dissipation by recalling that the first stage in the breather dynamics is the formation of a self-trapped exciton, which creates phonon excitations. In our simulations, phonon waves reflect from the molecular ends without vibrational energy loss. In real samples the phonons will be scattered on the defects that limit the conjugated segment, and only a fraction of their energy will return back to the breather. Interchain interaction could enhance these dissipative processes even further. Stronger defects will be less transparent for phonons. Therefore we expect that breathers will have longer lifetimes in the defected samples. To confirm this scenario we conducted simulations on an oligomer with two defects H(CH)₁₅  defect  (CH)₃₀  defect  (CH)₁₅H. A weak triple-bond (alkynyl) defect was used in one calculation, whereas a strong sp³ (CH₂ group) defect was used in a second calculation. To simulate dissipation we damped the motion between the chain end and the defect by setting the kinetic energy of nuclei to zero at every time step. We observe that in the case of weak defects the energy dissipates on the time scale of 70 fs (i.e., the breather amplitude drops to 10% of its initial maximum value), whereas in the chain with stronger defects the breather persists up to 200 fs.

In conclusion, strong nonlinear coupling among electronic and vibrational degrees of freedom leads to a complex spatially localized photoexcited vibronic dynamics, which is a distinctive property of conjugated polymers. We expect that this dynamical excitation self-trapping process is typical for conjugated materials featuring strong electron-phonon coupling, including disordered nanomolecules. Development of ultrafast subfemtosecond spectroscopic techniques (52) and efficient quantum-chemical approaches (42, 59, 60) allows monitoring of time-dependent vibronic processes, which opens the way to understand more deeply a number of photochemical, time-dependent spectroscopic, and carrier-transport processes in electronically active soft and biological nanoscale matter.