Fabrication and Characterization of Polypyrrole/Multi-Walled Carbon Nanotubes Thin Films Using Thermal Evaporation

Polypyrrole/multiwalled carbon nanotubes composites (PPy/MWCNTs) were produced in an acidic solution utilizing an in situ oxidative polymerization method using ferric chloride as an oxidizing agent and sodium dodecyl sulfate as a soft template. Thermal evaporation was used to fabricate thin films from polypyrrole/multiwalled carbon nanotube composites. The resulting composites were examined by different techniques to explore their morphology, structural and electrical characteristics. The surface morphology analysis revealed that polypyrrole structure is a two-dimensional film with impeded nanoparticles and the thickness of coated PPy around the MWCNTs decreases when increasing the amount of MWCNTs. XRD analysis revealed that the average crystallite size of the prepared composites is 62.26 nm. The direct energy gap for PPy is affected by a factor ranging from 2.41 eV to 1.47 eV depending on the contents of MWCNTs. The thin film’s optical properties were examined using experimental and TDDFT-DFT/DMOl3 simulation techniques. The optical constants and optical conductivity of the composites were calculated and correlated. The structural and optical characteristics of the simulated nanocomposites as single isolated molecules accord well with the experimental results. The nanocomposite thin films demonstrated promising results, making them a viable candidate for polymer solar cell demands. Under optimal circumstances, the constructed planar heterojunction solar cells with a 75 ± 3 nm layer of PPy/MWCNTs had a power conversion efficiency (PCE) of 6.86%.


Introduction
Conducting polymers (CPs) have recently attracted the interest of researchers due to their potential optoelectronic characteristics. Conjugated polymers with high absorption coefficients, high electrochemical activity, high conductivity, and strong chemical stability are known as CPs [1]. They have good optical and electrical qualities due to their unique physical and chemical structures and are utilized in a variety of applications such as electrochromic devices, emitting diodes, light photovoltaics, microwave shielding,

Fabrication of the Thin Films
Physical vapor deposition (PVD) was used to produce the thin films. The thin films were formed at an initial pressure of 5 × 10 −5 mbar onto an ITO/glass substrate and/or a single crystal of wafer, with inter-digitized electrodes spaced by 75 m. A quartz crystalline micro-balance with UNIVEX 250 Leybold, two tantalum boats at any location without vacuum breakdown, and a deposition rate of 3 Å/s are used to create the film thickness of approximately 200 nm [31,32]. Figure S2 presents the steps involved in the fabrication of the films.

Computational Study
According to DFT calculations, simulations utilizing CASTEP and DMol3 provided findings for the efficiency of molecular structure and frequency dimensions for [PPy] and [PPy/MWCNTs] as isolated states. The applications of DMol3 and CASTEP software were computed based on free molecules, GGA functional correlations, Perdew-Burke-Ernzerh exchange, the pseudo-conserving norm, and the DNP base set [33,34]. The plane-wave cutoff energy value for the structural matrix simulation computations was 220 eV. To demonstrate the structural and spectroscopic properties of isolated [PPy] and [PPy/MWCNTs] utilizing DMol3 and CASTEP frequency calculations at the GP, including an XRD system and optical characteristics. Functional Becke's non-local-based interchange correlation with functional B3LYP [35] and WBX97XD/6-311G were performed on doped [PPy] and [PPy/MWCNTs] as isolate state nanocomposites for optimized geometry and vibrational frequency (IR) measurements. Nanocomposite materials are evaluated for geometric parameters, vibration modes, optimal structure visualization, and energies using the GAUSSIAN 09W programmer (Carnegie Mellon University Gaussian, Inc). It has been discovered that DFT calculations are dependent on WBX97XD/6-311 G using the B3LYP method, which has resulted in numerous excellent findings for structural spectrum correlation, including some major experimental discoveries. The Gaussian Potential Approximation System (GAP) employs a range of descriptors, as well as the overall power and derivatives model, as well as the concurrent use of several unique uncertainty models, as well as the Gaussian and CASTEP models in the gas phases, to measure [PPy] and [PPy/MWCNTs].

Characterization
Several techniques were used to characterize the thin films, specifications are listed in Table 1.  Figure S3 shows FTIR spectra for PPY, MWCNTs, and/MWCNTs/PPy nanocomposites [12]. The 1550 and 1450 cm −1 bands are shown on this graph. These correspond to PPy's stretching pulses C=C and C-N. At 1170 cm −1 , the pyrrole ring's vibration is generated. C-H and N-H vibrated in-plane at 1040 cm −1 and out-of-plane at 860 cm −1 . The 1130 cm −1 band has a PPy (chloride ion doped) characteristic. The anticipated S=O expanded vibration at 1183 cm −1 could not be identified due to overlap with the pyrrole ring vibration at 1170 cm −1 . The FTIR spectrum of the PPy/MWCNT composites with PPy are comparable to that of clean MWCNTs and PPy. This validates the presence of all PPy and MWCNT characteristic bands in composites. The strength of the C-H band varied dramatically as the proportion of MWCNTs increased.

Surface Morphology
The scanning electron microscope (SEM) technique was used to examine the surface morphology of the synthesized neat PPy and PPy/MWCNTs composites. Figure 1a-c showed SEM images for PPy/MWCNTs composites. From Figure 1a, it can be noticed that the obtained polypyrrole structure with impeded carbon nanotubes. The shape and size of the obtained polypyrrole usually depend on the type of used surfactant during the polymerization process. The formation of PPy sheet/films has been synthesized in the presence of SDS surfactant using FeCl 3 as an oxidizing agent [36]. This difference depends on the synthesis conditions such as the variations in the molar ratios of the pyrrole and surfactants/oxidants, which might affect the formation of bonds between the resulting PPy chains. Figure 1b

X-ray Diffraction (XRD) Analysis
The XRD patterns indicate that the MWCNTs are embedded into the PPy nanocomposite. Furthermore, no significant peak shifts were observed in the XRD patterns. Figure  2 presents the pattern XRD for the experimental and simulated for PPy/MWCNTs composite. When spinning doped, the MWCNTs can be incorporated into the PPy crystal lattice. The XRD pattern obtained from the fabricated PPy/MWCNTs composites thin films. They were correlated to the isolated system matrix. The predicted crystallite size (D) and miller index (hkl) depend on the absolute values of full width at half maximum (FWHM) as given in Table 2. The data in database code amcsd 0020475 corresponds well with the interplanar distances d spacing [37]. The TDDFT-DFT and Crystal Sleuth Microsoft programs were used to designate peak lines calculated by diffraction that were close to the observed findings. The Debye-Scherrer formula was applied to assessed XRD for [PPy/MWCNTs] 0-3 thin films, the range of 5≤ 2θ ≤ 45 with 1 ℎ ⁄ = 0.06Å − 0.74Å , = 1.54 Å, / = 0.5, polarization = 0.5, and Pesedo-Voigt function. From Scherer's formula: where λ is the X-ray wavelength (1.54 Å). As presented in Table 2, for the fabricated [PPy/MWCNTs] composites thin films XRD data from the XRD pattern was used to examine factors and features such as FWHM, the crystallite size (D), hkl indices, d-spacing (d), and peak intensity. The crystallite size is = 62.26 was within the range of 32.49-144.87 nm. Polymorph calculated the theoretical X-ray diffraction models using content studio software (version 7.0) (See Figure 2). Inset Figure 2, the integrals were conducted on the Brillouin zone with 2 × 2 × 1 (Polymorph PPy/MWCNTs composites). For the corresponding experimental, a comparison was made between experimental X-rays structures and measured PXRD patterns for PPy/MWCNTs composites thin films. While both experimental and PXRD models the intensity and location of specific peaks vary only slightly, the emphasis here is mostly on their overall resemblance. Only the important comparison characteristics between the measured and the experimental data should therefore be evaluated. It is also known that instrumentation and data collection processes are only two of the many variables that can affect the experimental PXRD pattern. The simulated PPy/MWCNTs composites as isolated position in polycrystalline and provide a triclinic in the group P1 . For PPy/MWCNTs composites thin films, the main peaks at hkl (11 1), hkl (021), hkl (004), and hkl (1 1 2) at 2θ values of 17.54°, 18.30°, 20.04°, and 20.38°, respectively. A full assessment demonstrating a good agreement between the calculated PXRD patterns and the experimental patterns for PPy/MWCNTs composites thin films, validating the accuracy of the fabricated material's PXRD patterns. As presented in Figure  2, a combination of experimentally based diffraction and density functional theory calculations yields a great estimation of the atomic scale for PPy/MWCNTs composites thin films (2θ at hkl (11 1) .

X-ray Diffraction (XRD) Analysis
The XRD patterns indicate that the MWCNTs are embedded into the PPy nanocomposite. Furthermore, no significant peak shifts were observed in the XRD patterns. Figure 2 presents the pattern XRD for the experimental and simulated for PPy/MWCNTs composite. When spinning doped, the MWCNTs can be incorporated into the PPy crystal lattice. The XRD pattern obtained from the fabricated PPy/MWCNTs composites thin films. They were correlated to the isolated system matrix. The predicted crystallite size (D) and miller index (hkl) depend on the absolute values of full width at half maximum (FWHM) as given in Table 2. The data in database code amcsd 0020475 corresponds well with the interplanar distances d spacing [37]. The TDDFT-DFT and Crystal Sleuth Microsoft programs were used to designate peak lines calculated by diffraction that were close to the observed findings. The Debye-Scherrer formula was applied to assessed XRD for [PPy/MWCNTs] 0-3 thin films, the range of 5 ≤ 2θ ≤ 45 with 1/dhkl = 0.06Å −1 − 0.74Å −1 , λ = 1.54Å, I 2 /I 1 = 0.5, polarization = 0.5, and Pesedo-Voigt function. From Scherer's formula: where λ is the X-ray wavelength (1.54 Å). As presented in Table 2, for the fabricated [PPy/MWCNTs] composites thin films XRD data from the XRD pattern was used to examine factors and features such as FWHM, the crystallite size (D), hkl indices, d-spacing (d), and peak intensity. The crystallite size is D av = 62.26 nm was within the range of 32.49-144.87 nm. Polymorph calculated the theoretical X-ray diffraction models using content studio software (version 7.0) (See Figure 2). Inset Figure 2, the integrals were conducted on the Brillouin zone with 2 × 2 × 1 (Polymorph PPy/MWCNTs composites).
For the corresponding experimental, a comparison was made between experimental Xrays structures and measured PXRD patterns for PPy/MWCNTs composites thin films. While both experimental and PXRD models the intensity and location of specific peaks vary only slightly, the emphasis here is mostly on their overall resemblance. Only the important comparison characteristics between the measured and the experimental data should therefore be evaluated. It is also known that instrumentation and data collection processes are only two of the many variables that can affect the experimental PXRD pattern.  Figure 2, a combination of experimentally based diffraction and density functional theory calculations yields a great estimation of the atomic scale for PPy/MWCNTs composites thin films (2θ at hkl (111).

Geometric Study
For Figure 3a-f, the similarity of physical-chemical properties of the gaseous phase of PPy/MWCNTs composites was investigated using electrostatic potential and electron density [38,39]. DFT concepts Figure 3a,d, that use electron density as an essential operator for the assessment of isolated state of PPy/MWCNTs composites gas-phase electron systems, respectively. Figure 3b,e explains potential diagrams demonstrating significant potential growth for PPy/MWCNTs composites gas-phase potential, respectively. It supports the possibility of electron transfer to PPy/MWCNTs composites in the gas phase to compute the electrostatic potential (MEP) according to its surface density. Figure 3c,f showing 3D images of the active site of the MEP during the PPy/MWCNTs composites isolated molecules, respectively. The colors blue and red demonstrate the areas advantageous for nuclear and electrical attacks. In isolated molecule phase, the potential range [P]

Geometric Study
For Figure 3a-f, the similarity of physical-chemical properties of the gaseous phase of PPy/MWCNTs composites was investigated using electrostatic potential and electron density [38,39]. DFT concepts Figure 3a,d, that use electron density as an essential operator for the assessment of isolated state of PPy/MWCNTs composites gas-phase electron systems, respectively. Figure 3b,e explains potential diagrams demonstrating significant potential growth for PPy/MWCNTs composites gas-phase potential, respectively. It supports the possibility of electron transfer to PPy/MWCNTs composites in the gas phase to compute the electrostatic potential (MEP) according to its surface density. Figure 3c,f showing 3D images of the active site of the MEP during the PPy/MWCNTs composites isolated molecules, respectively. The colors blue and red demonstrate the areas advantageous for nuclear and  [40,41]. The MEP diagram shows potential negative areas of positive potential for hydrogen atoms. The blue hue represents the most attraction, while the red color represents the most repulsion [42,43]. The single electronegative atom pair was aligned with negative V regions (r). The oxygen atoms of the anion of nitrogen were detected in the examined molecule in negative regions, while maximal positive regions in pyrazine2-one are local to protonated groups (-N(H)-C=C=C-N(CH)-), which may be regarded as possible nucleonic attack sites with a maximum value of +3.87 a.u. According to the following computations, the MEP map showed negative potential locations for N-atoms and positive potentials surrounding the H-atoms. These websites offer a valuable understanding of the intermolecular interactions of the chemical. The presence of intermolecular hydrogen bonding is, therefore, confirmed in Figure 3c,f.  [40,41]. The MEP diagram shows potential negative areas of positive potential for hydrogen atoms. The blue hue represents the most attraction, while the red color represents the most repulsion [42,43]. The single electronegative atom pair was aligned with negative V regions (r). The oxygen atoms of the anion of nitrogen were detected in the examined molecule in negative regions, while maximal positive regions in pyrazine2-one are local to protonated groups (-N(H)-C=C=C-N(CH)-), which may be regarded as possible nucleonic attack sites with a maximum value of +3.87 a.u. According to the following computations, the MEP map showed negative potential locations for N-atoms and positive potentials surrounding the H-atoms. These websites offer a valuable understanding of the intermolecular interactions of the chemical. The presence of intermolecular hydrogen bonding is, therefore, confirmed in Figures 3c,f. The measured ∆ values were based on the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) theory discrepancy by utilizing the DMol 3 /DFT procedure as shown in insert Figure S5.
The boundary molecular orbits in the molecular orbit (FMOs), also known as HOMO and LUMO, are important parameters in quantum chemical simulations for complex analysis (insert Figure S4). The computed energy , and ∆ are presented in Table 3. The following equations were used to compute the reported values of the chemical potential, global hardness, electronegativity, global softness, global electrophilicity index, softness, and the maximum quantity of electronic charge [44].

The measured ∆E
Opt g values were based on the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) theory discrepancy by utilizing the DMol 3 /DFT procedure as shown in insert Figure S5.
The boundary molecular orbits in the molecular orbit (FMOs), also known as HOMO and LUMO, are important parameters in quantum chemical simulations for complex analysis (insert Figure S4). The computed energy E HOMO , E LUMO and ∆E Opt g are presented in Table 3. The following equations were used to compute the reported values of the chemical potential, global hardness, electronegativity, global softness, global electrophilicity index, softness, and the maximum quantity of electronic charge [44].   The negative values of E HOMO and E LUMO energies can be ascribed to product stability for the [PPy/MWCNTs] 0-3 matrix. Coordination position simulation was considered for the highest magnitude molecular orbital coefficients. The electrophilicity index (t), which evaluates energy stability when the device receives an additional electronic charge, is a critical quantum chemical feature [45].

Optical Properties
The observation of the absorption edge in the ultraviolet field will explain crystalline/noncrystalline materials' optical transitions and electronic band structures. Figure 4. demonstrates the PPy/MWCNTs composites thin films absorbance (a) which was computed for the normal incidence of light within a wavelength range 380 ≤ λ ≤ 1000 nm at room temperature (298 K). The relation between the absorption coefficient (α) and the transmittance (T) of a sample is α = 1/d ln(1/T), where d is the thickness (d ∼ = 75 ± 3 nm). The values of the absorption coefficient (α) for the [PPy/MWCNTs] composites are calculated by using the following equation: where A, E g, and n are a constant that is based on the transition probability, the band gap's width and an index characterize the optical absorption processes in the PPy/MWCNTs composites thin films, respectively [46]. Non-linear optical effects influence a new field that changes the point, λ(nm), amplitude, or other incident fields transmission as a result of the electromagnetic fields' relationships with different media. Non-linear optical has an important role in new applications including networking, optical links, and signal processing as the main roles of optical interference. The UV-Vis spectra, at the wavelength range (λ) 300 ≤ λ ≤ 1200 nm of the PPy/MWCNTs composites nanocomposite thin films are illustrated in Figure 4. For PPy/MWCNTs composites as-deposited thin films, in the first region, at a wavelength (λ) region 300 ≤ λ ≤ 705 nm the absorption bands are associated to n→π* electronic transition. The bands at 433, 57, 606, and 653 nm are associated with the electronic transitions from HOMO to LUMO for PPy/MWCNTs composite thin films, respectively. In the second region, at a wavelength range (λ) 705 ≤ λ ≤ 1200 nm, the absorption bands are associated with π→π* electronic transition. Likewise, the absorption bands at 775, 791, and 816 nm are associated with the electrons' transition from valence orbital to conduct orbital for PPy/MWCNTs composites thin film as experimental parts. The UV-Vis spectra, at the wavelength range (λ) 550 ≤ λ ≤ 2900 nm of the [PPy/MWCNTs] in an isolated gaseous state as simulated part (DMOl 3 /DFT) is illustrated in Figure S3. The results of the PPy/MWCNTs in an isolated gaseous state as simulated part (DMOl3/DFT) are in good agreement with PPy/MWCNTs composites thin film which is presented in Figure 4.  The absorption spectrum curve Abs. (λ) was used to measure the band-gap energy of the [PPy/MWCNTs] composites nanocomposite thin film. Tauc's interaction is used to evaluate the values of ( ℎ ) as from Abs. (λ) curve, where . and ℎ are the energy-independent constant and the incident photons energy, respectively [47]. the valence band to the conductive band. Adding the MWCNTs decreased the bandgap as shown in Figure 5. As shown in inset Figure S4, the value of = 2.888 was assessed by applying the DMol 3 process in DFT based on the discrepancy between HOMO and LUMO for free polypyrrole PPy. For Figure S4, the direct energy values are 2.648, 1.553, 1.503, and 1.48 eV for the [PPy/MWCNTs] composites, whereas the value of energy gaps for PPy/MWCNTs composites thin film decreased from 2.648 eV to 1.480 eV with increasing the molar ratio of MWCNTs. This decrease in . is because of increased dis- The absorption spectrum curve Abs. (λ) was used to measure the band-gap energy of the [PPy/MWCNTs] composites nanocomposite thin film. Tauc's interaction is used to evaluate the values of (αhν) A as from Abs. (λ) curve, where Ind. E and hν are the energy-independent constant and the incident photons energy, respectively [47]. The value A = 1 2 for indirect transitions (E Dir g ) and 2 for direct allowed (E Ind. g ) transitions. The expression α = Abs./d was utilized to evaluate the coefficient of absorption (α), where d is the thickness of the film. The E Dir g and E Ind. ruption, which allows the electrons transformation from the valence band to the conductive band. As shown in inset Figure 5, the value of = 1.58 was assessed by applying the DMol 3 process in DFT based on the discrepancy between HOMO and LUMO for PPy/MWCNTs composites in an isolated state. The results of the simulations by using DFT/DMOl 3 and experimental data (Figures 5 and 4s) are well agreed. The equation of Tauc`s and the obtained can be used to evaluate electrical and energy transfer methods effectively. The refractive index n(λ) is an important physical factor for a microscopic atomic interaction, which is important for polymer solar cells as well. (n) and (k) values can be calculated by applying reflectance (R) according to the following equations [48]: where (n) is the real part of the refractive index. Figure 6a-c reveals the dependence of n(λ) and k(λ)on photon energy (hν). As shown in Figure 6a   The refractive index n(λ) is an important physical factor for a microscopic atomic interaction, which is important for polymer solar cells as well. (n) and (k) values can be calculated by applying reflectance (R) according to the following equations [48]: where (n) is the real part of the refractive index. Figure 6a-c reveals the dependence of n(λ) and k(λ) on photon energy (hν). As shown in Figure 6a at ( Figure 6b, the intensity of four peaks observed is increased with an increase in the molar ratio of CNTs, respectively. From the behavior of the simulated nanocomposite PPy/MWCNTs composites as the isolated state in Figure 6c, the CASTEP/DFT computations were used to evaluate n(λ) and k(λ) values and compared to the experimental values, simulated values are close to those achieved by DFT with the CASTEP model. [PPy/MWCNTs] composites, respectively. From the behavior of PPy/MWCNTs composites in Figure 6b, the intensity of four peaks observed is increased with an increase in the molar ratio of CNTs, respectively. From the behavior of the simulated nanocomposite PPy/MWCNTs composites as the isolated state in Figure 6c, the CASTEP/DFT computations were used to evaluate n(λ) and k(λ) values and compared to the experimental values, simulated values are close to those achieved by DFT with the CASTEP model.  According to the single oscillator model, the dispersion curve can be described by dispersion ( ) and oscillating energy ( ) [49]; The ( − 1) is characterized as (ℎ ) (Figure 7a). Table 4 demonstrates typical dispersal parameters ( and ), from both the slope and the linear fit of the high-frequency field. The high-frequency dielectric constant is related to the contribution produced by the electronic polarization process to [PPy/MWCNTs] 0-3 total dielectric reaction. The oscillator strength can be computed = and demonstrated in Table 4. The high-frequency dielectric constant can be computed as follows: The relationship of vs. is demonstrated in Figure 7b. To measure both the ratio ( * ⁄ ) and , the slope and extrapolate can be employed in the resulting linear fittings. In addition, two significant periods like the one-oscillator parameters and depend on and as the following formula [50]: To identify a single oscillator approximation to the dielectric response of the material and find the average bond strength, these moments are estimated. These optical moments According to the single oscillator model, the dispersion curve can be described by dispersion (E d ) and oscillating energy (E 0 ) [49]; The n 2 − 1 −1 is characterized as (hν) (Figure 7a). Table 4 demonstrates typical dispersal parameters (E d and E o ), from both the slope and the linear fit of the high-frequency field. The high-frequency dielectric constant is related to the contribution produced by the electronic polarization process to [PPy/MWCNTs] 0-3 total dielectric reaction. The oscillator strength can be computed f = E 0 and demonstrated in Table 4. The high-frequency dielectric constant can be computed as follows: In te rc ep t = 22 .6 4, S lo pe = -9 . 8 9 In te rc ep t   The difference between ε1(λ) and ε2(λ) as a function of photon energy (hν) exposes various correlations with electrons and photons observed in PPy/MWCNTs composites as-deposited thin film. For PPy/MWCNTs composites as the deposited thin film, the n(λ) values are determined to achieve the lattice constant of a dielectric as the following formulas: and where the coefficient of extinction k(λ) is gained from the relation: The real ε1(λ) and ε2(λ) imaginary dielectric constants are used to illustrate the media's response spectra to electromagnetic radiation incidents [52]. Furthermore, for both status compositions, the computed ε1(λ) is higher than the computed ε2(λ). For PPy/MWCNTs composites, in the experimental part (Figure 8   The relationship of n 2 vs. λ 2 is demonstrated in Figure 7b. To measure both the ratio (N/m * ) and ε ∞ , the slope and extrapolate can be employed in the resulting linear fittings. In addition, two significant periods like the one-oscillator parameters M −1 and M −3 depend on E 0 and E d as the following formula [50]: and To identify a single oscillator approximation to the dielectric response of the material and find the average bond strength, these moments are estimated. These optical moments are also compared to macroscopic numbers as the effective dielectric constant, effective value of electrons in the substance studied [51].  Table 4. The low M −1 and M −3 values provide the signal for low polarization of the examined material. The averaged wavelength (λ 0 ) of the interband oscillator and the average intensity of the oscillator (S 0 ) can be estimated by a single Sellmeier oscillator at low energy: since λ 0 and S 0 values can be achieved from the slope and intercept of plotting (n 2 − 1) −1 against (λ) −1 curve as demonstrated in Figure 7c and tabulated in Table 4. The difference between ε 1 (λ) and ε 2 (λ) as a function of photon energy (hν) exposes various correlations with electrons and photons observed in PPy/MWCNTs composites as-deposited thin film. For PPy/MWCNTs composites as the deposited thin film, the n(λ) values are determined to achieve the lattice constant of a dielectric as the following formulas: and ε(λ) 2 = 2nk((λ) (19) where the coefficient of extinction k(λ) is gained from the relation: The real ε 1 (λ) and ε 2 (λ) imaginary dielectric constants are used to illustrate the media's response spectra to electromagnetic radiation incidents [52]. Furthermore, for both status compositions, the computed ε 1 (λ) is higher than the computed ε 2 (λ). For PPy/MWCNTs composites, in the experimental part (Figure 8 as-deposited thin film. For PPy/MWCNTs composites as the deposited thin film, the n(λ) values are determined to achieve the lattice constant of a dielectric as the following formulas: and where the coefficient of extinction k(λ) is gained from the relation: The real ε1(λ) and ε2(λ) imaginary dielectric constants are used to illustrate the media's response spectra to electromagnetic radiation incidents [52]. Furthermore, for both status compositions, the computed ε1(λ) is higher than the computed ε2(λ). For PPy/MWCNTs composites, in the experimental part (Figure 8    To obtain the spectrum of conductivity, ε1(λ) and ε2(λ) must be combined to the consequential relationship: * ( ) = ( ) + ( ) where the real part is and the imaginary part is fragments of the optical conductivity, (ω) and (εo) are the frequency of angular and the constant of dielectric free space, respectively. The optical conductivity real and imaginary parts of PPy/MWCNTs composites thin film are shown in Figure 9a. The σ2 has a greater value than that of σ1.  Using the CASTEP technique, the maximum value of ε 1 (λ) and ε 2 (λ) for PPy/MWCNTs composites in isolate state is ∼ = 7.50 eV at photon energy (eV) ∼ = 3.00 and 1.50, respectively as shown in Figure 8b. For average values ε 1 (λ) and ε 2 (λ), resulted from the experimental and simulation dimensions are found within the photon energy range values of 1-53 eV. As demonstrated in this figure, one peak was observed in the dielectric constant parts performance of PPy/MWCNTs composites. From the behavior of the simulated composite PPy/MWCNTs as an isolated state in Figure 8b, the CASTEP/DFT computations were used to evaluate ε 1 (λ) and ε 2 (λ)values and compared to the experimental values for PPy/MWCNTs composites thin films, simulated values are close to those achieved by DFT with the CASTEP model.
To obtain the spectrum of conductivity, ε 1 (λ) and ε 2 (λ) must be combined to the consequential relationship: where the real part is (σ 1 (λ) = ωε 2 (λ)ε 0 ) (22) and the imaginary part is (σ 2 (λ) = ωε 1 (λ)ε 0 ) fragments of the optical conductivity, (ω) and (ε 0 ) are the frequency of angular and the constant of dielectric free space, respectively. The optical conductivity real and imaginary parts of PPy/MWCNTs composites thin film are shown in Figure 9a. The σ 2 has a greater value than that of σ 1 . Photon energy (hv)(eV)   Figure S5 at different temperatures 290 ≤ T(K) ≤ 358. The nonlinear coefficient parameter (r) may come from the relationship: where R is the constant, and r is the slope of these curves. There are two areas in I-V curves, the r 1 values for smaller V and r 2 for larger V with r 2 > r 1 . The values of r 1 and r 2 , as recorded in Table 5, confirm the non-ohmic behavior of the I-V characteristics of the investigated films. The values of the non-linear coefficient parameters for Au/[PPy/MWCNTs] composites/n-Si/Al heterojunction diode are r 1 and r 2 , where the r 1 values are less than 2, but the values of r 1 and r 2 increase with increasing MWCNTs doping ratio. In addition, raising T(K) increases r 1 . On the other hand, the r 2 values are decreasing with increasing temperature. The nonlinear coefficient parameters (r) are utilized to statement the conduction mechanism in polymers [53]. We obtain ohmic behavior when r = 1. The dominant mechanism for r = 2 is trap-free space charge incomplete. Finally, if r > 2, the mechanism is defective in terms of trap charge [54]. When r is increased, traps becomes larger or deeper.
Higher V values undoubtedly enhance the rate of CNT aggregation, which is regarded as undesirable and poses a problem [55]. This problem necessitates decreasing the applied potential [56]. Aspect ratio that is high electrical percolation thresholds in MWCNT-based polymer nanocomposites (PNCs) are lower than in polymers loaded with carbon black, carbon fiber, or metals. The electrical percolation threshold is the concentration of filler at which electrical conductivity abruptly rises by many orders of magnitude [57]. The electrical characteristics of polymer nanocomposites are improved by the uniform dispersion of carbon nanotubes. The weak interfacial contact between CNTs and polymeric mixes, as well as the attraction between CNTs caused by the van der Waals force, which causes CNTs to aggregate, make achieving such a distribution difficult [58]. SEM scans indicated that an excellent dispersion of CNTs was obtained in this study. At a voltage of 20 V, DC is shown as a function of T(k). The interaction of PPy with MWCNTs boosts charge transfer inside nanocomposite films. If σ dc = 5.85 × 10 −5 S m −1 [59], PPy is classified as a semiconductor material. The σ dc of PPy/MWCNTs composites are 11.99 × 10 −5 , 34.18 × 10 −5 , and 56.08 × 10 −5 S m −1 , respectively, which are higher than the stated value for the free polypyrrole polymer (PPy) following E 0 g values (see Table 5). σ dc reaches 1.025 × 10 −5 , 1.955 × 10 −5 , 4.765 × 10 −5, and 6.60 × 10 −5 for PPy/MWCNTs composites films, nearly two orders of magnitude increase MWCNT loading was problematic due to (1) the blend's difficulty to absorb additional MWCNTs due to PPy's viscosity, and (2) MWCNTs' high surface energy and inclination to agglomerate. The MWCNTs [60] can be aligned by selecting the appropriate applied electrical field (E) to decrease agglomeration and expand networks from the negative electrode to the positive electrode. The following stages are recommended for generating conducting pathways: first, the CNTs are rotated to a specific angle due to applied E, which generates a dipole moment at the MWCNT edges, aligning them in the direction of E. Second, the CNTs attract one other until they make contact, resulting in the creation of three-dimensional networks. Third, MWCNTs move to and adhere to the negative electrode. In conclusion, the electronic conductivity of CNT is the primary cause of σ dc in nanocomposite films, whereas ionic conductivity is negligible. 10 −4 -10 −5 S m −1 are the values of σ dc in the semiconductor area. Raising T (i.e., semiconducting behavior) causes an increase in charge transfer, as seen in Figure 10. There are no 3D networks produced in the case of polymer and at lower CNT ratios, but the collected energy by charge carriers will activate them to leap potential barriers. Heating, on the other hand, will aid in the optimization of these routes, increasing in σ dc , following increasing MWCNT content and network development. decrease agglomeration and expand networks from the negative electrode to the positive electrode. The following stages are recommended for generating conducting pathways: first, the CNTs are rotated to a specific angle due to applied E, which generates a dipole moment at the MWCNT edges, aligning them in the direction of E. Second, the CNTs attract one other until they make contact, resulting in the creation of three-dimensional networks. Third, MWCNTs move to and adhere to the negative electrode. In conclusion, the electronic conductivity of CNT is the primary cause of in nanocomposite films, whereas ionic conductivity is negligible. 10 −4 -10 −5 S m −1 are the values of in the semiconductor area. Raising T (i.e., semiconducting behavior) causes an increase in charge transfer, as seen in Figure 10. There are no 3D networks produced in the case of polymer and at lower CNT ratios, but the collected energy by charge carriers will activate them to leap potential barriers. Heating, on the other hand, will aid in the optimization of these routes, increasing in , following increasing MWCNT content and network development. It was reported that polypyrrole (PPy)-coated multiwalled carbon nanotubes (MWCNTs) composites were synthesized by simple, cost-effective in situ oxidative polymerization method [61]. Gas sensors were prepared in pellet form, which would be robust, cheap and reasonably sensitive to ammonia vapor sensing. The gas sensitive characteristics of composite for different MWCNTs content and over wide range of NH3 vapor concentration were investigated at room temperature. Also, the effect of operating temperature, humidity, long term stability was also studied. Also, the temperature depend- It was reported that polypyrrole (PPy)-coated multiwalled carbon nanotubes (MWC-NTs) composites were synthesized by simple, cost-effective in situ oxidative polymerization method [61]. Gas sensors were prepared in pellet form, which would be robust, cheap and reasonably sensitive to ammonia vapor sensing. The gas sensitive characteristics of composite for different MWCNTs content and over wide range of NH 3 vapor concentration were investigated at room temperature. Also, the effect of operating temperature, humidity, long term stability was also studied. Also, the temperature dependence of DC conductivity of the PPy-MWCNT nanocomposites had been measured [62], the temperature deviation of DC conductivity was measured contained by a range 80 < T < 300 K. Resistivity decreases with the combination of MWCNTs in the PPy environment. PPy behaves as a good electron donor and MWCNTs are somewhat good electron acceptors. Therefore, there are some associations among the quinoid rings of PPy and MWCNTs which facilitate the charge transfer between the two components [63,64]. Figure 11 shows the current density-voltage (J-V) characteristics of the Au/[PPy/MWCNTs]/ n-Si/Al solar cell at various light intensities P in. In Table 6, the photovoltaic performance of the manufactured solar cell was tabulated, including the existence of current density (J sc ) at zero voltage and the existence of voltage (V oc ) at zero current density. When the intensity of the light is increased, both (J sc ) and (V oc ) values rise. The following equation is satisfied by increasing the J sc with the light intensity: (25) where A is a constant. By graphing (ln Jsc) vs. (ln P in ), the value of the exponent (γ) was found to be 1.28. Monomolecular and bimolecular recombination mechanisms are represented by the numbers 1.0 and 0.5 for (γ). The value of the exponent, on the other hand, falls between 0.5 and 1.0, indicating a continuous distribution of trapping centers. As a result, the obtained value of (γ) suggests that the cell under examination has a monomolecular recombination process. The current density J m and voltage V m correlate to the power density (maximum value PM). In Table 6, the matching current density and voltage were calculated and presented. The table shows that when the light intensity increases, the values of PM, J m, and V m rise. The device's fill factor (FF) and power conversion efficiency (η) were then determined using the formulae below [65]:

Photovoltaic Properties of PPy/MWCNTs Composites Films
In Table 6, the values of FF and η were provided. It can be seen that when the light intensity is increased, the efficiency improves. At 150 mW/cm 2 , the efficiency was determined to be 6.86%.
where A is a constant. By graphing (ln Jsc) vs. (ln Pin), the value of the exponent (γ) was found to be 1.28. Monomolecular and bimolecular recombination mechanisms are represented by the numbers 1.0 and 0.5 for (γ). The value of the exponent, on the other hand, falls between 0.5 and 1.0, indicating a continuous distribution of trapping centers. As a result, the obtained value of (γ) suggests that the cell under examination has a monomolecular recombination process.

Conclusions
Multi-walled carbon nanotubes were successfully coated by polypyrrole using in situ oxidative polymerization method. The polymerization process was carried out in an acidic medium by ferric chloride as an oxidizing agent in the presence of sodium dodecyl sulfate as a soft template. XRD data showed the experimental and simulated of PPy/MWCNTs composites had a triclinic crystal symmetry in the space group P1. The crystalline size is D av = 62.26 nm, was within the range of 32.49-144.87 nm. The direct energy gap for PPy varies by a factor of 2.413 eV to 1.475 eV depending on the MWCNTs loading. The TDDFT-DFT description of geometry parameters was determined by using both DMOl3 and CASTEP techniques. The film's I-V characteristics indicated non-ohmic behavior. They were discovered to be temperature and light-intensity-dependent. At low temperatures, the conduction mechanism in Au/n-[PPy/MWCNTs]TF/p-Si/Al was identified to be Schottky emission, but at high temperatures, the Poole-Frenkel effect was used in the [PPy/MWCNTs]TF. The current-voltage characteristics of the Au/n-[PPy/MWCNTs] thinfilm/p-Si/Al heterojunction device (thickness 75 ± 3 nm) were investigated under varying illumination intensities ranging from 30 to 150 W/m 2 . As a result of the production of charge carriers, the calculated current increased dramatically with increasing light intensity. At illumination intensities of 150 W/m 2 , the resultant average value of power conversion efficiency (PCE) of an Au/n-[[PPy/MWCNTs] thin film/p-Si/Al heterojunction solar cell was computed and found to be ∼ =6.86%. The nanocomposite thin films demonstrated promising results, making them a viable candidate for polymer solar cell demands.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.