First principles DFT study of polyethylene insulation containing chemical impurities - implementing counterpoise correction

Abstract: Density functional theory (DFT) calculations of polyethylene (PE) HVDC cable insulation have been performed for systems containing four different chemical impurities: acetophenone, cumene, $\alpha$-methyl styrene and $\alpha$-cumyl alcohol. Systems were generated by molecular dynamics (MD) equilibration at four different temperatures relevant for cable insulation applications: 277 K, 293 K, 343 K and 363 K. With the goal of gaining better measure of variations in hole and electron traps energies, four initial configurations were also stochastically generated at each temperature, which yielded four different final configurations after equilibration. The counterpoise correction scheme was implemented for DFT calculations, by distributing ghost atoms thought any empty pockets of space in between the PE chains. The PBE functional was selected for DFT simulations. The resulting band gaps were in agreement with those of earlier GGA-based studies, and thus lower by 3 eV than empirical band gaps. For all impurities, the first HOMO state and the first two LUMO states were generally located on the impurity molecule, forming one hole trap and two electron traps, but certain configurations generated increased electron trap numbers, or eliminated hole traps. No dependence could be derived between temperature and trap depth for either electron or hole traps. Mean electron trap energies were largely in agreement with results from earlier studies, they were deepest for acetophenone, and they varied by as much as 0.6 eV between different configurations. Hole traps are universally shallow and vary by up to 0.7 eV between configurations, and are similar in depth for all impurities. Results suggest that electron trap depths correlate with the presence of molecular features such as oxygen atoms and conjugated double bonds. The dependence of trap depth on the spatial configuration of the impurity molecule suggests that results could be improved by more precise quantum mechanical treatment of the dynamics of the impurity.