The role of π-linkers and electron acceptors in tuning the nonlinear optical properties of BODIPY-based zwitterionic molecules

Intramolecular charge transfer process can play a key role in developing strong nonlinear optical (NLO) response in a molecule for technological application. Herein, two series of boron dipyrromethene (BODIPY)-based push–pull systems have been designed with zwitterionic donor–acceptor groups, and their NLO properties have been evaluated using a density functional theory-based approach. Different π-conjugated linkers and electron acceptor groups were used to understand their roles in tuning the NLO properties. The molecules were analyzed through HOMO–LUMO gaps, frontier molecular orbitals, polarizabilities, hyperpolarizabilities, Δr indices, transition dipole moment densities, ionization potentials, electron affinities and reorganization energies for holes and electrons. These observations correlated well with the computed absorption spectra of the molecules. It is found that with the introduction of different π-linkers in the molecule, planarity is maintained and the HOMO–LUMO gap is systematically decreased, which leads to a large NLO response. It was noted that the electronic absorption wavelength maxima were found in the near-infrared region (934–1650 nm). The results show that compared to the pyridinium acceptor group, the imidazolium acceptor group in the BODIPY systems amplifies the NLO response to a larger extent. It is also observed that the BODIPY-based dye with an imidazolium acceptor and thienothiophene π-linker shows the highest first hyperpolarizability value of 3194 × 10−30 esu. Furthermore, the charge transfer occurs in the z-direction, as the z-component of the first hyperpolarizability is the dominant factor in this system. Here, the designed molecules show a characteristic reorganisation energy value, which is a deciding factor in the rate of hole/electron transport for favourable intermolecular coupling. As a whole, this theoretical work highlights that π-conjugated linkers and electron acceptor groups can be used judiciously to design new molecular systems for optoelectronic applications.


Introduction
Materials with high nonlinear optical (NLO) properties have gained much attention in recent years, as they offer potential applications in optoelectronic and photonic devices in the regulation of optical switching, micro fabrication and imaging, telecommunication, laser technology, data storage, etc. [1][2][3][4][5][6][7][8][9] Different types of materials of both inorganic 10 and organic [11][12][13] origin have been investigated for their interesting nonlinear optical properties. However, among these, organic materials with high NLO response have proven to be the optimum choice due to their tunability. 7,14 Most organic NLO materials are push-pull molecules containing p-conjugated systems linking an electron donor (D) to an electron acceptor (A) group. 15 The NLO properties in these compounds are developed due to high polarization of p electrons along the conjugated backbone. The rst hyperpolarizability (b) is associated with the intramolecular charge transfer (ICT) from the donor group to the electron acceptor group. Several studies conrmed that by optimizing the donor and acceptor groups, p-linker and ring twisting in the push-pull system, the ICT process can be amplied and hence the NLO response increases. [16][17][18] The main strategies for designing efficient NLO materials are based on the choice of suitable donor and acceptor groups and p-conjugated bridges. 19,20 Albert et al. suggested an approach in which by introducing zwitterionic behaviour in a conjugated molecule, one can enhance the NLO properties by providing low energy charge transfer. 21,22 Some researchers also revealed that zwitterionic D and A groups positioned at opposite ends of the conjugated systems of chromophores can ameliorate their NLO response. 12, 23 Xiong et al. synthesized a thermally stable zwitterionic picolinium(dicyano) esterquinodimethane chromophore with a b value of 1800 Â 10 À30 esu. 24 All these facts prompted us to design and investigate zwitterionic systems craed by different donors and acceptors with high NLO response.
The boron dipyrromethene (BODIPY)-based push-pull chromophore is an area of recent interest and has been used in the present investigation. BODIPY dyes are well known for their desirable photophysical properties, as they possess strong absorption bands in the UV-vis region, high uorescence quantum yields, 25,26 two-photon absorption properties, etc. 27 BODIPY dyes are widely used in photodynamic therapy, 28 as chemosensors, 29 uorosensors, 30 etc. Extensive research has been performed on the molecular structures and optical properties of BODIPY dyes; however, very few studies have reported BODIPY-based push-pull zwitterionic systems. In this work, the main strategy is the introduction of donor and acceptor groups to the opposite ends of the BODIPY core at the 2 and 6 positions. Imidazolium compounds are known for their signicant absorption in the entire UV region and for their uorescence properties. 31,32 It is reported that the imidazolium cation provides a remarkable contribution to enhanced NLO response. 33,34 To investigate the effects of different electron acceptors and p-conjugated linkers on the NLO properties of the BODIPY system, 10 new BODIPY-based D-p-A dyes were designed (Fig. 1). Pyridinium and imidazolium ions were taken as positively charged acceptor groups, and phenoxide group was taken as the negatively charged donor group. Four p-spacers, ethylene, cis-thiophene with respect to BODIPY, trans-thiophene with respect to BODIPY and thienothiophene, were used as the 1 st and 3 rd p-linkers, and BODIPY was used as the 2 nd p-linker.
This paper deals primarily with the theoretical calculation of the optoelectronic properties of newly designed molecules using density functional theory (DFT)-based methods. The focus of our work was to determine the inuence of different pconjugated linkers and p-acceptors on the NLO properties of molecules. DFT and TDDFT calculations were carried out to evaluate the molecular properties, including the frontier molecular orbitals, absorption spectra, polarizability and hyperpolarizability, ionization potential (IP), electron affinity (EA) and reorganization energy of the hole and electron. It is exhibited that with increasing p-electron delocalization and substitution of pyridinium cation by imidazolium cation in the BODIPY-based systems, the NLO response is amplied. The present work will provide a direction to researchers for the synthesis of novel NLO materials of BODIPY-based dyes for technological application.

Theoretical background and computational details
Geometry optimization of the designed molecules was performed in the gas phase using the Gaussian 09 program at the B3LYP/6-311++G(d,p) level of theory. 35 The molecules were optimized at the minimum energy level and no imaginary frequency was found, which signies the stationary point of the minima. All the obtained restricted B3LYP solutions were stable. The B3LYP and CAM-B3LYP functionals were found to be the most suitable for NLO calculations. Recently, Prakasam et al. investigated the second order hyperpolarizabilities and absorption properties of triphenylamine-based organic sensitizers with the B3LYP/6-311++G(d,p) method. 36 However, the B3LYP functional sometimes fails to predict accurate values in the case of polarizability of long chain molecules, charge transfer excitation, etc. Increasing the incorporated HF exchange fraction as well as the long range correction becomes signicant to improve the results, and many researchers have attempted to achieve this. Very recently, Avramopoulos and coworkers designed a photochromic material with switchable nonlinear optical properties by employing the density functional method with the CAM-B3LYP functional. 37 For this reason, we considered the B3LYP and CAM-B3LYP functionals for explanation of the NLO behaviour of the studied system. Hence, in order to analyze the NLO behaviour of the designed D-p-A systems, the dipole moment (m total ), average polarizability (a) and rst hyperpolarizability (b) were evaluated using the B3LYP functional as well as the CAM-B3LYP 38 functional with the 6-311++G(d,p) basis set to include a long range correction for improving the charge transfer characteristics in the Gaussian 09 soware. 39 The dipole moment was calculated using the following equation: where m x , m y , and m z are the components of the dipole moments in the x, y, and z directions, respectively. Average polarizability (hai) was calculated using the following equation: 40 where a xx , a yy and a zz are the polarizability tensor components.
The rst hyperpolarizability (b total ) can be expressed as 40 and b yyy are hyperpolarizability tensors along the x, y and z directions, respectively. The transition dipole moment (TDM) density has been evaluated on the basis of the dipole moment integral between the occupied and virtual molecular orbitals. The TDM density can be calculated and visualized using the Multiwfn wave function analyzer. 41 The charge transfer length during electron excitation in the D-p-A system is dened as the Dr index. This is dened by the following equation: 42 where the i and l indices run over all the occupied and vertical MOs, respectively, and F is the orbital wave function: where X i l and Y i l denote the conguration coefficients corresponding to i/l excitation and l/i de-excitation, respectively. The Dr indices were calculated using Multiwfn soware. The ionization potential (IP) denes the energy changes due to removing electrons from or adding holes to the neutral molecule. On the other hand, the electron affinity (EA) quan-ties the energy changes for adding electrons to or removing holes from the same. This can be expressed as: 43 E denotes the energy of the respective system. The hole transport reorganization energy (l hole ) and electron transport reorganization energy (l electron ) are calculated as: 44 Here, E 0 , E + and E À represent the energies of the neutral, cation and anion species, respectively, in their optimized geometries. E 0 +/À and E +/À 0 represent the energies of the cation/ anion with the optimized structure of the neutral species and the energy of the neutral species with the optimized structure of the cation/anion geometry, respectively.

Results and discussion
Electronic structure The electronic structures of molecules play a key role in regulating the nonlinear optical properties of the molecules, as they depend on the geometry of the molecular system. Here, Fig. 1 presents the BODIPY molecules with different p-linkers and acceptors and a phenoxide donor. The optimized structures of the molecules, dihedral angles and bond angles are given in the ESI (Fig. S1-S10 and Tables S1-S10 †). For a smooth electron transfer process from the donor unit to the acceptor unit, the dye must have co-planarity between the donor unit-bridging unit-acceptor unit (D-p-A). 45 It is found that most of the dyes maintain co-planarity, as the dihedral angles are close to 0 and 180 and the bond angles are close to 120 . It was also noted that dye 2d has perfect planar geometry, as the geometrical parameters are very close to ideal. Thus, we can say that the inclusion of different p-linkers favours planarity in the investigated systems. It is reported that upon attaching an unsubstituted phenyl ring to the BODIPY core, the uorescence quantum yield (4 f ) decreases compared to that of BODIPY alone. Meanwhile, incorporation of two methyl groups in the BODIPY system restricts the motion of the phenyl ring and increases the 4 f of the molecule. 46 Therefore, we can expect that incorporation of proper p-linkers will increase the NLO response.
The charge transfer process is highly dependent on the energy difference between the HOMO and LUMO (DE). Frontier Molecular Orbital (FMO) theory helps to predict the chemical stability of a molecule. 47,48 Usually, the LUMO denes the capacity of acceptance of electrons and the HOMO classies the electron donating ability. 49 A low DE value indicates a chemically so molecule, whereas a high DE value signies a chemically hard molecule. Polarizability can be enhanced with increasing soness in a molecule, which facilitates the NLO response.
The FMO energy levels along with the HOMO-LUMO gaps are represented in Fig. 2. Proper inclusion of acceptor and donor units at the 2 and 6 positions of BODIPY endow the molecule with zwitterionic character, and the HOMO-LUMO energy gap is found to be reduced. Currently, the compounds which are studied can be categorized into two sets, namely, BODIPY-based pyridinium acceptors (Series 1) and imidazolium acceptors (Series 2). Upon introduction of different p-linkers in the BODIPY framework in series 1 and 2, DE starts to diminish from the respective values of 0.88 eV and 1.08 eV in a very systematic way. Upon insertion of -CH]CH-(1a and 2a), cisthiophene (1b and 2b) and trans-thiophene (1c and 2c) groups in the BODIPY system, DE decreases successively. However, when thienothiophene groups are inserted (dye 1d and 2d), the lowest DE values of 0.45 eV and 0.50 eV are found. This clearly reveals that the thienothiophene p-linkers reduce the HOMO-LUMO gap more in the imidazolium acceptor system than in the pyridinium acceptor system. It is expected that by maintaining planarity, extended conjugation and a low DE gap, the molecules should achieve high NLO response. For some cases, the energy differences between the HOMO and LUMO of the designed dyes are very small. If we go along a series, the architecture of the molecule remains the same; only the conjugation length is increased, and as a consequence, the LUMO is stabilized. This could be the reason for these small differences. They also may be due to the limitations of the method used. Zhang et al. synthesized and performed DFTbased B3LYP calculations on a series of highly polarizable bis(N,N-diethyl) aniline-based chromophores. In this study, the reported energy differences of the compounds are also very low. 50 Therefore, these values can be taken as an explanation of the NLO behaviour of the investigated compounds. The orders of the HOMO-LUMO energy gaps of series 1 and 2 are 1d < 1c < 1b < 1a < 1 and 2d < 2c < 2b < 2a < 2. Therefore, it can be predicted that all these molecules will show higher absorption wavelengths and will be excellent candidates for NLO materials for various photoelectronic applications. The modication of the p-conjugated linker by tailoring the structure would be a strategy to obtain high NLO activity. In order to understand the distribution patterns of the HOMOs and LUMOs, we analyzed the frontier molecular orbitals of the investigated dyes, which are given in Fig. 3. From the FMO diagram, it is evident that for Series 1, the HOMOs and LUMOs of all the chromophores are delocalized throughout the molecules. Meanwhile, for Series 2, the HOMO is localized on the phenoxide donor part and p-conjugated part of the molecule, and the LUMO is localized mainly on the imidazolium acceptor part. This fact signies the electron donating ability of phenoxide ion and the electron acceptance nature of imidazolium cation. Overall, the calculated geometrical parameters indicate that planarity and delocalization in a molecule can be adjusted by varying the p-linkers, which are signicant characteristics for tuning optoelectronic properties.

NLO properties
To evaluate the NLO properties of the dyes, it was necessary to assess the magnitudes of the static dipole moment, linear response (polarizability) and nonlinear response (1 st hyperpolarizability). Benchmarking of the computational model for the NLO calculation was performed with some well-known push-pull systems given in the ESI 51 (Table S14 †). A comparison between the theoretically calculated results and experimental results of the hyperpolarizability show good agreement. The benchmarking exercise demonstrates that the hyperpolarizability value for the This journal is © The Royal Society of Chemistry 2020 RSC Adv., 2020, 10, 40300-40309 | 40303 B3LYP functional is closer to the experimental value compared to that for the CAM-B3LYP functional. Here, our main objective was to nd the effects of the zwitterionic donor-acceptor group and p-linkers on the hyperpolarizability value. Therefore, to obtain an explanation for the changes in the hyperpolarizability values with the structural modication in the framework of the donor-acceptor molecule in the gas phase, we kept the other tuning parameters unchanged. Hence, to determine the inuence of the p-conjugated linkers and electron acceptor groups on the NLO properties of the Series 1 and Series 2 compounds, these parameters were calculated with the B3LYP and rangeseparated CAM-B3LYP functionals with a standard basis set, 6-311++G(d,p), and the results are tabulated in Tables 1 and S12 (ESI †). In the following discussion, we illustrate the results obtained using the B3LYP functional. The average polarizability (Da) of the dyes increases in the order of a < b < c < d for both series, and the results found using the B3LYP and CAM-B3LYP functionals follow the same trend. Looking at Table 1, one can nd that the average polarizabilities of dye 1 and 2 are 103.37 Â 10 À24 esu and 83.72 Â 10 À24 esu, respectively. With the addition of ethylene p-linkers as the 1 st and 3 rd conjugators in the 1 and 2 dyes, Da increases by 61 Â 10 À24 esu and 65 Â 10 À24 esu, respectively. The polarizability further increases with the addition of the thiophene p-linker at the cis and trans positions of BODIPY. It was found that the average polarizabilities reached maxima and the values increased by 242 Â 10 À24 esu and 238 Â 10 À24 esu, respectively, when thienothiophene groups were incorporated in the parent compounds (Table 1). Therefore, Tables 1 and S12 † indicate that with modication by different p-conjugated linkers, the polarizability is enhanced systematically. The highest polarizability was found for dye 1d, 345.52 Â 10 À24 esu, with the smallest DE (Table 1). It is known that a small energy gap between the HOMO and LUMO inuences the polarizability of a molecule, and the results are in accordance with the calculated HOMO-LUMO gaps.
The rst hyperpolarizability is connected with the ICT process. This is due to the ow of electron density from D to A via a p-conjugated bridge, which is discussed in more detail later. Among the investigated molecules, the dipole moment was found to be highest for dye 2d, while the lowest value was found for dye 1. It is reported that a large dipole moment and hyperpolarizability play signicant roles in poled polymers. 52 Therefore, the high dipole moment of these molecules is likely to enhance the NLO response. The rst hyperpolarizability is graphically represented in Fig. 4. Calculation of the b total value using the B3LYP and CAM-B3LYP functionals follows the same trend of increasing order. An increase in the b total value of 557 Â 10 À30 esu is observed in the case of the imidazolium acceptor group (dye 2) compared to the pyridinium acceptor (dye 1), which proves the better electron acceptance power of imidazolium cation. 34 This is also found in the frontier molecular orbitals. Changes in the p-conjugation length due to the inclusion of the 1 st and 3 rd p-linkers affect the NLO properties of the studied systems. For Series 1 and Series 2, with the incorporation of the ethylene p-linker, the b total values are almost doubled. It was found that with the addition of thiolinkers in both series, the 1 st hyperpolarizability increased systematically. Compared to the cis-thiophene p-linkers, the trans-thiophene p-linkers improved the b total value by 60 Â 10 À30 esu in dye 1c and 24 Â 10 À30 esu in dye 2c. When  thienothiophene molecules were added to dyes 1 and 2, the b total value increased by 531 Â 10 À30 esu and 2293 Â 10 À30 esu, respectively. This can be correlated with the HOMO-LUMO gaps of the molecules, as shown in Fig. 2. It is known that the 1 st order hyperpolarizability increases with decreasing DE. 53,54 If we compare Fig. 2, 4 and Table 1, it can be found that for all cases, the DE gap is inversely proportional to the dipole moment, linear polarizability and 1 st order hyperpolarizability. The b total values of the two series decrease in the orders of 1d > 1c > 1b > 1a > 1 and 2d > 2c > 2b > 2a > 2 for the B3LYP as well as the CAM-B3LYP functional. It was observed that BODIPY-based dye 2d, which contains an imidazolium acceptor and thienothiophene p-linker, shows the highest b total and mb value of 3194 Â 10 À30 esu and 18.29 Â 10 À44 esu, respectively. This NLO response can be attributed to the effective charge transfer from the D to the A moiety. The computed b total value of 2d is 4110 times higher than the 1 st hyperpolarizability value of urea, which is used as a reference in organic systems. It should also be noted that in both series, the order of the b total value is in conformity with the average polarizability. For better insight, we report the components of the 1 st hyperpolarizability in Tables 2 and S13 † (ESI †).
Yu et al. argued that b y controls the b total of 2,3-naphto-15crown-5 ether for some metal cation complexes. 55 They conclude that during the polarization process, the charge ows in the y direction. The results in Tables 2 and S13 † indicate that for the investigated systems, b total is dominated by the b z component. The contribution of b x is negligible. Therefore, we can surmise that during the polarization process, the maximum charge is expected to transfer along the z direction. 56 To explain the origin and variation of the major transition energy of the NLO response, the Dr index, integral overlap of hole electron distribution (s) and distance between the centroid of the hole and electron (D) were tabulated (Table 3). From Table 3, we can nd that for all the cases, the Dr index is larger than 2.0; hence, they may be regarded as charge transfer (CT) mode. 42 If we compare Fig. 4 and Table 3, we can nd that both the b total and Dr indices increase, with the p-conjugated linkers showing the same trend of increment. The 2d compound shows the highest Dr indices and was found to show the highest b total values. If we compare these two series, it is exhibited that for   both series, the thienothiophene p-conjugated linkers show the highest Dr indices, and the imidazolium acceptor system proves to be better for charge transfer. This can be well correlated with the 1 st hyperpolarizability of the molecules. The integral overlap of the hole electron distribution (s) is a tool to measure the spatial separation of holes and electrons and the distance between the centroids of holes and electrons (D), which dene the charge transfer length. The higher the value of D, the greater the charge transfer length. For CT transition mode, D should be large and s should be small. The smallest s value of 0.015 and highest D value of 21.04 A were found for dye 2d, which elucidates the genesis of the highest NLO response. The transition dipole moment densities are plotted in Fig. 5.

UV-vis spectra of dyes
Time-dependent density functional theory (TDDFT) was used to calculate the vertical excitation energies of the designed dyes.
TDDFT computations were carried out using the B3LYP and CAMB3LYP functionals with the 6-311++G(d,p) basis set. The TDDFT method is a very popular choice for the calculation of vertical excitation energy. However, it has some limitations when predicting the transition energies for charge transferbased molecules. 57 Jacquemin and co-workers performed computational analysis of the spectral properties of BODIPYderived systems with different functionals using the TDDFT methodology 58 and compared the results with experimental data. In this paper, they also designed new aza-BODIPY molecules for near-infra red (NIR) applications with the same methodology. For the theoretical calculation of the vertical excitation energy of BODIPY-based molecules, TDDFT was found to be a suitable choice. Therefore, here, we used TDDFT to calculate the absorption spectra of the designed dyes absorbing light in the NIR region. A long-range correlated functional such as CAM-B3LYP is a better choice to calculate the vertical excitation energy to understand the charge transfer characteristics in molecules. 59,60 Benchmarking of the TDDFT methodology was performed with some common BODIPY systems, which are represented in Table 4. 25 The correlation of computed and experimental results induced us to choose the CAM-B3LYP functional for the estimation of the vertical excitation energy of our designed systems. The vertical excitation energies, maximum absorption wavelengths, oscillator strengths and nature of the transitions are reported in Tables 5 and S11 (ESI †). When we compare Fig. 2, Tables 5 and S11, † we can nd that the excitation energies in the case of the CAM-B3LYP functional are in consonance with the HOMO-LUMO gaps of the molecules. For most of the cases, the differences found in the excitation energy calculations with the B3LYP and CAM-B3LYP functionals varied by about $0.2 eV; this is due to the different functionals of DFT. In the following discussions, the results are illustrated by taking the CAM-B3LYP functional. In all cases, the S 0 / S 1 electronic transitions dominated with high oscillator strength, and for dyes 1b, 1c and 1d, the S 0 / S 2 electronic transitions were found to have moderate oscillator strength. The excitation energies of dye 1 and dye 2 are 1.173 eV and 1.327 eV, respectively. From these results, it was found that upon substitution of the pyridinium acceptor by the imidazolium acceptor in the BODIPY system, the absorption wavelengths blue shied due to the higher HOMO-LUMO gap compared to that of the pyridinium acceptor. Table 5 indicates that when an ethylene group is included in the D-p-A system in dyes 1 and 2, the absorption maxima are red shied by 29 nm and 88 nm, respectively. This is due to the increases of the p-conjugation length of the system. The excitation energy systematically decreases with addition of thiophene p-linkers. Compared to cis-thiophene, the trans-thiophene systems (1c and 2c) are red shied by 128 nm and 143 nm. The maximum absorption wavelengths were found for dyes 1d (1650 nm) and 2d (1471 nm) when thienothiophene groups were added as the 1 st and 3 rd plinker in the BODIPY donor-acceptor system, respectively. The structural modication in the BODIPY system pushes the chromophore to absorb light in the NIR region (933 nm-1471 nm). In both series, the absorption wavelengths red shied with the addition of p-linkers to the parent compounds. The increasing orders of l max are 1 < 1a < 1b < 1c < 1d and 2 < 2a < 2b < 2c < 2d. The highest oscillator strength was found for dye 1d, which is 1.4957. Overall, it can be concluded that with the addition of different p-linkers to the D-p-A system, the l max , oscillator strength, and ICT character can be systematically modied.

Ionization potential, electron affinity and reorganization energy
Isolated molecules with low reorganization energy are associated with high solid state charge carrier mobility. 61,62 It was also found that some N-containing conjugated organic compounds are promising n-type or p-type materials for the fabrication of optical light emitting diodes (OLEDs). 63,64 This arises due to their excellent optoelectronic properties, thermal stabilities, high electron mobilities, etc. As hyperpolarizability is associated with the charge transfer from the donor group to the acceptor group, determination of the ionization potential (IP), electron  Table 6. For all the molecules, the potential energy of the cationic state is higher than the potential energy of the neutral state; this provides a positive ionization potential. In the case of an anion, the potential energy is lower than the potential energy of the neutral state, which gives a negative EA. The IP and EA values of the studied dyes are in the ranges of 5.37 eV to 4.47 eV and 2.04 eV to 3.05 eV, respectively. It is known that molecules with higher EA values show higher electron transport abilities and molecules with lower IP values show higher hole transport ability. 65 Compound 2d has the lowest IP value of 4.47 eV, and 1d has the highest EA value of 3.05 eV. Therefore, molecules 1d and 2d will transport electrons and holes better than the other molecules. It is interesting to note that with the decrease of the HOMO-LUMO gap in both systems, the IP values gradually decrease and the EA values increase.
To examine the effects of various substitutions on the reorganization energy for hole/electron transport (l hole/electron ), l hole and l electron are tabulated in Table 6. To obtain high efficiency optical light-emitting diodes, it is exigent to optimize the charge carrier (electron/hole) injection process and the electron-hole recombination process. 66 Also, the charge injection and charge transport are crucial to ascertain the luminescence efficiency of OLEDs, and these can be obtained using a separate electron transport material and hole transport material. 67 According to Marcus theory, the reorganization energy is the deciding factor in the rate of charge transfer, and this energy is required for  structural changes associated with the charge transfer process. 68 It is well known that the smaller the hole transport reorganization energy (l hole ), the better the hole transport properties, and the smaller the electron transport reorganization energy (l electron ), the better the electron transport properties. As expected, with increasing p-conjugation length, the l hole and l electron values decreased for most of the cases. This implies that the hole/electron transport properties increase along the series. The smallest l hole and l electron values were found for the thienothiophene-substituted linker in both series, which is in accordance with the NLO properties of the molecules. This can be explained by the increased rigidity and planarity of the molecules. The smallest l hole and l electron values were found for dye 2d and dye 1d, which can be correlated with the obtained IP and EA values of the studied molecules. Most interestingly, it was found that for all the cases, l hole is smaller than l electron , which indicates that the molecules are better hole transporters than electron transporters. The differences between l hole and l electron of the compounds with p-spacers are less than 0.15 eV, which suggests that they can act as emitters with moderately high light-emitting efficiencies. 69 Thus, the studied molecules are promising candidates as hole transport materials.

Conclusions
In the present computational study, two series of novel zwitterionic BODIPY-based molecules with pyridinium and imidazolium electron acceptors, a phenoxide unit as the electron donor and various p-linkers were designed, and their nonlinear optical properties were investigated. DFT-based methods were employed to explore the electronic structures, dipole moments, polarizabilities, hyperpolarizabilities, absorption properties, IPs, EAs and reorganization energies of these molecules. Our results reveal that with the introduction of different p-linkers in the D-p-A systems, planarity was maintained and the HOMO-LUMO gap systematically decreased. It was observed that thienothiophene groups reduced DE more than the other p-linkers. This observation was also supported by the excitation energy calculated with the TDDFT approach. Small HOMO-LUMO gaps of the molecules prompted the systems to absorb light in the NIR region. The FMO analysis depicts that for the Series 1 dyes, the electrons are delocalized throughout the molecules, whereas for the Series 2 dyes, the donor unit largely stabilizes the HOMO along with the p-linker and the LUMO is localized on the acceptor unit. The intramolecular charge transfer from the phenoxide unit to the imidazolium unit through the pconjugated backbone plays a signicant role in obtaining large NLO responses in the Series 2 dyes. The computed b total values of the dyes were found to be 442 (for dye 1) to 4110 (for dye 2d) times greater than the value of urea molecule. Among the studied molecules, 2d exhibited the highest b total value. The imidazolium acceptor showed better NLO response than the pyridinium acceptor. It was also noted that the dipole moment, polarizability and hyperpolarizability are consistent with the HOMO-LUMO gaps of the molecules. The lowest electron transport reorganization energy (0.17 eV) and highest EA were found for dye 1d, and the smallest hole transport reorganization energy (0.06 eV) and lowest IP were found for dye 2d. Consequently, 1d and 2d appear to be highly efficient emitters with promising hole transport characteristics. Overall, all the investigated dyes show high NLO response. The rst hyperpolarizability of the molecules was found to respond dominantly in the z direction, which indicates the course of charge transfer. A good correlation was found between the Dr index, TDM density, hyperpolarizability value and hole/electron transport properties. As a whole, this work demonstrates that the structural modication of p-linkers and electron acceptors in designing D-p-A systems is a signicant approach to obtain high-performance NLO materials.

Conflicts of interest
There are no conicts to declare.