Methods to compute vibronic spectra of rigid systems are mature enough to become standard, however, when modelling complex systems is convenient to consider one part that represent the “true” system (i.e. chromophore in a protein), even though it may not be a minimum of the reduced system. This poses the problem of computing a vibronic spectrum of a system out-of-equilibrium. One possibility to solve this issue is to remove the degrees of freedom affected by the neglected part. Recently, we worked out a novel protocol in internal curvilinear coordinates to deal with these situations. We have selected polythiophenes (PT) as model because of the complexity of their multiple anharmonic inter-ring torsions. In spite of having been widely studied, the origin of their circular dichroism (ECD) spectrum is still up for debate. It has been shown that oligomer stacked chains forming a supramolecular helix can explain the shape of the ECD spectrum. By the other hand, experimental data from a single PT derivate molecule5 show similar results. In this work, we studied all-anti/all-syn single chain oligothiophenes with different sizes and torsion angles as model using quantum mechanical calculations based on the Density Functional Theory (DFT). Whereas all-anti chains show a strong monosignated ECD in the range of study - which correspond to the first excited state, S1 - all-syn chains possess another bright state slightly higher on energy (S2 in our model) with opposite sign to S1. The coupling of the contributions of the two states give rise to a ECD spectrum in good agreement with experiment.