Theoretical study of the photoisomerization of 1,2-bispyrazinyl-ethylene and the halogen ion salts of 1-Pyrazinyl-2-(4’-methylpyrazinyl)ethylene

Abstract

Abstract Context It has been reported that photoexcitation of azastilbene compounds like E-1,2-bispyrazinyl-ethylene (bpe) can undergo E → Z photoisomerization of its quaternary salts via the excited triplet state. However, experimentally it is possible to get low fluorescence and photoisomerisation quantum yields in a state with higher internal conversion than intersystem crossing. We modelled bpe and its methylated derivative (bpeMe), as well as its quaternary halogen salts (bpeMeX with X = F−, Cl−, Br− and I−) to study levels of fluorescence, phosphorescence and excited state potential energy surfaces (PES). Results support experimental observations of molecules where the anion of a salt is an efficient electron donor, that molecules with weak electron-donating anions like Cl− to give increased fluorescence and photoisomerization, as compared to molecules with stronger electron-donating anions like I−, which are dominated by competing electron transfer. The fluorescence of bpeMeF and bpeMeCl was found to be stronger than bpeMeBr and bpeMeI. A deep well in the triplet excited state of bpeMeI is considered responsible for the decreased photoisomerization, compared to what was experimentally observed for bpeMeCl. Uniquely, the bpeMeI molecule is characterised by near-zero splitting of the s1 and t1 excited states that can enhance charge transfer. The quaternary salt of bpeMe with stronger electron-donating Br− anion was observed to undergo fluorescence and phosphorescence at much lower energy compared to those with weak electron-donating F− and Cl− anions. This research shows how to control the excited state fluorescence, phosphorescence and isomerization of quaternary halogen salts of methyl derivatives of 1,2-bispyrazinyl-ethylene, which aids experimental design where excited state isomerization is considered. Methods Geometry optimization, molecular electrostatic potential (MESP), and time dependent density functional theory (TDDFT) calculations were conducted utilizing Gaussian 16 with the B3LYP functional and the 6–31 + G(d,p) basis set. The minimum energy path (MEP) for the E to Z isomerization of the molecules was established employing the Nudged-Elastic-Band (NEB) method, implemented in Orca 4.2. Precise energies of the E → Z isomerization reaction path were determined employing CASSCF and a more accurate multireference method, NEVPT2.