Intramolecular sp2-sp3 Disequalization of Chemically Identical Sulfonamide Nitrogen Atoms: Single Crystal X-Ray Diffraction Characterization, Hirshfeld Surface Analysis and DFT Calculations of N-Substituted Hexahydro-1,3,5-Triazines

In this manuscript, the synthesis and single crystal X-ray diffraction characterization of four N-substituted 1,3,5-triazinanes are reported along with a detailed analysis of the noncovalent interactions observed in the solid state architecture to these compounds, focusing on C–H···π and C–H···O H-bonding interactions. These noncovalent contacts have been characterized energetically by using DFT calculations and also by Hirshfeld surface analysis. In addition, the supramolecular assemblies have been characterized using the quantum theory of “atoms-in-molecules” (QTAIM) and molecular electrostatic potential (MEP) calculations. The XRD analysis revealed a never before observed feature of the crystalline structure of some molecules: symmetrically substituted 1,3,5-triazacyclohexanes possess two chemically identical sulfonamide nitrogen atoms in different sp2 and sp3-hybridizations.


Experimental Details
As it was mentioned above, the main objects of this work, N,N'-disulfamide substituted triazinanes 1-4, were prepared according to the method described in our preliminary communication [12] using the three-component Mg(ClO4)2 catalyzed condensation of arylsulfonamides with paraformaldehyde and n-or tert-butyl amine (Scheme 2, see the Supplementary Materials for detail of the experimental procedures and spectral data). The tert-butyl-and n-butylamines were chosen as the amino-components providing the highest yield of the target triazinanes. All obtained triazinanes are well-crystallized solids that allowed the growth of crystals suitable for XRD analysis.

Crystallographic Details
Single crystal X-ray diffraction experiments were performed at the Center for Shared Use of Physical Methods of Investigation at the Frumkin Institute of Physical Chemistry and Electrochemistry. The single crystal X-ray diffraction data for 1,3,5-triazacyclohexanes (1)(2)(3)(4) were collected on a Bruker Kappa Apex II automatic four-circle diffractometer (Bruker AXS, Madison, WI, USA) equipped with an area detector (Mo-Kα sealed-tube X-ray source, λ = 0.71073 Å, graphite monochromator) at 100 K for all compounds. The principal crystallographic data and structural refinements are summarized in Table 1. Atomic coordinates for compounds 1-4, have been deposited with the CCDC (number 1992667−1992670). The supplementary crystallographic data are available in the ESI section. The comparison of the crystal structure parameters with the analogous compounds were performed using ConQuest search in Cambridge Structural Database (CSD, Version 5.40). The histograms of angles values were obtained from a graphical search of Scheme 1. Compounds 1-4 studied in this work.

Experimental Details
As it was mentioned above, the main objects of this work, N,N -disulfamide substituted triazinanes 1-4, were prepared according to the method described in our preliminary communication [12] using the three-component Mg(ClO 4 ) 2 catalyzed condensation of arylsulfonamides with paraformaldehyde and nor tert-butyl amine (Scheme 2, see the Supplementary Materials for detail of the experimental procedures and spectral data). The tert-butyland n-butylamines were chosen as the amino-components providing the highest yield of the target triazinanes.

Experimental Details
As it was mentioned above, the main objects of this work, N,N'-disulfamide substituted triazinanes 1-4, were prepared according to the method described in our preliminary communication [12] using the three-component Mg(ClO4)2 catalyzed condensation of arylsulfonamides with paraformaldehyde and n-or tert-butyl amine (Scheme 2, see the Supplementary Materials for detail of the experimental procedures and spectral data). The tert-butyl-and n-butylamines were chosen as the amino-components providing the highest yield of the target triazinanes. All obtained triazinanes are well-crystallized solids that allowed the growth of crystals suitable for XRD analysis.

Crystallographic Details
Single crystal X-ray diffraction experiments were performed at the Center for Shared Use of Physical Methods of Investigation at the Frumkin Institute of Physical Chemistry and Electrochemistry. The single crystal X-ray diffraction data for 1,3,5-triazacyclohexanes (1-4) were collected on a Bruker Kappa Apex II automatic four-circle diffractometer (Bruker AXS, Madison, WI, USA) equipped with an area detector (Mo-Kα sealed-tube X-ray source, λ = 0.71073 Å, graphite monochromator) at 100 K for all compounds. The principal crystallographic data and structural refinements are summarized in Table 1. Atomic coordinates for compounds 1-4, have been deposited with the CCDC (number 1992667−1992670). The supplementary crystallographic data are available in the ESI section. The comparison of the crystal structure parameters with the analogous compounds were performed using ConQuest search in Cambridge Structural Database (CSD, Version 5.40). The histograms of angles values were obtained from a graphical search of sulfonamides (C-S(=O)2-NC2) with 3D parameters for angles. More than 7000 hits were analyzed. All obtained triazinanes are well-crystallized solids that allowed the growth of crystals suitable for XRD analysis.

Crystallographic Details
Single crystal X-ray diffraction experiments were performed at the Center for Shared Use of Physical Methods of Investigation at the Frumkin Institute of Physical Chemistry and Electrochemistry. The single crystal X-ray diffraction data for 1,3,5-triazacyclohexanes (1-4) were collected on a Bruker Kappa Apex II automatic four-circle diffractometer (Bruker AXS, Madison, WI, USA) equipped with an area detector (Mo-Kα sealed-tube X-ray source, λ = 0.71073 Å, graphite monochromator) at 100 K for all compounds. The principal crystallographic data and structural refinements are summarized in Table 1. Atomic coordinates for compounds 1-4, have been deposited with the CCDC Crystals 2020, 10, 369 3 of 14 (number 1992667−1992670). The supplementary crystallographic data are available in the ESI section. The comparison of the crystal structure parameters with the analogous compounds were performed using ConQuest search in Cambridge Structural Database (CSD, Version 5.40). The histograms of angles values were obtained from a graphical search of sulfonamides (C-S(=O) 2 -NC 2 ) with 3D parameters for angles. More than 7000 hits were analyzed.

Hirshfeld Surface Calculations
The Hirshfeld surface (HS) analysis [13][14][15] and their associated 2D fingerprint plots (full and decomposed) [16] were carried out employing the CrystalExplorer 17 program [17] in order to visualize and quantify various non-covalent interactions that stabilize the crystal packing. The HS was mapped over d norm property. The d norm property is a symmetric function of distances to the surface from nuclei inside and outside the Hirshfeld surface (di and de, respectively), relative to their respective van der Waals radii. The regions with red and blue color on the d norm represent the shorter and longer inter contacts while the white color indicates the contacts around the van der Waals radii. 2D fingerprint plots provide relevant information of intermolecular contacts in the crystal. The d norm surface was mapped with the color scale in the range −0.050 au (red) to 0.600 au (blue). 2D fingerprint plots (di vs. de) were displayed using the expanded 0.6-2.8 Å range.

Theoretical Methods
All DFT calculations were carried out using the  at the PBE1PBE-D3/def2-TZVP level of theory and using the crystallographic coordinates. The formation energies of the assemblies were evaluated by calculating the difference between the total energy of the assembly and the sum of the monomers that constitute the assembly, which have been maintained frozen. That is ∆E AB = E AB − E A − E B , where ∆E AB is the interaction energy; E AB is the energy of the dimer and E A and E B are the energy of the monomers. The BSSE has been used to correct the interaction energies by using the counterpoise =2 keyword in the Gaussian-19 program [18]. The molecular electrostatic potential was computed at the same level of theory and plotted onto the 0.001 a.u. isosurface. The Quantum Theory of Atoms-in-Molecules (QTAIM) [19] analysis was carried Crystals 2020, 10, 369 4 of 14 out at the same level of theory by means of the AIMAll program [20] to obtain the distribution of bond critical points (CPs) and bond paths [21].

Structural Description
According to the single crystal X-ray diffraction data, molecules 1-4 comprise the 1,3,5-triazacyclohexane ring bearing three substituents at the nitrogen atoms (see Figure 1).  Compound 1 crystalizes in the monoclinic crystal system in the space group P21/c. Selected bond lengths, angles and torsional angles are shown in Tables S2-S5. The S2-C21 and S1-N3 bond lengths of 1.758(1) and 1.642(1) Å respectively are in agreement with a single bond character of these bonds [24]. The S-O bond distances of the sulfamide moiety are in the range 1.432(1)-1.436(1) Å, which is typical for N-sulfamides. The SO2 group has typical angles if compare with CSD data for N-sulfamides ( Figure S1 in Supplementary Materials The crystal structure of this compound exhibits interesting assemblies in the solid state (see The asymmetric unit contains one molecule of each triazine. The general geometrical features of these systems are similar, which are the slightly distorted chair conformation of the six-membered ring, with the N-alkyl substituents (n-Bu or t-Bu) occupying the axial position and the two N-sulfamide fragments occupying the sterically favorable pseudo-equatorial orientation (see Table S2 for torsion angles of N-substituents). The atoms C2, C6, N3, and N5 of the triazinane cycle in all structures lie nearly in one plane (the deviation of one atom from the plane of the other three is less than 0.03 Å), while the deviation of N1 and C4 atoms from this plane range from −0.633(2) Å to −0.661(2) Å and from 0.643(14) Å to 0.712(2) Å, correspondingly, therefore the molecules have a classical chair-conformation of the central heterocycle (see Scheme 2 for atom numbering scheme). All CH 2 -N bond lengths and CH 2 -N-CH 2 bond angles are typical for 1,3,5-triazinanes and are listed in the corresponding tables in (Supplementary Materials Table S2-S15). The torsion angles CNCN in the triazacyclohexane ring are close to 60 • (see Table S2).
The most intriguing and distinguishing feature of the triazinanes under discussion is the unprecedented geometry of the N3 and N5 nitrogen atoms of the sulfonamide fragments in the N-butyl substituted heterocycles 1 and 2. As is generally known, the nitrogen atom in a sulfonamide group can adopt both sp 2 and sp 3 hybridization depending on substituents at the nitrogen atom [22][23][24]. However, according to the data of the CCDC, there are no known examples of 1,3,5-triazinanes or other saturated six-membered azaheterocycles simultaneously possessing two chemically identical sulfonamide nitrogen atoms in different hybridization. Analysis of the values of the sums of valence angles at nitrogen atoms in positions 3 and 5 allows to clearly identify atoms in sp 2 or in sp 3 hybridization ( Figure 1 and Table 2). The N3-atoms in compounds 1-3 are sp 2 -hybridizated and, as a result, they assume the flat trigonal configuration (the sum of the angles is close to 360 • ). N5-Atoms, chemically equivalent to N3-atoms, in the same molecules adopt the tetrahedral configuration (the sum of the angles lies in the diapason of 344-351 • ) and, therefore, are sp 3 -hybridizated. This is most clearly seen in the examples of compounds 1 and 2. The presence of the tert-butyl group at the N1 position in compounds 3 and 4, probably due to its high steric volume, symmetrizes the molecules, leveling the difference between both sulfonamide nitrogen atoms in a crystal. This is also observed in the equalization of the S-N distances of the sulfonamide groups in compounds 3 and 4 compared to 1 and 2, see Table 2. The largest difference between both S-N distances is observed in compound 2, i.e., 0.016 Å. In fact, the short S2-N5 distance is an indication of a partial double bond character, in agreement with the sp 2 -hybridization. Compound 1 crystalizes in the monoclinic crystal system in the space group P2 1 /c. Selected bond lengths, angles and torsional angles are shown in Tables S2-S5. The S2-C21 and S1-N3 bond lengths of 1.758(1) and 1.642(1) Å respectively are in agreement with a single bond character of these bonds [24]. The S-O bond distances of the sulfamide moiety are in the range 1.432(1)-1.436(1) Å, which is typical for N-sulfamides. The SO 2 group has typical angles if compare with CSD data for N-sulfamides ( Figure S1 in Supplementary Materials): O-S-O angle is around 120 • [O1-S1-O2 is 119.95(6) • ], while N-S-O angles are between 105 • and 115 • [N3-S1-O2 is 106.07(6) • , N3-S1-O1 is 106.63(5) • ].
The crystal structure of this compound exhibits interesting assemblies in the solid state (see Table S3 for H-bonds). For instance, Figure 2a shows a self-assembled dimer dominated by C-H···O interactions where the methyl group in para acts as a H-bond donor. The acidity of these protons is higher than usual for a methyl group due to the presence of the electron withdrawing sulfamide group. Quite remarkable is the ternary assembly shown in Figure 2b, where the aromatic π-cloud interacts simultaneously with the methyl group at one side and an aromatic C-H bond at the opposite side, thus forming a C-H···π···H-C assembly. It is worth noting that the C-H···π distances are very short (2.60 and 2.70 Å) thus, confirming their relevance in the solid state of this compound.
Compound 2 crystalizes in the monoclinic crystal system in the space group P2 1 /n and the main difference with compound 1 is the absence of methyl groups. In addition, the sulfonamide groups are attached to the 1,3,5-triazacyclohexane ring in different orientations as reflected by the C6-N5-S2-C21 and C2-N3-S1-C11 torsional angles of −105.0(1) • and −65.5(1) • , respectively. It also forms self-assembled dimers in the solid state, where both C-H···π and C-H···O (Table S3) interactions are established, as shown in Figure 3a. Moreover, it also forms infinite 1D supramolecular chains in the solid state promoted by C-H···O interactions involving the butyl chain and the sulfonamide group (see Figure 3b).  Compound 2 crystalizes in the monoclinic crystal system in the space group P21/n and the main difference with compound 1 is the absence of methyl groups. In addition, the sulfonamide groups are attached to the 1,3,5-triazacyclohexane ring in different orientations as reflected by the C6-N5-S2-C21 and C2-N3-S1-C11 torsional angles of −105.0(1)° and −65.5(1)°, respectively. It also forms self-assembled dimers in the solid state, where both C-H···π and C-H···O (Table S3) interactions are established, as shown in Figure 3a. Moreover, it also forms infinite 1D supramolecular chains in the solid state promoted by C-H···O interactions involving the butyl chain and the sulfonamide group (see Figure 3b).   Compound 2 crystalizes in the monoclinic crystal system in the space group P21/n and the main difference with compound 1 is the absence of methyl groups. In addition, the sulfonamide groups are attached to the 1,3,5-triazacyclohexane ring in different orientations as reflected by the C6-N5-S2-C21 and C2-N3-S1-C11 torsional angles of −105.0(1)° and −65.5(1)°, respectively. It also forms self-assembled dimers in the solid state, where both C-H···π and C-H···O (Table S3) interactions are established, as shown in Figure 3a. Moreover, it also forms infinite 1D supramolecular chains in the solid state promoted by C-H···O interactions involving the butyl chain and the sulfonamide group (see Figure 3b).  Compound 3 crystalizes in the monoclinic crystal system in the space group P2 1 /n and, similarly to compound 1, also forms ternary assemblies where the same aromatic ring establishes C-H···π interactions at both sides of the ring, thus forming a C-H···π···H-C assembly (see Figure 4). This compound also forms 1D supramolecular chains in the solid state where the H-atoms of the triazinane ring interact with the O-atoms of the sulfonamide groups (Table S3), as shown in Figure 4b. The formation of this assembly is facilitated by the relative orientation of the p-methyl-benzene-sulfonamide groups.
to compound 1, also forms ternary assemblies where the same aromatic ring establishes C-H···π interactions at both sides of the ring, thus forming a C-H···π···H-C assembly (see Figure 4). This compound also forms 1D supramolecular chains in the solid state where the H-atoms of the triazinane ring interact with the O-atoms of the sulfonamide groups (Table S3), as shown in Figure  4b. The formation of this assembly is facilitated by the relative orientation of the p-methyl-benzene-sulfonamide groups. Compound 4 crystallizes in the orthorhombic crystal system in the space group P212121 and also exhibits several motifs in the solid state that are mainly dominated by C-H···O interactions (Table  S3). As examples, two motifs are given in Figure 5, one corresponds to a discrete dimer where three H-bonds are formed and the other one to a 1D supramolecular polymer also governed by the formation of C-H···O bonds involving the aromatic H-atom located in para to the sulfonamide group (Figure 5b). Compound 4 crystallizes in the orthorhombic crystal system in the space group P2 1 2 1 2 1 and also exhibits several motifs in the solid state that are mainly dominated by C-H···O interactions (Table S3). As examples, two motifs are given in Figure 5, one corresponds to a discrete dimer where three H-bonds are formed and the other one to a 1D supramolecular polymer also governed by the formation of C-H···O bonds involving the aromatic H-atom located in para to the sulfonamide group (Figure 5b).

Hirshfeld Surfaces
The Hirshfeld surface analysis is a very convenient tool for analyzing intermolecular interactions. The HS surfaces mapped over dnorm property are displayed in Figure 6 highlighting the main intermolecular interactions and scheme of labels. The patterns of intermolecular interactions are similar in all structures, which prompted us to evaluate the contributions of the weak non-covalent contacts in the supramolecular assembly, as well as the importance of C-H···π interactions in stabilization of the crystal packing. The 2D fingerprint plots (Figure 7) of the

Hirshfeld Surfaces
The Hirshfeld surface analysis is a very convenient tool for analyzing intermolecular interactions. The HS surfaces mapped over d norm property are displayed in Figure 6 highlighting the main intermolecular interactions and scheme of labels. The patterns of intermolecular interactions are similar in all structures, which prompted us to evaluate the contributions of the weak non-covalent contacts in the supramolecular assembly, as well as the importance of C-H···π interactions in stabilization of the crystal packing. The 2D fingerprint plots (Figure 7) of the molecules illustrate significant differences between the intermolecular interaction patterns. The surfaces are shown as transparent to allow the visualization of the molecules. Contacts with distances equal to the sum of van der Waals (vdW) radii are represented as white regions and contacts with distances shorter than and longer than the vdW radii are shown as red and blue colors, respectively. Crystals 2020, 10, x FOR PEER REVIEW 9 of 15 HS area. A similar behavior was observed in the HS mapped over dnorm function for (4), which the six red spots observed (Figure 6d) are attributed to C26-H26···O4 (1), C16-H16···O4 (2), C24,H24···O2 (3), C12-12···O3 (4) and C6-H6A···O3 (5). Intermolecular interactions O···H/H···O are observed around 2.3 Å which is slightly shorter than those of other compounds with 30% contribution to the Hirshfeld contact surface. As was described previously, the structure of (1) is also stabilized by C-H···π interactions. The red area labeled 5 in the HS mapped over dnorm is attributed to C9-H9B···π. These C-H···π interactions are also evident from a pair of "wings" in the top left and bottom right region of the fingerprint plots for compounds 1-4 ( Figure 7). The shape of the wings and the sum of de and di show the importance of this interaction. The decomposition of the fingerprint plots shows that the C···H/H···C contributions comprising 17.4%, 17.2%, 18.8% and 18.7% of the total HS for each molecule of 1, 2, 3 and 4, respectively.

DFT Calculations
The DFT study is focused to analyze the supramolecular assemblies commented above in Figures 2-5, where combinations of C-H···π and C-H···O H-bonding networks are commonly formed in compounds 1-4. First of all, the molecular electrostatic potential (MEP) surfaces of compounds 1 and 2 have been computed in order to analyze the electron rich and electron poor regions of the molecules. The surfaces are represented in Figure 8 and it can be observed that the most negative regions correspond to the O-atoms of the sulfonamide group. The N-atoms of the triazinane ring are not good H-bond acceptors, likely because either the lone pair is delocalized into the SO2-group, in accordance with the S2-N5 and S1-N3 bond lengths of 1.632(1) and 1.642(1) Å, respectively indicating a double bond character of these bonds. The most positive region corresponds to the middle of the three axial H-atoms of the triazinane ring (+27 kcal/mol). The aromatic H-atoms and the CH3 substituents also present positive MEP values (+23 and +20 kcal/mol, The vdW forces (H···H contacts) have the largest contribution to the HS, and they are highlighted in the scattered middle points in the fingerprint plots with a minimum value of (d e + d 1 )~2.2 Å (Figure 7), which is the sum of the vdW radii. All red areas that are visible on the surfaces mapped over d norm function correspond to C-H···O contacts. For (1), the largest bright-red spot labeled 1 on the HS shows O···H/H···O contact associated with C17-H17C···O1 interaction, which constitutes the strongest among all interactions present in this compound. Two medium sized red spots labeled as 2 and 3 are associated with C15-H15A···O1 and C27-H27C···O4, respectively. These interactions are also visible as symmetrical sharp spikes centered at (de + di) 2.4 Å in the fingerprint plots (Figure 7a) with 25.2% contribution to the Hirshfeld contact surface. The intermolecular C25-H25A···N1 contact is visible in the HS as a red spot labelled 4, which comprises 2.3% of the total HS area. The HS of (2) mapped over d norm function (see Figure 6) shows four red spots, indicating the presence of C-H···O hydrogen bonds [C13-H13A···O1 (1), C26-H16A···O4 (2), C10-H10B···O1 (3), C14-H14A···O2 (4)]. The decomposed fingerprint plots (Figure 7b) show that intermolecular O···H/H···O contacts contribute 27.7% to the total HS area. The O···H/H···O contacts appeared as sharp spikes with (de + di) 2.35 Å. In the HS of (3), the O···H/H···O contacts ( Figure 6) are visible as six red spots attributed to C6-H6A···O3 (labeled 1), C10-H10B···O3 (labeled 1), C15-H15A···O4 (labeled 2) and C17-H17C···O4 (labeled 2). These interactions comprise the 24.1% of the HS area. A similar behavior was observed in the HS mapped over d norm function for (4), which the six red spots observed ( Figure 6) are attributed to C26-H26···O4 (1), C16-H16···O4 (2), C24,H24···O2 (3), C12-12···O3 (4) and C6-H6A···O3 (5). Intermolecular interactions O···H/H···O are observed around 2.3 Å which is slightly shorter than those of other compounds with 30% contribution to the Hirshfeld contact surface.
As was described previously, the structure of (1) is also stabilized by C-H···π interactions. The red area labeled 5 in the HS mapped over d norm is attributed to C9-H9B···π. These C-H···π interactions are also evident from a pair of "wings" in the top left and bottom right region of the fingerprint plots for compounds 1-4 ( Figure 7). The shape of the wings and the sum of d e and d i show the importance of this interaction. The decomposition of the fingerprint plots shows that the C···H/H···C contributions comprising 17.4%, 17.2%, 18.8% and 18.7% of the total HS for each molecule of 1, 2, 3 and 4, respectively.

DFT Calculations
The DFT study is focused to analyze the supramolecular assemblies commented above in  Figure 8 and it can be observed that the most negative regions correspond to the O-atoms of the sulfonamide group. The N-atoms of the triazinane ring are not good H-bond acceptors, likely because either the lone pair is delocalized into the SO 2 -group, in accordance with the S2-N5 and S1-N3 bond lengths of 1.632(1) and 1.642(1) Å, respectively indicating a double bond character of these bonds. The most positive region corresponds to the middle of the three axial H-atoms of the triazinane ring (+27 kcal/mol). The aromatic H-atoms and the CH 3 substituents also present positive MEP values (+23 and +20 kcal/mol, respectively). Finally, the MEP value over the aromatic rings is negative (−8 kcal/mol), thus adequate for establishing C-H···π interactions. The MEP analysis evidences that the most favored interactions from an electrostatic point of view are those involving the O-atoms as electron donors and either aromatic or triazinane protons as electron acceptors.  We have selected two supramolecular assemblies commented above in Figure 2 to analyze the energetic features of the H-bonds and C-H···π interactions in 1. The QTAIM distribution of critical points and bond paths are also given in Figure 9. The existence of a bond CP and bond path connecting two atoms is a good indicator of interaction [21]. For the self-assembled dimer (Figure  9a), in addition to the symmetrically related H-bonds (characterized by a bond CP (critical point) We have selected two supramolecular assemblies commented above in Figure 2 to analyze the energetic features of the H-bonds and C-H···π interactions in 1. The QTAIM distribution of critical points and bond paths are also given in Figure 9. The existence of a bond CP and bond path connecting two atoms is a good indicator of interaction [21]. For the self-assembled dimer (Figure 9a), in addition to the symmetrically related H-bonds (characterized by a bond CP (critical point) and bond path interconnecting the H and O-atoms), the QTAIM analysis reveals the existence of a π···π stacking interaction that further stabilizes the formation of the dimer. The dimerization energy is moderately strong (∆E 1 = −9.1 kcal/mol) due to the contribution of both H-bonds and the π-stacking. We have also analyzed the other motif, where C-H···π interactions are established. The interaction energy is very strong (∆E 2 = −20.4 kcal/mol), because in addition to the C-H···π contacts (two bond CPs and bond paths connect two aromatic H-atoms to two carbon atoms of the adjacent ring) an intricate network of H-bonds is established where six C-H···O and one C-H···N contacts are formed, which are highlighted in Figure 9b by yellow circles.  compounds 1 (a) and 2 (b). The energies at selected points of the surfaces are given in kcal/mol.
We have selected two supramolecular assemblies commented above in Figure 2 to analyze the energetic features of the H-bonds and C-H···π interactions in 1. The QTAIM distribution of critical points and bond paths are also given in Figure 9. The existence of a bond CP and bond path connecting two atoms is a good indicator of interaction [21]. For the self-assembled dimer ( Figure  9a), in addition to the symmetrically related H-bonds (characterized by a bond CP (critical point) and bond path interconnecting the H and O-atoms), the QTAIM analysis reveals the existence of a π···π stacking interaction that further stabilizes the formation of the dimer. The dimerization energy is moderately strong (ΔE1 = −9.1 kcal/mol) due to the contribution of both H-bonds and the π-stacking. We have also analyzed the other motif, where C-H···π interactions are established. The interaction energy is very strong (ΔE2 = −20.4 kcal/mol), because in addition to the C-H···π contacts (two bond CPs and bond paths connect two aromatic H-atoms to two carbon atoms of the adjacent ring) an intricate network of H-bonds is established where six C-H···O and one C-H···N contacts are formed, which are highlighted in Figure 9b by yellow circles.   Figure 10a shows the self-assembled dimer of compound 2 where up to eight C-H···O contacts are established between either aromatic or aliphatic H-atoms and the O-atoms of sulfonamide (each one characterized by a bond CP and bond path, see yellow circles in Figure 10a). Moreover, two symmetrically distributed C-H···π interactions are also present and characterized by a bond CP and bond path connecting the aliphatic H-atom to one C-atom of the aromatic ring. As a consequence of this combination of interactions, the dimerization energy is very large ∆E 3 = −16.7 kcal/mol, thus confirming the importance of this motif in the solid state of compound 2. Figure 10b shows a dimer extracted from the infinite 1D chain represented in Figure 3b. In this case, the interaction energy is modest (∆E 4 = −3.4 kcal/mol) because only one H-bond is established. The distribution of bond CPs and bond path also reveals the existence of van der Waals interactions between the alkyl chain and the aromatic ring. Figure 11 shows two dimers of compound 3, one corresponds to the C-H···π assembly commented above in Figure 4, where in addition to the C-H···π interaction (characterized by a bond CP and bond path) the assembly is further characterized by a C-H···O bond involving the methyl group. This assembly presents a modest interaction energy of ∆E 4 = −5.9 kcal/mol. In contrast, the dimer shown in Figure 11b, extracted from the infinite 1D assembly, exhibits a strong interaction energy (∆E 5 = −17.3 kcal/mol) due to the formation of four C-H···O contacts, which are characterized by a bond CP and bond path (see yellow circles in Figure 11b). The strong interaction energy agrees well with the MEP surface analysis commented above, since the H-bond donors belong to the triazinane ring that exhibit the most positive MEP values. Moreover, the H-bond acceptors are the O-atoms of the sulfonamide groups that present the most negative MEP values (see Figure 8).
one characterized by a bond CP and bond path, see yellow circles in Figure 10a). Moreover, two symmetrically distributed C-H···π interactions are also present and characterized by a bond CP and bond path connecting the aliphatic H-atom to one C-atom of the aromatic ring. As a consequence of this combination of interactions, the dimerization energy is very large ΔE3 = −16.7 kcal/mol, thus confirming the importance of this motif in the solid state of compound 2. Figure 10b shows a dimer extracted from the infinite 1D chain represented in Figure 3b. In this case, the interaction energy is modest (ΔE4 = -3.4 kcal/mol) because only one H-bond is established. The distribution of bond CPs and bond path also reveals the existence of van der Waals interactions between the alkyl chain and the aromatic ring. Figure 10. (a,b) Distribution of bond, ring and cage critical points (green, yellow and blue spheres, respectively) and bond paths in two dimers of complex 2. The C-H···O bonds are highlighted. Figure 11 shows two dimers of compound 3, one corresponds to the C-H···π assembly commented above in Figure 4, where in addition to the C-H···π interaction (characterized by a bond CP and bond path) the assembly is further characterized by a C-H···O bond involving the methyl group. This assembly presents a modest interaction energy of ΔE4 = −5.9 kcal/mol. In contrast, the dimer shown in Figure 11b, extracted from the infinite 1D assembly, exhibits a strong interaction energy (ΔE5 = −17.3 kcal/mol) due to the formation of four C-H···O contacts, which are characterized by a bond CP and bond path (see yellow circles in Figure 11b). The strong interaction energy agrees well with the MEP surface analysis commented above, since the H-bond donors belong to the triazinane ring that exhibit the most positive MEP values. Moreover, the H-bond acceptors are the O-atoms of the sulfonamide groups that present the most negative MEP values (see Figure 8). Finally, Figure 12 shows the dimeric motifs analyzed in compound 4. The dimer of Figure 12a presents an intricate combination of C-H···O bonds in addition to two C-H···π interactions involving the t-butyl group. As a consequence of the formation of six concurrent H-bonds, the dimerization energy is very large, ΔE7 = −20.2 kcal/mol. Figure 12b shows the dimer extracted from the 1D supramolecular polymer (see Figure 5b), which presents a modest interaction energy due to the formation of a single H-bond along with van der Waals contacts between the aromatic and aliphatic C-H bonds. Finally, Figure 12 shows the dimeric motifs analyzed in compound 4. The dimer of Figure 12a presents an intricate combination of C-H···O bonds in addition to two C-H···π interactions involving the t-butyl group. As a consequence of the formation of six concurrent H-bonds, the dimerization energy is very large, ∆E 7 = −20.2 kcal/mol. Figure 12b shows the dimer extracted from the 1D supramolecular polymer (see Figure 5b), which presents a modest interaction energy due to the formation of a single H-bond along with van der Waals contacts between the aromatic and aliphatic C-H bonds.
presents an intricate combination of C-H···O bonds in addition to two C-H···π interactions involving the t-butyl group. As a consequence of the formation of six concurrent H-bonds, the dimerization energy is very large, ΔE7 = −20.2 kcal/mol. Figure 12b shows the dimer extracted from the 1D supramolecular polymer (see Figure 5b), which presents a modest interaction energy due to the formation of a single H-bond along with van der Waals contacts between the aromatic and aliphatic C-H bonds.

Concluding Remarks
In summary, the synthesis and single crystal X-ray diffraction characterization of four N-substituted 1,3,5-triazinanes are reported along with a detailed analysis of the noncovalent interactions observed in the solid state. All complexes have in common the formation of several motifs characterized by a network of C-H···O interactions that exhibits very strong binding energies as a consequence of these cooperative H-bonds. Moreover, several structures also form interesting C-H···π···H-C ternary assemblies that have been described in detail. Besides, the MEP surfaces have been used to rationalize the noncovalent interactions and the QTAIM method to confirm the existence of the intricate combinations of H-bonds. Finally, the Hirshfeld surface analysis provides further evidence for the importance of C-H···π and C-H···O in the crystal packing of compounds 1-

Concluding Remarks
In summary, the synthesis and single crystal X-ray diffraction characterization of four N-substituted 1,3,5-triazinanes are reported along with a detailed analysis of the noncovalent interactions observed in the solid state. All complexes have in common the formation of several motifs characterized by a network of C-H···O interactions that exhibits very strong binding energies as a consequence of these cooperative H-bonds. Moreover, several structures also form interesting C-H···π···H-C ternary assemblies that have been described in detail. Besides, the MEP surfaces have been used to rationalize the noncovalent interactions and the QTAIM method to confirm the existence of the intricate combinations of H-bonds. Finally, the Hirshfeld surface analysis provides further evidence for the importance of C-H···π and C-H···O in the crystal packing of compounds 1-4. We assume, that observed case of intramolecular sp 2 -sp 3 disequalization makes the corresponding family of N-alkyl-N ,N"-substituted triazinanes an interesting object for research in the domain of local molecular disorder in organic crystals [25].  Table S1: Data and structure refinement for 1-4. Table S2: Selected distances (Å), angles ( • ) and torsion angles ( • ) for compounds 1-4.