Synthetic, structural, and computational investigations of N -alkyl benzo-2,1,3-selenadiazolium iodides and their supramolecular aggregates †

Despite their versatility, the application of telluradiazoles as supramolecular building blocks is considerably constrained by their sensitivity to moisture. Albeit more robust, their selenium analogues form weaker supramolecular interactions. These, however, are enhanced when one nitrogen atom is bonded to an alkyl group. Here we investigate general methods for the synthesis of such derivatives. Methyl, iso-propyl and tert -butyl benzo-2,1,3-selenadiazolium cations were prepared by direct alkylation or cyclo-conden-sation of the alkyl-phenylenediamine with selenous acid. While the former reaction only proceeds with the primary and tertiary alkyl iodides, the latter is very e ﬃ cient. Di ﬃ culties reported in earlier literature are attributable to the formation of adducts of benzoselenadiazole with its alkylated cations and side reactions initiated by aerobic oxidation of iodide. However, the cations themselves are resilient to oxidation and stable in acidic to neutral aqueous medium. X-ray crystallography was used in the identi ﬁ cation and characterization of the following compounds: [C 6 H 4 N 2 (R)Se] + X − , (R = CH(CH 3 ) 2 , C(CH 3 ) 3 ; X = I − , I 3 − ], [C 6 H 4 N 2 (CH 3 )Se] + I − , and [C 6 H 4 N 2 Se][C 6 H 4 N 2 (CH 3 )Se] 2 I 2 . Formation of Se ⋯ N secondary bonding interactions (chalcogen bonds) was only observed in the last structure as anion binding to selenium is a strong competitor. The relative strengths of those forces and the structural preferences they enforce were assessed with DFT-D3 calculations supplemented by AIM analysis of the electron density.


Introduction
While hydrogen bonding and the coordination of metal ions are preponderant in supramolecular chemistry, other interactions are receiving increasing interest in this field.A notable example is halogen bonding, 1-3 a special case of the phenomenon termed secondary bonding 4 that is recurrent in the structures of compounds of the heavy p-block elements.The potential of secondary bonding in supramolecular chemistry is well exemplified by the derivatives of the 1,2,5-telluradiazole ring (1, Chart 1).Their molecules usually associate by two antiparallel Te⋯N secondary bonding interactions forming the [Te-N] 2 supramolecular synthon, 5 a virtual four-membered ring (Scheme 1, E = Te). 6In the absence of substantial steric hindrance, each molecule can simultaneously form the [Te-N] 2 supramolecular synthon twice building in this way supramolecular ribbons. 7The driving force for association of these molecules is strong enough to overcome moderate steric repulsion by allowing structural distortion of the ribbons.This effect induces properties such as chromotropism 8 and second-harmonic generation 9 which are of interest for practical applications.2][13] In spite of their versatility, prospects of widespread application of telluradiazoles are limited by their sensitivity to moisture.
The lighter members of the chalcogenadiazole family (2, 3, Chart 1) are much more robust but SBIs centred on S and Se atoms are weaker than those formed by Te 14 and instances of supramolecular association are less frequent in such cases.However, crystallographic evidence suggests that attachment of Lewis acid (a transition metal ion 15,16 or a borane 17 ), on one nitrogen of the heterocycle strengthens the secondary interactions made by Se and Te, shortening the corresponding distances.Such effect is also strong in the cations formed by N-alkylation and provides an explanation for the peculiar trend in reduction potentials of the N-methyl-benzo-2,1,3-chalcogenadiazolium cations in acetonitrile. 18Such evidence of molecular association in solution through SBIs is consistent with the observation of [M 2 H] + supramolecular aggregates in the mass spectrum of benzo-2,1,3-telluradiazoles 19 as well as calculated interaction energies and AIM parameters. 12,18-Alkyl substituted selenadiazolium cationic rings could therefore be a convenient alternative to telluradiazoles in supramolecular architecture, provided that they are reasonably stable and the preparation methods are flexible enough to accommodate multiple variations in the structure.The most straightforward synthetic method would be the direct alkylation of the heterocycle in 4 (Chart 2).However, early reports of this approach describe very inefficient reactions plagued by multiple by-products with most common alkylating reagents. 20,21Clean and efficient alkylations have been reported rather recently using very reactive agents such as Me-OS(O) 2 -R′ (R′ = OMe, C 6 H 2 (NO 2 ) 2 -2,4), 20  ; 22 but these methods are applicable with only a handful of alkyl groups.Alternatively, the alkyl group could be attached to nitrogen before formation of the heterocycle.This approach was actually reported in an early publication 20,21 but its products were not structurally authenticated.A third route to N-alkyl selenadiazolium cations, also limited in flexibility, consists of the spontaneous de-alkylation of one nitrogen atom during the reaction of N,N′-bis(tert-butyl)diazabutadiene with SeCl 4 . 23 we are interested in using the N-alkylated selenadiazolium cations as supramolecular building blocks, in this report we compare the two most flexible synthetic methods for the preparation of the derivatives with primary, secondary and tertiary alkyl groups (5, Chart 2).We examine the crystal structures of their iodide salts, the affinity of the cations for halide anions, and spectroscopic properties useful for the characterization of these species in solution.

Materials and general procedures
Selenium dioxide, iodomethane were used as received from commercial suppliers.5][26] Benzo-2,1,3-selenadiazole was prepared by the reaction of 1,2-phenylenediamine with selenium dioxide according to literature. 27Anhydrous-grade solvents were used without further purification.
Instrumentation 1 H, 13 C and 77 Se NMR spectra were recorded on Bruker 200 and Bruker 600 spectrometers at ambient temperature.Chemical shifts are reported in δ values ( parts per million) with respect to the resonances of tetramethylsilane for 1 H and 13 C and with respect to the line of Me 2 Se for the chemical shift of 77 Se.In the last case, spectra were obtained using 1 H-77 Se HMBC unless otherwise stated.IR spectra were recorded on polyethylene pellets on a Nicolet 6700 FT-IR spectrometer with a resolution of 4 cm −1 .Raman spectra were obtained with a Renishaw inVia spectrometer exciting at 785 nm, 30 mW, averaging 10 × 10 s scans.Low-and highresolution electrospray spectral analyses were performed using Micromass GCT spectrometer and a Micromass Quattro Ultima for electron ionization spectra.In each case, the sample was introduced into the ionization chamber in a shortened borosilicate glass capillary on a probe rod.The samples were then heated and the temperature range that displayed the most intense parent ion peaks was used for spectral acquisition.A 70 eV electron stream ionized the sample and the positively charged ions were identified by a time-of-flight detector.A Cary 50 UV-vis spectrometer and its Cary Win-UV software package were used to acquire the spectra.Additional processing was carried out using the GramsAI suite (version 8.0).Melting points were measured on a Thomas-Hoover melting point apparatus and are reported uncorrected.Elemental analysis was performed by the Science Centre of London Metropolitan University, 29 Hornsey Road, London, UK, N7 7DD.
All single crystal X-ray diffraction samples were mounted with paratone oil at the tip of a glass fibre installed on a MiTeGen goniometer head, and kept under a cold stream of nitrogen while on the diffractometer.Data were collected at 100 K on a Bruker APEX2 diffractometer, using Mo Kα radi-ation (λ = 0.71073 Å) and outfitted with an Oxford cryostream low-temperature accessory.Ω and ϕ scans were collected in 0.5°steps with a crystal to detector distance of 4.954 cm.The preliminary unit cell parameters were determined using a minimum of 50 frames from three different orientations, and final cell refinement after integration in SAINT.Data were corrected for absorption and scaled using face-indices as well as redundant data in SADABS.Crystal structures were solved using SHELXL and most structures were refined by full-matrix least square of all F 2 values with the WinGX package.I] 2 (0.02 g, 0.109 mmol) was dissolved in methanol (1.6 mL) with stirring at 42-46 °C.The colour of the solution slowly changed from dark green to bright red.Cooling to room temperature produced red crystals which were recovered by filtration.Yield 17%.Method B. To a solution of N-methyl-phenylenediamine (0.11 mL, 1.023 mmol) and H 2 SeO 3 (0.1319 g, 1.023 mmol) in anhydrous ethanol, trifluoroacetic acid (0.07 mL, 1.023 mmol) was added dropwise under nitrogen.The reaction mixture was allowed to stir at room temperature for 20 min, time after which a solution of NaI (0.1533 g, 1.023 mmol) in ethanol (10 mL) was added with vigorous stirring, a dark red solid precipitate formed and was separated by filtration, washed with hexanes, dried under vacuum and stored under nitrogen.Yield 22%. 1 H NMR (D 2 O, 600 MHz): δ 4.60 (s, 3H, CH 3 ), 7.78, 7.96 (m, 4H, Ar-CH) ppm. 13 [

Syntheses
To a solution of N-isopropyl-benzene-1,2-diamine and H 2 SeO 3 in anhydrous ethanol, trifluoroacetic acid was added.The reaction mixture was allowed to stir at room temperature for 20 min under nitrogen.Upon the addition of NaI, the product immediately precipitated.The dark red solid was washed with hexane, dried and stored under a nitrogen atmosphere.Yield 30%. 1

Computational methods
All the structures considered in this study were fully optimized using the ADF DFT package (versions 2013.01-2014.01). 28,29he Adiabatic Local Density Approximation (ALDA) was used for the exchange-correlation kernel 30,31 and the differentiated static LDA expression was used with the Vosko-Wilk-Nusair parameterization. 32,33Calculation of model geometries was gradient-corrected with the exchange and correlation functionals of Perdew, Burke and Ernzerhof (PBE) 34 and applying the zeroth-order regular approximation (ZORA) [34][35][36][37][38] formalism with the specially adapted triple-ζ all-electron plus one-polarization-function basis sets.The contribution of dispersion was modelled with Grimme's correction. 39Analytical frequency calculations were performed to ensure that each geometry was at an energy minimum. 40,41The TD-DFT calculation was performed from the optimized geometry of [C 6 H 4 N(NCH 3 )Se] + using the statistical average of potentials model for exchange and correlation (SAOP) [42][43][44] and a basis set of quadruple-zeta quality plus polarization.

Syntheses by direct alkylation
Although benzo-2,1,3-selenadiazole is a weak base, it is readily protonated by strong acids. 45For this reason, the purported inefficiency 20,21 of its reaction with primary alkyl-iodides is unusual and merited reinvestigation using a more reactive tertiary iodide.Initial attempts at bench-top alkylation with (CH 3 ) 3 CI in a hot toluene solution yielded a complex mixture, as made evident by the 1 H NMR spectrum.Slow evaporation of the solution yielded a mixture of crystals with multiple morphologies; three of which were identified by single-crystal X-ray diffraction: The first species does contain the N-tert-butyl benzo-selenadiazolium cation (5c), which demonstrates that the alkylation reaction does proceed.However, the presence of the triiodide anion, plus protonated phenylenediamine and benzoselenadiazole point to the actual problem with this experiment: the iodide ion is likely oxidized by atmospheric oxygen to Interestingly, iso-propyl iodide was completely unreactive towards benzoselenadiazole under the same conditions.While the reaction of tert-butyl iodide is likely to proceed by a SN 1 mechanism, it is likely that the primary and secondary iodides prefer a SN 2 mechanism but steric hindrance prevents the reaction of the iso-propyl halide.

Condensation of selenous acid with N-alkylated phenylenediamines
Even in aqueous medium, benzoselenadiazole is readily formed by the reaction of H 2 SeO 3 and phenylenediamine.
N-substituted ortho-diamino benzenes were reported to undergo an analogous reaction in glacial acetic acid producing the benzoselenadiazolium cations, which are isolated by precipitation with sodium halide. 46,47This procedure works as expected but we found more convenient to carry out the reaction in a mixture of ethanol and trifluoroacetic acid in order to avoid using the less volatile acetic acid.Structural characterization by X-ray diffraction did confirm the identity of the products.We also found that exposure of the cations in solution to the atmosphere resulted in iodide oxidation and crystallization of the corresponding triiodide salts, which was verified by the strong scattering band at 111 cm −1 in the Raman spectrum 48 or crystallographic analysis.

X-ray crystallography
Crystallographic and refinement data for the salts of the three cations are presented in Table 1; selected distances and angles in each crystal structure are provided in Table 2.
The crystal structure of 4[5a] 2 [I] 2 consists of a molecule of 4 bound to two 5a, conforming a pseudo-trimer assembled by two asymmetric [Se-N] 2 supramolecular synthons (Fig. 1) with asymmetric SBI distances of 2.573(4) and 2.937(1) Å.This aggregate is analogous to the product formed by the reaction of benzoselenadiazole with [(CH 3 ) 3 O][BF 4 ], 22 in which the Se⋯N and Se⋯F secondary bonding distances are 2.573(2) and 2.966(6) Å and 2.970(4) Å respectively.In this arrangement, the molecules of 4 and 5a are nearly coplanar with a deviation angle of 5.06°.This structure also features Se-I short distances at 3.528(1) Å and 3.831(1) Å.Both iodides sit out of the average plane of the aggregate forming a virtual four-membered ring (i.e. the [Se I] 2 supramolecular synthon) with the Se atoms of the 5a molecules of two neighbouring pseudotrimers.In clear contrast with the structures of the pseudotrimers, 22 and the dimers observed for [5a][CF 3 SO 3 ] 18 and [C 2 H 2 N(NC(CH 3 ) 3 )Se][GaCl 4 ], 49 the crystal structures of the iodide and triiodide salts of the alkylated cations 5a-c feature no association of the heterocycles to each other.In all cases, strong Se-anion interactions prevent formation of the Se⋯N SBIs (Fig. 2 and 3).

Spectroscopy in solution
The iodide salts of 5a-c are stable in air.The compounds are slightly soluble in water at room temperature, their solubility decreases with the size of the R group.In fact, it is possible to recrystallize 5a from hot water.Aqueous solutions of the cations have a characteristic red-brown colour, their UV-vis absorption spectrum displays prominent maxima in the range of 330-343 nm and a weaker band between 440-460 nm.These spectra are in excellent agreement with the TD-DFT calculation for the 5a, which attributes the bands to the LUMO ← HOMO and LUMO ← HOMO−1 excitations (Fig. 4).The UV-vis spectra do not change from neutral to acidic pH, but hydrolysis ensues in basic medium.Given the modest solubility of the 5a-c iodides, their 77 Se NMR chemical shifts are best obtained through 1 H-77 Se heteronuclear multiple bond correlation (HMBC) experiments.The 77 Se δ values (1455-1490 ppm) appear at slightly lower frequency than the 1515 ppm reported for [C 6 H 4 N(NCH 3 )Se][CF 3 SO 3 ] in acetonitrile. 18For the pseudotrimer 4[5a] 2 [I] 2 , two resonances were observed with frequencies that correspond to those of the constituting neutral and cationic molecules, suggesting that the aggregate dissociates in solution.

DFT calculations
As the crystal structures show there is competition between the formation of the [Se-N] 2 supramolecular synthon and binding of the anion to selenium, the relative strengths of these interactions were evaluated using DFT and were extended to include the chloride and bromide salts.For computational expediency, GGA calculations were performed with the PBE functional, supplemented by dispersion and relativistic corrections.Model structures were optimized from the crystallographic coordinates and expanded to include hypothetical isomers and the analogues with chloride and bromide anions.Whenever an experimental structure was available, bond lengths were reproduced within 0.01 Å, the largest deviations correspond to the intramolecular Se-N distances.Seleniumcentred secondary bond distances were reproduced within 0.01 Å to nitrogen and 0.50 Å to iodide.The large deviation in the latter case is likely due to the additional interactions of the anion in the lattices.Interaction energies were evaluated using the transition-state method 50 that partitions the energy of interaction between two molecules or fragments in a hypothetical process in which the constituting units are first calculated individually, brought to their equilibrium positions without orbital mixing, and finally the electron density is relaxed by the interaction of fragment orbitals.The change from the first to the second step is regarded as the total steric interaction that results from the sum of the electrostatic (E Elstat ) and Paulirepulsion (E Pauli ).The third step gives the orbital interaction energy (E Orb ).Dispersion is treated as a separate contribution (E Disp ).The total interaction energy is given by the sum of all contributions (eqn (1)).
Three partitioning schemes (Fig. 5, i-iii) were used to evaluate the secondary interactions within the pseudo-trimer 4 [5a] 2 [I] 2 the results are summarized in Table 3.In this analysis the average interaction energy of the [Se-N] 2 supramolecular synthon was 56.4 kJ mol −1 .Formation of the interactions in iii and ii could be regarded as steps in the assembly of the pseudotrimer and the calculations show that the interaction energies are nearly additive.
In the case of the [5a][X] (X = Cl, Br, I) ion pairs, two geometries were considered, cis and trans to the alkylated nitrogen atom (Fig. 5, iv and v).The results, summarized in Table 4, indicate that in all three cases the trans structure is preferred.This is consistent with the results of previous calculations of the electrostatic potential maps of the borane adducts of chalcogenadiazoles which showed that the most prominent σ hole on the chalcogen is opposite to the substituted nitrogen. 12he dimerization energies the [5a][X] ion pairs were evaluated in two possible geometries (Fig. 5, vi and vii), the results The pseudo-trimer structures 4[5a] 2 [X] 2 (X = I, BF 4 ) are especially interesting because of their two distinct Se⋯N SBI distances.In order to assess the relative strength of each SBI, a bond-energy decomposition analysis was performed using the extended-transition-state method from the natural orbitals for chemical valence (ETS-NOCV).51 In this method the orbital component to the total bonding energy between fragments is calculated as in eqn (2).
The terms −F TS −k,−k and F TS k,k are diagonal transition-state Kohn-Sham matrix elements corresponding to NOCV's with eigenvalues −υ k and υ k respectively.This particular orbital interaction term provides both a qualitative and quantitative picture of the electronic density reorganization in the [Se-N] 2 supramolecular synthons.Each of the fragment interaction calculations (i-iii) shows that the shorter Se⋯N SBI, that opposite to the quaternary nitrogen, is stronger (−13.0 to −16.9 kcal mol −1 ) than the longer (−2.4 to −7.5 kcal mol −1 ).
The Se⋯N interactions in 4[5a] 2 [I] 2 were also examined through Bader's theory of atoms in molecules (AIM). 52revious studies 12 applied AIM to the examination the Lewis acidity and basicity of 1-3; their results highlighted the ambiphilic character of the chalcogenadiazoles molecules and showed that the attachment of a Lewis acid to one nitrogen of the heterocycle strengthens the [E-N] 2 supramolecular synthon formed by the other nitrogen.The calculated electron density (ρ(r)) and its Laplacian (∇ 2 ρ(r)) at the bond critical points (BCPs) of 4[5a] 2 [I] 2 are given in Table 6.Although the difference of Se⋯N distances is small, the AIM parameters suggest that they are very different in strength as there is indication that the energy of interaction scales linearly with the BCP density.
For comparison, the corresponding BCP parameters of the Se⋯I interaction in the pseudo-trimer are included in Table 6.It is interesting that while the BCP electron density is smaller, the its Laplacian is intermediate between the values for the Se⋯N SBIs.The comparison was further extended by calculating the BCP parameters of the dimers of [5a][X] (X = Cl, Br, I) in the two possible conformations considered under partition schemes vi and vii.The results are compiled in Table 7.All these values suggest that the Se⋯X and Se⋯N interactions are of similar strength.In some cases, the calculated BCP para-    meters approach those of the Te⋯N SBI in the neutral dimer of 1,2,5-telluradiazole (4.20 e Å −3 × 10 −2 and 10.9 e Å −5 × 10 −2 ). 12

Conclusions
As these investigations show, the most convenient method for the preparation of N-alkyl benzoselenadiazolium cations is the condensation of selenous acid with N-substituted phenylenediamines.Such intermediates can be conveniently obtained from o-fluoronitrobenzene and primary amines.As a wide variety of aromatic diamines are commercially available or readily prepared, this approach will certainly enable the preparation of a diversity of polycations that will be useful as building blocks for supramolecular structures, including species capable of forming large supramolecular aggregates.However, the N-alkyl benzoselenadiazolium cations have a strong affinity for halide anions, therefore the anions must be replaced for non-coordinating ions in order to enable the self-assembly of the building blocks through Se⋯N supramolecular interactions.In contrast to the telluradiazoles, and as a great advantage, the cations are not only stable in air but also in neutral to acidic aqueous media.

I 3 −
which then halogenates the solvent or benzoselenadiazole to generate HI.Crystals of [5c][I] were obtained when the reaction was carried out under an atmosphere of nitrogen, albeit in low yield.Similarly, the alkylation reaction proceeded with neat CH 3 I in anaerobic conditions, however, in this case the alkylation is incomplete due to the formation of the crystalline 4[5a] 2 [I] 2 , which was identified by X-ray diffraction.The salt [5a][I] was obtained from 4[5a] 2 [I] 2 by recrystallization from dilute solutions.

Fig. 3
Fig. 3 Packing arrangement in the crystal structures of the iodide salts of (a) 5a (b) 5b and (c) 5c.

Fig. 5
Fig. 5 Partition schemes used in the DFT-D3 fragment interaction analyses.Each fragment (region) is identified with a distinct colour.

Table 1
Crystallographic and refinement parameters for 4[5a] 2 [I] 2 and the salts of 5a-c

Table 2
Selected bond distances and angles for 4[5a] 2 [I] 2 and the salts of 5a-c 23 anions to Se has been observed in related species.For example, a virtual four-membered ring consisting of Se and Cl atoms is formed in the crystal structure of [C 2 H 2 NN(H)Se][Cl] (2.900(1) Å and (3.075 (1) Å).45The structure of [C 2 H 2 N(NC (CH 3 ) 3 )Se][Cl] also features a planar selenadiazole ring with one Se-Cl bond distance of 2.605 (1) Å.23In contrast, the GaCl 4-salt of the same selenadiazole heterocycle does not

Table 4
Contributions to the energy of interaction (kJ mol −1 ) between the halide anions and 5a in cis and trans geometries of [5a][X] (X = Cl, Br, I) ion pairs

Table 5
Calculated contributions to the binding energies for the dimers of the [5a][X] (X = Cl, Br, I) ion pairs

Table 7
Electron density and its Laplacian at the bond critical points of the Se⋯X interactions in the dimers of the [5a][X] (X = Cl, Br, I) ion pairs a Short (top) and long (bottom) Se⋯X interactions.

Table 6
Electron density and its Laplacian at the bond critical points of the Se⋯X and Se⋯N supramolecular interactions in 4[5a] 2 [I] 2 Open Access Article.Published on 06 January 2016.Downloaded on 10/06/2018 15:06:26.This article is licensed under a Creative Commons Attribution 3.0 Unported Licence.