Abstract

The Fourier transform infrared and Fourier transform Raman spectra of (E)-4-(benzylideneamino)-N-carbamimidoyl benzenesulfonamide were recorded and analyzed. Geometry and harmonic vibrational wavenumbers were calculated theoretically using Gaussian 03 set of quantum chemistry codes. Calculations were performed at the Hartree-Fock (HF) and density functional theory (DFT; B3PW91, B3LYP) levels of theory. The calculated wavenumbers (B3LYP) agree well with the observed wavenumbers. Potential energy distribution is done using GAR2PED program. The red shift of the N-H stretching band in the infrared spectrum from the computed wavenumber indicates the weakening of the N-H bond. The geometrical parameters of the title compound are in agreement with that of reported similar derivatives. The calculated first hyperpolarizability is comparable with the reported value of similar derivative and may be an attractive object for further studies of nonlinear optics. Potential energy surface scan studies have been carried out to understand the stability of planar and nonplanar structures of the molecule.

1. Introduction

Benzenesulfonamide derivatives find wide applications for the synthesis of pharmaceutical products which have bactericidal properties and of various bioactive agents, artificial fibers, dyes, and plasticizers, and the synthesis of high molecular weight substances. Sulfonamides are chemotherapeutics most commonly used in veterinary practices because of their inexpensiveness and wide spectrum antimicrobial activity [1]. Sulfonamides represent one of the classical chemotypes associated with potent CA inhibition [25]. The chemistry of sulfonamides has been known as synthons in the preparation of various valuable biologically active compounds [6, 7] used as an antibacterial [8], protease inhibitor [9], diuretic [10], antitumor [11], and hypoglycaemic [12]. Singh et al. [13] reported the synthesis and antimicrobial activity of Schiff’s and N-Mannich bases of Isatin and its derivatives with 4-amino-N-carbamimidoyl benzene sulfonamide. There has been growing interest in using organic materials for nonlinear optical (NLO) devices, functioning as second harmonic generators, frequency converters, electro-optical modulators, and so forth because of the large second order electric susceptibilities of organic materials. The organic compound showing high hyperpolarizability are those containing an electron-donating group and an electron-withdrawing group interacting through a system of conjugated double bonds. In the case of sulfonamides, the electron-withdrawing group is the sulfonyl group [14, 15]. To our knowledge, no theoretical HF or density functional theory (DFT) calculations or detailed vibrational infrared or Raman analyses have been performed on the title compound.

2. Experimental

All the chemicals were procured from Sigma-Aldrich, USA. 0.5 mg of sulphaguanidine and 0.3 mL of benzaldehyde in 20 mL ethyl alcohol are refluxed for 3 hrs. The white precipitate was filtered off, washed with ethanol, and dried. Elemental analysis: found/calculated (%): C 55.50/55.63; H 4.59/4.63; N 18.40/18.50; S 10.58/10.59. The FT-IR spectrum (Figure 1) was recorded on a DR/Jasco FT-IR 6300 spectrometer in KBr pellets, number of scans 16, resolution 2 cm−1. The FT-Raman spectrum (Figure 2) was obtained on a BRUKER RFS 100/S, Germany. For excitation of the spectrum the emission of an Nd:YAG laser was used, excitation wavelength 1064 nm, maximal power 150 mW, measurement of solid sample. One thousand scans were accumulated with a total registration time of about 30 min. The spectral resolution after apodization was 2 cm−1.

3. Computational Details

Calculations of the title compound were carried out with Gaussin 03 program [16] using the HF/6-31G*, B3PW91/6-31G* and B3LYP/6-31G*, basis sets to predict the molecular structure and vibrational wavenumbers. The wavenumber values computed contain known systematic errors [17] and we, therefore, have used the scaling factor values of 0.8929 and 0.9613 for HF and DFT basis sets. Parameters corresponding to optimized geometry (B3LYP) of the title compound (Figure 3) are given in Table 1. The absence of imaginary values of wavenumbers on the calculated vibrational spectrum confirms that the structure deduced corresponds to minimum energy. Potential energy distribution is done using GAR2PED program [18]. The assignment of the calculated wavenumbers is aided by the animation option of MOLEKEL program, which gives a visual presentation of the vibrational modes [19, 20]. Potential energy surface scan studies have been carried out to understand the stability of planar and non planar structures of the molecule. The profiles of potential energy surface for torsion angles C15–N14–C12–C3, N28–S25–C22–C20, N31–C30–N28–S25 and N32–C30–N28–S25 are given in Figures 47. The energy is minimum for −177.0 (−1309.42093 Hartree), 86.0 (−1309.41690 Hartree), 12.6 (−1309.42098 Hartree), and −169.9 (−1309.42070 Hartree) for the above torsion angles.

4. Results and Discussion

4.1. IR and Raman Spectra

The observed IR, Raman bands, and calculated (scaled) wavenumbers, and assignments are given in Table 2.

The asymmetric and symmetric stretching modes of SO2 group appear in the region 1360–1310 and 1165–1135 cm−1 [21]. The observed bands at 1358, 1130 cm−1 in the IR spectrum, 1360, 1133 cm−1 in the Raman spectrum, and 1371, 1148 cm−1 theoretically are assigned as SO2 stretching modes. These modes are not pure, but contain significant contributions from other modes also. Although the region of SO2 scissors (560±40 cm−1) and that of SO2 wagging (500±55 cm−1) partly overlap, the two vibrations appear separately [21]. These deformation bands of SO2 are assigned at 515, 469 cm−1 in the IR spectrum and at 510, 457 cm−1 theoretically. Chohan et al. [22] reported the SO2 stretching vibrations at 1345, 1110 cm−1 and SN and CS stretching modes at 833 cm−1 for sulfonamide derivatives. Hangen et al. [23] reported SO2 modes at 1314, 1308, 1274, 1157, 1147, and 1133 cm−1 and SN stretching modes at 917, 920, 932, and 948 cm−1 for sulfonamide derivatives. The twisting mode of SO2 is expected at around 350 cm−1 [21]. The calculated bands (B3LYP) at 365, 298 cm−1 are assigned as the twisting and rocking modes of SO2. The NS stretching vibration is expected [21] in the region 905±30 cm−1 and the band observed at 841 cm−1 in the IR spectrum and at 831 cm−1 theoretically are assigned as 𝜐SN. The C–S stretching vibration is assigned at 692 cm−1 (IR) and at 679 cm−1 (DFT) cm−1 [21].

The C=N stretching bands [24] are observed in the region 1550–1700 cm−1 and for the title compound, the B3LYP calculation give these bands at 1687(C12–N14), 1641 cm−1 (C30–N31). The CN stretching vibrations [21] are moderately to strongly active in the region 1275±55 cm−1. Primary aromatic amines with nitrogen directly on the ring absorb at 1330–1200 cm−1 because of the stretching of the phenyl carbon nitrogen bond [25]. For the title compound C15–N14 stretching mode is observed at 1281 cm−1 theoretically. Panicker et al. [26] reported the CN stretching modes at 1219, 1237 cm−1 (IR), 1222 cm−1 (Raman) and at 1292, 1234, 1200 cm−1 theoretically. The CN stretching modes C30–N28 and C30–N32 are assigned at 1100 and 942 cm−1 theoretically which is expected in the range 950–1115 cm−1 [21].

The NH2 stretching modes are expected in the region 3260–3390 cm−1 [21] and in the present case the bands observed at 3421 cm−1 in the IR spectrum and bands at 3501, 3399 cm−1 (DFT) are assigned as NH2 stretching modes. The bands corresponding to δNH2 vibrations are expected in the region 1610±30 cm−1 [21]. In the present case, the band at 1620 cm−1 (IR), 1625 cm−1 (Raman) and 1602 cm−1 (DFT) is assigned as δNH2. The rocking/twisting mode of NH2 is expected in the region 1195±90 cm−1 and this mode is assigned at 1182 cm−1 in IR, 1188 cm−1 in Raman spectrum at 1181 (DFT) cm−1. The wagging mode of NH2 is assigned at 778 cm−1 theoretically which is expected in the region 840±55 cm−1 [21]. The torsional NH2 mode is expected below 400 cm−1 [21] and these modes are identified and assigned (Table 2). For sulfonamide derivatives the NH2 modes are reported at 3390, 3395, at 3399 cm−1 [22].

The NH stretching vibrations generally give rise to bands at 3500–3300 cm−1 [27, 28]. In the present study, the N28–H29 stretching bands have split into a doublet, 3351, 3222 cm−1 in the IR spectrum owing to the Davydov coupling between neighboring units. The splitting of about 129 cm−1 in the IR spectrum is due to the strong intermolecular hydrogen bonding. Furthermore the NH stretching wavenumber is red shifted by 129 cm−1 in the IR spectrum with a strong intensity from the computed wavenumber which indicate the weakening of the NH bond [29]. The DFT calculations give the N31–H35 stretching mode at 3453 cm−1. In N-mono substituted amides, the in plane bending frequency and the resonance stiffened CN bands stretching frequency fall close together and therefore interact. The NH deformation band of guanidine structural motif is expected in the region 1395±25 cm−1 [21, 25]. The DFT calculations give these modes at 1442 and 1485 cm−1. The bands observed at 1500, 1450 cm−1 in the IR spectrum and at 1486, 1446 cm−1 in the Raman spectrum are assigned as the deformation bands of NH group. The out-of-plane NH deformation is expected in the region 650±50 cm−1 [21] and bands at 727 (IR), 729 (Raman) and 706, 608 cm−1 (DFT) are assigned as these modes.

Since the identification of all the normal modes of vibration of large molecules is not trivial, we tried to simplify the problem by considering each molecule as a substituted benzene. Such an idea has already been successfully utilized by several workers for the vibrational assignments of molecules containing homo- and heteroaromatic rings [3032]. In the following discussion, mono- and parasubstituted phenyl rings are designated as PhI and PhII, respectively. The modes in the phenyl rings will differ in wavenumber and the magnitude of splitting will depend on the strength of the interactions between different parts (internal coordinates) of the two rings. For some modes, this splitting is so small that they may be considered as quasidegenerate, and for some other modes a significant amount of splitting is observed. Such observations have already been reported [33, 34].

The phenyl CH stretching modes occur above 3000 cm−1 and are typically exhibited as multiplicity of weak to moderate bands compared with the aliphatic CH stretching [35]. In the present case, the DFT calculations give 𝜐CH modes of the phenyl rings in the range 3132–3057 cm−1. The bands observed at 3026, 3073 cm−1 in the IR spectrum and at 3066 cm−1 in the Raman spectrum were assigned as CH stretching modes of the phenyl rings.

The benzene ring possesses six ring stretching vibrations, of which the four with the highest wavenumbers (occurring near 1600, 1580, 1490, and 1440 cm−1) are good group vibrations. With heavy substituents, the bands tend to shift to somewhat lower wavenumbers. In the absence of ring conjugation, the band at 1580 cm−1 is usually weaker than that at 1600 cm−1. In the case of C=O substitution, the band near 1490 cm−1 can be very weak. The fifth ring stretching vibration is active near 1315±65 cm−1, a region that overlaps strongly with that of the CH in-plane deformation. The sixth ring stretching vibration, or the ring breathing mode, appears as a weak band near 1000 cm−1, in mono-, 1, 3-di-, and 1, 3, 5-trisubstituted benzenes. In the otherwise substituted benzenes, however, this vibration is substituent sensitive and difficult to distinguish from the ring in-plane deformation [21, 36].

The 𝜐Ph modes are expected in the region 1285–1610 cm−1 and 1280–1630 cm−1 for PhI and PhII rings, respectively [21]. The DFT calculations give the Ph stretching modes in the range 1290–1598 cm−1 and 1265–1578 cm−1 for PhI and PhII, respectively. For the title compound the phenyl ring stretching modes are observed at 1567, 1500 (IR), 1581, 1486 cm−1 (Raman) and 1587, 1528, 1403, 1130 cm−1 (IR), 1133 cm−1 (Raman), for PhI and PhII, respectively. Most of the modes are not pure but contain significant contributions from other modes also. The ring breathing mode of the parasubstituted benzenes with entirely different substituents [36] has been reported in the interval 780–880 cm−1. For the title compound, the ring breathing mode PhII is observed at 757 cm−1 in the IR spectrum which is supported by the computational result at 747 cm−1. The ring breathing mode of PhI is observed at 1001 cm−1 in the Raman spectrum and the theoretical calculated value is 990 cm−1, as expected [21]. The in-plane and out-of-plane CH deformation bands of the phenyl ring are expected above and below 1000 cm−1 [21]. The in-plane CH deformation bands are assigned at 1319, 1160, 1156, 1096, 1076, and 1013 cm−1 for PhI and at 1317, 1230, 1103, and 1044 cm−1 for PhII theoretically. The strong γCH occurring at 840±50 cm−1 is typical for 1, 4-disubstituted benzenes and the band observed at 804 cm−1 in the IR spectrum is assigned to this mode [21]. The in-plane and out-of-plane deformation modes of the phenyl ring are also identified and assigned (Table 2).

The IR bands in the region 2882–1915 cm−1 and their large broadening support the intramolecular hydrogen bonding [37].

4.2. Geometrical Parameters and First Hyperpolarizability

To the best of our knowledge, no X-ray crystallographic data of this molecule has yet been established. However, the theoretical results obtained are almost comparable with the reported structural parameters of the parent molecules.

Loughrey et al. [38] reported the bond lengths, S25–O27 = 1.4337, S25–O26 = 1.4256, S25–N28 = 1.6051, S25–C22 = 1.7737, C15–N14 = 1.4212, C12–N14 = 1.2712, whereas the corresponding values for the title compound are 1.4570, 1.4658, 1.7054, 1.7883, 1.4020, and 1.2816 Å. For the title compound, the DFT calculations give the bond angles, C22–S25–O27 = 108.9, C22–S25–O26 = 107.7, S25–N28–H29 = 112.2, O27–S25–O26 = 122.2, O27–S25–N28 = 108.0, O26–S25–N28 = 101.5, N28–S25–C22 = 107.7, C15–N14–C12 = 120.1, S25–C22–C18 = 118.7, S25–C22–C20 = 120.1, C22–C18–C16 = 119.3, C18–C16–C15 = 120.5, C16–C15–C17 = 119.2, N14–C15–C17 = 123.0, N14–C15–C16 = 117.8, C15–C17–C20 = 120.7, N14–C12–C3 = 122.8°, whereas the corresponding reported values are 106.5, 107.4, 110.0, 119.5, 106.1, 107.7, 109.3, 118.8, 120.5, 119.1, 120.4, 119.9, 120.3, 119.3, 122.4, 118.2, 120.4, and 124.4° [38].

Loughrey et al. [38] reported the torsion angles, O27–S25–C22–C18 = −142.0, O27–S25–C22–C20 = 38.5, O26–S25–C22–C18 = −12.8, O26–S25–C22–C20 = 167.7, N28–S25–C22–C18 = 103.8, N28–S25–C22–C20 = −75.7, C12–N14–C15–C16 = 143.9, C12–N14–C15–C17 = −39.2, C15–N14–C12–C3 = 117.6, S25–C22–C18–C16 = −179.9, C20–C22–C18–C16 = −0.4, S25–C22–C20–C17 = 179.2, C18–C16–C15–N14 = 178.1, N14–C15–C17–C20 = −178.6, N14–C12–C3–C2 = 8.1, N14–C12–C3–C4 = −172.6°. For the title compound, the corresponding torsion angles are 148.1, −30.8, 13.7, −165.2, −95.1, 86.0, −140.4, 42.1, −177.0, −179.4, −0.4, 178.5, −179.6, 178.6, 1.1, and −179.0°.

Petrov et al. [39] reported the molecular structure and conformations of benzenesulfonamide by gas electron diffraction and quantum chemical calculations and according to their results, the bond lengths, CS, SN, SO vary in the range 1.7756–1.7930, 1.6630–1.6925, 1.4284–1.4450 Å the bond angles, CSN, CSO, NSO, HNS, HNH vary in the range, 103.9–107.1, 107.6–107.8, 105.5–107.7, 111.0–113.7, 112.6–113.6 Å. These values are in agreement with the corresponding values for the title compound.

Lasibal et al. [40] reported the bond lengths SO = 1.4269–1.4291, SN = 1.6202, SC = 1.7582, N32–C30 = 1.4103, N31–C30 = 1.2723, N28–C30 = 1.3483 Å, whereas the corresponding values in the present case are 1.4658–1.4570, 1.7054, 1.7833, 1.3930, 1.2764, and 1.4033 Å. The values of bond angles O27–S25–O26 = 118.6, O26–S25–N28 = 108.9, O27–S25–N28 = 104.9, O27, 26–S25–C22 = 107.9–108.3, N28–S25–C22 = 107.9, C30–N28–S25 = 123.0° reported by Lasibal et al. [40] are in agreement with our values.

At C3 position, the bond angles C4–C3–C2, C4–C3–C12, and C2–C3–C12 are 119.2, 119.3, and 121.5° respectively. This asymmetry in angles reveals the interaction between azo-methane and the phenyl groups. At C15 position the angles C17–C15–N14 is increased by 3° and C16–C15–N14 is reduced by 2.2° from 120° which reveals the interaction between N14 and H21 atoms. At C22 position C20–C22–S25 = 120.1° and C20–C22–C18 = 121.2° which shows the interaction between SO2 group with H24 atom. At N28 position, S25–N28–H29 is reduced by 7.8° and S25–N28–C30 is increased by 6.8° from 120° which shows the interaction between H29 and O27. The C12 = N14 moiety is slightly tilted from monosubstituted phenyl ring as is evident from the torsion angles C5–C4–C3–C12 = −179.8, C4–C3–C12–N14 = −179.0, C1–C2–C3–C12 = 179.9, C2–C3–C12–N14 = 177.4° and is more tilted from the para substituted phenyl ring as is evident from the torsion angles C20–C17–C15–N14 = 178.6, C17–C15–C14–C12 = 42.1, C18–C16–C15–N14 = 179.6, and C16–C15–N14–C12 = −140.4°. The torsion angle S25–N28–C30–N32 = −169.9 and S25–N28–C30–N31 = 12.6°, which shows that the N31 and N32 atoms are in different planes.

Nonlinear optics deal with the interaction of applied electromagnetic fields in various materials to generate new electromagnetic fields, altered in wavenumber, phase, or other physical properties [41]. Organic molecules able to manipulate photonic signals efficiently are of importance in technologies such as optical communication, optical computing, and dynamic image processing [42, 43]. In this context, the dynamic first hyperpolarizability of the title compound is also calculated in the present study. The first hyperpolarizability (𝛽0) of this novel molecular system is calculated using B3LYP method, based on the finite field approach. In the presence of an applied electric field, the energy of a system is a function of the electric field. First hyperpolarizability is a third rank tensor that can be described by a 3 × 3 × 3 matrix. The 27 components of the 3D matrix can be reduced to 10 components due to the Kleinman symmetry [44]. The components of 𝛽 are defined as the coefficients in the Taylor series expansion of the energy in the external electric field. When the electric field is weak and homogeneous, this expansion becomes𝐸=𝐸0𝑖𝜇𝑖𝐹𝑖12𝑖𝑗𝛼𝑖𝑗𝐹𝑖𝐹𝑗16𝑖𝑗𝑘𝛽𝑖𝑗𝑘𝐹𝑖𝐹𝑗𝐹𝑘124𝑖𝑗𝑘𝑙𝛾𝑖𝑗𝑘𝑙𝐹𝑖𝐹𝑗𝐹𝑘𝐹𝑙+,(1) where 𝐸0 is the energy of the unperturbed molecule, 𝐹𝑖 is the field at the origin, 𝜇𝑖, 𝛼𝑖𝑗, 𝛽𝑖𝑗𝑘, and 𝛾𝑖𝑗𝑘𝑙 are the components of dipole moment, polarizability, the first hyperpolarizabilities, and second hyperpolarizabilities, respectively. The calculated first hyperpolarizability of the title compound is 6.25×1030esu, which is comparable with the reported values of similar derivatives, but experimental evaluation of this data is not readily available. Kucharski et al. [45] reported the first hyperpolarizability of certain sulfonamide amphiphiles by calculation and hyper-Rayleigh scattering in the range 0.2156-0.189×1030 esu. We conclude that the title compound is an attractive object for future studies of nonlinear optical properties.

5. Conclusion

The FT-IR and FT-Raman spectrum of the title compound were recorded and analyzed. The molecular geometry and vibrational wavenumbers were calculated using HF and DFT methods, and the optimized geometrical parameters (B3LYP) are in agreement with that of reported similar derivatives. The red shift of the N–H stretching band in the infrared spectrum from the computed wavenumber indicates the weakening of the N–H bond. The calculated first hyperpolarizability is comparable with the reported value of similar derivative and may be an attractive object for further studies of non linear optics. Potential energy surface scan studies have been carried out to understand the stability of planar and non planar structures of the molecule.