Trapping of excitons at chemical defects in polyethylene

Mount for spectroscopic analysis of samples under sustained tensile stress

Spectroscopic methods offer valuable insights into the molecular and structural changes induced by stress, but existing techniques are often unable to perform real-time measurements during deformation. A novel solid open mount design is presented that enables spectroscopic investigations of materials under sustained tensile stress while maintaining crucial alignment of the optical system. The mount design allows for sample movement in response to applied strain while maintaining the position of the sample plane, ensuring consistent and reliable spectroscopic measurements. The effectiveness of the mount design is demonstrated with vibrational sum-frequency generation measurements of an elastomer, cured hydroxyl-terminated polybutadiene, and a plastic, high-density polyethylene, taken before, during, and after tensile deformation. The application of this mount to other spectroscopic techniques is discussed. The ability to collect spectroscopic data during a stress event would provide valuable insights into the behavior of stressed materials.

Dynamics of Phonon-Assisted Holes Trapping and Transport over Chemical Defects in Polyethylene

Charge trapping and transport over chemical defects in polyethylene have significant impacts on its electrical and dielectric properties. However, the dynamics of this phenomenon and its underlying mechanisms remain unclear. To understand this fundamental aspect, we conducted a time-domain ab initio nonadiabatic molecular dynamics study of phonon-assisted holes dynamics in polyethylene over C═O and C-OH defect states. Our results suggest that the hole transfer and energy fluctuations substantially depend on temperature and local morphology. When the temperature decreases from 300 to 100 K, the hole transfer efficiency and the energy fluctuations are severely suppressed due to the weakened interactions between holes and phonons. Furthermore, amorphous polyethylene exhibits a severe suppression of the hole transfer process compared to crystalline polyethylene. An explanation for the influence of morphology on the hole transfer process can be found in the differences in the hole-phonon coupling and the electronic coupling between two chemical defect states in crystalline and amorphous polyethylene. Advancing the fundamental understanding of the dynamics of hole transfer over chemical effects in polymers is a key to improving their insulating properties for the next-generation high-voltage cables.

The ΔSCF method for non-adiabatic dynamics of systems in the liquid phase

Computational studies of ultrafast photoinduced processes give valuable insights into the photochemical mechanisms of a broad range of compounds. In order to accurately reproduce, interpret, and predict experimental results, which are typically obtained in a condensed phase, it is indispensable to include the condensed phase environment in the computational model. However, most studies are still performed in vacuum due to the high computational cost of state-of-the-art non-adiabatic molecular dynamics (NAMD) simulations. The quantum mechanical/molecular mechanical (QM/MM) solvation method has been a popular model to perform photodynamics in the liquid phase. Nevertheless, the currently used QM/MM embedding techniques cannot sufficiently capture all solute-solvent interactions. In this Perspective, we will discuss the efficient ΔSCF electronic structure method and its applications with respect to the NAMD of solvated compounds, with a particular focus on explicit quantum mechanical solvation. As more research is required for this method to reach its full potential, some challenges and possible directions for future research are presented as well.

Origin of the Failure of Density Functional Theories in Predicting Inverted Singlet-Triplet Gaps

Recent experimental and theoretical studies have shown several new organic molecules that violate Hund's rule and have the first singlet excited state lower in energy than the first triplet excited state. While many correlated single reference wave function methods have successfully predicted excited-state energetics of these low-lying states, conventional linear-response time-dependent density functional theory (TDDFT) fails to predict the correct excited-state energy ordering. In this article, we have explored the performance of combined DFT and wave function methods like doubles-corrected TDDFT and multiconfiguration pair-density functional theory for the calculation of inverted singlet-triplet gaps. We have also tested the performance of the excited-state DFT (eDFT) method for this problem. Our results have shown that it is possible to obtain inverted singlet-triplet gaps both by using doubles-corrected TDDFT with a proper choice of double-hybrid functionals or by using eDFT.

Using C-DFT to develop an e-ReaxFF force field for acetophenone radical anion

Increased electricity usage over the past several decades has accelerated the need for efficient high-voltage power transmission with reliable insulating materials. Cross-linked polyethylene (XLPE) prepared via dicumyl peroxide (DCP) cross-linking has emerged as the insulator of choice for modern power cables. Although DCP cross-linking generates the desired XLPE product in high yield, other by-products are also produced. One such by-product, acetophenone, is particularly intriguing due to its aromaticity and positive electron affinity. In this work, constrained density functional theory (C-DFT) was utilized to develop an e-ReaxFF force field suitable for describing the acetophenone radical anion. Initial parameters were taken from the 2021 Akbarian e-ReaxFF force field, which was developed to describe XLPE chemistry. Then, C-DFT geometry optimizations were performed wherein an excess electron was constrained to each atom of acetophenone. The resulting C-DFT energy values for the various electronic positions were added to the e-ReaxFF training set. Next, an analogous set of structures was energy-minimized using e-ReaxFF, and equilibrium mixture compositions for the two methods were compared at multiple temperatures. Iterative fitting against C-DFT energy data was performed until satisfactory agreement was achieved. To test force field performance, molecular dynamics simulations were performed in e-ReaxFF and the resulting electronic distributions were qualitatively compared to unconstrained-DFT spin density data. By expanding our e-ReaxFF force field for XLPE, namely, adding the capability to describe acetophenone and its interactions with an excess electron, we move one step closer to a comprehensive molecular understanding of XLPE chemistry in a high-voltage power cable.

Atomic Scale Mechanisms Controlling the Oxidation of Polyethylene: A First Principles Study

Understanding the degradation mechanisms of aliphatic polymers by thermal oxidation and radio-oxidation is very important in order to assess their lifetime in a variety of industrial applications. We focus here on polyethylene as a prototypical aliphatic polymer. Kinetic models describing the time evolution of the concentration of chain defects and radicals species in the material identify a relevant step in the formation and subsequent decomposition of transient hydroperoxides species, finally leading to carbonyl defects, in particular ketones. In this paper, we first summarize the most relevant mechanistic paths proposed in the literature for hydroperoxide formation and decomposition and, second, revisit them using first principles calculations based on Density Functional Theory (DFT). Our results partially confirm commonly accepted reaction energies, but also propose alternative, more favourable, reaction paths. We highlight the influence of the environment-crystalline or not-on the outcome of some of the studied chemical reactions. A remarkable result of our calculations is that hydroxyl radicals play an important role in the decomposition of hydroperoxides. Based on our findings, it should be possible to improve the set of equations and parameters used in current kinetic simulations of polyethylene radio-oxidation.

Approximating Quasiparticle and Excitation Energies from Ground State Generalized Kohn-Sham Calculations

Quasiparticle energies and fundamental band gaps in particular are critical properties of molecules and materials. It was rigorously established that the generalized Kohn-Sham HOMO and LUMO orbital energies are the chemical potentials of electron removal and addition and thus good approximations to band edges and fundamental gaps from a density functional approximation (DFA) with minimal delocalization error. For other quasiparticle energies, their connection to the generalized Kohn-Sham orbital energies has not been established but remains highly interesting. We provide the comparison of experimental quasiparticle energies for many finite systems with calculations from the GW Green function and localized orbitals scaling correction (LOSC), a recently developed correction to semilocal DFAs, which has minimal delocalization error. Extensive results with over 40 systems clearly show that LOSC orbital energies achieve slightly better accuracy than the GW calculations with little dependence on the semilocal DFA, supporting the use of LOSC DFA orbital energies to predict quasiparticle energies. This also leads to the calculations of excitation energies of the N-electron systems from the ground state DFA calculations of the ( N - 1)-electron systems. Results show good performance with accuracy similar to TDDFT and the delta SCF approach for valence excitations with commonly used DFAs with or without LOSC. For Rydberg states, good accuracy was obtained only with the use of LOSC DFA. This work highlights the pathway to quasiparticle and excitation energies from ground density functional calculations.

Surface Spectroscopic Signatures of Mechanical Deformation in High-Density Polyethylene (HDPE)

High-density polyethylene (HDPE) has been extensively studied, both as a model for semi-crystalline polymers and because of its own industrial utility. During cold drawing, crystalline regions of HDPE are known to break up and align with the direction of tensile load. Structural changes due to deformation should also manifest at the surface of the polymer, but until now, a detailed molecular understanding of how the surface responds to mechanical deformation has been lacking. This work establishes a precedent for using vibrational sum-frequency generation (VSFG) spectroscopy to investigate changes in the molecular-level structure of the surface of HDPE after cold drawing. X-ray diffraction (XRD) was used to confirm that the observed surface behavior corresponds to the expected bulk response. Before tensile loading, the VSFG spectra indicate that there is significant variability in the surface structure and tilt of the methylene groups away from the surface normal. After deformation, the VSFG spectroscopic signatures are notably different. These changes suggest that hydrocarbon chains at the surface of visibly necked HDPE are aligned with the direction of loading, while the associated methylene groups are oriented with the local C_2v symmetry axis roughly parallel to the surface normal. Small amounts of unaltered material are also found at the surface of necked HDPE, with the relative amount of unaltered material decreasing as the amount of deformation increases. Aspects of the nonresonant SFG response in the transition zone between necked and undeformed polymer provide additional insight into the deformation process and may provide the first indication of mechanical deformation. Nonlinear surface spectroscopy can thus be used as a noninvasive and nondestructive tool to probe the stress history of a HPDE sample in situations where X-ray techniques are not available or not applicable. Vibrational sum-frequency generation thus has great potential as a platform for material state awareness (MSA) and should be considered as part of a broader suite of tools for such applications.

Optical Properties of Saturated and Unsaturated Carbonyl Defects in Polyethylene

Polyethylene (PE), one of the simplest and most used aliphatic polymers, is generally provided with a number of additives, in particular antioxidants, because of its tendency to get oxidized. Carbonyl defects, a product of the oxidation of PE, are occurring in various forms, in particular saturated ones, known as ketones, where a C═O double bond substitutes a CH₂ group, and various unsaturated ones, i.e., with further missing hydrogens. Many experimental investigations of the optical properties in the visible/UV range mainly attribute the photoluminescence of PE to one specific kind of unsaturated carbonyls, following analogies to the emission spectra of similar small molecules. However, the reason why saturated carbonyls should not be optically detected is not clear. We investigated the optical properties of PE with and without carbonyl defects using perturbative GW and the Bethe-Salpeter equation in order to take into account excitonic effects. We discuss the calculated excitonic states in comparison with experimental absorption/emission energies and the stability of both saturated and unsaturated carbonyl defects. We conclude that the unsaturated defects are indeed the best candidate for the luminescence of oxidized PE, and the reason is mainly due to oscillator strengths.

Calculating excited state properties using Kohn-Sham density functional theory

The accuracy of excited states calculated with Kohn-Sham density functional theory using the maximum overlap method has been assessed for the calculation of adiabatic excitation energies, excited state structures, and excited state harmonic and anharmonic vibrational frequencies for open-shell singlet excited states. The computed Kohn-Sham adiabatic excitation energies are improved significantly by post self-consistent field spin-purification, but remain too low compared with experiment with a larger error than time-dependent density functional theory. Excited state structures and vibrational frequencies are also improved by spin-purification. The structures show a comparable accuracy to time-dependent density functional theory, while the harmonic vibrational frequencies are found to be more accurate for the majority of vibrational modes. The computed harmonic vibrational frequencies are also further improved by perturbative anharmonic corrections, suggesting a good description of the potential energy surface. Overall, excited state Kohn-Sham density functional theory is shown to provide an efficient method for the calculation of excited state structures and vibrational frequencies in open-shell singlet systems and provides a promising technique that can be applied to study large systems.