Physicochemical Properties and Atomic-Scale Interactions in Polyaniline (Emeraldine Base)/Starch Bio-Based Composites: Experimental and Computational Investigations

The processability of conductive polymers still represents a challenge. The use of potato starch as a steric stabilizer for the preparation of stable dispersions of polyaniline (emeraldine base, EB) is described in this paper. Biocomposites are obtained by oxidative polymerization of aniline in aqueous solutions containing different ratios of aniline and starch (% w/w). PANI-EB/Starch biocomposites are subjected to structural analysis (UV-Visible, RAMAN, ATR, XRD), thermal analysis (TGA, DSC), morphological analysis (SEM, Laser Granulometry), and electrochemical analysis using cyclic voltammetry. The samples were also tested for their solubility using various organic solvents. The results showed that, with respect to starch particles, PANI/starch biocomposites exhibit an overall decrease in particles size, which improves both their aqueous dispersion and solubility in organic solvents. Although X-ray diffraction and DSC analyses indicated a loss of crystallinity in biocomposites, the cyclic voltammetry tests revealed that all PANI-EB/Starch biocomposites possess improved redox exchange properties. Finally, the weak interactions at the atomic-level interactions between amylopectin–aniline and amylopectin–PANI were disclosed by the computational studies using DFT, COSMO-RS, and AIM methods.


Introduction
A variety of reports have emerged in recent decades documenting the use of sterically stabilized colloidal dispersions of conductive polymers, which are much more processable than traditional ones [1]. In the case of dispersion polymerization, and when the reaction medium contains a soluble stabilizer, the macroscopic aggregation and precipitation of the resulting conductive polymers can be avoided, leading to submicrometer well-dispersed particles [2]. In this context, several water-soluble polymers, such as cellulosic derivatives [3,4], poly(vinyl alcohol) [5], poly(ethylene oxide) [6], and poly(methyl vinyl ether) [7], were used as steric stabilizers to overcome the processing hardness of

Characterization
Structural characterization of samples was carried out by FTIR spectroscopy using a Nicolet 5700 spectrometer (Thermo Fisher, Waltham, MA, USA), equipped with a diamond ATR sampling accessory, on the powdered material, UV-Visible spectroscopy using a Cary 5000 UV−vis−NIR spectrometer (Varian, Palo Alto, CA, USA) on solutions prepared by dissolving 2 mg in 10 mL of DMF, and RAMAN spectroscopy using a DXR2 Raman spectrometer (Thermo Scientific, Madison, WI, USA), operating at 532 nm. The solubility tests were performed at different temperatures (20,40,50,60,80 • C) by dispersing the same weight of all materials in 10 mL of solvents, to obtain a dilute concentration (0.2% w/v). The solvents used are water, acetone, methanol, ethanol, glycerol, ethyl acetate, chloroform, toluene, hexane, cyclohexane, DMSO, and DMF. While colorless supernatants were obtained with insoluble mixtures, the presence of residual solids, in the case of partially soluble mixtures, was observed by a lens.
Homopolymers and biocomposites were subjected to laser granulometry analysis by using a CILAS 1190 analyzer (CPS Us, Inc., Madison, WI, USA). Thermal gravimetric analyses (with a TGA Q5000 instrument, TA Instruments, Hüllhorst, Germany) were performed until 700 • C (10 • C/min), while DSC (with a TA instrument Q20) was set from room temperature to 300 • C at 10 • C/min, under a nitrogen atmosphere; thermodynamic parameters were calculated by using the software TA Universal Analysis. X-ray diffraction analysis was performed using CuKα radiation (λ = 1.54 Å) (X'Pert3 Powder apparatus) from 2θ = 10 • to 60 • , at a scan rate of 0.02 s −1 ; the peak crystalline areas were determined by using the X'pert HighScore software. SEM images of the composites were obtained after metallization with sputter quorum Q 150R E settings using Quanta 200 apparatus (FEI, Hillsboro, OR, USA). Cyclic voltammetry tests (with PGZ 301 potentiometer, Radiometer Analytical SAS, Lyon, France) were examined in a three-electrode system constituted of a reference electrode (Ag/AgCl), an auxiliary platinum electrode, and a carbon electrode as the working electrode. The electrolyte was a KCl solution (1 M) and the voltammograms were recorded in a voltage window between −0.8 and +0.8 mV at various scan rates (10,50, and 100 mV s −1 ).
The HOMO (Highest Occupied Molecular Orbital) and LUMO (Lowest Unoccupied Molecular Orbital) frontier molecular orbitals (FMOs) were investigated [52]. The high value of the energy ξ HOMO justifies the inclination to transfer electrons to a suitable acceptor molecule. In contrast, the low value of the energy ξ LUMO justifies the molecule's capacity to accept electrons from donor molecules [53].
The electronic chemical potential (µ) and global hardness (η) were determined using the ξ HOMO and ξ LUMO energies [54][55][56]: The chemical potential µ is related to the electrophilicity index ω by the following relationship [57]: ω indicates an electrophile's capacity to gain an extra electrical charge [58]. The optimal fraction of transferred electrons (∆N) is given by [59]: COSMO-RS (Conductor-like Screening Model for Real Solvents) is a thermodynamicbased quantum chemistry technique for determining chemical potentials in solutions [60]. For each species in the solution, this technique can estimate sigma charge densities and chemical potentials. The computation is divided into two parts: first, the molecule was geometrically optimized; then, using the COSMOtherm program [61] (p. 2), sigma profiles and sigma potentials were computed using the acquired cosmo-files.
Intermolecular interactions may be calculated using the Atoms in Molecule (AIM) method. The main goal of AIM is to examine the character and strength of the bonding interaction in molecular systems by using the electron density ρ(r) of the molecules as a tool. Based on the second derivative, ∇^2 ρ(r) sign [62,63], the kind of chemical bonds may be identified at the bond critical point (BCP) where it is minimum. As a result, if the ρ value is large and ∇^2 ρ(r) is negative, the bond is covalent (polar). On the other hand, positive ∇^2 ρ(r) indicates that the kinetic energy G(r) is greater than the potential energy V(r) [64]. The |V|/G ratio indicates that the system's ability to group electrons, V(r), and its ability to dilute them through electronic mobility, G(r), compete within the system. As a result, the ratio |V|/G > 1 indicates that an excess of electric charges is utilized to generate the interaction, whereas |V|/G < 1 indicates the opposite. The total energy density, H(r) = V(r) + G(r), in the first situation, has a negative value, whereas in the second case, it has a positive value. Thus, the categorization of interactions is based on the values of that indicator: (i) pure interactions with closed layers, |V|/G < 1; (ii) interactions with closed layers, 1 < |V|/G < 2; and (iii) interactions with the shared layer, |V|/G > 2. The first two forms of interaction are ascribed to hydrogen bonding, whereas the third type is a covalent interaction. Furthermore, positive, ∇^2 ρ(r), and H values indicate electrostatic contact, whereas negative values indicate a covalent bond. A partly covalent bond is defined as a positive ∇^2 ρ(r) combined with a negative H value [65].
The Amsterdam Density Functional (ADF) program [66,67] was used to perform AIM calculations. For molecular structure optimization, the DFT method was employed. The AIM research [67,68] used the def-TZVP basis set to describe exchange correlation effects between these variables (Becke 3-Parameter, Lee, Yang, and Parr).

Physico-Chemical Properties
The scanning electron micrographs for polyaniline, starch, and polyaniline/starch biocomposites are shown in Figure 1a stabilizer. At the lowest concentration of aniline, the surface of the starch granules is rougher than that of pure starch, suggesting that aniline polymerization occurs on the surface of the granules of starch, forming an irregular coating over it [42,47]. However, as the concentration of aniline increases, more of the starch surface is covered and an overgrowth of polyaniline is observed over the starch particles [45]. pseudo-spherical granules with smooth and well-defined surfaces. Figure 1b-d represent SEM images of starch/polyaniline composites with different ratios of aniline to starch. The figure shows that, in the case of a composite, when aniline is polymerized in the presence of starch, the polyaniline chains produced are grafted or adsorbed on this steric stabilizer. At the lowest concentration of aniline, the surface of the starch granules is rougher than that of pure starch, suggesting that aniline polymerization occurs on the surface of the granules of starch, forming an irregular coating over it [42,47]. However, as the concentration of aniline increases, more of the starch surface is covered and an overgrowth of polyaniline is observed over the starch particles [45]. From the analysis of particle size by laser granulometry, it can be seen that the particles of the obtained PANI/starch materials exhibit a wide distribution of particle sizes revealing a high polydispersity of particles (relative mean variance of the particle size distribution), which may indicate the occurrence of some aggregation of the biocomposites particles [69]. This statement is consistent with the findings of other authors [70,71]. Furthermore, as shown in Table 2, homopolymers and biocomposites particles exhibit average diameters ranging from 28.8 µm to 68.9 µm. It is also worth noting that all composites possess particle sizes lower than that of polyaniline, whereby the composite with the highest aniline content exhibits the lowest average diameter. Thus, the growing polyaniline chains may affect the intramolecular forces along the polymer chains of amylose and amylopectin, leading to smaller particles [15].  From the analysis of particle size by laser granulometry, it can be seen that the particles of the obtained PANI/starch materials exhibit a wide distribution of particle sizes revealing a high polydispersity of particles (relative mean variance of the particle size distribution), which may indicate the occurrence of some aggregation of the biocomposites particles [69]. This statement is consistent with the findings of other authors [70,71]. Furthermore, as shown in Table 2, homopolymers and biocomposites particles exhibit average diameters ranging from 28.8 µm to 68.9 µm. It is also worth noting that all composites possess particle sizes lower than that of polyaniline, whereby the composite with the highest aniline content exhibits the lowest average diameter. Thus, the growing polyaniline chains may affect the intramolecular forces along the polymer chains of amylose and amylopectin, leading to smaller particles [15]. the majority of the solvents used even at high temperatures (80 • C). However, starch is soluble in DMSO but slightly soluble in water and DMF, whereas PANI-EB is slightly soluble in DMSO and soluble in DMF, and it is well known that it is not soluble in most common solvents [72]. The solubility behavior of PANI-EB/starch biocomposites differed from that of homopolymers where we noticed a good solubility in DMF and DMSO and an enhanced dispersion in water at 20-40 • C before more homogeneous mixtures were achieved at higher temperatures (60-80 • C). A further enhancement in the solubility was observed with glycerol and chloroform mainly at higher temperatures (>40 • C) with good dispersions at low temperatures. These findings indicate that the surfaces of the products have, at different temperatures, different affinities for each solvent. Nevertheless, the remaining solvents were not able to dissolve or even disperse PANI/starch biocomposites.

Structural Characterization
The UV-visible, FTIR, and RAMAN spectra of our products are shown in Figures 2-4, respectively. The UV-visible absorption spectra in Figure 2 display two broad absorption bands around 315-360 nm and 560-600 nm for all composites, which confirm the occurrence of the PANI-EB backbone structure. The first peak is assigned to π-π* transition occurring from the highest occupied molecular orbital to the lowest unoccupied molecular orbital of the amine benzenoid rings, while the second peak is ascribed to benzenoid-quinoid rings' charge transfer [73]. The difference in the relative intensities of the peaks observed can be due to the variation in the quantities of benzenoid and quinoid units present in the materials as the ratio of aniline to starch varies. Moreover, the red shift for the 3:1 sample may be due to a higher amount of longer chains as well as to a growing contribution of a less oxidized form of PANI [74,75]. Figure 3 shows the FTIR-ATR spectra of polyaniline, starch, and their biocomposites. For the native starch, the broadband characteristic for the O-H groups of starch appears at 3252 cm −1 and the peak around 2924 cm −1 is attributed to an asymmetric stretching vibration of the C-H bond in pyranoid rings. Several other adsorption bands between 854 cm −1 and 1415 cm −1 are attributed to the contribution of various functional groups such as C-O and C-O-C [47]. Moreover, the peaks around 991 to 994 cm −1 are ascribed to the "C-O" bonds resulting from "C-O-H" and those of the anhydrous glucose rings in the starch backbone [76]. The characteristic peaks of pure PANI have been reported in previous works [42,47]. The bands appearing at 1539 cm −1 and 1456 cm −1 can be assigned to the vibration of quinoid rings ("N=Q=N") and benzenoid rings ("N-B-N"), respectively. The bands at 1282 cm −1 and 1230 cm −1 correspond to the "C-H" stretching vibration within the aromatic conjugation. The presence of PANI is confirmed in all composites' Polymers 2022, 14, 1505 7 of 22 spectra in which PANI and starch are present at different concentrations. As PANI is more exposed on the surface of the particles, its signals dominate the spectra; moreover, the strong interaction between the chains may be responsible for the disappearance of some typical starch bands [77].   [47]. Moreover, the peaks around 991 to 994 cm −1 are ascribed to the "C-O" bonds resulting from "C-O-H" and those of the anhydrous glucose rings in the starch backbone [76]. The characteristic peaks of pure PANI have been reported in previous works [42,47]. The bands appearing at 1539 cm −1 and 1456 cm −1 can be assigned to the vibration of quinoid rings ("N=Q=N") and benzenoid rings ("N-B-N"), respectively. The bands at 1282 cm −1 and 1230 cm −1 correspond to the "C-H" stretching vibration within the aromatic conjugation. The presence of PANI is confirmed in all composites' spectra in which PANI and starch are present at different concentrations. As PANI is more exposed on the surface of the particles, its signals dominate the spectra; moreover, the strong interaction between the chains may be responsible for the disappearance of some typical starch bands [77].   Figure 3 shows the FTIR-ATR spectra of polyaniline, starch, and their biocomposites. For the native starch, the broadband characteristic for the O-H groups of starch appears at 3252 cm −1 and the peak around 2924 cm −1 is attributed to an asymmetric stretching vibration of the C-H bond in pyranoid rings. Several other adsorption bands between 854 cm −1 and 1415 cm −1 are attributed to the contribution of various functional groups such as C-O and C-O-C [47]. Moreover, the peaks around 991 to 994 cm −1 are ascribed to the "C-O" bonds resulting from "C-O-H" and those of the anhydrous glucose rings in the starch backbone [76]. The characteristic peaks of pure PANI have been reported in previous works [42,47]. The bands appearing at 1539 cm −1 and 1456 cm −1 can be assigned to the vibration of quinoid rings ("N=Q=N") and benzenoid rings ("N-B-N"), respectively. The bands at 1282 cm −1 and 1230 cm −1 correspond to the "C-H" stretching vibration within the aromatic conjugation. The presence of PANI is confirmed in all composites' spectra in which PANI and starch are present at different concentrations. As PANI is more exposed on the surface of the particles, its signals dominate the spectra; moreover, the strong interaction between the chains may be responsible for the disappearance of some typical starch bands [77].   The Raman spectra of PANI, starch, and their biocomposites, obtained using an exciting wavelength of 532 nm, are displayed in Figure 4. This method is very efficient for understanding the interaction among various components and provides information about the chemical structure of different domains in the composite systems.
From the spectrum of PANI, we can observe some characteristic bands rinsing from carbon-carbon, carbon-hydrogen, and carbon-nitrogen vibration modes. Benzenoid rings exhibited C-C stretching modes between 1520 cm −1 and 1650 cm −1 and C-H bending modes between 1100 cm −1 and 1210 cm −1 ; however, quinoid rings give rise to C=N, and C-C ring stretching vibrations, which appeared at 1382 cm −1 and 1595 cm −1 , respectively. Further bands are observed in the region 1210-1520 cm −1 for C-N stretching modes (amines, imines) and in the region 1100-1140 cm −1 characteristic of the PANI backbone. A weak band at 1170 cm −1 has been assigned for in-plane deformation of the C-C bond of the quinoid ring of PANI [46,78].
Similar to FTIR-ATR, the Raman curves of PANI/starch biocomposites display the typical peaks of PANI (Figure 4), with no meaningful shift. Moreover, knowing that Raman spectroscopy is a surface thickness-dependent tool [80,81], the response of PANI/starch biocomposites to laser excitation could be taken as a piece of further evidence that starch grains are covered by a thick PANI chain layer, which prevents the penetration of the laser beam to starch macromolecules and consequently impedes the appearance of starch characteristic peaks. These results concord with the SEM images in Figure 1.

X-ray Diffraction Pattern (XRD)
The X-ray diffraction patterns of starch, polyaniline, and PANI/starch biocomposites are presented in Figure 5. A semi-crystalline pattern can be observed for potato starch that exhibits strong reflection at 2θ = 17.2°, moderate reflections at 2θ = 5.7°, 15.04°, 22.22°, and 24.1°, in addition to weak reflections at 2θ = 15.2° and 19.6°. Such a pattern is typical for The observed bands in the Raman spectrum of starch, which are related to C-H and CH 2 vibrations (deformation, bending), appear in the region 1460-1340 cm −1 and at 862 cm −1 . The signals related to C-O-H bending vibrations arise at the wavenumbers 1335, 1123, and 1080 cm −1 . At 1257 cm −1 , the band related to vibrational modes of CH 2 OH appears, while those associated with the skeletal mode involving α-(1-4) linkage appear at 942 cm −1 and between 770-360 cm −1 [79].
Similar to FTIR-ATR, the Raman curves of PANI/starch biocomposites display the typical peaks of PANI ( Figure 4), with no meaningful shift. Moreover, knowing that Raman spectroscopy is a surface thickness-dependent tool [80,81], the response of PANI/starch biocomposites to laser excitation could be taken as a piece of further evidence that starch grains are covered by a thick PANI chain layer, which prevents the penetration of the laser beam to starch macromolecules and consequently impedes the appearance of starch characteristic peaks. These results concord with the SEM images in Figure 1.
The diffraction spacing and Miller indices (hkl) for our materials are displayed in Table 4. The values of the degree of crystallinity (Xc%) were obtained as a ratio between the peaks crystalline area and the total area of the XRD pattern under the curve. The following equation is employed [85]:  The values of the degree of crystallinity (Xc%) were obtained as a ratio between the peaks crystalline area and the total area of the XRD pattern under the curve. The following equation is employed [85]: where Hc and Ha are the intensities for the crystalline and amorphous profiles. The crystallinity degree values of polyaniline and starch, as well as of biocomposites are summarized in Table 5. The native starch shows a high crystallinity (37%) when compared with polyaniline (19%). For PANI-starch biocomposites, a contribution due to different starch content can be confirmed. Changes in XRD patterns could be noticed, which may reveal alterations in the ultrastructure of the microdomains of PANI and starch mainly due to changes in the crystallinity. Such an observation is supported by the disappearance of starch peaks at 2θ = 22.22 • and 26.58 • and all PANI-EB peaks. Thus, it appears that the existence of aniline in the reaction medium and the PANI chains in propagation, partially destroy the original crystalline lattice of starch structure, and this is particularly evident for M 21 and M 31 . On the other hand, it seems that polyaniline chains can form an organized structure, which is also reflected in the diffractograms, although with a lack of peaks attributed to the PANI-EB crystalline structure.

Thermal Characterization
The thermal behaviors of PANI, starch, and PANI/starch biocomposites were studied, under nitrogen, by TGA and DSC analyses. The corresponding thermogravimetry/differential thermogravimetry (TG/DTG) curves, depicted in Figure 6a-e, point out a weight loss within a broad range for polyaniline and as a sharp step for the degradation of starch. PANI-EB exhibits the first weight loss (Figure 6a) between 40 • C and 100 • C, which could be due to the evaporation of water and residual volatile molecules, followed by a small endothermic peak observed between 210 • C and 250 • C. A massive weight loss starts from 400 • C, owing to the degradation of the polyaniline chains [86]. The native starch undergoes 50% weight loss at 290 • C, while for polyaniline, a corresponding weight loss is noticed at 461 • C (Table 6). Furthermore, the DTG curves also show an overall decrease in the degradation rate for biocomposites with the rise in aniline amount, indicating a strong influence of polyaniline on the thermal stability of starch. As a matter of fact, the typical decomposition feature of PANI, with a strong DTG peak at 500 • C, is only present for M11, while M21 (not shown) and M31 show a quite regular weight loss till the limit of the explored range. Detailed data can be derived from the TGA thermograms, which are related to the thermal stability such as T5, T50, and T95, defined as the temperature at which 5, 50 and 95% of the mass is volatilized, respectively. As observed in Table 6, PANI-EB and biocomposites are remarkably more stable than starch. Table 6. Thermal characteristics of PANI, starch, and PANI/starch biocomposites (% weight loss at specific temperatures).  Table 6. Thermal characteristics of PANI, starch, and PANI/starch biocomposites (% weight loss at specific temperatures). Detailed data can be derived from the TGA thermograms, which are related to the thermal stability such as T5, T50, and T95, defined as the temperature at which 5, 50 and 95% of the mass is volatilized, respectively. As observed in Table 6, PANI-EB and biocomposites are remarkably more stable than starch. Figure 7 displays the DSC profiles of polyaniline, starch, and their biocomposites. Broad endothermic peaks are perceived in the temperature range of 60-120 • C, which may correspond to the evaporation of water as supported by the TGA and DTG data. Moreover, the PANI thermogram shows another endothermic step starting from 200 • C, which could be attributed to the glass transition temperature (Tg), where the movement of small chain segments starts [87,88]. On the other hand, the thermogram of starch reveals two narrow endothermic peaks centered at 266 • C and 278 • C that are attributed to the thermal decomposition of starch chains [89]. This result corroborates with the TGA obtained data. Regarding PANI/starch materials, we observe the presence of the endothermic peak related to water evaporation in the same temperature range, but the second endothermic peak, related to the variation in the Tg of polyaniline, shifts to a lower temperature, especially for M 11 , indicating structural modification in the biocomposites owing to the interactions between components. The results of DSC measurement are reported in Table 7.

Cyclic Voltammetry
The electroactivity of PANI and PANI-EB/starch composites was studied in the supporting aqueous electrolyte (KCl) by cyclic voltammetry, and the recorded voltammograms are shown in Figure 8. All the curves displayed a pair of oxidation and reduction peaks, arising from the redox behavior of PANI-EB. Indeed, the peaks appearing for each oxidation and reduction process for PANI and PANI-based composites can be assigned to leucoemeraldine/emeraldine transition states of polyaniline [44,90]. In most cases, voltammograms exhibit a well-defined oxidation peak at about 120 mV that appears to be associated with a less-evident reduction peak at about 110 mV, highlighting the occurrence of a redox electron transfer [91]. The species mobility within the composite was evaluated from the dependence of the redox peaks intensity on the scan rate, from 10 to 100 mV/s. As reported in Table 8, the current density increases with the increasing scan rate value, and the waves show an increase in their areas with increasing scan rates, revealing the good electroactivity as well as the electrochemical capacitive behavior of our materials [44,92]. In order to better understand the electronic behavior of PANI and its composites with starch, we have plotted the variations in the current density of the oxidation peak for each sample as a function of the square root of the scan rate √ v. We can observe in Figure 9 that the current density is linear with √ v for all the synthesized materials. The oxidation phenomenon is therefore limited by the ion diffusion, and the transport of charged species to the electrode surface is ensured by this transport mode. A = electrode area in cm 2 D = diffusion coefficient in cm 2 /s C = concentration in mol/cm 3 ν = scan rate in V/s Therefore, it is clear that the PANI synthesized by chemical oxidation polymerization is a conductor material. Note that according to our results, the M11 sample is the best conductor among these materials.

Computational Assessment
The main electronic properties calculated for the aniline, PANI, and starch are summarized in Table 9. As amylopectin is the most abundant component (85%) in the starch under investigation, it was considered a suitable model. According to the energy values, comparable values of the aniline isomers ξHOMO and ξLUMO indicate a quasi-equivalent tendency to donate or accept electrons. That said, isomers I1 and I2 have the highest values, showing that they are the more probable molecules participating in the polymerization process to produce the PANI chain (this is confirmed by the highest electron escaping tendency from an equilibrium system, µ = −3.027 eV). Both PANI and amylopectin are more reactive than the PANI parent species (isomers I1 or I2) (see Table 9). Regarding the gap values in Table 10, the amylopectin is the most stable in this set of polymers.
In conclusion, PANI has the highest ability to acquire an extra electron charge (highest ω value: 8.386 eV) and the lowest resistance to charge transfer (the lowest hardness η value, 0.895) (LUMO of PANI, see Table 10). It is thus concluded that amylopectin is the electron donor with the highest ξHOMO = −7.056 eV (see Table 2). The results for ΔN are In fact, the Randles Sevcik relation gives a straight line when pure diffusion controls the electrochemical reaction. I pa = 2.69 × 10 5 n 3/2 A D 1/2 C v 1/2 (7) ip = current maximum in amps n = number of electrons transferred in the redox event A = electrode area in cm 2 D = diffusion coefficient in cm 2 /s C = concentration in mol/cm 3 v = scan rate in V/s Therefore, it is clear that the PANI synthesized by chemical oxidation polymerization is a conductor material. Note that according to our results, the M11 sample is the best conductor among these materials.

Computational Assessment
The main electronic properties calculated for the aniline, PANI, and starch are summarized in Table 9. As amylopectin is the most abundant component (85%) in the starch under investigation, it was considered a suitable model. According to the energy values, comparable values of the aniline isomers ξ HOMO and ξ LUMO indicate a quasi-equivalent tendency to donate or accept electrons. That said, isomers I1 and I2 have the highest values, showing that they are the more probable molecules participating in the polymerization process to produce the PANI chain (this is confirmed by the highest electron escaping tendency from an equilibrium system, µ = −3.027 eV). Both PANI and amylopectin are more reactive than the PANI parent species (isomers I1 or I2) (see Table 9). Regarding the gap values in Table 10, the amylopectin is the most stable in this set of polymers.

An PANI Amylopectin
HOMO LUMO

An PANI Amylopectin
HOMO LUMO

An PANI Amylopectin
HOMO LUMO

An PANI Amylopectin
HOMO LUMO

An PANI Amylopectin
HOMO LUMO

An PANI Amylopectin
HOMO LUMO Table 11. Fukui analysis of aniline, polyaniline, and amylopectin: The nucleophilic term f_k^+; the electrophilic term f_k^-values of the Fukui function.

An PANI Amylopectin f+
In conclusion, PANI has the highest ability to acquire an extra electron charge (highest ω value: 8.386 eV) and the lowest resistance to charge transfer (the lowest hardness η value, 0.895) (LUMO of PANI, see Table 10). It is thus concluded that amylopectin is the electron donor with the highest ξ HOMO = −7.056 eV (see Table 2). The results for ∆N are 4.328 electrons (acquired) for PANI and 0.894 electrons (released) for amylopectin.
The Fukui results shown Table 11 Table 11. Fukui analysis of aniline, polyaniline, and amylopectin: The nucleophilic term f_kˆ+; the electrophilic term f_kˆ-values of the Fukui function.
An PANI Amylopectin f + LUMO Table 11. Fukui analysis of aniline, polyaniline, and amylopectin: The nucleophilic term f_k^+; the electrophilic term f_k^-values of the Fukui function. The sigma profile above is divided into three distinct regions:

An
LUMO Table 11. Fukui analysis of aniline, polyaniline, and amylopectin: The nucleophilic term f_k^+; the electrophilic term f_k^-values of the Fukui function. The sigma profile above is divided into three distinct regions:

An
LUMO Table 11. Fukui analysis of aniline, polyaniline, and amylopectin: The nucleophilic term f_k^+; the electrophilic term f_k^-values of the Fukui function. The sigma profile above is divided into three distinct regions: LUMO Table 11. Fukui analysis of aniline, polyaniline, and amylopectin: The nucleophilic term f_k^+; the electrophilic term f_k^-values of the Fukui function. The sigma profile above is divided into three distinct regions: LUMO Table 11. Fukui analysis of aniline, polyaniline, and amylopectin: The nucleophilic term f_k^+; the electrophilic term f_k^-values of the Fukui function. Positive sigma values represent negative polarities. Figure 10a shows that the highest peaks for all the selected molecules are in the non-polar region. They are showing the tremendous non-polar character of the molecule surfaces. Amylopectin has a net peak in the HBA region, while PANI is in the HBD region. In addition, PANI has the highest affinity to HBA, and amylopectin has the highest affinity to HBD (Figure 10b). This confirms the frontier molecular orbital study results.

An
The interaction energies reported in In addition, PANI has the highest affinity to HBA, and amylopectin has the highest affinity to HBD (Figure 10b). This confirms the frontier molecular orbital study results. The interaction energies reported in Table 12, reveal that the HB type contributes the most to the mixture's energy (−31.33496 and −40.29382 Ha). The PANI-amylopectin pair (ratio 1:1) achieves the highest value, while I2-amylopectin (ratio 3:1) has the lowest. The second contribution to the mixture's energy is Van der Waals (between −26.70916 and −27.83387 Ha), with the minimum value for the PANI-amylopectin couple (1:1) and the maximum value for I2-amylopectin (3:1). The mixture's energy ranges from −55.25390 for I2-amylopectin (3:1) to −63.51928 Ha for PANI-amylopectin (Emisfit). Table 12. Interaction energies of amylopectin with the various compounds, which are depicted in Figure 11.    Figure 11.  Figure 12 shows a molecular graph of the optimized PANI-amylopectin complex, while Tables 13 and 14 give the topological characteristics at the BCP of interaction contacts in the PANI-amylopectin complex. The computed values of ρ(r) in BCPs are in the range [2.21 × 10 −3 -7.83 × 10 −3 ] eV (Table 13), which are relatively low values, and ∇ 2 2ρ(r) is positive, indicating the existence of hydrogen bond interactions in BCPs 231, 236, 239, 323, 324, 361, and 363 (BCPs are given in Table 13). The positive values of ∇ 2 2ρ(r) and H demonstrate the H-bonding connection between the BCPs atoms (see Figure 12). Weak hydrogen bonding is characterized by E HB < 0. The interaction between PANI and amylopectin is thus considered to be physical (with no covalent bonding). Polymers 2022, 14, x FOR PEER REVIEW 17 of 22 Figure 11. Interactions in PANI/starch biocomposites using COSMO-RS. Figure 12 shows a molecular graph of the optimized PANI-amylopectin complex, while Tables 13 and 14 give the topological characteristics at the BCP of interaction contacts in the PANI-amylopectin complex. The computed values of ρ(r) in BCPs are in the range [2.21 × 10 −3 -7.83 × 10 −3 ] eV (Table 13), which are relatively low values, and ∇ 2 2ρ(r) is positive, indicating the existence of hydrogen bond interactions in BCPs 231, 236, 239, 323, 324, 361, and 363 (BCPs are given in Table 13). The positive values of ∇ 2 2ρ(r) and H demonstrate the H-bonding connection between the BCPs atoms (see Figure 12). Weak hydrogen bonding is characterized by EHB < 0. The interaction between PANI and amylopectin is thus considered to be physical (with no covalent bonding).      (Table 13), which are relatively low values, and ∇ 2 2ρ(r) is positive, indicating the existence of hydrogen bond interactions in BCPs 231, 236, 239, 323, 324, 361, and 363 (BCPs are given in Table 13). The positive values of ∇ 2 2ρ(r) and H demonstrate the H-bonding connection between the BCPs atoms (see Figure 12). Weak hydrogen bonding is characterized by EHB < 0. The interaction between PANI and amylopectin is thus considered to be physical (with no covalent bonding).

Conclusions
The results of this work demonstrate a facile method for the synthesis of polyaniline colloids. In particular, potato starch was employed as a steric stabilizer to obtain stable dispersions of polyaniline (emeraldine base). Polyaniline/starch biocomposites with various aniline/starch ratios were thus synthesized via an oxidative polymerization and characterized. To evaluate the effect of aniline/starch ratios variation, several tests (UV-vis, FTIR-ATR, Raman, XRD, TGA, DSC), as well as computational calculation (DFT, COSMO-RS, AIM), were performed. This study demonstrates the possibility of using starch as a steric stabilizer, and so to empower the processability of the polyaniline by improving the solubility in DMSO and DMF and the dispersion behavior in water, as well as in glycerol and chloroform, at higher temperatures.
The growth of polyaniline is observed on the starch particles' surfaces. The particle dimension and the degree of crystallinity appear to be strongly affected in the composite material, as demonstrated by the formation of smaller particles and the loss of the original starch crystalline lattice. Although less thermally stable, M11 appears to be the best in terms of morphology and retention of crystallinity degree. The cyclic voltammetry test shows better electro-activity in terms of the electron exchange and capacitive behavior of the composite materials, with the ratio of 1:1 being the best composition among these materials. Finally, the strong interactions between the PANI and amylopectin (the main starch component) were investigated by the DFT, COSMO-RS, and AIM methods, demonstrating their purely physical nature.
These results attest that the prepared composites as promising versatile materials, which are expected to broaden the electrochemical applications of PANI.