Synthesis and Spectroscopic Characterization of Selected Phenothiazines and Phenazines Rationalized Based on DFT Calculation

Two unique structures were isolated from the phosphorylation reaction of 10H-phenothiazine. The 5,5-dimethyl-2-(10H-phenothiazin-10-yl)-1,3,2-dioxaphosphinane 2-oxide (2a) illustrates the product of N-phosphorylation of phenothiazine. Moreover, a potential product of 2a instability, a thiophosphoric acid 2b, was successfully isolated and structurally characterized. Molecule 2a, similarly to sulfoxide derivative 3, possesses interesting phosphorescence properties due to the presence of d-pπ bonds. The X-ray, NMR, and DFT computational studies indicate that compound 2a exhibits an anomeric effect. Additionally, the syntheses of selected symmetrical and unsymmetrical pyridine-embedded phenazines were elaborated. To compare the influence of phosphorus and sulfur atoms on the structural characteristics of 10H-phenothiazine derivatives, the high-quality crystals of (4a,12a-dihydro-12H-benzo[5,6][1,4]thiazino[2,3-b]quinoxalin-12-yl)(phenyl)methanone (1) and selected phenazines 5,12-diisopropyl-3,10-dimethyldipyrido[3,2-a:3′,2′-h]phenazine (5) and 5-isopropyl-N,N,3-trimethylpyrido[3,2-a]phenazin-10-amine (6a) were obtained. The structures of molecules 1, 2a, 2-mercapto-5,5-dimethyl-1,3,2-dioxaphosphinane 2-oxide (2b), 3,7-dinitro-10H-phenothiazine 5-oxide (3), 5 and 6a were determined by single-crystal X-ray diffraction measurements.


Introduction
The phenothiazines and the phenazines are exceptional and versatile classes of electronrich nitrogen heterocycles with broad pharmaceutical profiles and rich materials sciences applications [1][2][3]. A seemingly small difference between these types of molecules, i.e., the presence of a six-membered ring with nitrogen and sulfur atoms for phenothiazines and two nitrogen atoms for phenazines, manifests itself profoundly through their physicochemical and structural properties. The unique characteristics of both systems are present in (4a,12a-dihydro-12H-benzo [5,6] [1,4]thiazino [2,3-b]quinoxalin-12-yl)(phenyl)methanone (1) and its analogs.
The phenazine class of N-heterocycles, numbering over 6000 compounds, is an important structural motif in functional materials, pharmaceuticals, agricultural, and natural products. Therefore, developing new efficient methods for preparing molecules with this privileged moiety is an essential endeavor in synthetic chemistry [4]. Their broad range of biological activity is thought to be due to their ability to undergo redox cycling in The phenothiazine molecule possesses a unique non-planar butterfly structure. Tian and co-workers drew attention to the derivatives having an S=O (or SO2) and P=O group, which contains the d-pπ bond that could hinder a possible intense intermolecular π-π stacking, which is beneficial for inhibiting excited triplet-triplet annihilation (TTA) [27]. Tian et al. presented some examples of molecules with a S=O bond (in the SO2 group), which shows a sharp increase in TTA by 19 times compared to that without this group.
Compound 1 crystallizes in the monoclinic system with the P21/n space group ( Figure  1) and demonstrates layered packing similar to the phenothiazine model system described by Tian and co-workers [27]. Intermolecular π-π stacking distances vary in the range of 3.670-3.736 Å (Figure 1c). Additionally, a rigid 2D network structure of 1 is stabilized via intermolecular C-H···N interactions with distances of 2.690 Å. The bending angles in the phenothiazine core of molecule 1 were defined as the ∠N12−S11−C10 and ∠N5−N6−C7, as shown in Figure 1a,b, and were found to be 130.88(5)° and 131.92(8)°. The distance between the centroids of the phenyl and azine rings is equal to 3.712(1)Å. Compound 1 crystallizes in the monoclinic system with the P2 1 /n space group ( Figure 1) and demonstrates layered packing similar to the phenothiazine model system described by Tian and co-workers [27]. Intermolecular π-π stacking distances vary in the range of 3.670-3.736 Å (Figure 1c). Additionally, a rigid 2D network structure of 1 is stabilized via intermolecular C-H· · · N interactions with distances of 2.690 Å. The bending angles in the phenothiazine core of molecule 1 were defined as the ∠N12−S11−C10 and ∠N5−N6−C7, as shown in Figure 1a,b, and were found to be 130.88 (5) • and 131.92 (8) • . The distance between the centroids of the phenyl and azine rings is equal to 3.712(1)Å. The phenothiazine molecule possesses a unique non-planar butterfly structure. Tian and co-workers drew attention to the derivatives having an S=O (or SO2) and P=O group, which contains the d-pπ bond that could hinder a possible intense intermolecular π-π stacking, which is beneficial for inhibiting excited triplet-triplet annihilation (TTA) [27]. Tian et al. presented some examples of molecules with a S=O bond (in the SO2 group), which shows a sharp increase in TTA by 19 times compared to that without this group.
The X-ray bond lengths and bond angles for molecule 1 are provided in detail in Supplementary Materials. The computational DFT results confirm a bent-core structure of compound 1 and are consistent with a crystal structure in bending angles and bond length presented in Supplementary Materials. The calculated bending angles (135.4° and 135.8°, respectively) are close to X-ray data, indicating that the bending structure is not Figure 1. ORTEP (Oak Ridge thermal ellipsoid plot) drawings of (4a,12a-dihydro-12H-benzo [5,6] [1,4] thiazino [2,3-b]quinoxalin-12-yl)(phenyl)methanone (1) with 50% probability. Top view (a), angle view (b) and packing (c) [23]. Hirshfeld surfaces for molecule 1 are mapped with d norm (d), shape index (e), and curvedness (f). The X-ray bond lengths and bond angles for molecule 1 are provided in detail in Supplementary Materials. The computational DFT results confirm a bent-core structure of compound 1 and are consistent with a crystal structure in bending angles and bond length presented in Supplementary Materials. The calculated bending angles (135.4 • and 135.8 • , respectively) are close to X-ray data, indicating that the bending structure is not primarily caused by lattice packing but is due to the steric hindrance inherent in phenothiazines.
The intermolecular interactions between atoms in the crystal of molecule 1 were also studied using Hirshfeld surfaces analysis [30,31] via the d norm , shape index, and curvedness (Figure 1d,e). The surfaces are shown as transparent to visualize the molecular skeleton. The d norm surface exhibits two intensive red spots, which correspond to short intermolecular CH . . . N interactions with distances of 2.690 Å (Figure 1d). The presence of π-π stacking interactions is indicated by red and blue triangles on the shape-index surface (brown ellipses in Figure 1e) and by a flat region on the curvedness (Figure 1f). A network of intermolecular CH . . . N and π-π stacking interactions hinder molecular movement in crystal 1, which helps to reduce possible energy losses through non-radiative relaxation channels.
Molecules 2a crystallize to a monoclinic system, P2 1 /n, with unit cell dimensions of a = 8.8248 (3) (Figure 2f,g). The shape index visualizes the complementarity of molecules in the crystal and shows that there are no π-π stacking interactions in the volume (Figure 2h,i).
Compound 3 crystallizes to a monoclinic system, P2 1 /n, with unit cell dimensions of a = 8.2753(3) Å, b = 19.9331(7) Å, c = 8.5875(3) Å, V= 1398.32(9) Å 3 , and Z = 4; the monoclinic angle (β) is 99.197(3) • , and the unit cell contains 4 molecules (Figure 3c). The oxygen atom from the -SO group is disordered over two positions with a refined occupancy ratio of 0.949(5):0.051 (5). Therefore, Hirshfeld surfaces analysis is not applicable for 3. Both molecules 2a and 3 crystallize in the monoclinic system and have no π-π interactions (for 2a) or are weak (for 3, Figure 3b). This can be explained by a reduction in the electron density on the benzene rings on both sides of phenothiazine upon the N-phosphorylation for compound 2a or the oxidation of sulfur (electron-donating group) to sulfoxides (electron-withdrawing group) for molecule 3. Experimental bending angles in the middle phenothiazine ring are 98.31 • for Θ S and 124.47 • for Θ N (Figure 3d). We can find some similarities in physicochemical properties among our molecule 3 and sulfonyl derivatives published by Tian et al. [27]. They also crystallize into the monoclinic system and have abundant intermolecular C-H· · · π and C-H· · · O interactions in the crystal. The multiple intermolecular interactions help to suppress molecular motions and thereby reduce the possible energy loss via non-radiative relaxation channels.   drawing group) for molecule 3. Experimental bending angles in the middle phenothiazine ring are 98.31°for ΘS and 124.47° for ΘN (Figure 3d). We can find some similarities in physicochemical properties among our molecule 3 and sulfonyl derivatives published by Tian et al. [27]. They also crystallize into the monoclinic system and have abundant intermolecular C-H···π and C-H···O interactions in the crystal. The multiple intermolecular interactions help to suppress molecular motions and thereby reduce the possible energy loss via non-radiative relaxation channels.  In the final step, we decided to replace heteroatoms of the phenothiazine core to establish the influence of such modification on the structural properties of studied compounds. According to the literature, changing the nitrogen atom into a sulfur atom, which leads to thianthrene derivatives, would preserve the bent spatial arrangement of the atoms, hindering π-π stacking interactions. Such a configuration is abundant in molecules already described in previous paragraphs. Aiming to increase the number of π-stacks in the crystal structure, we chose more planar phenazines as synthetic targets, where nitrogen is placed instead of a sulfur atom. There are already published examples of synthesizing this class of compounds in the literature. Therefore, further experiments were conducted to explore some synthetic routes that gave the aforementioned compounds. Additionally, we tried to explain the mechanism of producing symmetrical and unsymmetrical phenazines.
We have recently reported the direct oxidative condensation to introduce the carbonnitrogen bond by a one-step method [35]. To synthesize the symmetrical phenazine 5 and the unsymmetrical phenazines 6, we applied two other types of chemical transformations. One led to symmetrical molecules and was based on the oxidative condensation of two identical molecules, i.e., 8-isopropyl-2-methylquinolin-5-amine (4a) (Scheme 3). In the second route, we chose the oxidative condensation between two electron-rich amines (ERA) N 1 ,N 1 -dimethylbenzene-1,4-diamine and selected 8-(alkyl)-2-methylquinolin-5-amines 4 (Scheme 4). We chose silver carbonate adsorbed on celite as an oxidizing agent in both methods. The transformation occurred with a good yield for symmetrical phenazine 5 (Scheme 3) and a poor yield for the unsymmetrical phenazines 6 (Scheme 4). The deposition of metallic silver during the reaction progress suggests the involvement of electron transfer processes from ERA to silver(I) salt. This assumption is grounded on the preparation of the 1,2,4,5-tetrakis(dimethylamino)benzene dication iodide-triiodide through iodination of 1,2,4,5-tetrakis(dimethylamino)benzene patented by Staab et al. [36].
(Scheme 4). We chose silver carbonate adsorbed on celite as an oxidizing agent in both methods. The transformation occurred with a good yield for symmetrical phenazine 5 (Scheme 3) and a poor yield for the unsymmetrical phenazines 6 (Scheme 4). The deposition of metallic silver during the reaction progress suggests the involvement of electron transfer processes from ERA to silver(I) salt. This assumption is grounded on the preparation of the 1,2,4,5-tetrakis(dimethylamino)benzene dication iodide-triiodide through iodination of 1,2,4,5-tetrakis(dimethylamino)benzene patented by Staab et al. [36].   Seth and co-workers reported the synthesis of symmetrical phenazines on non-radical C-H activation mode of C-N coupling during the oxidative self-coupling of anilines catalyzed by Pd-Ag and Cu-Ag nanoclusters [37,38]. They show the superiority of Pd-Ag over Cu-Ag nanoclusters during the synthesis of phenazines, with yields up to 87%. As another product, they isolated azoaryls with good yields (up to 89%) by applying the same nanoclusters, with superiority of Cu-Ag. In order to discern the character of the reaction mechanism, ionic or radical, they used radical scavenger (2,2,6,6-tetramethylpiperidin-1yl)oxyl (TEMPO). Due to the lack of evidence of the impact of TEMPO, they excluded the radical mechanism of phenazine formation. Inspired by the work of Seth and co-workers [37,38], we employed only silver carbonate as a mildly basic oxidizing agent. We observed a similar double C-N bond formation involving two different amine molecules. However, the synthesis of unsymmetrical phenazines 6 is more complicated. Both used amine molecules 4 and N 1 ,N 1 -dimethylbenzene-1,4-diamine could participate in electron transfer to silver carbonate, however, at a different rate, due to their redox potentials (Scheme 5). EPR techniques do not readily detect free radicals 4a' derived from ERA [39,40]. The evidence of their presence may be an azoarene of the 7 types (Scheme 5), the formation of which accompanied the synthesis of phenazines, similar to Seth and co-workers' studies [37,38]. For any synthesis of unsymmetrical phenazine 6, we identified compounds 6e and 7, which were confirmed by the original sample addition. This observation partially explains the much lower yields of unsymmetrical phenazines 6 (Scheme 4). On the basis of the literature and our present results, we propose a plausible mechanism for the formation of phenazines (symmetrical 5 and unsymmetrical 6) under applied reaction conditions. For simplicity, the postulated pathways are based on amine 4e (Scheme 5). The presented Seth and co-workers reported the synthesis of symmetrical phenazines on non-radical C-H activation mode of C-N coupling during the oxidative self-coupling of anilines catalyzed by Pd-Ag and Cu-Ag nanoclusters [37,38]. They show the superiority of Pd-Ag over Cu-Ag nanoclusters during the synthesis of phenazines, with yields up to 87%. As another product, they isolated azoaryls with good yields (up to 89%) by applying the same nanoclusters, with superiority of Cu-Ag. In order to discern the character of the reaction mechanism, ionic or radical, they used radical scavenger (2,2,6,6-tetramethylpiperidin-1-yl)oxyl (TEMPO). Due to the lack of evidence of the impact of TEMPO, they excluded the radical mechanism of phenazine formation. Inspired by the work of Seth and co-workers [37,38], we employed only silver carbonate as a mildly basic oxidizing agent. We observed a similar double C-N bond formation involving two different amine molecules. However, the synthesis of unsymmetrical phenazines 6 is more complicated. Both used amine molecules 4 and N 1 ,N 1 -dimethylbenzene-1,4-diamine could participate in electron transfer to silver carbonate, however, at a different rate, due to their redox potentials (Scheme 5). EPR techniques do not readily detect free radicals 4a' derived from ERA [39,40]. The evidence of their presence may be an azoarene of the 7 types (Scheme 5), the formation of which accompanied the synthesis of phenazines, similar to Seth and co-workers' studies [37,38]. For any synthesis of unsymmetrical phenazine 6, we identified compounds 6e and 7, which were confirmed by the original sample addition. This observation partially explains the much lower yields of unsymmetrical phenazines 6 (Scheme 4). On the basis of the literature and our present results, we propose a plausible mechanism for the formation of phenazines (symmetrical 5 and unsymmetrical 6) under applied reaction conditions. For simplicity, the postulated pathways are based on amine 4e (Scheme 5). The presented mechanism illustrates the influence of the silver(I) cation as an oxidizing agent and oxygen from the air. During transformation, molecules of water are created.  Figure 4 presents the molecular (a) and crystal (b) structure of phenazine 5. Figure 4a-c shows that the planar molecular skeleton of molecule 5 promotes the formation of π-stacks with intermolecular distances of 3.652-3.690 Å. Hirshfeld's dnorm surface indicates these regions as white or light blue colors in Figure 4d. More clearly, these π-stacking arrangements are identified in the shape index map with red and blue triangles (brown ellipses in Figure 4e) and the curvedness map as large green flat areas along the aromatic structural fragments of compound 5 (Figure 4f).  Figure 4 presents the molecular (a) and crystal (b) structure of phenazine 5. Figure 4a-c shows that the planar molecular skeleton of molecule 5 promotes the formation of π-stacks with intermolecular distances of 3.652-3.690 Å. Hirshfeld's d norm surface indicates these regions as white or light blue colors in Figure 4d. More clearly, these π-stacking arrangements are identified in the shape index map with red and blue triangles (brown ellipses in Figure 4e) and the curvedness map as large green flat areas along the aromatic structural fragments of compound 5 (Figure 4f).
The Hirshfeld d norm surface for molecule 6a illustrates the prominent red spot, which corresponds to the intermolecular CH· · · N interaction with a distance of 2.892 Å (Figure 5b,c). The intermolecular H . . . H and C· · · H contacts (Figure 5b,c) appear as pale red areas in the d norm . The latter belongs to the C-H· · · π interactions. The shape index (Figure 5d) confirms the presence of π-π stacking interactions as closely spaced red and blue triangles encircled by brown ellipses. Figure 5e shows π-π stacking with a flat zone in a curvedness map.  Figure 4 presents the molecular (a) and crystal (b) structure of phenazine 5. Figure 4a-c shows that the planar molecular skeleton of molecule 5 promotes the formation of π-stacks with intermolecular distances of 3.652-3.690 Å. Hirshfeld's dnorm surface indicates these regions as white or light blue colors in Figure 4d. More clearly, these π-stacking arrangements are identified in the shape index map with red and blue triangles (brown ellipses in Figure 4e) and the curvedness map as large green flat areas along the aromatic structural fragments of compound 5 (Figure 4f). Compound 6a crystallizes in the triclinic system (P-1 space group) with a = 9.0731(6) Å, b = 9.3169(7) Å, and c = 11.5962(8) Å. Figure 5 presents the molecular (a) and packing (b) structure of phenazine 6a. The distance between the centroids of the azine rings of neighboring molecules is equal to 3.628(1) Å.   Compound 6a crystallizes in the triclinic system (P-1 space group) with a = 9.0731(6) Å, b = 9.3169(7) Å, and c = 11.5962(8) Å. Figure 5 presents the molecular (a) and packing (b) structure of phenazine 6a. The distance between the centroids of the azine rings of neighboring molecules is equal to 3.628(1) Å. The Hirshfeld dnorm surface for molecule 6a illustrates the prominent red spot, which corresponds to the intermolecular CH···N interaction with a distance of 2.892 Å ( Figure  5b,c). The intermolecular H…H and C···H contacts (Figure 5b,c) appear as pale red areas in the dnorm. The latter belongs to the C-H···π interactions. The shape index (Figure 5d) confirms the presence of π-π stacking interactions as closely spaced red and blue triangles encircled by brown ellipses. Figure 5e shows π-π stacking with a flat zone in a curvedness map.
Mortzfeld et al. [41] reported the alkylation methodology of phenazine. As alkylation reagents, they chose methyl iodide, trimethyloxonium tetrafluoroborate, and dimethyl sulfate. They received stable mono N-alkylated products with good to poor yield: 86% for dimethyl sulfate, 25% for trimethyloxonium tetrafluoroborate, and only 19% for methyl iodine

Anomeric Effect of Compound 2a
The fundamental impact of stereochemical properties on a molecule's biological activity makes the conformational aspect of research essential in many fields of organic chemistry. One of the crucial phenomena altering conformational behavior and, consequently, the reactivity of saturated heterocyclic systems is the anomeric effect, also known as the Edward-Lemieux effect [42,43]. It characterizes the propensity of heteroatomic substituents neighboring a heteroatom embedded in the cyclohexane ring to favor the axial position instead of the less sterically hindered equatorial orientation. In the IUPAC Compendium of Chemical Terminology, this is "a special case of a general preference for synclinal (gauche) conformations about the bond C-Y in the system X-C-Y-C where X and Y are heteroatoms having non-bonding electron pairs, commonly at least one of which is nitrogen, oxygen or fluorine" [44].
In the 1970s, Bentrude et al. initiated conformation studies concerning saturated 1,3,2-dioxaphosphinane rings [45,46], which structures mirror features of the carbohydrate molecules. As inferred, the preferential axial position of substituents is the one adjacent to oxygen. Hence, the anomeric effect was determined to be the most dominant conformation-controlling factor for compounds with 1,3,2-dioxaphosphinane rings.
Our results show that compound 2a also undergoes the anomeric effect. Both forms of chair conformers were conveniently identified by 31 P{ 1 H}-NMR spectroscopy and X-ray diffraction measurements (Figures 2, 6 and 7). For molecule 2a, the non-bonding p oxygen orbitals interact with phosphorus, and the charge is transferred through the nitrogen atom to the thiazine-type ring, which possesses empty π-accepting orbitals. Moreover, the energies of the two possible chair conformers differ by about 1.45 kcal/mol or 6.07 kJ/mol (gas-phase B3LYP/6-31g(d,p) calculation, Figure 6), which makes the existence of an equilibrium between these two forms at room temperature possible. As shown in the 31 P{ 1 H}-NMR spectrum (Figure 7), compound 2a has two separated conformers, δ1 = −6.91 and δ2 = −7.02 ppm. The relative signal intensity of both molecules 2a is in a ratio (ca. 72:1). Mortzfeld et al. [41] reported the alkylation methodology of phenazine. As alkylation reagents, they chose methyl iodide, trimethyloxonium tetrafluoroborate, and dimethyl sulfate. They received stable mono N-alkylated products with good to poor yield: 86% for dimethyl sulfate, 25% for trimethyloxonium tetrafluoroborate, and only 19% for methyl iodine

Anomeric Effect of Compound 2a
The fundamental impact of stereochemical properties on a molecule's biological activity makes the conformational aspect of research essential in many fields of organic chemistry. One of the crucial phenomena altering conformational behavior and, consequently, the reactivity of saturated heterocyclic systems is the anomeric effect, also known as the Edward-Lemieux effect [42,43]. It characterizes the propensity of heteroatomic substituents neighboring a heteroatom embedded in the cyclohexane ring to favor the axial position instead of the less sterically hindered equatorial orientation. In the IUPAC Compendium of Chemical Terminology, this is "a special case of a general preference for synclinal (gauche) conformations about the bond C-Y in the system X-C-Y-C where X and Y are heteroatoms having non-bonding electron pairs, commonly at least one of which is nitrogen, oxygen or fluorine" [44].
In the 1970s, Bentrude et al. initiated conformation studies concerning saturated 1,3,2-dioxaphosphinane rings [45,46], which structures mirror features of the carbohydrate molecules. As inferred, the preferential axial position of substituents is the one adjacent to oxygen. Hence, the anomeric effect was determined to be the most dominant conformationcontrolling factor for compounds with 1,3,2-dioxaphosphinane rings.
Our results show that compound 2a also undergoes the anomeric effect. Both forms of chair conformers were conveniently identified by 31 P{ 1 H}-NMR spectroscopy and X-ray diffraction measurements (Figures 2, 6 and 7). For molecule 2a, the non-bonding p oxygen orbitals interact with phosphorus, and the charge is transferred through the nitrogen atom to the thiazine-type ring, which possesses empty π-accepting orbitals. Moreover, the energies of the two possible chair conformers differ by about 1.45 kcal/mol or 6.07 kJ/mol (gas-phase B3LYP/6-31g(d,p) calculation, Figure 6), which makes the existence of an equilibrium between these two forms at room temperature possible. As shown in the 31 P{ 1 H}-NMR spectrum (Figure 7), compound 2a has two separated conformers, δ 1 = −6.91 and δ 2 = −7.02 ppm. The relative signal intensity of both molecules 2a is in a ratio (ca. 72:1).
The nature of the anomeric effect is still a matter of controversy, and it is interesting to note that this challenging situation has attracted the attention of chemists and physicists for over half a century. It has been shown that steric, electrostatic, and stereoelectronic factors are all important forces that influence the anomeric effect. In the literature, there are several proposed explanations for the validity of the anomeric effect. One could be a stabilizing hyperconjugation contribution between the unshared electron pair on the oxygen atom in a 1,3,2-dioxaphosphinane ring and the σ* orbital for the P-C bond in the axial position ( Figure 6). For conformer 2, the stabilizing hyperconjugation effect is weaker. Therefore, this form is less stable and was not isolated but only identified in the 31      The occurrence of the anomeric phenomenon entails certain geometric features that are consistent with the above-noted hyperconjugative delocalization: some bonds are lengthened, others are shortened, and the angles are widened or reduced [47]. The assignments of conformation were based on NMR, X-ray, literature analysis, and our previous results [28,32,45]. According to the Karplus equation and the couplings constant, the order of chemical shifts for CH 2 protons, δ Haxial > δ Hequatorial , was determined. The methylene protons of compound 2a showed two distinct double doublets in the 1 H-NMR spectrum at δ = 3.60 and 3.74 ppm. The expected long-range couplings ( 3 J P,Haxial = 3.7 Hz and 3 J P,Hequatorial = 22.9 Hz) between phosphorus and one of the bridged methylene protons are observed. Phenothiazine 2a has two absorption bands located at 280 and 234 nm, which can be ascribed to the n-π* and π-π* transitions, respectively, typical of many phenothiazines [48]. In addition, we can observe a band below 200 nm, similar to other phenothiazines [49]. In comparison with 10H-phenothiazine (λ max = 252 nm and 316 nm [50]), the presence of N-phosphoryl substituent in compound 2a caused a hypsochromic effect. Non-planar structure of 2a makes the n-π assignment questionable. ygen atom in a 1,3,2-dioxaphosphinane ring and the σ* orbital for the P-C bond in the axial position ( Figure 6). For conformer 2, the stabilizing hyperconjugation effect is weaker. Therefore, this form is less stable and was not isolated but only identified in the 31 P{ 1 H}-NMR spectrum (Figure 7).
The occurrence of the anomeric phenomenon entails certain geometric features that are consistent with the above-noted hyperconjugative delocalization: some bonds are lengthened, others are shortened, and the angles are widened or reduced [47]. The assignments of conformation were based on NMR, X-ray, literature analysis, and our previous results [28,32,45]. According to the Karplus equation and the couplings constant, the order of chemical shifts for CH2 protons, δHaxial > δHequatorial, was determined. The methylene protons of compound 2a showed two distinct double doublets in the 1 H-NMR spectrum at δ = 3.60 and 3.74 ppm. The expected long-range couplings ( 3 JP,Haxial = 3.7 Hz and 3 JP, Hequatorial = 22.9 Hz) between phosphorus and one of the bridged methylene protons are observed.

Photophysical Properties of the Phenothiazine Derivatives 2a and 3
Figure 8a depicts the absorption spectra of molecule 2a in MeCN in the 600-190 nm region. Phenothiazine 2a has two absorption bands located at 280 and 234 nm, which can be ascribed to the n-π* and π-π* transitions, respectively, typical of many phenothiazines [48]. In addition, we can observe a band below 200 nm, similar to other phenothiazines [49]. In comparison with 10H-phenothiazine (λmax = 252 nm and 316 nm [50]), the presence of N-phosphoryl substituent in compound 2a caused a hypsochromic effect. Non-planar structure of 2a makes the n-π assignment questionable.  According to quantum chemical calculations by the TDDFT/B3LYP/6-31G(d,p) method, the weak absorption in the range 300-275 nm (absorption band 280 nm in Figure 8) in the experimental spectrum of molecule 2a is manifested as a broadening on the right side of the spectrum and corresponds to a series of low-intensity transitions S 0 → S 1 , S 0 → S 2 and S 0 → S 3 at 286 (f = 0.0148), 273 (f = 0.0141), and 265 nm (f = 0.0401), respectively. The electronic transitions and frontier molecular orbitals (HOMO-highest occupied molecular orbital; LUMO-lowest unoccupied molecular orbital) are presented in Table S1 and Figures 8b, S9 and S10 in Supplementary Materials. All these three transitions into S 1 , S 2 , and S 3 states of compound 2a possess a similar orbital nature, S 0 → S 1 (HOMO → LUMO), S 0 → S 2 (HOMO → LUMO + 1), and S 0 → S 3 (HOMO → LUMO + 2, HOMO-1→ LUMO), and are due to the charge transfer from the phenothiazine core to benzene fragments (Figure 8b). The same orbital nature of the first absorption band is also characteristic of N-methyl phenothiazine [48]. The next intense absorption band (234 nm) was calculated in the range of 250-220 nm with a maximum at 238 nm ( Figure S7), and it can be assigned to the electronic transitions into the S 4 , S 5 , and S 6 excited states (Table S1 in Supplementary Materials).
The fluorescent band of compound 2a (Figure 9) is located at 374 nm, and it is also blueshifted in comparison with the emission band of 10H-phenothiazine (λ em = 454 nm [50]) and N-methyl-phenothiazine (λ em = 450 nm [51]). The TDDFT calculations for molecule 2a clearly indicate the observed emission band at 374 nm (Figure 9), which corresponds to the vertical S 1 →S 0 transition calculated at 386 nm (Table S2, Supplementary Materials). and are due to the charge transfer from the phenothiazine core to benzene fragments (Figure 8b). The same orbital nature of the first absorption band is also characteristic of Nmethyl phenothiazine [48]. The next intense absorption band (234 nm) was calculated in the range of 250-220 nm with a maximum at 238 nm ( Figure S7), and it can be assigned to the electronic transitions into the S4, S5, and S6 excited states (Table S1 in Supplementary Materials).
The fluorescent band of compound 2a (Figure 9) is located at 374 nm, and it is also blue-shifted in comparison with the emission band of 10H-phenothiazine (λem = 454 nm [50]) and N-methyl-phenothiazine (λem = 450 nm [51]). The TDDFT calculations for molecule 2a clearly indicate the observed emission band at 374 nm (Figure 9), which corresponds to the vertical S1→S0 transition calculated at 386 nm (Table S2, Supplementary Materials). We also investigated the emission behavior of compound 2a at 77 K ( Figure 10). When excited at 291 nm, molecule 2a shows an emission band with two peaks at 434 nm and 457 nm. Its profile of phosphorescence is similar to the phosphorescence spectra of pure 10H-phenothiazine (λ 500 nm [50], λ 535 nm [52]). Similar to absorption and fluorescence, it is clear that the phosphoryl moiety in compound 2a causes a large hypsochromic effect (c.a. 70 nm) in the phosphorescence spectrum. Phosphorescence bands of molecule 2a are located in the same region as the phosphorescence of two N-acylphenothiazine We also investigated the emission behavior of compound 2a at 77 K ( Figure 10). When excited at 291 nm, molecule 2a shows an emission band with two peaks at 434 nm and 457 nm. Its profile of phosphorescence is similar to the phosphorescence spectra of pure 10H-phenothiazine (λ 500 nm [50], λ 535 nm [52]). Similar to absorption and fluorescence, it is clear that the phosphoryl moiety in compound 2a causes a large hypsochromic effect (c.a. 70 nm) in the phosphorescence spectrum. Phosphorescence bands of molecule 2a are located in the same region as the phosphorescence of two N-acylphenothiazine derivatives, namely 10H-phenothiazin-10-yl)(phenyl)methanone and 1-(10H-phenothiazin-10-yl)ethan-1-one) [27]. derivatives, namely 10H-phenothiazin-10-yl)(phenyl)methanone and 1-(10H-phenothiazin-10-yl)ethan-1-one) [27]. We also recorded the phosphorescence spectrum of compound 2a in crystals at 298 K ( Figure 11). Molecule 2a exhibits a relatively low-intensity fluorescence peak at 373 nm and an intense phosphorescence peak centered at 540 nm. Phosphorescence characteristics of N-phosphorylated phenothiazine 2a are similar to phosphorescence spectra of 10Hphenothiazin-10-yl)(phenyl)methanone and 1-(10H-phenothiazin-10-yl)ethan-1-one) that We also recorded the phosphorescence spectrum of compound 2a in crystals at 298 K ( Figure 11). Molecule 2a exhibits a relatively low-intensity fluorescence peak at 373 nm and an intense phosphorescence peak centered at 540 nm. Phosphorescence characteristics of N-phosphorylated phenothiazine 2a are similar to phosphorescence spectra of 10Hphenothiazin-10-yl)(phenyl)methanone and 1-(10H-phenothiazin-10-yl)ethan-1-one) that have phosphorescence emission at 460-600 nm, with the maximum located at 530 nm and 525 nm, respectively [8]. The TDDFT results for compound 2a show that the lowest triplet (T 1 ) state lies at 2.46 eV and involves mixed electronic configurations HOMO → LUMO, HOMO-1 → LUMO + 2, HOMO-2→ LUMO + 3, HOMO-3 → LUMO + 1. Similar to other phenothiazines [52,53], for molecule 2a, S 1 , S 2 , S 3 , and T 1 possess charge transfer characters with a certain degree of locally excited nature. Thus, a thermally activated reverse singletto-triplet intersystem crossing (RISC) can probably occur from the upper triplet to singlet levels providing phosphorescence in compound 2a since the singlet-triplet (S 1 -T 1 ) energy gap is large (0.76 eV).  The absence of π-π stacking contacts in molecule 2a and strong charge-transfer character of the T 1 state (Figure 8b) provide efficient σ-π orbital mixing, which is responsible for the increase in spin-orbit coupling and intersystem crossing in this compound affecting high phosphorescence [54]. Figure 12 shows the steady-state emission spectra of the phenothiazine derivative 2a under different conditions. Dinitrophenothiazine S-oxide 3 has four absorption bands located at 271 nm, 359 nm, 398 nm, and 528 nm (Figure 13a). In comparison with phenothiazine S-oxide (λ max 229 nm, 271 nm, 303 nm, 336 nm [55]), the two additional nitro groups presented in the phenothiazine S-oxide skeleton caused a large bathochromic effect. TDDFT calculations predict that the first high-intensity absorption maximum in molecule 3 ( Figure 13) belongs to the S 0 →S 1 transition calculated at 376 nm (f = 0.6989) in DMSO solution (Table S1). The S 0 →S 1 transition is provided by the HOMO→LUMO configuration and possesses a charge transfer (CT) nature; this is mainly CT from sulfur to nitro groups (Figure 13b). For comparison, the first S 0 →S 1 transition for phenothiazine S-oxide was calculated at 305 nm, which confirms the presence of a bathochromic effect due to nitro groups. The difference in the calculated and experimental values of the wavelengths can be explained by an incomplete allowance for solvation effects within the framework of the TDDFT/B3LYP/6-31G(d,p) level of theory. Figure 13b indicates that the HOMO is an antibonding orbital with respect to C-N and C-S bonds in the central ring, showing that there is a strong C-C bonding character, whereas the LUMO is mostly a non-bonding orbital in the center. Thus, the HOMO-LUMO excitation leads to much shorter and stronger C-N and C-S bonds in the S 1 excited state, but C-C bonds are becoming weaker ( Figure S4).   (Table S1). The S0→S1 transition is provided by the HOMO→LUMO configuration and possesses a charge transfer (CT) nature; this is mainly CT from sulfur to nitro groups (Figure 13b). For comparison, the first S0→S1 transition for phenothiazine S-oxide was calculated at 305 nm, which confirms the presence of a bathochromic effect due to nitro groups. The difference in the calculated and experimental values of the wavelengths can be explained by an incomplete allowance for solvation effects within the framework of the TDDFT/B3LYP/6-31G(d,p) level of theory. Figure 13b indicates that the HOMO is an antibonding orbital with respect to C-N and C-S bonds in the central ring, showing that there is a strong C-C bonding character, whereas the LUMO is mostly a non-bonding orbital in the center. Thus, the HOMO-LUMO excitation leads to much shorter and stronger C-N and C-S bonds in the S1 excited state, but C-C bonds are becoming weaker ( Figure S4). Compound 3 shows a fluorescent band located at 647 nm ( Figure 14). The phenothiazine S-oxide or its N-alkyl analogs exhibit fluorescence in the 350-450 nm region [56,57]. For molecule 3, TDDFT calculations indicate the observed emission band at 647 nm as the S1 → S0 transition calculated at 568 nm (Table S2). Two electron-withdrawing nitro groups remove the electron density from the sulfur atom and decrease the LUMO energy, and the emission band of dinitrophenothiazine S-oxide is red-shifted by 200 nm. Moreover, molecule 3 shows a large Stokes shift (122 nm), which indicates that the geometry of the singlet excited state differs from the geometry of the ground state. In the S0 state of molecule 3, the bending angles in the middle phenothiazine ring denoted as ΘS and ΘN were calculated to be 95.08° and 122.78°, respectively, while in S1, these angles are 104.50 and 127.27 ( Figure S4, SI). Much larger distortions are calculated in the C-N and C-S bond lengths upon the S1 → S0 transition ( Figure S4) according to the strong differences in the Compound 3 shows a fluorescent band located at 647 nm ( Figure 14). The phenothiazine S-oxide or its N-alkyl analogs exhibit fluorescence in the 350-450 nm region [56,57]. For molecule 3, TDDFT calculations indicate the observed emission band at 647 nm as the S 1 → S 0 transition calculated at 568 nm (Table S2). Two electron-withdrawing nitro groups remove the electron density from the sulfur atom and decrease the LUMO energy, and the emission band of dinitrophenothiazine S-oxide is red-shifted by 200 nm. More-over, molecule 3 shows a large Stokes shift (122 nm), which indicates that the geometry of the singlet excited state differs from the geometry of the ground state. In the S 0 state of molecule 3, the bending angles in the middle phenothiazine ring denoted as Θ S and Θ N were calculated to be 95.08 • and 122.78 • , respectively, while in S 1 , these angles are 104.50 and 127.27 ( Figure S4, SI). Much larger distortions are calculated in the C-N and C-S bond lengths upon the S 1 → S 0 transition ( Figure S4) according to the strong differences in the bonding-nonbonding characters of the frontier orbitals (Figure 13b). Typical distortions also occur in the nitro groups since LUMO is bonding for the C-N links and anti-bonding for the N-O chemical bonds. In contrast to molecule 2a and 10H-phenothiazine, the 5,5-dioxide-based push-pull chromophore [27] compound 3 shows no room-temperature phosphorescence, which is typical for nitroaromatic compounds. Asymmetric NO2 stretching and hindered rotation of nitro-groups represent accepting and promoting modes, respectively, for efficient phosphorescence quenching of compound 3.

Instrumentation
NMR spectra were obtained with Avance 400 and 500 spectrometers (Bruker, Billerica, MA, USA) operating at 500.2 or 400.2 MHz ( 1 H) and 125.8 or 100.6 MHz ( 13 C) and 202.47 ( 31 P) at 21 °C. Chemical shifts referenced to ext. TMS (tetramethylsilane) ( 1 H, 13 C) or 85% H3PO4 ( 31 P), or using the residual CHCl3 signal (δH 7.26 ppm) and CDCl3 (δC 77.1 ppm) as internal references for 1 H and 13 C-NMR, respectively. Coupling constants are given in Hz. The LCMS-IT-TOF analysis was performed on an Agilent 1200 Series binary In contrast to molecule 2a and 10H-phenothiazine, the 5,5-dioxide-based push-pull chromophore [27] compound 3 shows no room-temperature phosphorescence, which is typical for nitroaromatic compounds. Asymmetric NO 2 stretching and hindered rotation of nitro-groups represent accepting and promoting modes, respectively, for efficient phosphorescence quenching of compound 3.

Instrumentation
NMR spectra were obtained with Avance 400 and 500 spectrometers (Bruker, Billerica, MA, USA) operating at 500.2 or 400.2 MHz ( 1 H) and 125.8 or 100.6 MHz ( 13 C) and 202.47 ( 31 P) at 21 • C. Chemical shifts referenced to ext. TMS (tetramethylsilane) ( 1 H, 13 C) or 85% H 3 PO 4 ( 31 P), or using the residual CHCl 3 signal (δ H 7.26 ppm) and CDCl 3 (δ C 77.1 ppm) as internal references for 1 H and 13 C-NMR, respectively. Coupling constants are given in Hz. The LCMS-IT-TOF analysis was performed on an Agilent 1200 Series binary LC system coupled to a micrOTOF-Q system mass spectrometer (BrukerDaltonics, Bremen, Germany). High-resolution mass spectrometry (HRMS) measurements were performed using a Synapt G2-Si mass spectrometer (Waters, New Castle, DE, USA) equipped with an ESI source and quadrupole-time-of-flight mass analyzer. To ensure accurate mass measurements, data were collected in centroid mode, and mass was corrected during acquisition using leucine enkephalin solution as an external reference (Lock-Spray) (Waters, New Castle, DE, USA). The measurement results were processed using the MassLynx 4.1 software (Waters, Milford, MA, USA) incorporated within the instrument. A Nicolet iS50 FTIR spectrometer was used for recording spectra in the IR range of 4000-400 cm −1 . FTIR spectra were recorded on a Perkin Elmer (Schwerzenbach, Switzerland) spectrophotometer in the spectral range of 4000-450 cm −1 with the samples in the form of KBr pellets. Elementary analysis was performed using a Vario EL III apparatus (Elementar, Langenselbold, Germany). Differential scanning calorimetry (DSC) measurements were performed using a Q2000 calorimeter (TA Instruments, New Castle, DE, USA) in a nitrogen stream at a scanning rate of 10 • C/min. Samples were analyzed in aluminum pans in the temperature range of 50 to 350 • C. Melting points were determined on an MPA100 OptiMelt melting point apparatus (Stanford Research Systems, Sunnyvale, CA, USA) and are uncorrected. UV/Vis absorption spectra were measured on a SHIMADZU UV-VIS-NIR spectrophotometer using quartz (Suprasil) cuvettes (10 mm path length). Fluorescence spectra were recorded on an FLS-900 spectrofluorimeter (Edinburgh Instruments, Edinburgh, UK) and were measured in quartz (Suprasil) cuvettes (10 mm path length) at 25 • C. The phosphorescence spectrum was recorded on an FLS-900 spectrofluorimeter (Edinburgh Instruments, Edinburgh, UK) at 77 K using a 3 mm (inner diameter) quartz tube inside a quartz liquid nitrogen Dewar flask.
To a suspension of NaH (0.96 g, 40.0 mmol) in THF (100 mL), 10H-phenothiazine (4.00 g, 20.1 mmol) was added and stirred until the evolution of H 2 ceased. Reagents were stirred under reflux for 30 min. under argon. 2-Chloro-5,5-dimethyl-1,3,2-dioxaphosphinane 2-oxide (3.70 g, 20.1 mmol) was then added to the reaction mixture, which was refluxed overnight. After evaporation of the solvent to give a solid, water (20 mL) and chloroform (100 mL) were added. The organic layer was separated, and the aqueous layer was extracted four times with chloroform. The combined organic layers were dried over MgSO 4 . After the solvent evaporated, the crude product was purified by column chromatography on silica gel using methanol/dichloromethane as eluent to afford a crude solid, and finally, crystallization from a mixture of CH 2 Cl 2 and hexane occurred to yield solids as follows: 10H-Phenothiazine (10.0 g, 50.2 mmol) was added in one portion to 100 mL of 65% nitric acid stirred in a wide beaker at room temperature. The immediate evolution of brown vapor was observed, and the reaction mixture's temperature increased to 45 • C.
After stirring for 1h, the deep red solution was filtered through a sintered glass funnel under reduced pressure, and the obtained solid was flushed with 200 mL of water and 100 mL of methanol. The crude product was boiled for 5 min in 300 mL of THF, filtered, and dried under reduced pressure. To obtain the crystals, the saturated solution of the purified product in CHCl 3 was kept at −20 • C for a week.
Synthesis of 2,8-dimethylquinoline. Toluene (50 mL) and crotonaldehyde (2.6 mL, 2.2 g, 31.4 mmol) were added to the solution of o-toluidine (1.7 g, 15.7 mmol) in aqueous 6 M HCl (200 mL) and were heated under reflux for 16 h. The mixture was allowed to cool down to room temperature. The aqueous layer was separated and neutralized with an aqueous solution of K 2 CO 3 . After extraction with CH 2 Cl 2 (3 × 50 mL), the organic layer was separated and dried over MgSO 4 , then filtered and distilled bp 100-110 • C/3 mmHg. The liquid mixture was dissolved in concentrated 36% HCl (100 mL) at 5 • C, and ZnCl 2 (2.7 g, 20.0 mmol) was added with vigorous stirring for one hour. The precipitate was filtered and washed with cold 3 M aq. HCl and dried in air. The solid was washed with iPrOH and dried. The received white solid was added to a 10% ammonia solution and extraction with Et 2 O (3 × 50 mL). The organic layer was separated and dried over MgSO 4 4 and HNO 3 (4.5 and 10.5 mL, respectively) at 5 • C. After stirring for 1 h at room temperature, no evolution of gas was observed, so the reaction mixture was heated up to 70 • C and stirred overnight. After this time, the reaction mixture was poured down into a beaker containing 25 g of ice and 25 mL of water, and the precipitated solid was filtered off, washed with 10 mL of cold water, and dried on air, giving: 2,8-Dimethyl-5-nitroquinoline as a beige solid, 1.  The solution of 8-isopropyl-2-methylquinolin-5-amine (4a) (400 mg, 2.00 mmol) in toluene (50 mL) was brought to gentle reflux, and silver carbonate on celite (2.5 g, 50% w/w, 4.50 mmol Ag 2 CO 3 ) was added in small portions over 30 min. The reaction mixture was stirred under reflux for a further 4 h and filtered under reduced pressure without cooling. The solid residue was flushed with ca. 50 mL of ethyl acetate until a clean, colorless filtrate was observed. The obtained solution was left to stand in the refrigerator overnight. The next day, the formed light yellow precipitate was collected by filtration, flushed with warm methanol, and recrystallized from chloroform to afford yellow crystals:

Synthesis of Unsymmetrical Phenazines (6a-d)
The solution of appropriate quinoline (4a-c) or isoquinoline 4d (2.00 mmol) and N 1 ,N 1 -dimethylbenzene-1,4-diamine (286 mg, 2.10 mmol) in toluene (50 mL) was brought to a gentle reflux, and silver carbonate on celite (2.5 g, 50% w/w, 4.50 mmol Ag 2 CO 3 ) was added in small portions over 30 min. The reaction mixture was stirred under reflux for a further 4 h and filtered under reduced pressure without cooling. The solid residue was 3.9. Synthesis of 6e and 7 The solution of N 1 ,N 1 -dimethylbenzene-1,4-diamine (286 mg, 2.10 mmol) in toluene (50 mL) was brought to a gentle reflux, and silver carbonate on celite (2.5 g, 50% w/w, 4.50 mmol Ag 2 CO 3 ) was added in small portions over 30 min. The reaction mixture was stirred under reflux for 4 h and filtered under reduced pressure without cooling. The solid residue was flushed with ca. 50 mL of ethyl acetate until a clean, colorless filtrate was observed. The obtained solution was concentrated and purified by column chromatography on silica gel using chloroform and/or ethyl acetate/hexane eluent system to afford: N 2 ,N 2 ,N 7 ,N 7 -tetramethylphenazine-2,7-diamine (6e) as a red powder, 90 mg (0. 34

Hirshfeld Surface Calculations
Hirshfeld surfaces were calculated using Crystal Explorer 21.5 [31] to analyze the close contacts in the phenothiazine and phenazine compounds. The normalized contact distance (d norm ) values are mapped with red-blue-white colors, where the closer contacts with d norm < 0 are represented in red areas, while the blue region indicates longer contacts with d norm > 0, and white areas correspond to contacts with distances equal to the van der Waals radii (d norm = 0). The shape index visualizes complementary molecules in the crystal with red (hollows) and blue (bumps) colors, where molecular surfaces touch one another. Curvedness maps show planar regions on the surface in green, which is separated by blue edges corresponding to large positive curvature.

Density Functional Theory (DFT) Computations
The molecular structures of the phenothiazines and phenazines in the ground singlet state (S 0 ) were optimized with the DFT/B3LYP/6-31G(d,p) approach [60][61][62][63]. For the calculated molecules, all vibrational wavenumbers were found to be real, indicating the location of the true minimum on the hypersurface of the total energy. The first excited singlet (S 1 ) and triplet (T 1 ) states, as well as vertical S 0 -S n transitions, were calculated within the time-dependent (TD) DFT approximation [64] with the same B3LYP/6-31G(d,p) method [61][62][63]. The solvent effects on the absorption and fluorescence spectra were taken into account using the polarized-continuum model [65]. The Gaussian 16 program package [66] was used for all the DFT and TDDFT calculations.

X-ray Diffraction Experiments
The data for 1, 2a, 3, 5, and 6a were collected using a SuperNova diffractometer (Agilent Technologies, currently Rigaku Oxford Diffraction). Accurate cell parameters were determined and refined using the CrysAlis Pro program [67]. Additionally, integration of the collected data was performed with this program. The data for 2b were collected using an Xcaliburdiffractometer (Oxford Diffraction, currently Rigaku Oxford Diffraction). Determination and refinement of the cell parameters, as well as the integration of the collected data, were performed using the CrysAlis Pro program [68]. All the structures were solved using direct methods with the SHELXS97 program and then refined using the SHELXL-2018/3 program [69]. Nonhydrogen atoms were refined with anisotropic displacement parameters. The hydrogen atoms were fixed at calculated distances and allowed to ride on the parent atoms.

Conclusions
A series of N-heteroaromatics have been strategically designed and synthesized to compare their spectroscopic properties and rationalize them on the basis of theoretical calculations (Scheme 1). We described and analyzed the crystal structure of the N-benzoylated phenothiazine-phenazine hybrid molecule 1. Compounds of this type possess rich future chemistry potential for photoluminescent materials [27] and polymerization photoinitiating systems [23]. The presented synthesis protocols allowed obtaining the symmetrical phenazines with good yields (molecules 5 and 6e) and with low yields for unsymmetrical representatives 6a-d, giving six novel compounds. The crystal structures of phenazines 5 and 6a were determined by single-crystal X-ray diffraction measurements. Interestingly, molecules 6a-d are structurally related to the well-known biological stain neutral red, the analogs of which exhibit promising photosensitizing [70] and electrochemical properties [71].
Our results show that phosphorylation of phenothiazine proceeds exclusively to one new product 2a, with a newly formed P-N bond together with rearrangement leading to a sulfurylated thiophosphoric acid 2b. The S 1 and T 1 excited states of compound 2a possess hybridized local and charge-transfer characters and show excellent luminescent properties. This benefit makes molecule 2a a potentially highly luminous molecule. The N-phosphoryl moiety does not participate directly in the low-energy excitations of molecule 2a (it does not contribute to frontier molecular orbitals, Figure 8b) but provides a hypsochromic shift. According to the energy gap law [72], this blue shift leads to hindered non-radiative quenching of the T 1 state and to sustainable phosphorescence.
Comparing the molecular orbitals involved in excitations for 2a and 3, it should be noted that for both phenothiazines, these MOs are localized in the phenothiazine ring, but their shapes and charge transfer nature differ significantly. The presence of phosphorus in molecule 2a provides a strong stabilizing hyperconjugation effect, and the charge is transferred through the nitrogen atom to the thiazine-type ring. Additionally, the absence of π-π stacking interactions in molecule 2a and the strong charge-transfer character of the T 1 state provide efficient σ-π orbital mixing responsible for the increase in spin-orbit coupling and intersystem crossing in this compound, which contributes to high phosphorescence. Phenothiazine 3 has no stabilizing hyperconjugation effect; the first S 0 →S 1 transition possesses a charge transfer character mainly from sulfur to nitro groups. In contrast to molecule 2a, compound 3 does not exhibit phosphorescence at room temperature, which is typical for nitroaromatic compounds.
Supplementary Materials: The following supporting information can be downloaded at https: //www.mdpi.com/article/10.3390/molecules27217519/s1: Tables S1-S16 and Figures S1-S18. The following are available online: CCDC 2177507 for 1, CCDC 2177508 for 2a, CCDC 2178722 for 2b, CCDC 2177509 for 3, CCDC 2177510 for 5 and CCDC 2177511 for 6a contains the supplementary crystallographic data for this paper, accessed on 13 June 2022. These data can be obtained free of charge from http://www.ccdc.cam.ac.uk/conts/retrieving.html (accessed on 13 June 2022) or from the Cambridge Crystallographic Data Centre, 12 Union Road, Cambridge CB2 1EZ, UK; Fax: (+44) 1223-336-033; or e-mail: deposit@ccdc.cam.ac.uk. Calculations have been carried out at the Wroclaw Centre for Networking and Supercomputing (http://www.wcss.wroc.pl (accessed on 13 June 2022)).  Data Availability Statement: The data set presented in this study is available in this article.