Charge carrier dynamics in conducting polymer PEDOT using ab initio molecular dynamics simulations

As conducting polymers become increasingly important in electronic devices, understanding their charge transport is essential for material and device development. Various semi-empirical approaches have been used to describe temporal charge carrier dynamics in these materials, but there have yet to be any theoretical approaches utilizing ab initio molecular dynamics. In this work, we develop a computational technique based on ab initio Car-Parrinello molecular dynamics to trace charge carrier temporal motion in archetypical conducting polymer poly(3,4-ethylenedioxythiophene) (PEDOT). Particularly, we analyze charge dynamics in a single PEDOT chain and in two coupled chains with different degrees of coupling and study the effect of temperature. In our model we first initiate a positively charged polaron (compensated by a negative counterion) at one end of the chain, and subsequently displace the counterion to the other end of the chain and trace polaron dynamics in the system by monitoring bond length alternation in the PEDOT backbone and charge density distribution. We find that at low temperature (T = 1 K) the polaron distortion gradually disappears from its initial location and reappears near the new position of the counterion. At the room temperature (T = 300 K), we find that the distortions induced by polaron, and atomic vibrations are of the same magnitude, which makes tracking the polaron distortion challenging because it is hidden behind the temperature-induced vibrations. The novel approach developed in this work can be used to study polaron mobility along and between the chains, investigate charge transport in highly doped polymers, and explore other flexible polymers, including n-doped ones.

A computational study of cellulose regeneration: Coarse-grained molecular dynamics simulations

Understanding the microscopic mechanisms of regeneration of cellulose is prerequisite for engineering and controlling its material properties. In this paper, we performed coarse-grained Martini 3 molecular dynamics simulations of cellulose regeneration at a scale comparable to the experiments. The X-ray diffraction (XRD) curves were monitored to follow the structural changes of regenerated cellulose and trace formation of cellulose sheets and crystallites. The calculated coarse-grained morphologies of regenerated cellulose were backmapped to atomistic ones. After the backmapping we find that the regenerated coarse-grained cellulose structures calculated for both topology parameters of cellulose I_β and cellulose II/III, are transformed to cellulose II, where the calculated XRD curves exhibit the main peak at approximately 20-21 degrees, corresponding to the (110)/(020) planes of cellulose II. This result is in good quantitative agreement with the available experimental observations.

Exploiting mixed conducting polymers in organic and bioelectronic devices

Efficient transport of both ionic and electronic charges in conjugated polymers (CPs) has enabled a wide range of novel electrochemical devices spanning applications from energy storage to bioelectronic devices. In this Perspective, we provide an overview of the fundamental physical processes which underlie the operation of mixed conducting polymer (MCP) devices. While charge injection and transport have been studied extensively in both ionic and electronic conductors, translating these principles to mixed conducting systems proves challenging due to the complex relationships among the individual materials properties. We break down the process of electrochemical (de)doping, the basic feature exploited in mixed conducting devices, into its key steps, highlighting recent advances in the study of these physical processes in the context of MCPs. Furthermore, we identify remaining challenges in further extending fundamental understanding of MCP-based device operation. Ultimately, a deeper understanding of the elementary processes governing operation in MCPs will drive the advancement in both materials design and device performance.

Tuning of the elastic modulus of a soft polythiophene through molecular doping

Molecular doping of a polythiophene with oligoethylene glycol side chains is found to strongly modulate not only the electrical but also the mechanical properties of the polymer. An oxidation level of up to 18% results in an electrical conductivity of more than 52 S cm^-1 and at the same time significantly enhances the elastic modulus from 8 to more than 200 MPa and toughness from 0.5 to 5.1 MJ m^-3. These changes arise because molecular doping strongly influences the glass transition temperature T_g and the degree of π-stacking of the polymer, as indicated by both X-ray diffraction and molecular dynamics simulations. Surprisingly, a comparison of doped materials containing mono- or dianions reveals that - for a comparable oxidation level - the presence of multivalent counterions has little effect on the stiffness. Evidently, molecular doping is a powerful tool that can be used for the design of mechanically robust conducting materials, which may find use within the field of flexible and stretchable electronics.

Faradaic Pixels for Precise Hydrogen Peroxide Delivery to Control M-Type Voltage-Gated Potassium Channels

H2 O2 plays a significant role in a range of physiological processes where it performs vital tasks in redox signaling. The sensitivity of many biological pathways to H2 O2 opens up a unique direction in the development of bioelectronics devices to control levels of reactive-oxygen species (ROS). Here a microfabricated ROS modulation device that relies on controlled faradaic reactions is presented. A concentric pixel arrangement of a peroxide-evolving cathode surrounded by an anode ring which decomposes the peroxide, resulting in localized peroxide delivery is reported. The conducting polymer (poly(3,4-ethylenedioxythiophene) (PEDOT), is exploited as the cathode. PEDOT selectively catalyzes the oxygen reduction reaction resulting in the production of hydrogen peroxide (H2 O2 ). Using electrochemical and optical assays, combined with modeling, the performance of the devices is benchmarked. The concentric pixels generate tunable gradients of peroxide and oxygen concentrations. The faradaic devices are prototyped by modulating human H2 O2 -sensitive Kv7.2/7.3 (M-type) channels expressed in a single-cell model (Xenopus laevis oocytes). The Kv7 ion channel family is responsible for regulating neuronal excitability in the heart, brain, and smooth muscles, making it an ideal platform for faradaic ROS stimulation. The results demonstrate the potential of PEDOT to act as an H2 O2 delivery system, paving the way to ROS-based organic bioelectronics.

Ion Diffusion through Nanocellulose Membranes: Molecular Dynamics Study

One of the most promising applications of nanocellulose is for membranes for energy storage devices including supercapacitors, batteries, and fuel cells. Several recent studies reported the fabrication of cellulose-based membranes where ionic conductivity was confined to channels. So far, theoretical understanding of the effect of the nanoconfinement and surface charged groups on the diffusion coefficient of ions in cellulose nanochannels is missing. In the present study, we perform atomistic molecular dynamics simulations to provide this theoretical understanding and unravel mechanisms affecting the ionic diffusion in nanochannels. We demonstrate that the diffusion coefficient of ions in cellulose nanochannels is reduced in comparison to its bulk value. The change of the diffusion coefficient depends on the density of charged surface groups in nanochannels and the channel height, and it is primarily caused by the Coulomb interaction between the ions and the surface. We believe that our results reveal an important structure/property relationship in cellulose nanochannels, and they show that accounting for the dependence of the diffusion coefficient on the charge of the surface groups and channel height can be important for the Nernst-Plank-Poisson modeling of the ion conductivity in nanomembranes as well as for accurate fitting the experimental data to extract the material parameters.

How Long are Polymer Chains in Poly(3,4-ethylenedioxythiophene):Tosylate Films? An Insight from Molecular Dynamics Simulations

Poly(3,4-ethylenedioxythiophene) (PEDOT) is one of the most important conductive polymers utilized in a variety of applications in organic electronics and bioelectronics and energy storage. PEDOT chains are believed to be rather short, but detailed knowledge of their length is missing because of the challenges in its experimental determination due to insolubility of PEDOT films. Here, we report a molecular dynamics (MD) study of in situ oxidative chemical polymerization and simultaneous crystallization of molecularly doped PEDOT focusing on the determination of its chain lengths at different polymerization temperatures. We find the average chain length to be 6, 7, and 11 monomers for 298, 323 and 373 K, respectively. At the same time, the length distribution is rather broad, for example, between 2 and 16 monomer units for T = 323 K. We demonstrate that the limiting factor determining the chain length is the diffusivity of the reactants (PEDOT monomers and oligomers). We also study the polymer film formation during solvent evaporation, and we find that although crystallization starts and proceeds already during the polymerization and doping phases, it mostly occurs during the evaporation phase. Finally, we believe that our results providing the oligomer chain length and polymerization and crystallization mechanisms obtained by means of MD "computational microscopy" provide an important insight into the morphology of PEDOT that cannot be obtained by other means.

PEDOT:PSS nano-particles in aqueous media: A comparative experimental and molecular dynamics study of particle size, morphology and z-potential

PEDOT

PSS is the most widely used conducting polymer in organic and printed electronics.

PEDOT

PSS films have been extensively studied to understand the morphology, ionic and electronic conductivity of the polymer. However, the polymer dispersion, which is used to cast or spin coat the films, is not well characterized and not well understood theoretically. Here, we study in detail the particle morphology, size, charge density and zeta potential (z-potential) by coarse-grained MD simulations and dynamic light scattering (DLS) measurements, for different pH levels and ionic strengths. The PEDOT:PSS particles were found to be 12 nm-19 nm in diameter and had a z-potential of -30 mV to -50 mV when pH was changed from 1.7 to 9, at an added NaCl concentration of 1 mM, as measured by DLS. These values changed significantly with changing pH and ionic strength of the solution. The charge density of PEDOT:PSS particles was also found to be dependent on pH and ionic strength. Besides, the distribution of different ions (PSS^-, PEDOT⁺, Na⁺, Cl^-) present in the solution is simulated to understand the particle morphology and molecular origin of z-potential in PEDOT:PSS dispersion. The trend in change of particle size, charge density and z- potential with changing pH and ionic strength are in good agreement between the simulations and experiments. Our results show that the molecular model developed in this work represents very well the PEDOT:PSS nano-particles in aqueous dispersion. With this study, we hope to provide new insight and an in-depth understanding of the morphology and z-potential evolution in PEDOT:PSS dispersion.

UV-to-IR Absorption of Molecularly p-Doped Polythiophenes with Alkyl and Oligoether Side Chains: Experiment and Interpretation Based on Density Functional Theory

The UV-to-IR transitions in p-doped poly(3-hexylthiophene) (P3HT) with alkyl side chains and polar polythiophene with tetraethylene glycol side chains are studied experimentally by means of the absorption spectroscopy and computationally using density functional theory (DFT) and tight-binding DFT. The evolution of electronic structure is calculated as the doping level is varied, while the roles of dopant ions, chain twisting, and π-π stacking are also considered, each of these having the effect of broadening the absorption peaks while not significantly changing their positions. The calculated spectra are found to be in good agreement with experimental spectra obtained for the polymers doped with a molybdenum dithiolene complex. As in other DFT studies of doped conjugated polymers, the electronic structure and assignment of optical transitions that emerge are qualitatively different from those obtained through earlier "traditional" approaches. In particular, the two prominent bands seen for the p-doped materials are present for both polarons and bipolarons/polaron pairs. The lowest energy of these transitions is due to excitation from the valence band to a spin-resolved orbitals located in the gap between the bands. The higher-energy band is a superposition of excitation from the valence band to a spin-resolved orbitals in the gap and an excitation between bands.

Density Functional Theory Mechanistic Study on H2O2 Production Using an Organic Semiconductor Epindolidione

Organic semiconductors have recently emerged as promising catalytic materials for oxygen reduction to hydrogen peroxide, H₂O₂, a chemical of great importance in industry as well as biology. While examples of organic semiconductor-mediated photocatalytic and electrocatalytic processes for H₂O₂ production become more numerous and improve in performance, fundamental understanding of the reaction mechanisms at play have been explored far less. The aim of the present work is to computationally test hypotheses of how selective oxygen reduction to H₂O₂ generally occurs on carbonyl dyes and pigments. As an example material, we consider epindolidione (EPI), an industrial pigment with demonstrated semiconductor properties, which photocatalytic activity in oxygen reduction reaction (ORR) and thereby producing hydrogen peroxide (H₂O₂) in low pH environment has been recently experimentally demonstrated. In this work, the ability of the reduced form of EPI, viz. EPI-2H (which was formed after a photoinduced 2e^–/2H⁺ process), to reduce molecular triplet oxygen to peroxide and the possible mechanism of this reaction are computationally investigated using density functional theory. In the main reaction pathway, the reduction of O₂ to H₂O₂ reaction occurs via abstraction of one of the hydrogen atoms of EPI-2H by triplet dioxygen to produce an intermediate complex consisting of the radicals of hydrogen peroxide (HOO^•) and EPI-H^• at the initial stage. HOO^• thus released can abstract another hydrogen atom from EPI-H^• to produce H₂O₂ and regenerates EPI; otherwise, it can enter another pathway to abstract hydrogen from a neighboring EPI-2H to form EPI-H^• and H₂O₂. EPI, after reduction, thus plays in ORR the role of hydrogen atom transfer (HAT) agent via its OH group, similar to anthraquinone in the industrial process, while HAT from its amino hydrogen is found unfavorable.

Side Chain Redistribution as a Strategy to Boost Organic Electrochemical Transistor Performance and Stability

A series of glycolated polythiophenes for use in organic electrochemical transistors (OECTs) is designed and synthesized, differing in the distribution of their ethylene glycol chains that are tethered to the conjugated backbone. While side chain redistribution does not have a significant impact on the optoelectronic properties of the polymers, this molecular engineering strategy strongly impacts the water uptake achieved in the polymers. By careful optimization of the water uptake in the polymer films, OECTs with unprecedented steady-state performances in terms of [μC^* ] and current retentions up to 98% over 700 electrochemical switching cycles are developed.

Reversible Electronic Solid-Gel Switching of a Conjugated Polymer

Conjugated polymers exhibit electrically driven volume changes when included in electrochemical devices via the exchange of ions and solvent. So far, this volumetric change is limited to 40% and 100% for reversible and irreversible systems, respectively, thus restricting potential applications of this technology. A conjugated polymer that reversibly expands by about 300% upon addressing, relative to its previous contracted state, while the first irreversible actuation can achieve values ranging from 1000-10 000%, depending on the voltage applied is reported. From experimental and theoretical studies, it is found that this large and reversible volumetric switching is due to reorganization of the polymer during swelling as it transforms between a solid-state phase and a gel, while maintaining percolation for conductivity. The polymer is utilized as an electroactive cladding to reduce the void sizes of a porous carbon filter electrode by 85%.

Structural Control of Charge Storage Capacity to Achieve 100% Doping in Vapor Phase-Polymerized PEDOT/Tosylate

Vapor phase polymerization (VPP) is used to fabricate a series of tosylate-doped poly(3,4-ethylenedioxythiophene) (PEDOT) electrodes on carbon paper. The series of VPP PEDOT/tosylate coatings has varying levels of crystallinity and electrical conductivity because of the use (or not) of nonionic triblock copolymers in the oxidant solution during synthesis. As a result, the impact of the structure on charge storage capacity is investigated using tetra-n-butylammonium hexafluorophosphate (0.1 M in acetonitrile). The ability to insert anions, and hence store charge, of the VPP PEDOT/tosylate is inversely related to its electrical conductivity. In the case of no nonionic triblock copolymer employed, the VPP PEDOT/tosylate achieves electrochemical doping levels of 1.0 charge per monomer or greater (≥100% doping level). Such high doping levels are demonstrated to be plausible by molecular dynamics simulations and density functional theory calculations. Experiments show that this high doping level is attainable when the PEDOT structure is weakly crystalline with (relatively) large crystallite domains.

Ion Electron-Coupled Functionality in Materials and Devices Based on Conjugated Polymers

The coupling between charge accumulation in a conjugated polymer and the ionic charge compensation, provided from an electrolyte, defines the mode of operation in a vast array of different organic electrochemical devices. The most explored mixed organic ion-electron conductor, serving as the active electrode in these devices, is poly(3,4-ethyelenedioxythiophene) doped with polystyrelensulfonate (PEDOT:PSS). In this progress report, scientists of the Laboratory of Organic Electronics at Linköping University review some of the achievements derived over the last two decades in the field of organic electrochemical devices, in particular including PEDOT:PSS as the active material. The recently established understanding of the volumetric capacitance and the mixed ion-electron charge transport properties of PEDOT are described along with examples of various devices and phenomena utilizing this ion-electron coupling, such as the organic electrochemical transistor, ionic-electronic thermodiffusion, electrochromic devices, surface switches, and more. One of the pioneers in this exciting research field is Prof. Olle Inganäs and the authors of this progress report wish to celebrate and acknowledge all the fantastic achievements and inspiration accomplished by Prof. Inganäs all since 1981.

Computational microscopy study of the granular structure and pH dependence of PEDOT:PSS

Computational microscopy based on Martini coarse grained molecular dynamics (MD) simulations of a doped conducting polymer poly(3,4-ethylenedioxythiophene)polystyrene sulfonate (best known as PEDOT:PSS) was performed focussing on the formation of the granular structure and PEDOT crystallites, and the effect of pH on the material morphology. The PEDOT:PSS morphology is shown to be sensitive to the initial distribution of PEDOT and PSS in the solution, and the results of the modelling suggest that the experimentally observed granular structure of PEDOT:PSS can be only obtained if the PEDOT/PSS solution is in the dispersive state in the initial crystallization stages. Variation of the pH is demonstrated to strongly affect the morphology of PEDOT:PSS films, altering their structure between granular-type and homogeneous. It also affects the size of crystallites and the relative arrangement of PEDOT and PSS chains. It is shown that the crystallites in PEDOT:PSS are smaller than those in PEDOT with molecular counterions such as PEDOT:tosylate, which is consistent with the available experimental data. The predicted changes of the PEDOT:PSS morphology with variation of the pH can be tested experimentally, and the calculated atomistic picture of PEDOT:PSS films (not accessible by conventional experimental techniques) is instrumental for understanding the material structure and building realistic models of PEDOT:PSS morphology.

The intrinsic volumetric capacitance of conducting polymers: pseudo-capacitors or double-layer supercapacitors?

The capacitance of conducting polymers represents one of the most important material parameters that in many cases determines the device and material performances. Despite a vast number of experimental studies, the theoretical understanding of the origin of the capacitance in conducting polymers remains unsatisfactory and appears even controversial. Here, we present a theoretical method, based on first principle capacitance calculations using density functional theory (DFT), and apply it to calculate the volumetric capacitance of two archetypical conducting polymers: poly(3,4-ethylene dioxythiophene) (PEDOT) and polypyrrole (PPy). Our aim is to achieve a quantitate description of the volumetric capacitance and to provide a qualitative understanding of its nature at the atomistic level. We find that the volumetric capacitance of PEDOT and PPy is ≈100 F cm⁻³ and ≈300 F cm⁻³, respectively, which is within the range of the corresponding reported experimental results. We demonstrate that the capacitance of conducting polymers originates from charges stored in atomistic Stern layers formed by counterions and doped polymeric chains. The Stern layers have a purely electrostatic origin, since the counterions do not form any bonds with the atoms of the polymeric chains, and no charge transfer between the counterions and conducting polymer takes place. This classifies the conducting polymers as double-layer supercapacitors rather than pseudo-capacitors. Further, we analyze contributions to the total capacitance originating from the classical capacitance C_C and the quantum capacitance C_Q, respectively, and find that the latter provides a dominant contribution. The method of calculations of the capacitance developed in the present paper is rather general and opens up the way for engineering and optimizing the capacitive response of the conducting polymers.

Using the density functional theory, the intrinsic volumetric capacitance of conducting polymers is calculated. It is shown that conducting polymers operate as double-layer supercapacitors rather than pseudo-capacitors.

Morphology and ion diffusion in PEDOT:Tos. A coarse grained molecular dynamics simulation

A Martini coarse-grained Molecular Dynamics (MD) model for the doped conducting polymer poly(3,4-ethylenedioxythiophene) (PEDOT) is developed. It is shown that the diffusion coefficients decrease exponentially as the hydration level is reduced.

Charge transport and structure in semimetallic polymers

Owing to changes in their chemistry and structure, polymers can be fabricated to demonstrate vastly different electrical conductivities over many orders of magnitude. At the high end of conductivity is the class of conducting polymers, which are ideal candidates for many applications in low‐cost electronics. Here, we report the influence of the nature of the doping anion at high doping levels within the semi‐metallic conducting polymer poly(3,4‐ethylenedioxythiophene) (PEDOT) on its electronic transport properties. Hall effect measurements on a variety of PEDOT samples show that the choice of doping anion can lead to an order of magnitude enhancement in the charge carrier mobility > 3 cm²/Vs at conductivities approaching 3000 S/cm under ambient conditions. Grazing Incidence Wide Angle X‐ray Scattering, Density Functional Theory calculations, and Molecular Dynamics simulations indicate that the chosen doping anion modifies the way PEDOT chains stack together. This link between structure and specific anion doping at high doping levels has ramifications for the fabrication of conducting polymer‐based devices. © 2017 The Authors. Journal of Polymer Science Part B: Polymer Physics Published by Wiley Periodicals, Inc. J. Polym. Sci., Part B: Polym. Phys. 2018, 56, 97–104

pH Dependence of γ-Aminobutyric Acid Iontronic Transport

The organic electronic ion pump (OEIP) has been developed as an "iontronic" tool for delivery of biological signaling compounds. OEIPs rely on electrophoretically "pumping" charged compounds, either at neutral or shifted pH, through an ion-selective channel. Significant shifts in pH lead to an abundance of H⁺ or OH^-, which are delivered along with the intended substance. While this method has been used to transport various neurotransmitters, the role of pH has not been explored. Here we present an investigation of the role of pH on OEIP transport efficiency using the neurotransmitter γ-aminobutyric acid (GABA) as the model cationic delivery substance. GABA transport is evaluated at various pHs using electrical and chemical characterization and compared to molecular dynamics simulations, all of which agree that pH 3 is ideal for GABA transport. These results demonstrate a useful method for optimizing transport of other substances and thus broadening OEIP applications.

Semi-metallic polymers

Polymers are lightweight, flexible, solution-processable materials that are promising for low-cost printed electronics as well as for mass-produced and large-area applications. Previous studies demonstrated that they can possess insulating, semiconducting or metallic properties; here we report that polymers can also be semi-metallic. Semi-metals, exemplified by bismuth, graphite and telluride alloys, have no energy bandgap and a very low density of states at the Fermi level. Furthermore, they typically have a higher Seebeck coefficient and lower thermal conductivities compared with metals, thus being suitable for thermoelectric applications. We measure the thermoelectric properties of various poly(3,4-ethylenedioxythiophene) samples, and observe a marked increase in the Seebeck coefficient when the electrical conductivity is enhanced through molecular organization. This initiates the transition from a Fermi glass to a semi-metal. The high Seebeck value, the metallic conductivity at room temperature and the absence of unpaired electron spins makes polymer semi-metals attractive for thermoelectrics and spintronics.