Semiconducting conjugated polymers (CPs) are pivotal in advancing organic electronics, offering tunable properties for solar cells and field-effect transistors. Here, we carry out first-principle calculations to study individual cis-polyacetylene (cis-PA) oligomers and their ensembles. The ground electronic structures are obtained using density functional theory (DFT), and excited state dynamics are explored by computing nonadiabatic couplings (NACs) between electronic and nuclear degrees of freedom. We compute the nonradiative relaxation of charge carriers and photoluminescence (PL) using the Redfield theory. Our findings show that electrons relax faster than holes. The ensemble of oligomers shows faster relaxation compared to the single oligomer. The calculated PL spectra show features from both interband and intraband transitions. The ensemble shows broader line widths, redshift of transition energies, and lower intensities compared to the single oligomer. This comparative study suggests that the dispersion forces and orbital hybridizations between chains are the leading contributors to the variation in PL. It provides insights into the fundamental behaviors of CPs and the molecular-level understanding for the design of more efficient optoelectronic devices.

In this work, we implemented an approximate algorithm for calculating nonadiabatic coupling matrix elements (NACMEs) of a polyatomic system with ab initio methods and machine learning (ML) models. Utilizing this algorithm, one can calculate NACMEs using only the information of potential energy surfaces (PESs), i.e., energies, and gradients as well as Hessian matrix elements. We used a realistic system, namely CH2NH, to compare NACMEs calculated by this approximate PES-based algorithm and the accurate wavefunction-based algorithm. Our results show that this approximate PES-based algorithm can give very accurate results comparable to the wavefunction-based algorithm except at energetically degenerate points, i.e., conical intersections. We also tested a machine learning (ML)-trained model with this approximate PES-based algorithm, which also supplied similarly accurate NACMEs but more efficiently. The advantage of this PES-based algorithm is its significant potential to combine with electronic structure methods that do not implement wavefunction-based algorithms, low-scaling energy-based fragment methods, etc., and in particular efficient ML models, to compute NACMEs. The present work could encourage further research on nonadiabatic processes of large systems simulated by ab initio nonadiabatic dynamics simulation methods in which NACMEs are always required.

Nonadiabatic effects are ubiquitous in physics, chemistry, and biology. They are strongly amplified by conical intersections (CIs), which are degeneracies between electronic states of triatomic or larger molecules. A few years ago it was revealed that CIs in molecular systems can be formed by laser light, even in diatomics. Because of the prevailing strong nonadiabatic couplings, the existence of such laser-induced conical intersections (LICIs) may considerably change the dynamical behavior of molecular systems. By analyzing the photodissociation process of the D2+ molecule carefully, we found a robust effect in the angular distribution of the photofragments that serves as a direct signature of the LICI, providing undoubted evidence of its existence.