Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.
Scientific Reports volume 13, Article number: 20186 (2023 ) Cite this article Biodegradable Equipment
The donor–acceptor (D-A) type of conjugated polymers has emerged as the paradigm of the third generation of electronically conducting polymers demonstrating improved infrared activity and intrinsic electronic conductivity. Judicious selection of donor (D) and acceptor (A) monomers for copolymerization can further fine-tune these properties. Notably, for such refinement, natural compounds provide many conjugated molecules with various functional groups. Berberine cation (Ber+) found in Coscinium fenestratum has extensive conjugation and contains both an electron deficient isoquinolium A moiety and electron-rich D-type methylenedioxy and methoxy groups. The incorporation of natural products in electronic materials is a novel area of research which opens a wide scope for future electronic and optoelectronic devices. Investigation of their fundamental properties via computer simulations is therefore important. In this study, quantum chemical calculations are performed using density functional theory (DFT) to investigate the electronic and optical properties of oligomers of Ber+ and 3,4-ethylenedioxythiophene (EDOT) and to explore the possibilities for homo-polymerization of Ber+ and its copolymerization with EDOT. It has been revealed that homo-polymerization is not favoured but copolymerization with EDOT is possible. As such, Ber+ was copolymerized with EDOT and the copolymers formed by electro-polymerization are extensively characterised and the D-A behaviour of the copolymers verified. Furthermore, the theoretical predictions have been compared with the experimental data.
Abdullah F. Alahmadi, Jingyao Zuo & Frieder Jeckle
Bing Zheng, Lijun Huo & Yongfang Li
Ryo Asakawa, Takuji Seo, … Yutaka Ie
The study of donor–acceptor (D-A) types of electronically conducting polymers has garnered significant attention in recent years due to their potential applications in optoelectronics, solar cells, and other electronic devices1,2,3,4,5,6,7. The D-A behaviour in polymers refers to the interactions between electron-rich and electron-poor units within a polymer. The electron-rich units, known as donors (Ds), may transfer electrons to the electron-poor units, known as acceptors (A)s4. This D-A behaviour plays a crucial role in determining the electronic and optical properties of the polymer in its intrinsic state without deliberate doping. Incorporation of alternating D and A units into a conjugated polymer allows control of the bandgap since the LUMO is mainly localised on the A unit and the HOMO on the D unit. Hence, judicious variation of D and A units controls the HOMO–LUMO separation, a key factor in the optical and redox properties of the materials.The electron-rich units typically consist of conjugated molecules, including homologous polyaromatic compounds and those containing heteroatoms such as nitrogen, sulphur, and oxygen. Conversely, electron-poor units are typically based on molecules with electron-withdrawing groups such as carbonyl, nitrile, quaternary ammonium, and cyano groups.
Berberine, a natural compound found in various plants, has been recognized as a potential donor–acceptor (D-A) material8. It contains a quaternary ammonium ion serving as the acceptor (A) moiety, with methylenedioxy and methoxy groups attached to phenyl rings acting as donor (D) moieties. Berberine exhibits several medicinally important properties, including antimicrobial, anti-inflammatory, and anti-cancer activities9,10,11. While the medicinal applications of berberine have been extensively studied, its utilization in the fabrication of electronically conducting polymers has not been explored thus far. Theoretically, Ber+ has the potential to undergo homo-polymerization, resulting in a poly(Ber+) homopolymer. However, oxidative polymerization may require a high positive potential due to the challenging oxidation of such a complex conjugated ion12. Alternatively, Ber+ can be copolymerized with easily oxidizable conjugated monomers such as EDOT or terthiophene (T3), following similar procedures employed in the synthesis of various complex D-A polymers13,14,15,16,17. Therefore, it is crucial to conduct theoretical investigations to explore the feasibility of copolymerization between Ber+ and EDOT.
However, to the best of our knowledge, there are no reports on either homo-polymerization of Ber+ or its copolymerization with easily oxidizable monomers such as EDOT. In the present context, it is crucial to investigate the donor–acceptor (D-A) behaviour of Ber+, along with its co-oligomers and co-polymers with EDOT. Figure 1A illustrates the D (highlighted in blue) and A (highlighted in green) segments, as well as the extended conjugation (highlighted in red) within the Ber+ cation. Figure 1B shows the electron-rich nature of EDOT.
Chemical structures of (A) berberine cation (Ber+) showing electron rich (blue circled) and electron deficient (green circled) moieties and (B) ethylenedioxythiophene (EDOT) showing lone pairs of electrons capable of increasing the electron density of thiophene moiety.
The strength of D-A interactions can be influenced by several factors, including the polymer's chemical structure, degree of conjugation, nature of the interface between electron-rich and electron-poor units, and steric hindrance. The precise control of these interactions holds significant value in the design and advancement of polymers for diverse applications. Previously we have conducted investigations into the electrochemistry and optical spectroscopy of monomeric berberine cation with a view towards its eventual use in D-A conductive polymers8. In this work, we report DFT calculations on Ber+/EDOT copolymers and an experimental electrochemical investigation of the copolymerization and charge transport in a Ber+/EDOT copolymer.
Computer simulations were executed to analyse the electrochemical and spectroscopic attributes of Ber+ and EDOT, employing the Gaussian 0918 software. The optimized geometries of Ber+/EDOT were determined using DFT with the B3LYP functional and 6-31G(d) basis set. Excited state energies were estimated using TD-DFT at the same level of theory. We have previously ascertained that the B3LYP functional in conjunction with the implicit solvent model (default PCM in Gaussian 09) is sufficient to describe the optical and electrochemical behaviour of monomeric berberine cation8.
Tetrabutylammonium hexafluorophosphate (TBAPF6), procured from Sigma Aldrich, and acetonitrile, obtained from Honeywell—Riedel-de-Haen, were used. The extraction and isolation of Ber+ from Coscinium fenestratum (tree turmeric) and its spectroscopic and electrochemical characterization were performed as previously reported8.
The UV–visible absorption spectra of the samples in acetonitrile and aqueous solution were recorded using a Shimadzu UV-1800 spectrometer (Serial No: A11635305394 CD). To perform scanning electron Microscopy (SEM) and energy dispersive X-ray analysis (EDX), the samples were deposited on fluoride-doped tin oxide (FTO) transparent conducting glass substrate surfaces. The samples were then cleaned with acetone to remove any residual solvent and air dried. The SEM images were obtained with a FEI Quanta 3D Dual Beam Microscope available in the high vacuum mode. The samples were analysed under analytical mode, with a working distance of 10 mm. The EDX and elemental maps were obtained with an X-MaxN 50 spectrometer (Oxford Instruments) mounted on the SEM machine.
Cyclic voltametric studies of Ber+ and EDOT were conducted using a Metrohm PGSTAT204 electrochemical workstation (Serial No: AUT50184), following the methodology outlined in our previous work8. For each experiment, a solution was prepared by dissolving the respective monomer mixtures of Ber+ and EDOT in 25.0 mL of acetonitrile, resulting in concentrations of 5 mmol dm–3 in Ber+ and 10 mmol dm–3 in EDOT. To this solution, a tetrabutylammonium hexafluorophosphate (TBAPF6) background electrolyte (BGE) was added at a concentration of 0.10 mol dm–3. The electro-polymerization and cyclic voltametric (CV) measurements were performed using a one-compartment cell equipped with three electrodes. The working electrode (WE) was either a glassy carbon (GC) or an FTO, the reference electrode (RE) was a saturated calomel electrode (SCE), and a Pt-wire served as the counter electrode (CE). Before each experiment, the solution was degassed by purging with high purity nitrogen gas for 20 min, and a slow flow of nitrogen was maintained above the solution to prevent the re-entry of air. Electro-copolymerization was carried out by performing 25 repetitive CV cycles within the potential range of 0.0 V to + 1.8 V. The resulting polymer films on the WE surface were washed with acetone to remove any residual organic matter, and the CVs were recorded in the neat BGE. The CV data were recorded in the potential range of − 2.0 V to + 2.0 V to encompass all redox peaks of the materials, and a scan rate of 100 mV s–1 was employed, unless otherwise specified.
The Nyquist and Bode plots portraying alternating current (AC) impedance properties were generated over a frequency spectrum spanning 0.1 Hz to 1 MHz. These measurements were conducted at distinct direct current (DC) potential biases, selectively chosen, within the potential range of − 2.0 V to 2.0 V. This range comprehensively encompassed the pertinent information regarding the investigated electrochemical system.
Both DFT and TD-DFT calculations were performed to investigate the electronic structure of berberine dimers as an intermediate in the homopolymerisation of berberine. The structure of the dimer was optimised, and then vertical excitation energies were calculated to the lowest excited singlet (S1) and triplet (T1) states; the data is summarised in Table 1 alongside a comparison of the results for the monomeric berberine cation8. The HOMO and LUMO energies are lower in the dimer, owing to the electrostatic effect, but overall, the HOMO–LUMO gap increased by about 0.095 eV. Nevertheless, the energy of the S1 state is decreased in the dimer by 0.42 eV and the S1-S0 gap in the monomer has been shown to correlate well with the electrochemical gap8. The importance of this lowering lies in the mechanism of electropolymerisation, which for a typical conjugated monomer requires that the oligomers are oxidised at the same or a lower potential than the monomer19,20. However, this is insufficient on its own because it is necessary to consider the geometry of the dimer; the optimised geometry of the berberine dimer is shown in Fig. 2.
(A) Cyclic voltammetry of berberine in acetonitrile application of − 2.0 V to + 2.0 V with 100 mV s–1 scan rate. The working electrode was a 0.50 cm diameter glassy carbon disc, the reference was SCE and the concentration of berberine was 3.0 mM. The solution was purged with N2 to remove dissolved oxygen. The direction of change with scan number is to a decrease in the 1.65 V and 1.86 V peaks. (B) Ber+–Ber+ optimized using the B3LYP/6-31G(d) model chemistry. Planes of two Ber+ moieties in the Ber+-dimer are shown in purple and green.
Our attempts to homopolymerize Ber+ in an acetonitrile medium, with TBAPF6 as the BGE, did not yield any deposition on the surface of the working electrode and the currents in consecutive CVs decreased. Figure 2A shows the CV of Ber+Cl– in the potential range from – 2.0 to + 2.0 V.
The CV of Ber+ in acetonitrile containing TBAPF6 given in Fig. 2A shows three distinctive oxidation peaks centred at + 1.21 V, + 1.65 V, and + 1.86 V and the latter two are due to the oxidation of Ber+ at high positive potentials resulting in compounds formed as shown in Fig. 3A. The first peak is characteristic of \({\mathrm{Cl}}^{-}\) to Cl2 oxidation as shown in our previous publication8. The repeated CVs did not yield any polymer deposited on the WE surface, and the current densities declined in the consecutive repeated scans. As such, it is clear that the oxidation of Ber+ does not result in electro-polymerization to yield a conductive poly(Ber+). The steric hindrance between two Ber+ cations prevent the normal coplanar arrangement of rings in conductive polymers as shown in Fig. 2B; we argue that this factor is partly responsible for the lack of success in achieving homopolymerisation of berberine. We must also consider whether the one-electron oxidation product of berberine undergoes other chemical reactions in competition with oligomerisation.
(A) Possible oxidized products (OX2 and OX3) of Ber+ in the potential range from 1.6 V to 2.0 V region. (B) Calculated IR spectra of Ber+ (blue), Oxidized product 2 (OX2 Cal–brown), Oxidized product 3 (OX3 Cal–red) and experimental FT-IR spectra of Oxidized product 2 (OX2 Exp–purple), Oxidized product 3 (OX3–green). IR frequencies of two oxidized products were calculated using acetonitrile solvent model and 6-31G(d) basis set. (a.u.–arbitrary units).
Based on FT-IR data, we uncover evidence supporting the structural resemblance between these compounds in our comparative analyses of calculated OX221 and OX3 out of 12 structures and (Supplementary information SI Figs. 03, 04) and experimental OX2, OX3 oxidation products12. Notably, the FT-IR bands for several functional groups are nearly identical, indicating that they share a similar bonding nature. For the second oxidation product OX2, wavenumbers align for C–O stretching (–C–O–CH3) (OX2 Cal: 1181 cm–1, OX2 Exp: 1183 cm–1) 22, –C–N– stretching (OX2 Cal: 1359 cm–1, OX2 Exp: 1334 cm–1), –C–H bending (OX2 Cal: 1404 cm–1, OX2 Exp: 1392 cm–1)23, =C–H stretching (OX2 Exp: 1472 cm–1, OX2 Cal: 1486 cm–1)24, and –C–H3 stretching vibrations (OX2 Cal: 3063 cm–1, 3071 cm–1, and OX3: 3076 cm–1), further affirming structural similarity with an additional band centred at 1778 cm–1 indicating the presence of a carbonyl group at some position in the oxidised product OX 2 The most likely position for oxidation to form a carbinyl group is the methylenedioxy carbon atom that has also been suggested by the Referee 1 also. We, therefore, propose the structure OX2 given in Fig. 3A as the most likely oxidation product appearing at + 1.65 V wrt Ag/AgCl/KCL(sat.) electrode. This may have been formed by the reaction with traces of oxygen present in the solution although it was purged with nitrogen gas for about 30 min. This is further evidenced by the observation of the increased peak current at + 1.65 V when the solution is purged with air. For the third product, OX3, an additional band centred at 3627 cm-1 in the experimental spectrum is due to –OH group stretching at which appears at 3416 cm-1 for the OX3 in the calculated spectrum22,25, OX 3 could be formed when OX2 is reacted with traces of water present in the solution. Therefore, the possible structure for the oxidation product appearing at + 1.68 V wrt Ag/AgCl/KCl(sat.) is the OX 3 in Fig. 3A.
We calculated the HOMO and LUMO energy levels of Ber+ and EDOT in three distinct media: vacuum, water, and acetonitrile. The primary aim of the theoretical study was to determine the solvent medium that yields the smallest HOMO–LUMO energy difference for Ber+ and EDOT in combination. The obtained results serve to provide valuable insights into the electronic properties of Ber+ and EDOT within different solvents, aiding in the comprehension of its behaviour across various media. The outcomes of these calculations are presented in Table 2.
The data in Table 2 indicate that although the HOMO energies in Ber+ and EDOT are similar in polar media, the LUMO energy of Ber+ is much lower by about 2.5 eV than that of EDOT. The calculated electrostatic potential (ESP) for the ground state of EDOT was calculated at the B3LYP/6-31G(d) level. Figure 4A depicts the analysis of high electron density contours, which show that carbon atoms numbered 2 and 5 exhibit a greater inclination to donate electrons to A molecules like in Fig. 4B. This can be attributed to the presence of two lone pairs on the S atom that are bonded to these carbon atoms.
(A) Distribution of electrostatic potential map (ESP) map on the EDOT (positive contours in blue and negative contours in red). (B) Spin density diagram of berberine (Ber2+) cation.
In our previous studies, we found that if 1:1 stoichiometry of the two monomers A and B is used, the electro-co-polymerization leads to the formation of an alternating copolymer, (A-B)n. Additionally, if x:y stoichiometric ratio is used the resulting polymer is a copolymer [Ax-By]n4,14,26,27. However, it is important to carefully control reaction conditions and other relevant parameters to understand the behaviour of the cations and the reaction mechanism. It is well established that electro-polymerization of EDOT leading to PEDOT undergoes via 2 and 5 positions (two ortho positions to S in the thiophene unit). However, several possible regioisomers of Ber+-EDOT co-oligomers are possible (Fig. 5A). Below we use A, B and C to denote the pattern of coupling of Ber+ and EDOT irrespective of the EDOT/Ber+ ratio n.
(A) The proposed monomer structures of EDOT-(Ber+)-EDOT when Ber+ and EDOT reacted in 1:2 molar ratio. (B) Simulated absorption spectra of three possible (EDOT-Ber+-EDOT) monomers. Model chemistry B3LYP/6-31G(d) was used with vacuum (dash line) and acetonitrile (solid line) solvent models.
To find out energetically the most favourable regioisomers, DFT calculations were performed to investigate the interaction between Ber+ and EDOT and the data are summarised in Tables 3 and 4. The removal of a hydrogen atom from carbon number 5 in isomer A is impeded by the presence of the electron-rich oxygen atom and electron cloud in the O-CH3 group results in steric hindrance, which complicates the combination of the electron rich EDOT molecule with carbon 5. The simulation data shows that the isomer C has the lowest overall electronic energy, with B lying 0.19 eV higher and isomer A at 0.70 eV. O this basis isomer C would be favoured, but other factors must be considered.
The simulated optical spectra of isomers A, B, and C are compared in Fig. 5B. When the acetonitrile solvent model is considered, the S1 absorption bands are observed centred at 513 nm (A), 458 nm (B), and 466 nm (C) in acetonitrile. This suggests isomer A may be most easily oxidised assuming a correlation between electrochemical and optical gap holds.
Figure 6B shows changes in the orbital energies of the 10 highest occupied and 10 lowest unoccupied orbitals of (EDOT)n(Ber+)(EDOT)n for n = 1, 2 and 3. As expected, it highlights that the difference between the HOMO and LUMO energies decreases with increasing monomer count, resulting in a reduction of the energy difference between the HOMO and LUMO levels. This data does not however clearly favour any regioisomer and the details of the mechanism require consideration.
(A) The spin densities (blue/green) of structure (a), (b) and (c) for (EDOT)n(Ber+)(EDOT)n for n = 2; n = 3 and n = 4. All calculations were performed at the geometry optimized structure using B3LYP/6-31G(d) with implicit solvation by acetonitrile. (B) The orbital energies of the 10 highest occupied and 10 lowest unoccupied orbitals of (EDOT)n(Ber+)(EDOT)n for n = 1, 2 and 3. (C) Proposed mechanism of electro-copolymerization. Top. The reaction is shown as a coupling of EDOT cation radicals, one already part of a growing chain anchored to the electrode surface on the left of the diagram. Middle. Reoxidation of the newly added EDOT. Bottom. Further addition of radical cations is possible if spin density is significant on the end of the growing chain.
The HOMO, LUMO energies and the gaps for (EDOT)n(Ber+)(EDOT)n species for the three structures A, B and C, for n = 1, 2, 3 and 4 are given in Table 5. The maps of spin density in these structures are given in Fig. 6A. These simulations can be interpreted in the context of Fig. 6C to provide evidence that the most promising mole ratio for copolymerization is 2:1, EDOT:Ber+ in isomer B.
Figure 6C shows a schematic illustration of the copolymerization of Ber+ and EDOT at the end of a growing polymer chain that is anchored (on the left of the diagram) to the electrode surface. Coupling of new radicals to the growing chain is feasible if there is significant spin density at the end of the chain far from the electrode. Figure 6A shows the spin densities for the one-electron oxidation products of species (EDOT)n(Ber+)(EDOT)n for n = 2, 3 and 4 and isomers A, B and C of Fig. 5A In all cases, except one, the spin density is localized on one side of the Ber+ unit. This suggests that these systems will not readily polymerise because the growing ends of the polymer chain remote from the electrode will not acquire spin density to couple to further radical cations. However, isomer B with n = 2 does show the spin density on both sides of the Ber+ unit. We interpret this effect in terms of a mismatch between the energies of the EDOT and Ber+ orbitals – as the EDOT chain lengthens, the orbital energies change. This series of calculations suggests that the ratio EDOT:Ber+ should not be too large for successful copolymerisation and that isomer B is the most likely product.
Figure 7 illustrates the CVs acquired during the attempted electro-polymerization process of Ber+ with EDOT at 1:2 molar ratio. Figure 7A exhibits the repetitive CVs that depict the successive electro-polymerization leading to the formation of a polymer film on the surfaces of GC electrodes. The increasing current observed in consecutive CVs in Fig. 7A indicates the deposition of conducting polymer films onto the GC working electrode surface. In Fig. 7B, the CV obtained for the suspected polymer in the BGE without monomers is displayed. The CVs of the polymers shown in Fig. 7B provide clear evidence of the incorporation of the EDOT and Ber+ in the copolymer where the presence of EDOT is indicated by an oxidation peak at approximately + 0.80 V during the forward scan, followed by a corresponding reduction peak at around + 0.70 V during the reverse scan. Additionally, an oxidation peak at approximately + 1.8 V is observed for the first oxidation step of Ber+.
(A) Cyclic voltammograms (CVs) of 5.0 mM berberine cation (Ber+) and 10.0 mM EDOT in 0.10 M TBAPF6 solution between 0 and 2.0 V potential range, and (B) CV of the polymer in 0.10 TBAPF6 solution (neat background electrolyte without the monomers) between − 2.0 V–2.0 V potential range. (C) CVs of 10.0 mM PEDOT in 0.10 M TBAPF6 in acetonitrile solution between 0 V–2.0 V potential range. (D) CV of the polymer in 0.10 TBAPF6 acetonitrile solution between − 2.0 V–2.0 V potential range. (E) UV–Visible absorption spectra of Ber+-EDOT co-polymer and PEDOT. (F) Spectra in (E) in the range from 270 to 450 nm showing π \(\to \) π* and n \(\to \) π* transitions.
The CVs of the PEDOT in Fig. 7C and D show their typical capacitive behaviour if recorded within 0 V to + 1.0 V with typical oxidation peak centred at + 0.25 V in the forward scan with its reduction peak at − 0.41 V in the reverse scan. There are hardly any currents in the negative potential range from – 1.0 V to – 2.0 V wrt Ag/AgCl/KCl(sat.) electrode. PEDOT gets over-oxidized above + 1.50 V as we see from the sharp rise in currents above + 1.53 V. This CV is typical of PEDOT as it is conducting at positive potentials due to the formation of polarons and biplarons. When the potential is reversed, it takes time to neutralise the charge carriers formed and currents are observed down to – 0.7 V beyond which it is not conducting as there are no appreciable currents.
The absorption spectra, shown in Fig. 7E, have distinct differences. Both show absorption in the entire region 270 nm to 900 nm range that we could measure indicating that the materials are not simple molecules or oligomers. The fact that both materials show high absorptions in the entire range indicates that they are electronically conducting polymers. PEDOT shows π \(\to \) π* transition at 282 nm while the same of the copolymer is centred at 299 nm with a 17 nm bathochromic shift. The lowering of band gap due to the localisation of HOMO and LUMO levels at D and A units, respectively, gives this bathochromic shift of the π \(\to \) π* transition as shown in Fig. 7F. Interestingly, the n \(\to \) π* transition of the copolymer appearing at 351 nm shows a small hypsochromic shift of 6 nm from that of PEDOT (357 nm) though the maximum of the absorption band is difficult to locate. If there is such a small hypsochrmic shift it may be due to the decrease of n-electron concentration in the conjugated chain due to the presence of some Ber+ units instead of EDOT units in the copolymer. However, note that EDOT has 8 non-bonded electrons in its ethylenedioxy moiety while Ber+ has 16 of them in methylenedioxy and methoxy groups. However, the conjugated chain has electron deficient N+ making the n-electron density calculation difficult. The overall effect could be the reason for small difference observed in \(n \to \) π* transitions of the two materials. Additionally, the low energy polaron and bipolaron bands are clearly observed in both materials.
Electrochemical analysis provided initial insights into the polymer composition, which were further supported by SEM–EDX analysis. SEM images (Fig. 8) reveal the morphology of the polymer films on FTO glass. In comparison, SEM images of bare FTO showed grain structures of FTO, which facilitated epitaxial polymer growth and resulted in a thermodynamically favoured cloud-like structure. The atomic ratios obtained from SEM–EDX analysis were complex; however, elemental mapping at the macroscopic scale exhibited acceptable levels of carbon, nitrogen, oxygen, and sulphur, providing evidence in support of copolymerization.
Scanning electron microscopy images (SEM) of (A) Berberine with magnification 1300x (B) the suspected polymer with magnification 7000x.
Tables 6 and 7 present the theoretically calculated and experimentally obtained atomic percentages of berberine cation and the (EDOT-Ber+-EDOT)n copolymer. These tables provide detailed information regarding the atomic composition of the copolymer, further supporting the successful formation of the desired polymer structure.
Data in Table 6 show some differences between the theoretical and empirical atomic percentages of Ber+Cl-. This is because the SEM-EDAX is incapable of measuring element H and hence calculating atomic percentages with respect to the measured elements. When the atomic percentages were calculated excluding H atoms a somewhat closer agreement can be observed. However, empirical results show good agreement for N and Cl atomic percentages when H is included.
As can be seen from the data, the presence of C, N, O, S, P and F suggests that both EDOT and Ber+PF6– are incorporated in the polymer. Since the polymer was grown from a solution containing large excess of TBAPF6, it is possible to incorporate PF6– ion as the counter ion and the chloride ion remaining in the solution. As explained, SEM-EDAX is incapable of measuring element H and hence the atomic percentages determined are based on those measured. When we consider the possibility for oxidation of Ber+ forming structures given in Fig. 3, at positive potentials, the polymer may contain these species also thus increasing atomic percentage of O. Since the polymer is electronically conducting there may be more than one counter ions ([PF6]–) in the repeat unit as EDOT may exist as EDOT.+ in it. Considering these possibilities theoretical atomic percentages were calculated for all three cases with and without including H atoms. The best correlation is obtained for Theoretical (EDOT-Ber+PF6–EDOT)n excluding H. It is, therefore, possible that the 1:2 stoichiometry is favoured in the polymer because EDOT is the one that is easiest to oxidize forming cation radical EDOT.+ (at around + 0.8 V), and it can then combine with a radical on the growing polymer chain. However, Ber+ is difficult to oxidize to form Ber2+.radical. The possible way forward for propagation of the reaction is, therefore, to react EDOT.+-Ber+ with another EDOT.+ forming EDOT.+- -Ber+- EDOT.+ as the repeat unit giving a copolymer with EDOT:Ber+ stoichiometry of 1:2 as suggested in Fig. 6.
Impedance spectra were recorded at different applied potentials to evaluate the electronic and ionic conductivity of poly(Ber+-EDOT). Nyquist [–Im(Z) vs Re(Z)] plots obtained at selected potentials of − 0.8 V and 0.0 V for the 2:1 co-polymer is given in Fig. 9. The full dataset is shown in the SI. The impedance spectra are presented with logarithmic axes scales in order to display clearly their main features; a high frequency semi-circle and a low frequency region in which the impedance is either capacitive (Fig. 9A) or has the form of a flattened semi-circle (Fig. 9B). We analysed these spectra using the dual-rail transmission line model28 shunted by a capacitor (C2) to describe the low frequency region and a parallel RctC1 circuit to describe the high frequency region (Fig. 9A).
Impedance spectra for poly(Ber+-EDOT) films over the frequency range 1 MHz–1 Hz. The electrolyte was 0.1 mol dm–3 TBAPF6. (A) dc potential of 0.0 V; (B) dc potential of − 0.8 V. The symbols are experimental data, and the solid line is the least squares fit of the impedance of the equivalent circuit of Fig. 10A. (C) Nyquist plots of range − 1.8 V–1.8 V. (D) Experimental dc potentials of range − 1.8 V–0.2 V.
The impedance spectrum defined by the equivalent circuit of Fig. 10A was fitted to the experimental data by the complex least squares method. Three main regions are observed in these spectra and denoted (i), (ii) and (iii) in the discussion. (i) a semicircle at large frequency (R about 200 Ω); (ii) an incomplete semicircle at lower frequency and (iii) in the equivalent circuit (blue) the impedance eventually becomes capacitive at the very lowest frequency, but the data does not.
(A) Equivalent circuit used to fit the impedance spectra of poly(Ber+-EDOT). DRTL indicates a dual-rail transmission line with resistive rails (common resistance, R) representing the electronic and ionic conduction pathways and distributed capacitance (value not shown). (B) lnR against dc potential. Ru is the uncompensated solution resistance and was found to be of the order of 50 Ohm.
We observed that at potentials above 0.2 V, the quality of the fit is poor; this is the region of the cyclic voltammogram (Fig. 7B) in which the polymer is p-doped, and its impedance is small. At positive potentials, there is also the possibility of oxidation of berberine8 and the simple equivalent circuit of Fig. 10A cannot model such a faradaic process. The fit of this model to the data for E < 0.2 V is satisfactory except at very low frequencies (iii) where the impedance spectrum is strongly influenced by background faradaic processes and deviates from the predicted pseudo-capacitive behaviour of any model based on a finite thickness polymer film and no background faradaic process. Nevertheless, we can fit the two main features of interest (i), (ii) and extract an effective charge transfer resistance (Rct) for the high-frequency semicircle and the resistance of the rails of the transmission line (R) for the low-frequency data.
The resistance Rct extracted from the high frequency semi-circle is relatively insensitive to the applied potential, with no clear trend and a value of 210 ± 63 Ohm almost independent of potential. Rct is sensitive to the background electrolyte concentration and at -0.6 V it has a value of 207 Ohm, 327 Ohm and 416 Ohm in concentrations of 0.1, 0.05 and 0.03 mol dm-3 TBAPF6. On this basis, we do not assign Rct as a true charge transfer resistance but suggest that it is an effective resistance associated with the high-frequency response of the cell and electronics – the associated RctC1 time constant is of the order of 0.5 μs.
The second, lower frequency, incomplete semi-circle (region (ii)) contains information on charge transport in the polymer, which we represent by a transmission line element (Fig. 10A). The dual-rail transmission line is characterized by the resistances of the two rails, representing the electronic and ionic conduction and the distributed capacitance29. We found that a common value of resistance (R) for both rails was sufficient to fit the experimental data. The apparent equality of the resistance of electronic and ionic resistances has also been observed by previous workers30. Figure 10B shows the variation of R with applied potential over the range 0.2 V to − 1.8 V. There is a clear increase in R from a value of about 40 Ohm to one of about 3 kOhm as the potential becomes more negative than − 0.2 V. This is the region of the cyclic voltammogram (Fig. 7B) at which the polymer is nominally undoped. We estimate the transition potential to be about − 0.3 ± 0.05 V. In the case of pure poly(EDOT)31, the impedance spectrum shows a large semicircle in the undoped region whose diameter increases by several orders of magnitude as the potential becomes more negative. In contrast, although R increases below − 0.2 V in poly(Ber+-EDOT) it reaches an approximately constant value of the order of 3 kOhm from − 0.4 V down to the most negative potentials studied. We interpret this non-zero conductivity in the absence of oxidative doping as evidence of some self-doping of poly(Ber+-EDOT) because of the presence of its electron-accepting isoquinolinium moiety.
We find that direct polymerization of berberine does not occur, in agreement with calculations and previous works on the oxidation chemistry of berberine cation12. However, copolymerization of berberine from solutions containing mixtures of EDOT and berberine is possible. The cyclic voltammetry of the polymerization media shows the expected growth in current as the polymer film forms on the electrode surface. Voltammetry of the film in the absence of monomer shows both the broad waves expected of EDOT doping/undoping and additional peaks at positive and negative potentials which are characteristic of berberine8. DFT calculations suggest that the most likely structure of the copolymer is that formed by attaching EDOT units via their 2 or 5 positions to C atoms 15 and 18 of the berberine cation for (EDOT)n-Ber+-(EDOT)n type polymers and C atoms 1 and 15 for (EDOT-Ber+-EDOT)n type polymers. Impedance spectroscopy provides evidence of some residual conductivity (“self-doping”) in the potential region in which the polymer is nominally undoped according to the cyclic voltammetry.
The datasets used and/or analyzed during the current study are available from BRH, HMNPG upon reasonable request.
Liu, H. et al. Electrically conductive coordination polymers for electronic and optoelectronic device applications. J. Phys. Chem. Lett. 12, 1612–1630 (2021).
Article CAS PubMed Google Scholar
Wu, X., Fu, W. & Chen, H. Conductive polymers for flexible and stretchable organic optoelectronic applications. ACS Appl. Polym. Mater. 4, 4609–4623 (2022).
Dauzon, E. et al. Pushing the limits of flexibility and stretchability of solar cells: A review. Adv. Mater. 33, 2101469 (2021).
Rajapakse, R. M. G. et al. Implementing the donor–acceptor approach in electronically conducting copolymers via electropolymerization. RSC Adv. 12, 12089–12115 (2022).
Article ADS CAS PubMed PubMed Central Google Scholar
Murad, A. R., Iraqi, A., Aziz, S. B., Abdullah, S. N. & Brza, M. A. Conducting polymers for optoelectronic devices and organic solar cells: A review. MDPI Polym. 12, 2627 (2020).
Moliton, A. & Hiorns, R. C. Review of electronic and optical properties of semiconducting π-conjugated polymers: Applications in optoelectronics. Polym. Int. 53, 1397–1412 (2004).
Khokhar, D., Jadoun, S., Arif, R. & Jabin, S. Functionalization of conducting polymers and their applications in optoelectronics. Polym. Plast. Technol. Mater. 60, 463–485 (2020).
Rajapakse, R. M. G. et al. Berberine isolation from Coscinium fenestratum: Optical, electrochemical, and computational studies. RSC Adv. 13, 17062–17073 (2023).
Article ADS CAS PubMed PubMed Central Google Scholar
Khan, S. et al. A review of the berberine natural polysaccharide nanostructures as potential anticancer and antibacterial agents. Biomed. Pharmacother. 146, 112531 (2022).
Article CAS PubMed Google Scholar
Kumar, A. et al. Current knowledge and pharmacological profile of berberine: An update. Eur. J. Pharmacol. 761, 288–297 (2015).
Article ADS CAS PubMed Google Scholar
Malhotra, B. et al. Recent advances on Berberis aristata emphasizing berberine alkaloid including phytochemistry, pharmacology and drug delivery system. J. Herb. Med. 27, 100433 (2021).
Skopalová, J. et al. Electrochemical oxidation of berberine and mass spectrometric identification of its oxidation products. Bioelectrochemistry 87, 15–20 (2012).
Karunathilaka, D. et al. Correlation of solid-state order to optoelectronic behavior in heterocyclic oligomers. CrystEngComm 24, 4564–4572 (2022).
Sparks, N. E. et al. Electrochemical copolymerization of isoindigo-based donor-acceptor polymers with intrinsically enhanced conductivity and near-infrared-II activity. ChemElectroChem 7, 3752–3760 (2020).
Ranathunge, T. A. et al. Designing hierarchical structures of complex electronically conducting organic polymers via one-step electro-polymerization. J. Mater. Chem. C Mater. 8, 5934–5940 (2020).
Rajapakse, R. M. G. et al. Advances in electro-copolymerization of NIR emitting and electronically conducting block copolymers. J. Mater. Chem. C Mater. 7, 3168–3172 (2019).
Ranathunge, T. A. et al. Radically accessing D-A type ambipolar copolymeric materials with intrinsic electrical conductivity and visible-near infrared absorption via electro-copolymerization. Macromol. Chem. Phys. 220, 1970042 (2019).
Gaussian 09, Revision D.01., Frisch MJ, GW Trucks, HB Schlegel, Scuseria GE, Robb MA, Cheeseman JR, Scalmani G, Barone V, Petersson GA, Nakatsuji H, Li, X, M, Caricato, 1999; [ PubMed ] Marenich A, BG Janesko, R Gomperts, B Mennucci, HP Hratchian, JV Ortiz, AF Izmaylov, JL Sonnenberg, D Williams-Young, F Ding, F Lipparini, F Egidi, J Goings, B. Peng, A. Petrone, T. Henderson, D. Ranasinghe, VG Zakrzewski, J. Gao, N. Rega, G. Zheng, W. Liang, M. Hada, M. Ehara, K. Toyota, , Fukuda R , Hasegawa J , Ishida M , Nakajima T , Honda Y , Kitao O , Nakai H , Vreven T , Throssell K , Montgomery JA , Peralta JE , Ogliaro F F , M Bearpark , JJ , Heyd , E Brothers , Kudin KN , Staroverov VN , Keith T , Kobayashi R , Normand K , Raghavachari K , Rendell A , Burant JC , Iyengar SS , Tomasi J , Cossi M , . Millam JM, Klene M, Adamo C, Cammi R, Ochterski JW, Martin RL, Morokuma K, Farkas O, Foresman JB, and DJ Fox, Gaussian, Inc., Wallingford CT (2016).
Genies, E. M., Bidan, G. & Diaz, A. F. Spectroelectrochemical study of polypyrrole films. J. Electroanal. Chem. Interfacial Electrochem. 149, 101–113 (1983).
Waltman, R. J. & Bargon, J. Electrically conducting polymers: A review of the electropolymerization reaction, of the effects of chemical structure on polymer film properties, and of applications towards technology. Can. J. Chem. 64, 76–95 (2011).
Diculescu, V. C., Enache, T. A., Oliveira, P. J. & Oliveira-Brett, A. M. Electrochemical oxidation of berberine and of its oxidation products at a glassy carbon electrode. Electroanalysis 21, 1027–1034 (2009).
Bashmakova, N., Kutovyy, S., Zhurakivsky, R., Hovorun, D. & Yashchuk, V. Vibrational spectra of berberine and their interpretation by means of DFT quantum-mechanical calculations. Ukr. J. Phys. 56, 130–137 (2011).
Thusa, R. & Mulmi, S. Analysis of phytoconstituents and biological activities of different parts of Mahonia nepalensis and Berberis aristata. Nepal J. Biotechnol. 5, 5–13 (2017).
Comincini, S. et al. Development of Berberine-loaded nanoparticles for astrocytoma cells administration and photodynamic therapy stimulation. Pharmaceutics 15, 1078 (2023).
Article CAS PubMed PubMed Central Google Scholar
Leona, M. & Lombardi, J. R. Identification of berberine in ancient and historical textiles by surface-enhanced Raman scattering. J. Raman Spectrosc. 38, 853–858 (2007).
Article ADS CAS Google Scholar
Higgins, S. J., Lovell, K. V., Rajapakse, R. M. G. & Walsby, N. M. Grafting and electrochemical characterisation of poly-(3,4-ethylenedioxythiophene) films, on Nafion and on radiation-grafted polystyrenesulfonate–polyvinylidene fluoride composite surfaces. J. Mater. Chem. 13, 2485–2489 (2003).
Velauthamurty, K. et al. Synthesis and characterization of monomeric and polymeric Pd(II) and Pt(II) complexes of 3,4-ethylenedioxythiophene-functionalized phosphine ligands. J. Mater. Chem. 19, 1850–1858 (2009).
Albery, W. J. & Mount, A. R. Dual transmission line with charge-transfer resistance for conducting polymers. J. Chem. Soc. Faraday Trans. 90, 1115–1119 (1994).
Albery, W. J. et al. Spectroscopic and electrochemical studies of charge transfer in modified electrodes. Faraday Discuss. Chem. Soc. 88, 247–259 (1989).
Albery, W. J. & Mount, A. R. A second transmission line model for conducting polymers. J. Electroanal. Chem. Interfacial Electrochem. 305, 3–18 (1991).
Sundfors, F., Gustafsson, H., Ivaska, A. & Kvarnström, C. Characterisation of the aluminium-electropolymerised poly(3,4- ethylenedioxythiophene) system. J. Solid State Electrochem. 14, 1185–1195 (2010).
This material is based upon work supported by the United States Air Force Office of Scientific Research/AOARD funding under award number FA2386-21-1-4096. Computational facility was provided by the Department of Chemistry, Faculty of Science, University of Peradeniya, Sri Lanka and Newcastle University, UK. This research also made use of the Rocket High Performance Computing service at Newcastle University.
Department of Chemistry, University of Peradeniya, Peradeniya, 20400, Sri Lanka
RM Gamini Rajapakse, HMNP Gunarathna, AU Malikaramage, MGSAMEWDDK Egodawele, WHMRNK Herath, Lahiru Sandakelum, VMYSU Bandara, WVNS Bowatta, JM Susanthi Jayasinghe, VN Seneviratne, Udayana Ranatunga & LLK Perera
School of Natural and Environmental Sciences, Newcastle University, Newcastle Upon Tyne, NE1 7RU, UK
Department of Decision Sciences, University of Moratuwa, Katubedda, Moratuwa, Sri Lanka
Department of Physical Science and Technology, Faculty of Applied Sciences, Sabaragamuwa University of Sri Lanka, Belihuloya, Sri Lanka
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
R.M.G.R. and B.R.H. supervised the work, raised funding and designed the study. The electrochemical studies, impedance, IR and optical characterisation was carried out by A.U.M., H.M.N.P.G., M.G.S.A.M.E.W.D.D.K.E., W.H.M.R.N.K.H., L.S., V.M.Y.S.U.B., W.V.N.S.B., J.M.S.J., U.R., V.N.S., S.M.D., L.L.K.P. and C.P.U. The DFT and TD-DFT computational studies were performed by B.R.H. and H.M.N.P.G. The manuscript was prepared by R.M.G.R., B.R.H. and H.M.N.P.G. All coauthors edited and approved the manuscript.
Correspondence to R. M. Gamini Rajapakse.
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Gamini Rajapakse, R.M., Horrocks, B.R., Gunarathna, H.M.N.P. et al. Computational analysis and experimental verification of donor–acceptor behaviour of berberine, and its co-oligomers and co-polymers with ethylenedıoxythıophene. Sci Rep 13, 20186 (2023). https://doi.org/10.1038/s41598-023-47541-7
DOI: https://doi.org/10.1038/s41598-023-47541-7
Anyone you share the following link with will be able to read this content:
Sorry, a shareable link is not currently available for this article.
Provided by the Springer Nature SharedIt content-sharing initiative
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.
Scientific Reports (Sci Rep) ISSN 2045-2322 (online)
Cracking System Small Test Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.