Marcos A. P.
Martins
*a,
Paulo R. S.
Salbego
a,
Guilherme A.
de Moraes
b,
Caroline R.
Bender
a,
Priscilla J.
Zambiazi
b,
Tainára
Orlando
a,
Anderson B.
Pagliari
a,
Clarissa P.
Frizzo
a and
Manfredo
Hörner
*b
aNúcleo de Química de Heterociclos (NUQUIMHE), Department of Chemistry, Federal University of Santa Maria (UFSM), 97105-900, Santa Maria, RS, Brazil. E-mail: marcos.nuquimhe@gmail.com
bNúcleo de Investigação de Triazenos e Complexos (NITRICO), Department of Chemistry, Federal University of Santa Maria (UFSM), 97105-900, Santa Maria, RS, Brazil. E-mail: hoerner.manfredo@gmail.com
First published on 22nd November 2017
The crystallization of a series of 1-(4-halophenyl)-3-phenyltriazenide N1-oxides (1–4) and 1-(phenyl)-3-phenyltriazenide N1-oxide (5) was evaluated using the supramolecular cluster approach. This method is an efficient tool to assess the crystallization mechanism of compounds and, consequently, the steps involved in crystal formation. Compounds 1 and 4 show crystallization in two main steps while compounds 2 and 3 present three main steps, in which column formation occurs in the first step. The crystallization process for 5 occurs in 3 main steps, starting from a robust dimer formation (−16.82 kcal mol−1). Two new parameters – NCG% (topological and energetic contribution percentage) and NG/NC (energetic parameter/topological parameter ratio) – assisted in the interpretation of crystal growth. Compounds 1–4 showed NCG% = 50 in the first step while compound 5 reached only 50% of the contribution in the second step. The differences in NCG% were attributed to strong hydrogen bonds in the non-halogenated compound. The dominant parameter in each step of the crystallization process was indicated by the NG/NC parameter. The crystallization mechanism in all compounds was initially driven by an energetic process followed by a topological process. The existence of X⋯π interactions was shown and was observed to be a consequence of a topological process and without any major contributions to crystal formation. Thermal analysis and UV-vis spectral data were also discussed regarding the properties of these compounds.
According to Kitaigorodskii,3 the assembly route of molecules to crystals is governed by molecules approaching each other to form the first close-pack, which are then close-packed by other layers or chains of molecules. Because of this, the crystallization process can be viewed as a stepwise progression in which molecule association increases system complexity, leading to the final 3D crystal. This progression occurs analogously to Darwinian evolution,4 where the process of self-organization blends competition, cooperation, and complementarity. The Darwinian process improves individual organization, which leads to the emergence of the crystalline arrangement.
Following this idea, our research group proposed a method in a previous study5 where crystallization mechanisms for isoxazoles were presented in order to better understand crystalline formation. This proposal was based on the observation of the energetic contributions of different dimers that formed the crystal and used the supramolecular cluster as part of the study. Our findings show that the supramolecular cluster has an energetic hierarchy related to the formation of interactions; in this case, the hydrogen bond was the strongest interaction and the first to form. The π-interactions, on the other hand, were weaker and were observed in other steps of the crystallization.
Triazene N1-oxide derivatives are an important class of molecules due to the hard basic character of the coordinating atoms O and N. This characteristic imposes electronic properties that allow molecular self-organization via intermolecular secondary interactions.6 These species have different applications in chemistry and biology, such as in sensors, electronic displays,7,8 and antileukemic and anticancer drugs.9,10 Some of these applications are a result of non-covalent intermolecular interactions and the packing behavior of these compounds.
Therefore, our aim is to assess the comprehension of the system by studying the analysis of crystal formation while considering the topological and energetic normalized data of the supramolecular cluster as a tool. In pursuance of this, a series of five N1-oxide triazenes was used in addition to the proposal of crystallization mechanisms for each system in order to observe the contrasts between different substituents. Important information regarding the contribution of each crystallization step will be provided. This will be done by observing the dominant parameter of the process in addition to observing the influence of halogen⋯π interactions on crystal formation. Geometric and QTAIM data were used to answer the questions on the role of halogen⋯π interactions in crystal formation. Additionally, thermal analyses (TGA and DSC) and UV-vis spectroscopy were carried out to determine the potential applications of these molecules.
Fig. 1 Projection of the molecular structures of compounds 1–5. H atoms are represented with arbitrary radii. Non-hydrogen atoms are represented with thermal ellipsoids at the 50% probability level. |
All compounds discussed in this study belong to the class of triazene N1-oxides, such as derivatives from the classical triazenes, in which the double-bonded terminal N atom in the nitrogen triad is substituted by an “oxo” group in the solid state. In solution, these compounds are in chemical equilibrium with the respective tautomer form in which the oxo group is protonated, resulting in a hydroxyl functional group. The diazo(oxide)amine [–NN(→O)–N(H)–] functional group can be stabilized by substitution of the terminal N atoms with alkyl, phenyl, or aryl substituents. In the present study, 1-phenyl-3-para-halophenyl (p-F, p-Cl, p-Br and p-I) substitution was selected in view of the comparison of its crystallographic features with that of the unsubstituted 1,3-bis(phenyl)triazene N1-oxide.
Compound 1 belongs to a monoclinic crystal system in a primitive three-dimensional lattice corresponding to Bravais type P and chiral space group P21 (Z = 2; Z′ = 1). Given that compound 1 does not present a heavy atom with anomalous scattering, we were unable to determine its absolute configuration. Thus, the selected molecular configuration was based on the lowest Flack value observed, which was 0.1(3). An attempt to solve the structure considering the centrosymmetric space group P21/m in comparison with the acentric space group P21 was fruitless. Compounds 2 and 3 are isostructural and belong to the centrosymmetric monoclinic space group P21/c (Z = 4; Z′ = 1). Compound 4 belongs to the orthorhombic and chiral P212121 space group (Z = 4; Z′ = 1), and its molecular configuration was selected based on the lowest Flack value observed, which was 0.048(18). Compound 5 crystallizes in the centrosymmetric monoclinic space group P21/n (Z = 4; Z′ = 1). All molecules present the expected trans stereochemistry regarding the NN double bond in the [–NN(→O)–N(H)–] diazoamine N1-oxide moiety and an almost planar shape. Molecules 1–5 deviate slightly from planarity [r.m.s. deviation from the least-squares plane through the non-hydrogen atoms: 0.0857 Å (1), 0.0353 Å (2), 0.0318 Å (3), 0.0466 Å (4), 0.0468 Å (5)]. The N13–N12–N11–O dihedral angle in the molecules {[−0.3(2)° (1)], [−0.9(4)° (2)], [−1.1(4)° (3)], [1.4(12)° (4)] and [0.5(3)° (5)]} is similar to the angle [N13–N12–N11–O1 = 0.7(5)°] found in 3-phenyl-1-methyltriazene N1-oxide and contributes to the near-planar conformation of the molecules.11 The average value of 1.295 Å for the N11–O bond distance observed in molecules 1–5 (ESI,† Table S1) is in agreement with the N11–O1 bond distance of 1.246(4) Å observed in a similar compound, 3-phenyl-1-methyltriazene N1-oxide, which indicates that the molecules adopt the N-oxide tautomeric form in the solid state.11 The angles between the planes of the rings presented values of 10.32, 3.01, 2.16, 5.93, and 3.53 for 1–5, respectively. These values are a result of the influence of the crystal packing on the different intermolecular interactions of each structure.
The neighboring molecules around M1 represent the molecular coordination number (MCN). These dimers are represented by the interaction between the M1 and MN molecules (e.g., M1⋯M2, M1⋯M3, …, M1⋯M15). All molecules which share the surface contact area of the van der Waals radii with the reference molecule M1 are considered part of the cluster. The supramolecular cluster with the surface contact area of each dim was obtained using the Voronoi–Dirichlet polyhedra (VDP) method with the ToposPro 4.0 software.13 This method has already been reported5,14–16 as an efficient tool to assess crystalline understanding, especially regarding the steps involved in crystal formation, i.e. the crystallization mechanism.5
The different supramolecular clusters can be compared by the normalization12 of raw topological and energetic data in order to observe the distinct behavior of the dimers in the formation of each crystalline structure. Normalization is a tool that places all the different data in the same scale, allowing precise comparisons and interpretations of different systems. Hence, all N molecules in the cluster contribute in some way to the whole system, and it is necessary to assess the contribution of each one in relation to a standard, in our case, in relation to the number of molecules present in the cluster (around reference molecule M1).
For example, in a cluster with an MCN of 16, all dimers must have an ideal contribution of 1 for all considered parameters (e.g. energy), which leads to a 100% contribution of 16. However, this is not the reality for most molecules where distinct interactions occur, leading to different dimeric contributions. Therefore, all raw data should be reduced to the same scale using the MCN as reference for this normalization.
From normalized values, intermolecular interactions can be divided into four types.12 Type I represents a molecular dimer with a high interaction energy on a small contact surface (e.g. hydrogen bonds). Type II corresponds to a high interaction energy on a large contact surface (e.g., π⋯π interactions). Type III corresponds to both small interaction energy and contact surfaces with a maximum difference of ±0.5 between both parameters, and type IV is a contribution based on large contact surfaces with low interaction energy.
Topological and energetic normalized data of each dimer from the supramolecular cluster of compounds 1, 2, and 5 (for 3 and 4, see the ESI,† Fig. S1) with the respective dimer classification are described in Fig. 2.
Fig. 2 Topological (black) and energetic (red) normalized data of each dimer from the supramolecular cluster of compounds 1, 2, and 5. |
Fig. 2 demonstrates that compound 1 has two type I dimers and two type II dimers, which correspond to dimers involved in hydrogen interactions and stacking, respectively. The remaining dimers exhibited a type III behavior with low area and energy values. Notably, data from compound 1 are similar to data from compound 4 (see the ESI,† Fig. S1). Compound 2 (similar for compound 3, see the ESI,† Fig. S1) had only one dimer with type I behavior and one dimer with type II behavior, which is different from what was observed in compounds 1 and 4. These data provide important information on the formation of the crystalline structure, especially concerning the first dimers formed (the most energetic ones). Finally, compound 5 displayed a distinct characteristic when compared to compounds 1–4. The supramolecular cluster presented only one type I dimer of the highest energetic contribution of all systems and two type II dimers with lower contribution regarding stacking interactions when compared to 1–4. This normalized data approach provides significant data of the characteristics of the dimers and enables the understanding of the crystal formation, which is discussed in the following section.
Firstly, individual dimer data regarding stabilization energy were used to define the highest and lowest dimer hierarchies in order to assess the steps of crystal formation. Then, the main steps were defined, where in some cases more than one hierarchy is present in the same crystallization step. The proposed crystallization mechanism for compounds 1, 2, and 5 are presented in Fig. 3, 4 and 5, respectively.
Fig. 3 Proposed crystallization mechanism of compound 1. NG = normalized energy. The shaded area represents the portion in the previous step. |
Fig. 4 Proposed crystallization mechanism of compound 2. NG = normalized energy. The shaded area represents the portion in the previous step. |
Fig. 5 Crystallization mechanism of compound 5. NG = normalized energy. The shaded area represents the portion in the previous step. |
Once the main steps were determined, we returned to the data of the normalized dimer in order to obtain two main parameters to assess the elucidation of the process:
• The topological and energetic contribution for each step of the crystal growth. The percentage of contribution (NCG%) from each given step can be found by the sum of the normalized topological and energetic data (NC + NG) of each step, in which the total NC + NG of the cluster is 100%. As a result, we are able to establish the steps with the highest contribution for both parameters, i.e. surface area and energy. The step contribution can indicate where the definitive step towards the considered crystalline structure is.
• The ratio between NG and NC to find the parameter that rules each step. The NG/NC ratio provides important information on the governing parameter in each step, in which values over 1.0 indicate a higher energetic contribution and values under 1.0 indicate a dominant topologic parameter. Values near 1.0 indicate a balance between parameters. Thus, it is possible to inform where the changes between the dominant parameters occur.
Fig. 3 presents the steps of the crystallization mechanism for compound 1. Compound 4 also presents the same general steps, and an illustration can be found in the ESI† (Fig. S8). The cluster molecules involved in the steps were depicted and their normalized energy values were highlighted.
The proposed crystallization mechanism for structures 1 and 4 shows two major steps (Fig. 3). The first step occurs with the formation of interconnected columns by hydrogen interactions with a stabilization energy of −9.78 kcal mol−1 (NG = 2.99) for each of the two dimers involved and an energy value of −8.85 kcal mol−1 (NG = 2.71) for the molecule stacking of each dimer. The dimers involved in this step present the highest portion of surface contact area and stabilization energy in relation to the whole supramolecular cluster (see Fig. 2). This step contributes 53.4% to crystal formation when both stabilization energy (NG) and surface contact area (NC) are considered, namely the topological and energetic parameters are taken into consideration in the final contribution (NCG%).
Step 2 presents the final crystal formation (3D growth) with the interaction between the columns formed in step 1. These dimer interactions occur with stabilization energy below −2.7 kcal mol−1 (NG < 0.83). All these dimers represent the remaining 46.6% of topological and energetic contribution to form the 3D network.
A value of 1.46 was found for the first step, regarding the second parameter (NG/NC), which indicates a strong impact from the energetic parameter in this step of crystal growth, i.e. almost 50% more influence than the topological factor. Moreover, a value of 0.65 was found for the second step, demonstrating a change in the dominant parameter in this step of the crystallization where topology plays a more relevant role than stabilization energy.
Compounds 2 and 3 present similar clusters with isostructural behavior, which leads to identical proposals of crystallization mechanism. Fig. 4 depicts the mechanism of compound 2 (for compound 3, see the ESI,† Fig. S9).
There are three main steps in the proposed crystallization mechanism for compounds 2 and 3. In the first step, the two main energetic dimers form the displaced columns: the dimer involved in a hydrogen interaction of −15.17 kcal mol−1 (NG = 3.31) and the stacking dimer of −13.65 kcal mol−1 (NG = 2.98). Two dimers from the supramolecular cluster with a stabilization energy of −4.48 kcal mol−1 (NG = 1.96) appear in this step due to the formation of the displaced columns. The dimers considered in the first step had an NCG% of 54.6%, indicating that this first step is a point of no return for the crystal since over 50% of all stabilization energy and surface contact area of the crystal rely on this first step.
The second step is the interaction between the columns (formed in step 1) with an interaction energy of −6.34 kcal mol−1, which leads to the 2D growth of the crystal. In this step, crystallization is increased by a NCG% of 9.6%. Finally, the third step is responsible for the formation of the 3D crystal network. In this final step, the 2D columns interact through dimers with a lower stabilization energy below −3.2 kcal mol−1 (NG < 0.69). This final step represents the remaining 35.8% of the topological and energetic contributions. The NG/NG parameter presented values of 1.17, 1.06, and 0.77 for the first, second, and third steps, respectively. These data demonstrate the evolution of the dominant parameter through the crystallization process with a slight dominance of the stabilization energy in the first step and the subsequent rise of topological governance in the other two steps.
Fig. 5 demonstrates the proposed crystallization mechanism of compound 5, which differs from the other supramolecular structures.
The first step for compound 5 (Fig. 5) differs from those of the previous compounds due to the formation of a dimer with a robust stabilization energy of −16.82 kcal mol−1. The normalized data (NG = 4.49) show the magnitude of this dimer with a value four and a half times greater than expected, which represents a strong hydrogen interaction. This first step corresponds to 19.2% of the contribution of the whole crystal, even though it is only one dimer. The second step is the stacking interaction of these dimers in order to form columns with a stabilization energy of −9.44 kcal mol−1 (NG = 2.52) and −3.89 kcal mol−1 (NG = 1.04). At this moment, there is an increase of 37% in the contribution. The third step is similar to that of compounds 1 and 4, where the columns interact with each other with energies below −2.6 kcal mol−1 (NG < 0.69) in order to form the 3D crystal lattice. The last step contributes with an NCG% of 43.8. The crystallization behavior of this structure differs from that of the others in the contribution of each step, since there is a unique first step with a strong hydrogen interaction between the two molecules. The NG/NC ratio showed values of 2.72, 1.05, and 0.63 for steps 1, 2, and 3, respectively. These data clearly demonstrate the transition between the dominant parameter in each step. The first step with high energy in a smaller surface area had a value almost three times higher, thus highlighting the importance of the energetic parameter. The transition in the dominant parameters in the second step indicated a balance between the surface and energy parameters. The last step, due to the presence of lower stabilization energies in column approximation, was highly influenced by the topological parameter.
For better visualization, the three proposed mechanisms of crystallization for compounds 1–5 are summarized in Fig. 6.
Fig. 6 Summary of crystallization steps for compounds 1–5 with the topologic and energetic contribution (NCG%). The NG/NC parameter is also demonstrated for each stage. |
The data presented in Fig. 6 show the main event in each step for compounds 1–5 until 3D crystal growth. In each step, the combined contribution considering topological and energetic parameters is described in order to have a comprehensive view of the steps in all the studied supramolecular structures. Compounds 1–4 had contributions of over 50% in the first step. Compound 5, which formed a strong first dimer, reached a contribution of 50% only when the first and second stages were considered. However, it is possible to state that the first dimer is formed irreversibly, since the hydrogen interaction formed corresponds to almost 20% of the entire interaction energy of the crystal.
The NG/NC parameter allows us to observe the transition between the dominant parameters during the crystallization process. Compounds 1–5 behaved similarly during the steps, in other words, the first dominance of the stabilization energy of the interactions involved in the considered dimers, its subsequent decrease, and the increase of the topological relevance to the crystal growth. The final step of the crystallization process relies almost exclusively on a topologically dependent process.
Thus, based on such findings, differences between the angles of the rings in each structure mentioned in the molecular description section can be better understood. Compounds 1 and 4 have the highest angles, which is unexpected due to the differences between F and I. However, the influence on these torsion angles relies on crystal packing, and it was possible to observe that 1 and 4 have similar mechanisms of crystallization, achieving a similar 3D crystal network. Structures 2 and 3 are isostructural and consequently present similar torsion values. The structure without a substituent (5) has an intermediate value, which is distinct from the others and consistent with its different crystal packing.
This discussion will focus mainly on halogen interactions and its importance to crystal packing. These intermolecular interactions were identified according to the standardization of the Platon program.17 The arrangements involving halogen⋯π interactions for compounds 1–4 are shown in the ESI† (Fig. S40), where a supposed prominence of these intermolecular interactions is demonstrated. Fig. 7 presents the orientation of the halogen atoms for compounds 1–4 in relation to the plane of the aromatic ring forming the halogen⋯π interactions. The geometric parameters for all halogen⋯π interactions are presented in Table 1.
Fig. 7 Orientation of the halogen atoms in relation to the plane of the aromatic ring forming the halogen⋯π type interactions for compounds 1–4. |
Parameter | 1 | 2 | 3 | 4 |
---|---|---|---|---|
α = X⋯centroid–plane. β = C–X⋯centroid. | ||||
Dist. ra (Å) | 3.2654(5) | 3.5379(4) | 3.5332(1) | 3.6053(1) |
Dist. rb (Å) | 3.4127(5) | 3.4023(4) | 3.4054(1) | 3.9080(2) |
Dist. rr (Å) | 4.1018(7) | 3.7591(5) | 3.7499(1) | 3.9080(2) |
β (°) | 91.443(5) | 154.064(7) | 154.457(2) | 170.601(1) |
α (°) | 44.757(4) | 70.169(8) | 70.344(2) | 49.518(1) |
Symmetry operators | 2 − x, −½ + y, −z | 1 + x, −½ − y, ½ + z | 1 + x, ½ − y, ½ + z | 1 − x, −½ + y, 1½ − z |
The interactions presented in 1 and 4 are classified as localized and type (L).1,2 This is due to their parameter ra (distance from the halogen to a specific carbon) as the shortest distance when compared to the parameters rr (distance from the halogen to the centroid) and rb (distance from the halogen to the bond between two carbons of the aromatic ring). Table 4 shows ra values equal to 3.2654(5) and 3.6053(1) Å for compounds 1 and 4. The distance of the van der Waals radii for F⋯C and I⋯C are 3.17 and 3.68 Å, respectively, indicating a strong halogen⋯π interaction for compound 4, whereas values for distance ra are lower than the sum of the van der Waals I⋯C radii.
Compounds 2 and 3 are classified as semi-localized and type (SL). Their distances of lower values registered for parameter rb are 3.4023(4) and 3.4054(1) Å, respectively, where the distances of the interactions Cl⋯π and Br⋯π are 3.4023(4) and 3.4054(1) Å, respectively. In this sense, values below the sum of the van der Waals radii, which are 3.45 and 3.55 Å, respectively, for Cl⋯C and Br⋯C, indicate a significant interaction. The halogen⋯π interaction nomenclature for compounds 1–4 can also be attributed to the geometric criteria. For compound 1, the F⋯π interaction shows angle β(C–F⋯centroid) as 91.443(5)°, a value very close to the orthogonal one, thus indicating that it is a lone pair⋯π type interaction according to Mooibroek et al.18 For compounds 2–4, the halogen⋯π interactions have a β angle of 154.064(7), 154.457(2), and 170.601(1), respectively, and present values close to linearity (180°). This indicates that these are σ-hole type interactions.
In this geometric approach, it is possible to assess some magnitude of the X⋯π interactions using only the atom–atom distances, in accordance with the established literature. However, how can the existence of the X⋯π interactions be ensured only by means of raw geometric data? And what about the strength of those interactions, how can it be assessed without considering an energetic approach? How can it be properly assessed that these types of interactions are relevant to crystal formation?
Fig. 8 Representation of X⋯π (1–4) and C–H⋯π (1 and 5) interactions through the BCPs observed by QTAIM analysis. The energy involved in each interaction is highlighted (in kcal mol−1). |
The first question regarding the existence of the intermolecular interactions is then answered: each interaction involved in the considered dimer has been identified by the BCPs in the QTAIM. The second and third questions are still pendent and, thus, it is necessary to go beyond the energy calculation of the dimer and the observed interaction. A referential is necessary to assess the relevance of this dimer/interaction, and the supramolecular cluster can aid us in acquiring this information. Therefore, the energy involved in each interaction was obtained through the fragmentation of the energy from the considered dimer using a previously reported method,5,14 which is highlighted in Fig. 8.
The QTAIM data regarding the interactions involved in these dimers are presented in Table 2. These interactions observed by the BCPs, in a separate form, can induce the observer to conclude that X⋯π interactions are relevant to crystal formation. However, when considering the energy of the involved dimers that range from −1.76 to −2.69 kcal mol−1, it becomes clear that these dimers have a low energetic value when considering total crystal energy. Consequently, this information helps answer the two remaining questions.
Comp. | Dimer | G M1⋯MN (kcal mol−1) | Interaction | ρ INT (a. u.) | G AI (x⋯Y) (kcal mol−1) | %d |
---|---|---|---|---|---|---|
a Stabilization energy obtained by DFT calculation. b Obtained using QTAIM. c Interaction energy involved in each interaction obtained by the fragmentation. d Interaction contribution regarding the total energy cluster. | ||||||
1 | M1⋯M12 | −2.69 | F⋯π | 0.004047 | −1.10 | 3.7 |
M1⋯M13 | H⋯π | 0.005811 | −1.59 | 5.4 | ||
2 | M1⋯M9 | −1.76 | Cl⋯π | 0.006136 | −1.76 | 5.5 |
M1⋯M14 | ||||||
3 | M1⋯M9 | −2.61 | Br⋯π | 0.007093 | −2.61 | 7.9 |
M1⋯M14 | ||||||
4 | M1⋯M16 | −2.56 | I⋯π | 0.005375 | −2.56 | 7.0 |
M1⋯M19 | ||||||
5 | M1⋯M9 | −2.33 | H⋯π | 0.002655 | −0.83 | 2.8 |
M1⋯M16 | H⋯π | 0.004812 | −1.50 | 5.0 |
The fragmentation of the interaction energy of the dimer in the different interactions provides important additional data. The contribution of X⋯π interactions is discernible when the whole system is considered. Moreover, X⋯π interactions are responsible for less than 7.9% of the entire stabilization energy of the supramolecular cluster, contributing 3.7%, 5.5%, 7.9%, and 7.0% for compounds 1 to 4, respectively. These data provide an important insight into X⋯π interactions and allow us to answer the previous question of the importance of this kind of interaction for crystal formation. This broad analysis has shown us that in these series of crystal molecules, halogen⋯π interactions are a topological consequence of the crystallization process because they are present in the last steps of the crystallization process in the approximation of the columns. The complex approach, which is not solely dependent on geometric parameters, enabled us to confidently assess that halogen⋯π interactions do not drive the crystallization process.
Using the molecular electrostatic potential (MEP), it is possible to observe the regions with negative and positive potential in the molecules.20 In the case of compounds 2–4, it is possible to observe the positive region in the halogen atom (σ-hole) (Fig. 9), which allows the σ-hole⋯π interaction observed in the previous section to occur. Compound 1 with the fluorine atom, on the other hand, does not present the σ-hole, which leads to a lone-pair⋯π interaction.20 Additional MEPs regarding the involved dimers are presented in the ESI† (Fig. S10).
DSC and TG analyses were carried out to elucidate the thermal performance of these structures in a broad range of temperatures. Firstly, TGA screening was carried out for compounds 1–4 at a temperature range of 25 to 900 °C in order to observe the decomposition steps and decomposition temperature of the structures as a function of the halogen p-substituted in the aromatic ring (see the ESI;† Fig. S11–S15). The decomposition percentage indicates that the decomposition of the samples is almost complete (>90%) in the evaluated temperature range. Nevertheless, Fig. 10 shows the TGA thermograms of compounds 1–4 in a range from 100 to 400 °C, in which the main changes in weight loss versus temperature were observed. For a better visualization, the TGA thermogram of 5 was omitted from Fig. 10 and depicted in the ESI.†
From Fig. 10, it is possible to note that these structures decompose in two main steps, the first from 90 to 220 °C and the second from 150 to 450 °C, approximately. Furthermore, by peak differences in the derivative weight as a function of temperature plots, it can be noted that the halogens influence the decomposition pattern of triazene N-oxides. Table 3 shows data regarding the initial decomposition temperature (Ti), final decomposition temperature (Tf), decomposition percentage (%), and temperature of maximum decomposition (Td). Ti is the temperature at which the compound starts to lose mass and indicates that no thermal degradation based on mass loss occurs at temperatures below Ti. Tf indicates the temperature in which mass loss is complete. Td corresponds to the temperature of the maximum peaks in the derivative weight (%) as a function of temperature (T) in the TGA thermogram (Fig. 10).
Compound | Step | T i (°C) | T f (°C) | T d (°C) | %d |
---|---|---|---|---|---|
a Initial decomposition temperature. b Final decomposition temperature. c Temperature of maximum decomposition. d Decomposition percentage. | |||||
1 | 1 | 129.36 | 155.16 | 143.97 | 52.50 |
2 | 155.16 | 396.20 | 177.54 | 41.47 | |
2 | 1 | 137.24 | 173.59 | 159.96 | 67.77 |
2 | 173.59 | 400.05 | 241.27 | 23.68 | |
3 | 1 | 94.92 | 207.79 | 160.92 | 63.34 |
2 | 207.79 | 449.78 | 254.66 | 27.58 | |
4 | 1 | 121.70 | 217.35 | 162.83 | 66.61 |
2 | 217.35 | 455.52 | 282.40 | 20.45 | |
5 | 1 | 128.40 | 154.26 | 138.76 | 60.49 |
2 | 154.26 | 312.88 | 166.60 | 35.84 |
Halogens, in general, tend to increase the stability of triazene compounds as observed in Table 3. An order for decomposition temperatures can be estimated for the compounds in this study, in which 5 < 1 < 2 < 3 < 4. This order in relation to compounds in which X = halogen corresponds to F < Cl < Br < I.
In the literature, halogen compounds exhibit a pattern of decomposition behavior generally related to the strengths of the carbon–halogen (C–X) bonds. The greater the C–X bond dissociation enthalpies, the more thermally stable the halogenated compounds.23 For alkyl halides, for example, the bond dissociation enthalpies are 106, 81, 69, and 54 kcal mol−1 for C–F, C–Cl, C–Br and C–I, respectively, in which the stability order is I < Br < Cl < F.24 This order represents the opposite trend found in this study. We suppose that the order found for triazene N-oxides is a result of two mutual factors: (i) halogen electronegativity and (ii) atomic radius. An increase in decomposition temperature of the compounds as the halogen atomic radius increases and electronegativity decreases can be highlighted. Unexpectedly, this indicates that triazene N-oxides substituted in the para position by hard substituents (lower atomic radius and greater electronegativity, i.e. F atom) decomposed at lower temperatures than the same compounds with soft substituents in the same position (greater atomic radius and lower electronegativity, i.e. I atom) that possess a more polarizable electronic cloud.
Compounds were submitted to DSC analysis in order to determine their thermal events (e.g. melting and crystallization temperatures). Each sample was subjected to three cycles of heating and cooling at a safe temperature range varied according to the decomposition temperature of the compounds present in the sample. This was done to avoid decomposition within the DSC equipment. However, no thermal events were observed for 1–4 in all cycles within the temperature ranges evaluated (see the ESI;† Fig. S16–S19). This leads us to believe that the melting point of the compounds is likely very close to the point at which they decompose. This assumption is in accordance with the study reported by Serebryanskaya et al.25 The authors observed thermal decompositions for 1,3-bis(2-alkyltetrazol-5-yl)triazenes and their Fe(II), Co(II), and Ni(II) complexes using DSC in addition to finding that the triazene complexes decomposed without melting at temperatures above 240 °C.
Fig. 11 (a) Absorbance versus wavelength (λ) and (b) absorbance versus concentration for 1 at 200 nm in ethanol at 25 °C. |
By using the Lambert–Beer law, the angular coefficient of the straight line observed in Fig. 11(b) corresponds to the molar absorptivity (ε) of compound 1 in ethanol at 200 nm. The absorbance versus λ for 2–4 and the absorbance versus concentration plots for 1 at 234 and 348 nm and 2–4 in all λ of absorption are depicted in ESI,† Fig. S28–S31. The straight lines of all experimental data were obtained with an excellent determination coefficient. Molar absorptivity (ε) measures the ability of a mole of substance to attenuate light incident on a given wavelength. It is an intrinsic property of substances and provides an estimation of how strongly a substance absorbs radiation of a given frequency. Table 4 shows the λ1, λ2, and λ3 of absorption peaks and the molar absorptivity (ε) in each λ for compounds 1–4.
Compound | λ 1 (nm) | λ 2 (nm) | λ 3 (nm) | ε 1 (L mol−1 cm−1) | ε 2 (L mol−1 cm−1) | ε 3 (L mol−1 cm−1) |
---|---|---|---|---|---|---|
a ε 1, ε2 and ε3 correspond to molar absorptivity in λ1, λ2, λ3, respectively. | ||||||
1 | 348 | 234 | 200 | 23098 | 15289 | 26757 |
2 | 360 | 240 | 200 | 12562 | 6424 | 11531 |
3 | 352 | 238 | 200 | 21238 | 12257 | 23963 |
4 | 360 | 241 | 203 | 25083 | 12874 | 21899 |
From Table 4, λ1, λ2, and λ3 were assigned to n → π*, n → σ* and π → π* electron transitions, respectively. The triazene chemical groups probably involved in these transitions are NN, –NN–N, and CC, respectively. The absorption transition of these triazene N1-oxides are in agreement with data previously reported in the literature for triazene structures without coordinated oxygen.26,27 Furthermore, the ε of 1–4 demonstrates that all compounds possess a strong absorbance capacity, which shows promising optical properties.
Domingues et al.26 reported an ortho-substituted triazene with a fluorine atom possessing a greater ε value than that of the ortho-substituted analogue compound with a chlorine atom. In this study, compound 1 with an aromatic ring p-F substituted shows significantly larger ε1, ε2, and ε3 values than the p-Cl substituted analogue compound (2). Nevertheless, the ε1, ε2, and ε3 values for compounds 2 and 3 were similar to values found for compound 1. This indicates that the electron transitions associated with the absorbance phenomena in triazene N1-oxides is not highly dependent on electronic effects caused by different halogen bonds in the aromatic ring.
Compound | 1 | 2 | 3 | 4 | 5 |
---|---|---|---|---|---|
Moiety formula | C12H10FN3O | C12H10ClN3O | C12H10BrN3O | C12H10IN3O | C12H11N3O |
Sum formula | C12H10FN3O | C12H10ClN3O | C12H10BrN3O | C12H10IN3O | C12H11N3O |
Mw | 231.23 | 247.68 | 292.14 | 339.13 | 213.24 |
Crystal system | Monoclinic | Monoclinic | Monoclinic | Orthorhombic | Monoclinic |
Space group | P21 | P21/c | P21/c | P212121 | P21/n |
a [Å] | 9.6487(15) | 6.1680(9) | 6.2398(3) | 4.7465(3) | 5.4917(5) |
b [Å] | 5.2866(7) | 9.5740(16) | 9.7185(5) | 9.9491(6) | 16.1419(13) |
c [Å] | 11.1382(19) | 19.477(3) | 19.5527(10) | 26.6177(14) | 12.4976(10) |
α [°] | 90 | 90 | 90 | 90 | 90 |
β [°] | 110.361(5) | 96.539(9) | 96.462(2) | 90 | 100.382(7) |
γ [°] | 90 | 90 | 90 | 90 | 90 |
V [Å3] | 532.65(14) | 1142.7(3) | 1178.17(10) | 1256.98(13) | 1090.20(16) |
Z | 2 | 4 | 4 | 4 | 4 |
D calcd.[g cm−3] | 1.442 | 1.440 | 1.647 | 1.792 | 1.299 |
μ [mm−1] | 0.108 (Mo-Kα) | 0.320 (Mo-Kα) | 3.475 (Mo-Kα) | 2.535 (Mo-Kα) | 0.087 (Mo-Kα) |
F(000) | 240 | 512 | 584 | 656 | 448 |
T [K] | 293(2) | 293(2) | 293(2) | 293(2) | 295(2) |
Reflections collected | 19621 | 9904 | 11495 | 14579 | 11977 |
Independent reflections | 2392 | 2379 | 2602 | 2719 | 3451 |
R int | 0.0533 | 0.0638 | 0.0475 | 0.0470 | 0.0587 |
Observed reflections [ I > 4σ(F)] | 2125 | 1489 | 1813 | 1758 | 1177 |
Parameters refined | 154 | 154 | 155 | 154 | 145 |
Flack x parameter | 0.1(3) | — | — | 0.048(18) | — |
R(|F|) reflections [ I > 4σ(F)] | 0.0351 | 0.0660 | 0.0337 | 0.0466 | 0.0687 |
wR(F2) | 0.0870 | 0.1660 | 0.0688 | 0.1145 | 0.1464 |
S | 1.058 | 1.015 | 1.007 | 1.037 | 0.998 |
Δρmin [e Å−3] | +0.235 | +0.506 | +0.228 | +0.373 | +0.164 |
Δρmax [e Å−3] | −0.248 | −0.664 | −0.220 | −0.1.031 | −0.148 |
Footnote |
† Electronic supplementary information (ESI) available: 759893, 1570873 and 1570875–1570877. For ESI and crystallographic data in CIF or other electronic format see DOI: 10.1039/c7ce02015e |
This journal is © The Royal Society of Chemistry 2018 |