Electron Transport and Nonlinear Optical Properties of Substituted Aryldimesityl Boranes: A DFT Study

A comprehensive theoretical study was carried out on a series of aryldimesityl borane (DMB) derivatives using Density Functional theory. Optimized geometries and electronic parameters like electron affinity, reorganization energy, frontiers molecular contours, polarizability and hyperpolarizability have been calculated by employing B3PW91/6-311++G (d, p) level of theory. Our results show that the Hammett function and geometrical parameters correlates well with the reorganization energies and hyperpolarizability for the series of DMB derivatives studied in this work. The orbital energy study reveals that the electron releasing substituents increase the LUMO energies and electron withdrawing substituents decrease the LUMO energies, reflecting the electron transport character of aryldimesityl borane derivatives. From frontier molecular orbitals diagram it is evident that mesityl rings act as the donor, while the phenylene and Boron atom appear as acceptors in these systems. The calculated hyperpolarizability of secondary amine derivative of DMB is 40 times higher than DMB (1). The electronic excitation contributions to the hyperpolarizability studied by using TDDFT calculation shows that hyperpolarizability correlates well with dipole moment in ground and excited state and excitation energy in terms of the two-level model. Thus the results of these calculations can be helpful in designing the DMB derivatives for efficient electron transport and nonlinear optical material by appropriate substitution with electron releasing or withdrawing substituents on phenyl ring of DMB system.


Introduction
A great deal of interest has recently been aroused in organic materials for their active role in opto-electronic devices such as photovoltaic cells and organic lightemitting diodes (OLEDs) [1]. The practically useful OLEDs require good luminescence, efficient charge transport and fast energy transfer properties. The simplest OLED consist of anode, an electron transport material (ETM) layer, a hole transport material (HTM) layer and a cathode. Organic material used as ETM lag behind due to obstacles such as scarce solubility, difficult processibility, and instability in air [2]. Thus, one of the key challenges in developing OLEDs is the design of stable ETM with high efficiency of electron ejection and adequate electron mobility.
The photochemical investigation of organoboron compounds reveals that triphenyl boranes in presence of donor solvents undergo photodegradation, in contrast to trimesityl boranes, which were found photochemically inert in all these solvents. It was observed that two mesityl groups were significant in providing an extra degree of stability, which attracted considerable research interest in carrying out photophysical studies of p-substituted aryldimesityl boranes (DMB) [2][3][4][5][6]. The three coordinated boron is electron deficient with its vacant pz orbital. Due to strong p-acceptor character, a significant delocalization occurs when conjugated with an adjacent organic p-system. The acceptor character of the dimesitylboryl group was found similar to cyano group than the nitro group. But unlike other p-acceptors, it is a weak s-electron acceptor as electronegativity of boron is less than that of carbon. Thus, when boron bonded to organic systems it may even function as a s-electron donor [4][5][6][7][8][9]. The presence of vacant pz orbital make boron compounds prone to nucleophiles, resulting in either bond cleavage or the formation of a four-coordinate borate species that causes interruption in conjugation with adjacent p-systems. However, the organic system like mesityl (2,4,6-trimethylphenyl), forms a propeller like structure around the borane, and a cage is formed by methyl groups on aryl ring around the vacant pz-orbital [3][4][5][6][7], which enhances inertness of the borane compound. It has been found that these kinds of materials are applicable in a wide range of fields, such as nonlinear optics, single and biphotonic absorption materials, electron transporting and emissive materials in organic light emitting devices as well as chemosensors [7][8][9][10][11][12][13][14][15][16][17][18][19][20][21]. Shirota and others have reported several three-coordinate organoboron compounds as electron transport/hole-blocking materials for OLEDs [21][22][23][24][25][26][27].
The development of high-performance charge transport materials is a key issue for the fabrication of high-performance organic light emitting devices. The charge transport in molecular organic materials can be viewed in light of two theories, the band theory [28] and the hopping model [29]. In the band theory, transport of charge is an activationless process occurring through bands formed by the overlapping MOs between neighboring molecules. However in the hopping model charge transfer occurs between small coupled neighboring molecules. According to Marcus electron-transfer theory the rate of electron or hole transfer k et is given by following equation [30,31] Where l and DH ab are the reorganization energy for the intramolecular electron transfer and the electronic coupling integral between donor-acceptor pair, respectively, and DG 0 is the Gibbs free energy change of the process. Efforts for several years in theoretical investigations of charge transport properties for several OLED materials confirm that both the reorganization energies (l) and electron coupling parameter (DH ab ) play an important role in the electron/hole transport process [32,33] and can be evaluated by Quantum chemical calculations. The most well-known electron transport material used in OLEDs is Alq3, in which three 8-hydroxyquinolato (q) ligands chelate to the Al 3+ ion in an octahedral environment. Tris(8-hydroxyquinolinato)aluminum (Alq3) has been widely used in organic light-emitting diodes (OLEDs)both as electron transport and light-emitting materials. According to the Marcus theory electron transport mainly depends on electronic coupling between donor and acceptor and the reorganization parameters for the two states upon electron transfer. Lin et al, calculated the reorganization energies of Alq3 for hole and electron transport. However, the slightly smaller value of reorganization energy for hole transport (0.242 eV) as compared to reorganization energy for electron transport (0.276 eV) suggest that Alq3 posses better hole transport than electron transport property, apparently contradicting the experimental observations [34]. Therefore, it was concluded that the electronic coupling factor (DH ab )may be a dominant factor in influencing charge transport in a system [35]. The first-principle bandstructure calculations using density-functional theory with generalized gradient approximation by Yang et al., reveals that the inter band gaps within the unoccupied bands are generally smaller than those for the occupied bands, which indicate that the electron can hop from one band to another much easier than the hole, through electron-phonon coupling for instance, therefore, explaining the larger mobility for the electron than for the hole [36]. In the present work, we have adopted quantum chemical calculations to analyze several substituted DMB derivatives, at electronic structure level, for their potential character as electron transport materials. The choice of substituents has been made in such a way that they range from strong electron releasing groups (-(CH 3 ) 2 N) to strong electron withdrawing ones (-NO 2 ). Our investigation on boron compounds is motivated by the scope of this material as potential conducting material having application in OLEDs.
Photo-physical properties of DMB make them suitable candidates having potential use in the design and development of materials exhibiting large secondorder nonlinear optical (NLO) properties. NLO response of molecular system is governed by the parameters depending on the extend of the conjugation and the strength of substituent groups in terms of their donor or acceptor nature [37]. Earlier studies of the structure-property relationships of such molecules indicate that the hyperpolarizability (b) of the molecule changes with a change in donor and acceptor strength and with increase in the extent of pconjugation [38,39]. At molecular level, nonlinear optical susceptibilities are determined by the first and second hyperpolarizabilities. The present work is intended for the optimization of the molecular structure-NLO properties relationships of DMB with conventional electron donor or acceptor substituent groups on the phenyl ring. Investigations have been made to study the substitution effect in the DMB derivatives and its influence on their energy levels, nonlinear optics and electrontransport properties. In particular, the molecular geometries, HOMO and LUMO levels, ionization potential, electron affinity, reorganization energy, and molecular orbitals shapes have been investigated as they are directly related to such OLED's characteristics like efficiency and durability.

Theoretical Calculations
In this work all theoretical calculations were carried out at the density functional theory (DFT) level of theory, performed by using Gaussian 03 computational package [40]. Geometry optimizations were performed on isolated entities in gaseous phase, employing Becke's 3 Perdew Wang 91 (B3PW91) [41] exchangecorrelation functional using 6-311++G (d, p) basis set [42]. All the geometries were characterized as minima with zero imaginary frequencies. The neutral molecules were treated as closed-shell systems; while for the radical anion openshell system optimizations were carried out using a spin unrestricted wave function (UB3PW91 procedure). Vertical electronic excitation energies were computed by using the time-dependent density functional level of theory (TD-DFT) [43]. The TDDFT/B3PW91 level of theory has been found quite efficient, for performing significantly single-excitation calculations for the low-lying valence excited states of both closed shell and open-shell molecules [44]. The static polarizability and hyperpolarizabilities were calculated using derivative method which relates different derivatives of energy or dipole moment to various coefficients of power series expansions. The total static dipole moment (m), average linear polarizability (a) and first-order hyperpolarizability (b) using the x, y, z components from Gaussian 03 W [40] output were calculated by the equation 2, 3 and 4 [45][46][47].

Results and Discussion
Ab initio calculations were carried out on a series of DMB monomers with donor (2-7) and acceptor (8)(9)(10)(11)(12)(13) substituents as shown in Figure 1. The substituents were chosen on the basis of their ability to affect electron density of the phenyl ring, covering possibly wide range of the Hammett parameter (s) [48] from strong electron releasing groups to strong electron withdrawing ones. According to Hammett equation, the electron releasing substituents have negative and the electron withdrawing groups have positive value of Hammett parameter. In this work, we present a comparative effect of substitution by correlating theoretically calculated values of interest with Hammett parameter and the analysis of the observed relationships. In the first section, we have studied the electronic properties and their correlation with Hammett parameter and geometrical parameters. In second section, we have studied the nonlinear optical properties of DMB in terms of polarizability and hyperpolarizability for all systems. The calculations have been carried out for neutral molecules, radical-anion as well as radical-cation.

Optimized Geometry
The optimized geometrical structures of the studied systems 1-13 are shown in Figure S1, the coordinates of all optimized DMB derivatives are given in Table S1 and some selected bond lengths, bond angles of parent DMB are listed in Table  S2. We found a reasonable agreement between the calculated geometrical parameters (bond lengths and bond angles) and the experimental values obtained from the X-ray diffraction data [49] for the DMB(1) molecules (Table-S2). The slight deviations observed may be attributed to the fact that theoretical calculations are done for single isolated molecule in gaseous phase and the experimental values are those of the crystalline solid phase. This indicates that the adopted density functional and basis set are feasible for the studied systems. The un-substituted DMB(1) is composed of two mesityl rings and a phenyl ring bonded to a central boron atom. The central BC 3 moiety is planar and forms an equilateral triangle. The two mesityl and phenyl rings are twisted from that plane so the whole molecule can be visualized as a slightly distorted propeller. The C-B bond lengths for two mesityl rings are equal to1.570 Å , and C-B length for phenyl ring is 1.581 Å and the torsion angle of the phenyl ring with the mesityl amounts to about 21.52 0 in the neutral state. The optimized geometric parameters of DMB(1) calculated in this study are presented in Table 1. We studied in detail different derivatives of DMB, choosing carefully different functional groups in varying order of their electron releasing (N(CH 3 ) 2 , NH 2 , OH, OCH 3 , C 2 H 5 , CH 3 ) and electron withdrawing(F, Cl, Br, COOH, CN, NO 2 ) character at para-position (phenyl) of DMB, as shown in Figure 1. From the optimized geometries, we found that upon substitution of hydrogen with electron releasing groups, C 18 -B bond length increase as [2,3,455,657] compared to parent molecule (1), and   ]. The effect of substitution on the aromaticity was analyzed by bond length alternation (BLA). We calculated a local BLA associated with the C 23 -C 20 , C 18 -B, C 20 5C 18 lengths according to the definition given by Fu et al. [50]. Thus, the BLA is calculated by subtracting the arithmetic mean of the two central C-C bonds and the central C5C bond length The increase in electron density due to electron donating groups on parent DMB show an increased double bond character of C 23 -C 20 and C 18 -B bonds and decreased double bond character of C 20 5C 18 bond in the order of 7,6,5,4,3,2. The BLA values increases from 0.056 (2) ,0.059(3) ,0.061(4) 50.061(5) ,0.063(6) ,0.064 (7), indicating that the degree of conjugation is increasing upon the substitution of the parent DMB with electron donating groups. However, there is no good correlation between BLA and the electron withdrawing capacity of substituents. The DMB (12) and DMB(13) have maximum electron withdrawing capacity, with highest BLA values 0.083 and 0.085 respectively, but in case of 8, 9,10 electron withdrawing capacity decreases as 8.9.10, but the BLA values does not differ much (0.070(8) 50.070 (9) ,0.071 (10).

Electron Transport Properties
To study electron transport property of DMB systems, we optimized the anion geometry to understand the effects of charge injection on the molecular conformational stability, in terms of structural distortions and changes in electronic structure, by varying the nature of substituents. The calculations show that neutral molecules show an overall negligible change in geometrical parameters upon the formation of radical anion in all the derivatives (Table 1). There is a change in dihedral angle ,CBCC between the mesityl ring and phenyl ring from 22.15 0 to 24.05 0 and the bond length of C 18 -B bond slightly changes from 1.581 Å to 1.594 Å in case of DMB(1). The structural reorganization from neutral to anionic species, reflected by the change in C 18 -B bond length upon electron injection, shows the following trend; (-N(CH 3 ) 2 (2) ,-NH 2 (3) ,-OH(5) ,-OCH 3 (4) ,-C 2 H 5 (6) ,-CH 3 (7) ,-F(8) ,-Cl (9) ,-Br (10)  reflects that the contribution to activation barrier from geometrical reorganization for electron transport is small in case of electron releasing derivatives as compared to electron withdrawing derivatives. The magnitude of structural modifications increases on moving from the electron releasing substituents to the electron withdrawing substituents. This indicate the existence of strong correlation between Hammett parameter (s) and the geometrical parameters, as is evident from the plots of geometrical parameters verus s, shown in Figure 2 and Figure  S2. The enhanced change in dihedral angle due to electron withdrawing groups represents the greater tendency of these molecules to undergo structural change upon radical-anion formation from their neutral states, as compared to molecules with electron releasing groups. The main influence of the injection of one negative charge is manifested in the increase in bond length of B-C 18 substituted DMBs radical anion, which may be responsible for decrease in its rigidity. However, we found that, in general, the optimized anionic geometries show small structural changes comparative to the neutral molecules, the structural changes being close to Alq 3 [35]. Thus, the injection of one negative charge does not affect the structural stability of DMBs derivatives and this property may allow tuning them into better electron transport systems.
We calculated electron affinities EA, ionization potential (IP) and reorganization energies (l 1 , l 2 , and l i ) for electron transport using the Figure 3. According to Marcus Model, the rate of intramolecular electron transfer depends mainly on the reorganization energy and the electronic coupling between the donor-acceptor pair, in addition to the overall exergonicity of the process. For a self-exchange electron transfer process, the Gibbs free energy change (DG) is zero and the rate of electron transfer only depends on intrinsic activation barrier, signified by (inner and outer) reorganization energy (l), and the electronic coupling parameter(DH ab ). The inner intramolecular reorganization energy (l i ) for selfexchange process has two contributions, arising from the geometry relaxation along inter-nuclear coordinate upon moving from neutral-state geometry to the charged-state geometry and vice versa; l i~l1 zl 2 ð6Þ l 1 and l 2 can be calculated directly from the adiabatic potential energy surface as; Where E 0 (G 0 ) and E -(G -) are the ground-state energies of the neutral and anionic states, respectively, and E 0 (G -) and E -(G 0 ) are the energy of the neutral molecule at the optimal anionic geometry and the energy of the anionic molecule at the optimal geometry of the neutral molecule.
The calculated reorganization energies for electron transport l 1 , l 2 and l i are shown in Table 2 From Table 2, it is clear that the value of l i decrease both for DMB derivatives with electron releasing (2-7) and electron withdrawing (7)(8)(9)(10)(11)(12)(13) substituents as compared to non-substituted DMB (1), with exception of compound 8 and 9, where it slightly increases. The calculated values of inner sphere reorganization parameter (l i )for different derivatives largely correlates well with the order of structural changes from neutral to anionic species, reflected by the change in C 18 -B bond length, upon electron injection. The order of change in reorganization energy is much more uniform in case of electron releasing substituents, as compared to electron withdrawing substituents. In case of electron withdrawing group the reorganization energies follow the trend: _CN ,_NO 2 ,_Br ,_COOH ,_Cl ,_F, with maximum for fluorine (0.214 eV) derivative which is, however, lower than to Alq 3 (0.276 eV) [35,51]. The smaller reorganisation energy value of DMB derivatives as compared Alq 3 will benefit the charge carrier transport. For the DMB with the electron releasing group, the reorganization energies has the following order: _CH 3 . Figure 4 shows a correlation of reorganization energies l with Hammett parameter. In case of electron withdrawing groups the trend is irregular, floro shows large geometrical changes as compared to strong electron withdrawing substituents, such as _NO 2 , _COOH, and _CN comparable to those of un-substituted DMB(1), whereas electron releasing substituents show clear trend; the stronger (_NH 2 ,_N(CH 3 ) 2 ) cause the DMB derivatives to undergo lesser geometry changes compared to weak electron releasing substituents.   Another important factor which characterizes the electron transport properties of such materials is their electron affinity. In an electron transport material, the electron affinity must be high enough to allow an efficient injection of electrons into low-lying unoccupied molecular orbitals. We calculated vertical and adiabatic electron affinities by using the equations 9 and 10 at B3PW91/6-311G (d,p)++ level of theory and the results are shown in Table 2.
The theoretically calculated AEA for the molecular oxygen at the B3LYP/6-31+G is 0.59 eV [52], while the experimental value, determined by photoelectron spectroscopy, is reported as 0.448¡0.006 eV [53]. Except compounds (2-5), EAs calculated for the whole set of derivatives are higher than the values reported for the molecular oxygen. So, these compounds in anionic form seem to be quite stable towards quenching caused by molecular oxygen present in the atmosphere.
Electrical properties of molecular systems depend on the energies of levels allowed for electrons/holes and energy gap between these levels. The energies of the HOMO and LUMO levels used for the calculation of the energy gap were calculated using DFT employing B3PW91/6-311++G(d, p) level performed for isolated molecules (Table 3). It is clear from Table 3 that the electron releasing substituents(2-7) increase the LUMO energy of the DMB with increasing order of Electron Transport and Nonlinear Optical Response their electron donating behavior and, therefore, electron affinity of these DMB derivatives decreases in that order ( Figure 5(a)). Whereas, the electron withdrawing groups (8-13) decrease the LUMO energy of the DMB with increasing order of their electron withdrawing character and, as a consequence, electron affinity of these DMB derivatives increases with increase in the electron withdrawing character of the substituents ( Figure 5(b)). A similar trend is observed for the HOMO energies and ionization potential. Simply this behavior is explained by build of enhanced effect of p-conjugation due to electron-releasing environment created by electron rich substituents and decrease of electron densities by electron withdrawing substituents. However, it is observed that with electron releasing groups as substituents (2-7), the changes in the HOMO energies are larger than those of LUMO energies, whereas with electron withdrawing groups as substituents, the changes in LUMO energies are greater than HOMO energies. One important observation is that the HOMO-LUMO gap, DE(H-L), of the studied compounds decrease as compared to un-substituted DMB and this gap goes on decreasing with increase in the electron donating as well as electron withdrawing character of the substituents. Therefore, it is possible to tailor the HOMO-LUMO gap in such electron transport materials by rational variation of substitutions on the basis their position on the Hammett parameter (s) scale. It was not possible to compare the calculated IP and EA with the experimental values due to lack of experimental data, however the plot of IP and EA correlated quite well with s, which points out an obvious substitution effect on energy levels in the studied systems ( Figure 6). The frontier molecular orbitals (FMOs) of all the systems are shown in Figure 7. It was observed that the highest  occupied molecular orbitals (HOMOs) of the neutral molecules delocalize primarily over the two mesityl rings, whereas the lowest unoccupied molecular orbitals (LUMOs) localize largely on the Boron (B) atom and the phenylene of the DMB moiety, with a small contribution from electron releasing or withdrawing substituent. Thus, mesityl rings constitute the donor state, while the phenylene   and B atom as the acceptors state in these systems. Compared to Alq3 [51] the LUMO is more delocalized than HOMO in case of DMB derivatives, indicating that the electron can easily move among the DMB molecule which will facilities enhanced electron transportation, particular in case of DMB derivatives with  electron donating substituents. The HOMO-LUMO electronic transitions can be attributed to charge transfer from the two mesityl rings to the phenylene and the vacant p z of Boron atom. Thus, the minimization of the barrier for electron injection and transport can be achieved by proper adjustment of EA of given derivative.
The theoretical data obtained for electron transport of DMB derivatives can be used to model an OLED of desired properties. We used a single OLED based on ITO(Indium tin oxide) as an anode, 9,9-dioctyfluorene as a hole transport layer, Ca-Ag as cathode and some studied TBA derivatives as an electron transport/ emitting layer(EL) (Figure 8). The first parameter that was considered in analysis was electron injection from Ca-Ag into EL. Efficiency of this process depends upon the energy barrier between the Ca-Ag work function and LUMO level of EL. Lower the barrier, the higher the injection efficiency. On comparing, we see that the DMB derivatives with electron donating substituents have lower energy barrier. Another parameter is efficiency of hole transfer from HTL to EL which is significant in case of derivatives 2 and 3, having least energy barrier between HOMO of ITO and the derivatives.

Nonlinear Response
According to Oudar and Chelma, first hyperpolarizability (b) is related to lowlying charge-transfer transition through a two level model [54] on the basis of the complex sum over states (SOS) expression as; where (m ee -m gg ) is the difference between the dipole moments of the mth excited state and ground state respectively, m ge is the transition dipole, and E ge is the transition energy.
In the two-level expression, the first hyperpolarizability (b), is proportional to the square of transition dipole and is inversely proportional to the square of the transition energy. For a molecule, the transition energy is therefore the decisive factor for large b values. According to two-level model materials with wellperforming NLO property should possess a low lying Charge Transfer excited states with good oscillator strength [55][56][57]. The calculated values of b, m ge and E ge employing TDDFT using B3LYP/6311++G (d, p) level of theory are given in Table 4. From the Table 4, it is clear that there is an enhancement of b values with increasing donor/acceptor abilities by substitutions as compared to DMB (1). Thus, on account of above equation the value of b increase with the increasing m ge 2 and decreasing E ge values. The results show that in case of electron releasing groups the observed trend for m ge 2 /E ge 2 values is 2.3.455.7.6 and is in agreement with the hyperpolarizability values. On comparing the -OCH 3 and -N(CH 3 ) 2 substituted molecules (5 and 7) with the methyl substituted molecule (6), it is found that compound 6 possess smaller m ge 2 value. Thus, the hyperpolarizability of the -N(CH 3 ) 2 substituted molecule is five times that of the methyl substituted molecules. Comparing the value of m ge 2 /E ge 2 of electron releasing groups with electron withdrawing groups, the former have higher values of m ge 2 /E ge 2 than the latter ones, in general. The value of m ge 2 /E ge 2 for methoxy substituted compounds is lower than the halogen derivatives (8,9,10), however the hyperpolarizability of the halogen substituted molecules is higher than the methoxy substituted molecules. The probable reason is large electronegativity of halogen substituents as compared to methoxy group result in higher dipole moment along the X-direction cause enhance in b value of halogen derivatives [53]. The transition energy of methoxy derivative (4.369 eV) was found to be higher than the halogen DMB derivatives. Thus, the contribution of the square form (m ge 2 /E ge 2 ) has major effect on hyperpolarizability in the two-level model. The values of b for the two extreme such as nitro (strongest electron withdrawing group) and dimethylamine (strong electron releasing group) are completely according to two model equation and depend upon the transition dipole and transition energy of the group. Comparing the bromo and carboxy derivatives of DMB, the former has higher value of m ge 2 /E ge 2 but the latter has higher value of polarizability. Thus, the value of polarizability is also effected by other parameters like size, electronegativity of substituents etc. The first electronic transition is primarily composed of the HOMO-LUMO transition. As the carboxy group is replaced by the bromo, the HOMO energy does not change significantly; however, the LUMO is lowered. The electron density of the HOMO orbital of carboxy derivative would favor the charge transfer to the LUMO orbital, confined toward the acceptor part of a molecule due to reduction in HOMO-LUMO gap. Similarly, on comparing the dimethylamine and amine derivates of DMB, it was observed that the HOMO of the dimethylamine-DMB rise and the electron density shift, resulting in intense charge transfer compared to amine derivative and gives high b value for dimethylamine DMB derivative. From this observation, we may conclude that the first transition energy is an important measure of the hyperpolarizability of such molecules.
We compare the selected geometry parameters with values of b of DMB derivatives. In molecules, with the releasing groups, the B-C 18 bond length was observed to be minimum for dimethylamine derivative and maximum for ethyl derivative and follows the order 2,3,455,6,7, and the B-C 2 and B-C 10 bond length slightly decrease from dimethylamine derivative to ethyl derivative. In case of electron withdrawing groups, all the three B-C 18 , B-C 2 and B-C 10 bond lengths decreases but the maximum decreases appears in B-C 18 bond [859. 10.11512.13], as substituted phenlyene ring is linked to boron through B-C 18 bond. Probable reason may be extended-p conjugation improved by electron releasing groups. The well-known resonance and inductive effect [58,59] may appropriately explain the relationship between donor-acceptor groups and the geometry. The increase in pconjugation due to +I (inductive) effect of electron releasing groups cause enhance in extent of resonance, results in geometrical modification of derivative on account of attachment of electron releasing substituent. The Hammett constant were correlated with hyperpolarizability. The plots of the b of DMB derivatives verses s in figure 9 reveals that hyperpolarizability of DMB derivatives with electron releasing group decrease with increasing s value. While as, in case of DMB derivatives with electron withdrawing groups, the b value increase with increase in Hammett value. The results fit a straight line with a correlation coefficient of 0.92 (or higher) for electron releasing groups, the other one has a correlation coefficient of 0.81 (or lower) for electron withdrawing groups. According to Ulman, inductive and resonance effects have contribution in Hammett substituent constant and both the effects are important to the molecules [58]. The donor group will easily push electrons towards the boron and remain itself slightly electron deficient. This situation affects the magnitude of dipole moments as well as hyperpolarizability. Thus, due to the presence of variety of functional groups (either having electron releasing or electron withdrawing abilities) affect the magnitude of static hyperpolarizability of DMBs.
To highlight the effect of change in bond length and bond angle upon substitution on hyperpolarizability, we plotted the dipole moment and hyperpolarizability of DMB derivatives versus change in key bond lengths and bond angle (see Figure S3 and S4. The observed plots revealed a close correlation between them for electron releasing groups. In DMB derivatives with electron releasing groups, the hyperpolarizability decrease as from 2.3.4.5.6.7  groups from methyl to dimethylamine have a pronounced impact on the hyperpolarisabilities. In general the trend is similar for other bond lengths for DMB derivatives with electron releasing derivatives. However, there is hardly a correlation between the hyperpolarisabilities and the change in geometrical parameters when the attached groups are electron withdrawing. In case of electron withdrawing groups, the DMB derivative with nitro group show maximum change in geometrical parameters and hyperpolarizability lower than that of fluorine DMB derivatives. We have correlated the Key geometrical parameter BLA (bond length alternation) calculated by the equation 5 with the hyperpolarizability. In case of DMB derivatives with electron donating groups, the hyperpolarizability decreases from 2.3.4.5.6.7 and the BLA values increases from 0.056(2) ,0.059(3) ,0.061(4) 50.061(5) ,0.063(6) ,0.064 (7). It shows that an increase in the strength of donor groups from methyl to dimethylamino have a clear impact on the b values. Thus, the minimization of the bond length alternation by varying the electron releasing groups is an important manifestation of change in hyperpolarisability. However, there is no good correlation between the hyperpolarizability and the BLA when the substituents are electron withdrawing. The DMB (12) and DMB(13) have lowest b values, with highest BLA values 0.083 and 0.085 respectively, but in case of 8, 9, 10 hyperpolarizability values decreases as 8.9.10, but the BLA values does not differ much (0.070(8) 50.070(9) ,0.071 (10). This indicates that the withdrawing groups have little effect on the BLA value and do not contribute to the change of hyperpolarisabilities. It can be concluded that changing the electron donor groups have a better effect on the BLA value and the hyperpolarizability in such electron transport systems.

Conclusion
In this work, density functional theory has been employed to analyze the electron transport and nonlinear optical properties of DMB derivatives. The geometrical parameters optimized at B3PW91/6-311++G (d, p) level are in agreement with the experimental ones. From the optimized geometries, we have calculated HOMO and LUMO energy levels, electron affinity, reorganization energy, frontiers molecular contours, polarizability and hyperpolarizability in addition to innersphere reorganization energies. The absorption energies have been calculated with time-dependent density-functional theory (TDDFT) at the optimized geometries. Replacing the hydrogen of DMB(1) ranging from -N(CH 3 ) 2 to NO 2 groups can cause change in extent of electron transport and nonlinear optical properties of DMB systems. In this study we have correlated the electron transport and nonlinear optical properties with different parameters including geometrical and Hammett parameters. There occurs an enhanced change in dihedral angle due to electron withdrawing substituents, reflecting increased tendency to change in radical-anion as compared to electron releasing groups. A strong correlation was observed between Hammett parameter and the geometrical parameters as the amplitude of the structural modifications increase on moving from electron releasing substituents to electron withdrawing substituents. The effect on reorganization energy was much more prominent in the case of electron withdrawing substituent, with the maximum for fluorine as compared to electron releasing substituents. It was observed that the highest occupied molecular orbitals (HOMOs) of the neutral molecules delocalize primarily over the two mesityl rings, whereas the lowest unoccupied molecular orbitals (LUMOs) localize largely on the B atom and the phenylene of the DMB moiety, with a small contribution from electron releasing or withdrawing substituent. Using two-level model, we observed that for determining the value of hyperpolarizability, the first transition energy is an important factor. In DMB derivatives with electron releasing groups, the hyperpolarizability decrease as from 26889.28610 230 esu.15924.7610 230 esu.10813.67610 230 esu.8321.161610 230 esu.3156.49610 230 esu.2512.611610 230 esu and the key BLA value (B-C 18 ) 0.0004,0.0005,0.000850.0008,0.0009, respectively. We hope that our theoretical study can give some insight to experimentalists for devising a wide variety of materials with extensive electron transport and nonlinear optical properties.