Figure shows the linear relationship between the calculated (∆𝑮 ‡ ) and σ para for 1-3 and 7-13 (top), and (1-3)-OMe and (7-13)-OMe (bottom). We used a broad range of EWG-and EDG-substituents (σ para ranges from -0.66 -0.81). The R 2 values are 0.89 and 0.78 for the top and bottom plots, respectively. The linear correlation suggests there is a moderately strong relationship between the ability of a substituent to delocalize the π NN -electrons and the ∆𝑮 ‡ and corresponding t 1/2 . The slopes in Figures and We present our computed results for all 26 hemi-azothiophenes in the scatterplot in Figure . We compared our computed λ max and t 1/2 to those of hemi-azothiophene reported by Wegner in 2020 (λ m ax of 365 nm and t 1/2 of 17.7 hours). The scatterplot is organized into four quadrants; The red quadrant (bottom left) contains photoswitches where both λ max and t 1/2 are lower relative to the past system, while the yellow quadrants (top left and bottom right) have a red-shifted λ max and shorter t 1/2 (topleft) or a longer t 1/2 and a lower λ max (bottom right) compared to the past system. The green quadrant (top right) contains eight photoswitches from our set with red-shifted λ max and longer t 1/2 than hemi-azothiophene. plotted against the E-isomer λ max (4a). We divided the scatterplot into four quadrants. The eight azothiophenes in the green quadrant are shown in Figure . The λ max values were computed using ωb97-XD/aug-cc-pVTZ and IEFPCM water , and the t 1/2 were calculated using PBE0-D3BJ/6-31+G(d,p) and IEFPCM water .