We simulate the interchain polaron recombination process in conjugated polymer systems using a nonadiabatic molecular dynamics method, which allows for the coupled evolution of the nuclear degrees of freedom and multiconfigurational electronic wavefunctions. Within the method, the appropriate spin symmetry of the electronic wavefunction is taken into account, thus allowing us to distinguish between singlet and triplet excited states. It is found that the incident polarons can form an exciton, form a bound interchain polaron pair, or pass each other, depending on the interchain interaction strength and the strength of an external electric field. Most importantly, we found that the formation of singlet excitons is considerably easier than triplet excitons. This shows that in real organic light emitting devices, the electroluminescence quantum efficiency can exceed the statistical limitation value of 25%, in agreement with experiments.