Density matrix renormalisation group calculations of a suitably parametrised model of long polyenes (polyacetylene oligomers), which incorporates both long range Coulomb interactions and adiabatic lattice relaxation, are presented. The 1 3 B + u and 2 1 A + g states are found to have a 2-soliton and 4-soliton form, respectively, both with large relaxation energies. The 1 1 B − u state forms an exciton-polaron and has a very small relaxation energy. The relaxed energy of the 2 1 A + g state lies below that of the 1 1 B − u state. The soliton/anti-soliton pairs are bound.Electronic interactions in polyenes and polyacetylene (PA) induce strong spin-density-wave correlations in the ground state, resulting in low energy spin-flip (or covalent) triplet ( 3 B + u ) excitations. These combine to form even-parity (dipole-forbidden) singlet ( 1 A + g ) excitations. Optical (dipole-allowed) transitions to the odd-parity singlet state ( 1 B − u ) are essentially ionic in character, resulting in charge transfer from one site to another. In the non-interacting limit the 1 3 B + u and 1 1 B − u states are degenerate, and the 2 1 A + g state always lies higher in energy. However, electron correlations can lead to a reversal of the energetic ordering of the 1 1 B − u and 2 1 A + g states. Electron-electron correlations in π conjugated systems, such as PA, are conveniently modelled by the one-band Pariser-Parr-Pople (P-P-P) model, which includes long range Coulomb interactions.Electron-phonon interactions in the non-interacting limit are described by the SSH model. In the adiabatic limit it predicts a wealth of non-linear excitations, including charged/spinless (S ± ) and neutral/spin 1/2 (S σ ) solitons. It is the inter-play of both electron-electron and electron-phonon interactions in PA which leads to an extremely rich variety of excitations. To describe these excitations we employ the density matrix renormalisation group (DMRG) [1] method to solve the P-P-P-SSH model, and utilise the Hellmann-Feynman (H-F) theorem to calculate the low-lying excited states and the lattice relaxation associated with them.Earlier work on the solitonic structure of the low-lying excitations include, a renormalisation group calculation of the Hubbard-SSH model of up to 16 sites [2]; a meanfield study of the Heisenberg-Peierls model [3]; an exact diagonalisation of a 12 site extended Hubbard-SSH model [4]; and a strong coupling and perturbation calculation of the Hubbard-SSH model [5]. The DMRG method has recently been used by and Yaron et al. [6] and Fano et al. [7] to solve the P-P-P model for linear and cyclic polyenes, respectively. Jeckelmann [8] studied the metal-insulator transition in doped PA by solving the extended Hubbard-SSH with the DMRG method. Likewise, Kuwabara et al.[9] used the DMRG method to study the relative stability of bipolarons using the same model.The P-P-P-SHH Hamiltonian is defined aswhere, t i = t 0 + ∆iis the bond order operator of the ith. bond. We use the Ohno function for the Coulomb interaction: V ij = U/ 1 + βr 2 i...