Abstract

The study under consideration represents the computational calculations of Azo-based direct dye named p-(dimethylamino)azobenzene (DMAB) under the effect of solvents with different relative permittivities. A density functional theory (DFT) method at the B3LYP level with 6-311G++ was applied for the spectroscopic and structural analysis of the title compound. Calculations of geometric parameters (bond orders, bond lengths, and dihedral angles), electron densities, thermodynamic parameters, and orbital energies were performed for the title compound. Mulliken population analysis (MPA) as well as natural population analysis (NPA) was also performed at the B3LYP level with different solvents for finding solvent effects. In order to predict the reactivity of DMAB, molecular electrostatic potential (MESP) calculations were carried out for it. For vibrational analysis, the infrared (IR) spectra were computed for the title compound at the B3LYP/6-311G++ level in the gas phase and in different solvents with good agreement to the experimental FT-IR spectrum. The different modes of vibrations were assigned using potential energy distribution (PED). The computed Raman spectra also showed appreciable agreement with the experimental recorded Raman spectrum. The electronic absorption spectra of the title compound have been computed employing DFT/B3LYP with the 6-311G++ basis set in the gas phase and in four different solvents, that is, DMSO, ethanol, acetonitrile, and water which were compared with the experimental spectra with appreciable agreement. NBO analysis was carried out for understanding the intramolecular and intermolecular bonding of the compound and the density transfer from completely filled to unfilled orbital was found. The HOMO-LUMO energies were determined for analyzing the mechanism of intramolecular charge transfer.

1. Introduction

The increasing interest in the azo compounds is due to their versatile properties which make them favorable for use in fiber dyeing, biological reactions (such as in tumor growth inhibition), and in analytical chemistry. Azo dyes are greatly important in the field of pigments, dyes, printing systems, optical storage, and advanced materials. More than 50% colorants with the azo group are most widely employed in areas of dyeing such as that of textile fibers and in coloration of various materials. Their applications in advanced organic synthesis and biological and medical fields are well recognized [13]. Antibacterial activity of azo compounds and their metal chelates have also been reported [4]. The photosensitivity of azo dyes and their ability of bonding with N=N groups [5] make them favorable for different applications. Brilliant color and chromophoric strength have been observed with aromatic azo dyes, especially known for their excellent sublimation fastness and light-based properties [6]. Their large color range, dyeing facility, and good stabilities make them highly demanding in industry. They are highly versatile for use in different fields and in technologies like chemosensing, nonlinear optics, and for optical data storage [7]. Azobenzenes are one of the most prominent classes of photo chromic molecules with two sis and trans isomers [811]. The introduction of multiple azobenzenes in oligonucleotides assisted in regulation of DNA duplex formation [12]. The role of azobenzenes as convenient actinometers has also been reported [13].

In the present work, azo direct dye, named p-(dimethylamino)azobenzene (DMAB), was studied theoretically in the gas phase as well as in different solvents and its vibrational and structural properties have been reported in comparison with some experimental data. Multidisciplinary studies involving experimental as well as theoretical methods proved to be very useful for understanding the chemistry of the title compound and provide insight into molecular analysis [14, 15]. This manuscript mainly deals with the structural and spectroscopic (IR, Raman, and UV-Vis absorption spectra) properties, NBO analysis, and transition states of the title compound and its behavior in different solvents. The energetic behavior has been examined in gas and different solvents (water, DMSO, ethanol, and acetonitrile).

2. Materials and Methods

2.1. Experimental Detail

The title compound, p-(dimethylamino)azobenzene (DMAB), was purchased from Sigma–Aldrich, having greater than 99% purity and no additional purification was made. An IFS 66V spectrophotometer using KBr pellets was employed for the recording of the FTIR spectrum in the 4000–400 cm−1 frequency range at room temperature. The scanning speed of 30 cm−1 min−1 was utilized with spectral resolution of 2.0 cm−1. The electronic absorption spectrum was recorded on the Hitachi U-2800 spectrophotometer [16]. The FRA 106 Raman module in the range of 3500–50 cm−1 was used for the recording of the FT-Raman spectrum.

2.2. Computational Details

Computational calculations were performed by employing the Gaussian 09 software [17]. The ground-state and excited-state geometries have been optimized in analytical gradients in Gaussian 09 using the DFT method in employing B3LYP correction with the 6-311G++ basis set [18]. The optimized parameters (bond length, bond angle, and dihedral angles) of DMAB have been obtained by using the same levels of theory [19]. The electronic spectrum of the compound is calculated in gas, water, and some other solvents (water, ethanol, DMSO, and acetonitrile) at DFT/B3LYP/6-311G++. The minima of the potential energy hypersurfaces are considered to be the stationary points and were confirmed from the absence of any imaginary frequency. Moltran software was employed in computing thermodynamic parameters of DMAB [20]. The VEDA 4 program was used for the detailed vibrational analysis taking into account the computed potential energy distribution (PED). The computation of anharmonic frequencies for the fundamental vibrational modes with the corresponding IR intensities and the associated overtones was carried out [21]. A complete bond analysis was carried out using NBO analysis. The excitation energies and oscillator strengths for the lowest singlet–singlet transitions at the optimized geometry in the ground state were obtained as well by DFT calculations with the same basis set. Absorption spectra were taken with the same computational method. The spectral shifts obtained with different sets of samples were identical in most of the cases.

3. Result and Discussion

3.1. Structural Description

The compound under study is p-(dimethylamino)azobenzene. Molecular formula of the title compound is C14H15N3. The preoptimized geometry at the Hartree-Fock level was used to obtain optimized geometry at DFT/B3LYP with the 6-311G++ basis set. The optimized structure with atom numbering is shown in Figure 1(b). The optimized structural parameters with DFT/B3LYP/6-311G++ in the gas phase MP2/6-311G++ results are listed in Table 1. Geometric parameters were computed in various solvents (water and ethanol) which showed great accordance with the computed values at a higher level calculation of MP2/6-311G++. These optimized structures are helpful in computing a variety of parameters further. The optimized bond lengths of C14-H16, C14-H17, C14-H18, C15-H19, C15-H20, and C15-H20 are the same with the value of 1.09 Å. The dihedral angel between the (C22-N1-N2-C3) of two rings of DMAB is 180°. Computed values of most of the bond lengths and bond angles have shown very close agreement with the MP2 level of calculations. The actual substance exist in the solid phase, hence weak Van der Waals interactions might be developed in molecule; therefore, small differences between computed and experimental parameters occur [22].

3.2. Vibrational Studies

DMAB is consisted of 32 atoms, and therefore, it shows 90 different modes of vibration. The IR spectra of DMAB computed at the DFT/B3LYP/6-311G++ level in different solvents including gas, water, DMSO, acetonitrile, and ethanol and scaled by 0.9688 [23] are given in Figure 2. A comparison of the calculated and observed IR spectra shows good agreement. Thus, more reliable assignment of IR active modes of vibration can be carried out by comparing experimental FTIR intensities with the computed wavenumbers. The computed (scaled) frequencies as well as the PEDs of selected vibrational modes of DMAB along with the experimentally observed FTIR/FT-Raman frequencies are listed in Table 2 [24, 25]. The experimental and calculated scaled vibrational (IR) spectra are shown in Figure 2. In Figure 3, the computed Raman spectrum at the DFT/B3LYP/6-311G++ level of theory scaled by 0.9688 [23] is compared with the recorded FT-Raman spectrum which reflect the good agreement.

The carbon–hydrogen symmetric stretching vibrations of the ring were observed in the frequency region of 3194, 3160, 3139, 3112, 3099, 3065, and 2997 cm−1. Generally, in the aromatic system, stretching vibrations of C-H bond occur in the form of bands of variable intensity at 3073, 3065, 3020, and 3017 cm−1 [25]. The C-C aromatic stretching and CCH bending occur at 1575, 1564, and 1400 cm−1 [26]. The C=C bond of aromatic system gave stretching vibrations in the region of 1695, 1620, 1607, and 1566 cm−1. C-H in-plane bending was shown in the ranges of 1556, 1553, 1539, 1525, and 1485 cm−1 [27]. Asymmetric stretching vibration in C-C bond was exhibited at the frequency of 1356 cm−1. The sharp band in the region of 1146, 1335, and 1085 cm−1 is assigned as stretching vibration of CN bond. CN and NN bond stretching and bending in NCH fall in frequency ranges of 1092, 1085, 1078, 1057, and 1044 cm−1 [28]. Wagging of the CH3 group occurred at frequency of 963 and 936 cm−1. Out-of-plane bending in CH bond was appeared at 866 and 860 cm−1 [29]. The frequency ranges of 773, 739, 691, and 671 cm−1 depicted the rocking and out-of-plane bending of the CH group. In-plane bending in CCC bond angle of aromatic ring was observed at frequency ranges of 630, 563, 549, and 522 cm−1. In-plane bending in NCC bond angle was experienced at 522 and 461 cm−1 frequency ranges.

3.3. Natural Bond Orbital (NBO) Analysis

The basic objective of NBO analysis is to get chemical insight in the compound under discussion. Such insights are related to bonding concepts such as the Lewis structure and bonding type and atomic charges. Important NBO-computed parameters exhibiting the formation of Lewis and non-Lewis orbitals by valence hybrids are shown in Table 3. In order to estimate the donor-acceptor interaction, the second-order perturbation theory analysis of the Fock matrix was carried out. NBO analysis was performed on title molecule at the DFT/B3LYP/6-311G++ level in order to explain the rehybridization, intramolecular interactions, and delocalization of electron density within the molecule. The results of Fock matrix analysis are presented in Table 4. The larger the E(2) value, the stronger will be the interaction between electron donors and electron acceptors. The greater the donating tendency of electron donors, the greater will be the extent of conjugation of the whole system. The most prominent interactions in the gas phase of the title compound are due to interactions between donor BD (N1-N2) and acceptor RY (C3), exhibiting the greatest interaction energies E(2) value, that is, 4.1 kcal/mol [30].

The strong charge-transfer interaction through p-conjugated bridge results in considerate mixing of ground-state donor and acceptor orbitals, and thus, the charge-transfer band appears in the electronic absorption spectrum. Resultantly, from the donor side of the p-conjugated system, an electron density transfer is observed to its electron-withdrawing part [31].

3.4. Electronic Absorption Spectra

In order to understand the electronic transitions of DMAB, TD-SCF/DFT/B3LYP6-311G++ calculations were used. The absorption spectra were recorded in gas and other solvents (water, DMSO, ethanol, and acetonitrile) which are represented in Figure 4 in comparison with the experimental one [32]. The maximum absorption in gas took place at 499.78 nm, but the absorption peaks in all other solvents exhibited small variations than in gas possibly due to the solvent effect. The UV absorption was found to be much stronger than that in the visible region. It was found that the calculated line shape and relative strengths of peaks are in good agreement with those of the experimental results [33]. Table 5 shows the calculated and experimental wavelengths from absorption, % contribution from each transition, transition energies, and oscillator strength computed at the B3LYP/6-311G++ level for DMAB in gas and different solvents (DMSO, ethanol, acetonitrile, and water).

3.5. HOMO-LUMO Energy and Thermodynamic Parameter Studies

The HOMO and LUMO orbitals computed at the B3LYP6-311G++ level are shown in Figure 5. Their energy gaps were observed to be the lowest-energy transitions as depicted in Table 6. The difference of energy between HOMO (the highest occupied molecular orbital) and LUMO (lowest unoccupied molecular orbital) was determined to elucidate excitation energy [34] and was calculated by equation,

In a reaction, these orbitals (HOMO and LUMO) act as the main chemical participants. The HOMO orbital is an electron donor, on the other hand, LUMO has the capacity to act as an acceptor, and their energy gap predicts the reactivity as well as chemical stability of the compound. The lower energy gap between the orbitals indicates high reactivity and low stability of the compound [35]. Most probable electronic transition is from π-π. Slight variations in HOMO-LUMO energy gaps of all the solvents are shown in Table 6 along with some thermodynamic parameters like thermal correction to enthalpy (Hcorr) and thermal correction to Gibbs free energy (Gcorr).

Thermodynamic functions including internal energy, enthalpy, entropy, Gibbs free energy, and heat capacity for DMAB in the gas phase in a temperature range of 0 to 500 K are presented in Figure 6. The parameters were calculated using Moltran software. Thermodynamic parameters were found to be varied with varying temperature. The computed values of entropy were observed to increase with the increasing temperature [36]. Different thermodynamic parameters like enthalpy, Gibbs free energy, and entropy are used to describe the state as well as direction of a chemical reaction. These parameters also examine the spontaneity of a reaction, its energy profile, that is, either endothermic or exothermic, and describe temperature effects on different thermodynamic properties [37].

3.6. Mulliken Population Analysis

The Mulliken atomic charges were computed by the determination of electron population of each atom. The atomic charges present a major role in the application of quantum mechanical calculations in a molecular system. These charges are observed to affect the properties like dipole moment, electronic parameters, polarizability, and refractivity [38]. The Mulliken charge distribution of the title compound in gas and other solvents (water, DMSO, acetonitrile, and ethanol) computed by employing the DFT/B3LYP/6-311G++ level is presented in Table 7. Figure 7 exhibits the comparison of MPA with NPA. Natural population analysis better describes the electronic charge distribution over the system, and numerical stability is also higher as compared to MPA. As charge distribution seems to be method sensitive that is why we have made its comparison within three types of analyses. The table represents that all hydrogen atoms have positive charges. The carbon atoms in the ring carry both positive and negative charges. The substituted nitrogen has negative charge [39].

3.7. Molecular Electrostatic Potential (MESP)

The main objective of molecular electrostatic potential calculations is to analyze electron density distribution over the molecule. It is much useful for the prediction of reactivity of the molecule towards nucleophilic and electrophilic attacks [40]. It also guesses which types of intermolecular interactions are possible [41]. In the case of electrophilic attack in the molecule, it will preferably go to the most negative portion of the molecule; that is, the site with most dominant electron effect.

Figure 8 shows the molecular electrostatic potential map of DMAB created on the isodensity surface. The portion with red color is representative of negative charge with higher electron density, while the region in blue is electropositive with lower electron density. The figure depicts that maximum electron density is concentrated near nitrogen atom so this will be the best site for electrophilic attack. The positive electrostatic potential is concentrated near carbon atoms so it will be a good site for attack of nucleophile.

4. Conclusion

This manuscript presents thorough analysis of the geometrical parameters, Mulliken population analysis in comparison with natural population analysis electronic absorption spectra, NBO analyses, HOMO–LUMO energy orbitals, MESP isodensity surface computation and the vibrational spectra of the title compound, and DMAB in the gas phase as well as in four different solvents (water, DMSO, acetonitrile, and ethanol). This compound was investigated by employing the DFT/B3LYP/6-311G++ method. All the significant vibrations observed in the experimental spectra (FT-IR/FT-Raman) and computed frequency (IR/Raman) spectra of the compound are assigned to the various modes of vibration, and most of the modes have expected wavenumber range. Small discrepancies in vibrational wavenumbers could be due to different media in the recorded and computed spectra. The complete vibrational assignment is made on the basis of potential energy distribution (PED). Both vibrational spectra (IR/Raman) and electronic spectra help to predict diverse properties of the title molecule. The excellent agreement found in the observed and calculated vibrational spectra discloses the advantages of a higher basis set for quantum chemical calculation. NBO analysis provides an efficient method to explain accurate calculations into chemical insights. Experimental and calculated values of most of the bond lengths and bond angles of the title compound are found very close. The thermodynamic parameters were studied to explain the thermal behavior of the title compound. Mulliken charge investigation provides insight into chemical reactivity of our studied compound.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

Authors are highly thankful to Higher Education Commission (HEC), Pakistan for the financial support for the current study (NRPU-2016-17/Project no. 5613).