1. Introduction
In recent years, the study of conjugated polymers (CPs) has attracted significant interest, particularly because of their photoluminescent characteristics. The unique arrangement, featuring a sequence of alternating single and double bonds, is a critical factor contributing to the distinct properties of CPs [
1,
2]. These conjugations enable the electrons to delocalize through the polymeric chain. As a result of their properties, CPs have been used in many important applications such as organic photovoltaic [
3], organic thermo-electrics [
4], supercapacitors [
5], electrochromic [
6], biosensing [
7], photocatalysis [
8], optoelectronic device technology [
9,
10], such as light emitting diodes (LEDs) [
2,
11], transistors (thin-film transistors, TFTs) [
11], and solar cells [
9,
11]. Currently, these compounds are among the most promising materials, due to their distinctive characteristics, such as processability, affordability, and thermal stability in film formation [
12,
13]. A considerable amount of experimental and theoretical work has been dedicated to elucidating the electronic and optical characteristics of CPs [
14].
CPs often have structures that are amorphous or semi-crystalline [
15]. They exhibit short-range order without long-range regularity, which complicates their solubility and results in polydispersity. This makes it challenging for computational models to accurately represent these polymers with a variety of sizes and shapes. Consequently, the properties of shorter oligomeric chains have been studied, from which one can infer the characteristics of much longer polymeric chains. However, accurate electronic properties of polymers cannot simply be extrapolated from oligomers, as revealed by Hutchison et al. [
16]. Studies have shown that, through meticulous chemical design involving the strategic placement of fluorine atoms within the polymer’s backbone, aqueous CP nanoparticles exhibit elevated photoluminescence (PL) quantum yields, particularly in the far-red and near-infrared (NIR) regions [
17]. The functionalization of CPs with flexible side chains, including alkyl groups, leads to the enhancement of solubility and device performance. These side chains do not directly participate in charge transport, but they can influence the overall device performance by affecting the molecular packing and physical robustness of the films [
18].
Among the diverse CPs, polyacetylene (PA) is distinguished by its prototypical role and noteworthy photoelectric qualities [
19]. Monosubstituted PA exhibits a regular helical structure that can be fine-tuned and polymerized under non-stringent conditions. These attributes allow for precise control over their luminescence properties, making them suitable for various PL applications [
20,
21]. The singular dimensionality of CPs facilitates soliton formation, which impacts the charge separation and transport within these materials [
22,
23]. Photogenerated excitations in trans-PA tend to relax into soliton states, significantly reducing PL in the visible range while primarily confining it to the infrared range. This characteristic, stemming from the unique topological and electronic structures of trans-PA, inherently limits its effectiveness in applications requiring visible-range PL. In contrast, cis-PA does not exhibit this soliton-induced reduction in PL, making it suitable for broader applications in optoelectronic devices, where efficient and strong light emission is essential [
22,
24].
The application of density functional theory (DFT) is essential for studying the electronic structures of semiconductors, providing a comprehensive understanding of their ground states and mechanisms involved in charge transfer at material interfaces [
25]. The computational methodologies divide the description of charge transfer into two distinct segments: the initial structural modeling and electronic state analysis, and the subsequent categorization and dynamic assessment of excited states. The interplay between dissipative quantum dynamics and molecular trajectories, informed by DFT and refined through approaches such as surface hop** [
26,
27,
28,
29,
30,
31] and Redfield [
32,
33,
34], remains a rich area of research to unravel the complexities of charge-transfer excitations and their temporal evolution.
In previous work [
35], we focused on the ground state properties of ensembles of cis-PA oligomers. We have also studied the nonradiative relaxation of charge carriers and PL of individual cis-PA oligomers [
36]. In this study, we examine the excited state dynamics of ensembles of oligomers, comparing the results with those of individual oligomers. Specifically, the ground electronic structures are obtained from DFT calculations. We then compute nonadiabatic couplings (NACs) between electronic and nuclear degrees of freedom along adiabatic molecular dynamic (MD) trajectories. The excited state dynamics is computed using the reduced density matrix (RDM) [
37,
38,
39] formalism within the Redfield theory. This method provides a deeper understanding of the complex behavior of charge carriers, enabling us to map out the possible relaxation processes for both individual and ensembles of oligomers.
3. Results
Figure 2a describes the energy distributions of the frontier orbitals for the single oligomer. We utilize dual notations for these orbitals:
for the unoccupied orbitals and
for occupied ones. The energy of the lowest unoccupied molecular orbital, i.e., LU (
), is at
and that of the highest occupied molecular orbital, i.e., HO
, is at
. When examining higher
values, there is an increase in the energy levels of the unoccupied orbitals. This notation is according to the quantum confinement paradigm. A way to count orbitals of electrons upward from the conduction band minimum (CBM) and holes downward from the valence band maximum (VBM) is an efficient way to establish an analogy with momentum dispersion in direct gap infinite polymers;
;
. Note the negative size of dispersion for holes. For finite size oligomers, both
follow the trend of a particle in a box, reflected in the special pattern of subgaps. These background concepts provide a foundation for a more organized analysis of electronic structures.
Our findings from the single oligomer model indicate an energy gap
of 1.2 eV between the HO and LU orbitals. Furthermore,
Figure 2b illustrates that the gap between HO and LU is more pronounced than the subgaps noted between HO and HO
1, as well as LU and LU+1 orbitals for the single oligomer. The symbols A, B, C, D, and M are used for the single oligomer, and A′ and M′ are used for the ensemble of oligomers. A−D represent transitions from HO to LU (A), HO−1 to LU (B), HO to LU+1 (C), and HO
1 to LU+1 (D). A′ shows the transition from HO to LU (
). A−D and A′ correspond to interband transitions, whereas M and M′ correspond to intraband transitions that may arise in PL.
Although the ensemble of cis-PA oligomers follows the pattern of a particle in a box, there is a slight difference from the single oligomer. There are interactions for the ensemble of oligomers. Due to inter-oligomer interactions, the degenerate states of individual oligomers are lifted in the ensemble. The spatial proximity of the oligomers induces the interaction between their individual eigenstates. Generally, an interaction between any two quantum levels results in an increase in their energy offset. If the offset is zero, indicating degeneracy, the interaction lifts this degeneracy, making the offset nonzero. Compared to a single oligomer, an ensemble of oligomers has a greater chance of hybridization between the electronic states.
In our earlier work [
36], we did show that TDDFT calculations with a hybrid functional provide gap values comparable to experimental results. By comparing DFT and TDDFT data, one may extract a systematic correction. For a single oligomer,
and
. The difference is
. Thus, each transition energy in our calculations may be corrected as follows:
An alternative way is to scale the energies by a factor of
, where
. In this approach, all transition energies in our calculations may be corrected as follows:
.
Figure 2c illustrates the spatial iso-contours of the frontier orbitals for a single oligomer of cis-PA. Each orbital demonstrates two properties: first, a repeating pattern in the vicinity of each carbon atom, and second, a slowly changing envelope function that modulates the density near each atom as a function of position along the oligomer. Notably, the LU orbital demonstrates the maximum amount of the charge density peaks within its central domain, whereas the edges display reduced densities. In analogy to the particle in a box model, the maximum charge density is located centrally in the box, with diminished densities at the edges. For
, the charge density’s peak shifts from its initial position, indicative of the eigenstate’s index increment. This results in the orbital splitting into two segments, with the highest charge density concentrated in the central regions and negligible density at the edges, as depicted for the ensemble of cis-PA oligomers in
Figure S1. The ensemble of oligomers exhibits a charge distribution pattern akin to the particle in a box model.
The calculated absorption spectra for the single and ensemble of oligomers are reported in
Figure 3. The peaks are labeled in the order of ascending transition energy. The initial peak A′ and A are observed at 0.8 eV for the ensemble and at 1.2 eV for the individual oligomer, respectively, indicative of transitions from HO to LU. The other peaks in the spectra are representative of various electronic transitions occurring among different orbitals. Peaks for an ensemble of oligomers are redshifted, broader, and of lower intensity compared to peaks for a single oligomer. The broadening is due to inter-oligomer interactions.
Table 1 provides representative transitions based on peak labels from
Figure 3, which are used as initial conditions for excited state dynamics.
In the next step, we explore nonradiative transitions between electronic states. The RDM approach and Redfield formalism provide components of the Redfield tensor after processing the autocorrelation function of NACs. The Redfield tensors provide the electronic level-to-level transition rates.
Figure 4a–d demonstrate the autocorrelation function of NACs and the elements of the Redfield tensor, respectively, for the single and ensemble of oligomers. Using the autocorrelation function, we can explain the average nonadiabatic interaction in the electron–lattice interaction. This autocorrelation function is the correlation function, which provides information about the intensity of dissipative electronic transitions for the given indices
with several amplitudes. Interestingly, for any tested combination of indices, the autocorrelation function decays abruptly within 5 fs for the single oligomer and within 2 fs for the ensemble. The ensemble decays at a faster rate than the single oligomer, since inter-chain interactions affect electronic transitions.
Dissipative excited state dynamics are provided in
Figure 5. There are two distinct sectors, each marked by a different color. The yellow sections signify an electron charge surplus
, while the blue sections denote a shortage in hole charge
, compared to the neutral charge distribution
.
Figure 5a visualizes the nonradiative relaxation of the photoexcited electrons and holes for a single oligomer, which can take ~15 ps and ~150 ps, respectively. The energy dissipation dynamics are tracked through the expected values of energy for both electrons and holes. Between roughly ~10 ps and ~15 ps, one observes the electron relaxation from LU
1 to LU. Following a period of around 150 ps, there is an observed energy shift for the holes. This increase marks the migration of the hole population from HO
1 to HO.
Figure 5b delves into the progression of the non-equilibrium charge density distribution along the z-axis. The moments when electrons and holes achieve maximum charge density transfer are marked by vertical lines for electrons
and for holes
, respectively. The change in spatial distribution of charge can be interpreted in terms of the pattern of a particle in a box, observed for progressions of occupied and unoccupied orbitals. The central and edge regions of the box lack yellow lines in the range of 0~15 ps due to the specific pattern of the HO
1 and LU
1 orbitals as depicted in
Figure 2. After 15 ps, the yellow lines have already condensed in the center, and the blue areas are still available. This corresponds to LU with maximal density in the central area. We observe that at a time between 15 ps and 150 ps, the electrons and holes stay in different spatial areas of the single oligomer. At ~150 ps, the blue lines also condense in the center of the oligomer as a signature of HO getting populated.
Figure 5d,e illustrates the dynamic results for the ensemble of oligomers. The initial photoexcitation condition is HO
6
LU+4. At the initial time, the localization of holes in oligomers 3 and 4 is evident, while electrons are primarily localized in oligomers 1 and 6. Then, one sees the hot charge carriers undertaking stepwise changes until the electrons and holes reach frontier orbitals, and in the process dissipating excess energy transforms to heat. Along the excited state trajectories, the ensemble of oligomers experiences the subsequent population of several intermediate electronic states. At 1~5 ps, electrons and holes stay in different areas of the space. Specifically, electrons stay at oligomers 1 and 6, whereas holes stay at oligomers 2, 3, and 4. Finally, after ~5 ps, electrons are at oligomers 5 and 6, whereas holes are at oligomers 1, 2, and 3.
Figure S1 depicts HO
1, HO, LU, and LU
1 orbitals, which clearly shows the localization pattern.
Table 2 and
Table 3 summarize the charge carrier relaxation rates for the single oligomer and ensemble of oligomers under different initial photoexcitation conditions. It is found that electrons relax faster than holes in both single and ensemble of oligomers. This is due to smaller subgaps and more pronounced phonon interactions for electrons than holes. Additionally, the ensemble of oligomers shows faster relaxation compared to the single oligomer. This originates from the higher density of states and a larger number of relaxation pathways in the ensemble compared to the single oligomer. In
Figure 6, we explore the dependence of relaxations rates of charge carriers on dissipated energies. For both models, one finds the decrease in relaxation rates for electrons and holes with increasing dissipated energies in accordance with the so-called energy gap law [
57].
Figure 7 depicts the dynamics of radiative relaxation for the single oligomer and ensemble of oligomers. The initial conditions are the same as in
Figure 5. Panels (a, d) show the dynamics of excitation energy dissipation. In these spectra, one finds features A–C for the single oligomer and features A′–C′ for the ensemble of oligomers. These features are due to the cascading down of excitation energy from the initial excitation to the lowest excitation. Panels (b, e) show the time-resolved emission spectra. For the ensemble, one finds several emission features with vanishing intensities. Peaks A′ and C′ are more pronounced than peak B′. When the population cascades, it may visit the dark states. The main features are peaks M and M′, which correspond to intraband emission features below the band gap. These intraband transitions appear in the emission spectra but not in the absorption spectra. Panels (c, f) show the time-integrated emission spectra. Peaks in these spectra are consistent with peaks from the time-resolved emission spectra. However, some emission characteristics such as B′, which are quite faint and almost imperceptible in time-resolved emission spectra, become apparent in the time-integrated emission spectra. The onset of these spectral features occurs later in the dynamic process. This prolonged activity contributes significantly to the time-integrated emission spectra.
The emission spectra of both models demonstrate both interband and intraband emission. Peaks C and C′ corresponding to the initial excitations will vanish for lower energy excitations. At sufficiently high energy excitations, these peaks will be visible. The current initial conditions are intentionally chosen to simplify analysis. The ensemble of oligomers shows three effects: broader line widths, redshift of transition energies, and lower intensities compared to the single oligomer.
Figure 8 illustrates the dynamics and spatial distribution of electrons and holes for the single oligomer and ensemble of oligomers. It also showcases the varying nature of molecular orbitals, which can be observed as either localized or delocalized. These patterns are directly associated with the spatial distribution of electrons and holes of each system. The initial conditions are the same as in
Figure 5. For the single oligomer, the delocalization patterns for both electrons and holes are visible at the initial time. After relaxation, electrons are more localized and holes are more delocalized. Comparing electrons with holes throughout the relaxation pathway, electrons remain more localized than holes. For the ensemble, electrons are delocalized over four oligomers, and holes are localized over two oligomers at the initial time. After relaxation, both electrons and holes are more localized. Comparing electrons with holes throughout the relaxation pathway, electrons are generally less localized than holes, except in the intermediate states, where both electrons and holes show greater delocalized patterns.
4. Discussion
A systematic study on the dependence of chain length is known in the literature and is ongoing as a continuation of this work. An increase in chain length is expected to provide a mild redshift of transition energies. Our exploration into models with shorter-length oligomers suggests a decrease in computed conductivity. This decrease is attributed to quantum confinement effects in shorter oligomers, which increase band gaps and subgaps.
In the context of CPs, oligomers demonstrate the formation of bound excitons of finite sizes, primarily due to electron–hole interactions. Brey and Burghardt [
58] elucidate a coherent transient localization mechanism for exciton transport in regioregular poly(3-hexylthiophene), highlighting the interplay of exciton polarons and interchain vibrational modes. In shorter oligomers, the confinement is more pronounced, affecting the electronic properties significantly. In contrast, longer oligomers exhibit transition energies that stabilize and are defined predominantly by the size of the exciton and polaron rather than by further increases in length. This behavior aligns with observations in the literature, suggesting a general trend for exciton behavior in CPs. For instance, Kraner et al. [
59] find that exciton sizes increase uniformly with molecular π-system sizes ranging from 10 to 40 Å. Beyond this range, no significant increase in exciton size is observed. These findings guide our decision to utilize shorter oligomers.
We focus on the influences and consequences of inter-chain interaction on nonradiative relaxation and spectral line widths. The fundamental photophysical properties such as charge carrier dynamics and PL behavior are inherent to cis-PA. These properties demonstrate consistency across the different molecular weights we have studied, suggesting that the trends we observe could extend to cis-PA products with broader molecular weights. As soon as the chain length exceeds the size of the polaron, its influence on such observables becomes negligible. There is ongoing research to confirm this hypothesis. Adding additional defects, such as kinks, non-planar structures, and twists, will provide additional effects. This approach aligns with the comprehensive design strategies outlined in the literature, providing a foundation for our continued research into develo** more efficient and application-specific CPs. In the current study, dispersion forces and orbital hybridizations between chains are the major effect.
The incorporation of methyl groups as terminal substituents in our oligomer configurations offers significant benefits. Methyl groups enhance molecular stability by providing steric protection, which reduces the susceptibility of the polymers to oxidative degradation. This protection is crucial for maintaining the structural integrity of the oligomers during experimental procedures and analyses. Additionally, these methyl groups significantly influence the electronic structures of the oligomers. They adjust the electron density along the polymer backbone, thereby affecting the band gap and electronic transitions crucial for the PL properties under study. The strategic inclusion of methyl groups modifies the energy states of HO and LU, enhancing the optical properties of the oligomers. Therefore, the choice to use methyl groups instead of hydrogen atoms is not merely a structural decision but a strategic one aimed at optimizing the optoelectronic properties of the oligomers.
The current investigation is designed to deepen the theoretical understanding and enhance the accuracy of computational models within a defined theoretical framework. While experimental validation, partially provided in ref. [
60], remains critical for the broader application of these computational study, the present study is limited to comparative analyses between different computational models. The computation-to-experiment comparison, while highly valuable, falls outside the scope of the current study. We expect that adding the scissor operator or scaling correction would lead to accurate transition energies and bring the spectra in computational results closer to experimental values.