The ‘phonon-glass electron-crystal’ concept has triggered most of the progress that has been achieved in inorganic thermoelectrics in the past two decades. Organic thermoelectric materials, unlike their inorganic counterparts, exhibit molecular diversity, flexible mechanical properties and easy fabrication, and are mostly ‘phonon glasses’. However, the thermoelectric performances of these organic materials are largely limited by low molecular order and they are therefore far from being ‘electron crystals’. Here, we report a molecularly n-doped fullerene derivative with meticulous design of the side chain that approaches an organic ‘PGEC’ thermoelectric material. This thermoelectric material exhibits an excellent electrical conductivity of >10 S cm−1 and an ultralow thermal conductivity of <0.1 Wm−1K−1, leading to the best figure of merit ZT = 0.34 (at 120 °C) among all reported single-host n-type organic thermoelectric materials. The key factor to achieving the record performance is to use ‘arm-shaped’ double-triethylene-glycol-type side chains, which not only offer excellent doping efficiency (~60%) but also induce a disorder-to-order transition upon thermal annealing. This study illustrates the vast potential of organic semiconductors as thermoelectric materials.

polymers have also emerged as potential candidates for organic thermoelectrics, [7,8] potentially delivering flexible, large-area, and low-cost energy generation or heatingcooling devices for appealing applications, e.g., wearable energy harvesting, that are currently not possible for traditional brittle and usually either toxic or rare inorganic crystalline thermoelectric materials. Thermoelectric materials are evaluated by the dimensionless figure of merit ZT = S 2 σT/κ, where S, σ, T and κ represent the Seebeck coefficient, electrical conductivity, absolute temperature and thermal conductivity, respectively. Most conjugated polymers are characterized with low κ values, intrinsically contributing to a high ZT. This point has been verified by recent extensive thermoelectric studies based on p-type conjugated polymers such as poly(3,4-ethylenedioxythiophene) (PEDOT) with ZT > 0.25. [9,10] The performance of p-type and n-type thermoelectric materials should pair with each other ahead of any practical applications. However, n-type conjugated polymer-based thermoelectric devices are still far inferior to their p-type counterparts in terms of power factor (S 2 σ). [11,12] Therefore, the development of efficient The ORCID identification number(s) for the author(s) of this article can be found under

No abstract

Incorporating ethylene glycols (EGs) into organic semiconductors has become the prominent strategy to increase their dielectric constant. However, EG’s contribution to the dielectric constant is due to nuclear relaxations, and therefore, its relevance for various organic electronic applications depends on the time scale of these relaxations, which remains unknown. In this work, by means of a new computational protocol based on polarizable molecular dynamics simulations, the time- and frequency-dependent dielectric constant of a representative fullerene derivative with EG side chains is predicted, the origin of its unusually high dielectric constant is explained, and design suggestions are made to further increase it. Finally, a dielectric relaxation time of ∼1 ns is extracted which suggests that EGs may be too slow to reduce the Coulombic screening in organic photovoltaics but are definitely fast enough for organic thermoelectrics with much lower charge carrier velocities.

Organic semiconductors consisting of molecules bearing polar side chains have been proposed as potential candidates to overcome the limitations of organic photovoltaics owing to their enhanced dielectric constant. However, introducing such polar molecules in photovoltaic devices has not yet resulted in higher efficiencies. A microscopic understanding of the impact of polar side chains on electronic and structural properties of organic semiconductors is paramount to rationalize their effect. Here, the impact of such side chains on bulk heterojunction overall morphology, molecular configurations at donor-acceptor (DA) interfaces, and charge carrier energy levels is investigated. The multiscale modeling approach used allows to resolve DA interfaces with atomistic resolution while taking into account the large-scale self-organization process which takes place during the processing of an organic thin film. The polar fullerene-based blends are compared to the well-studied reference system, poly(3-hexyl-thiophene) (P3HT):phenyl-C 61butyric acid methyl ester (PCBM). Introduction of polar side chains on a similar molecular scaffold does not affect molecular orientations at the DA interfaces; such orientations are, however, found to be affected by processing conditions and polymer molecular weight. Polar side chains, instead, are found to impact considerably the charge carrier energy levels of the organic blend, causing electrostatic-induced broadening of these levels.

The quality of molecular dynamics simulations strongly depends on the accuracy of the underlying force fields (FFs) that determine all intra- and intermolecular interactions of the system. Commonly, transferable FF parameters are determined based on a representative set of small molecules. However, such an approach sacrifices accuracy in favor of generality. In this work, an open-source and automated toolkit named Q-Force is presented, which augments these transferable FFs with molecule-specific bonded parameters and atomic charges that are derived from quantum mechanical (QM) calculations. The molecular fragmentation procedure allows treatment of large molecules (>200 atoms) with a low computational cost. The generated Q-Force FFs can be used at the same computational cost as transferable FFs, but with improved accuracy: We demonstrate this for the vibrational properties on a set of small molecules and for the potential energy surface on a complex molecule (186 atoms) with photovoltaic applications. Overall, the accuracy, user-friendliness, and minimal computational overhead of the Q-Force protocol make it widely applicable for atomistic molecular dynamics simulations.

The low efficiency of organic photovoltaic (OPV) devices has often been attributed to the strong Coulombic interactions between the electron and hole, impeding the charge separation process. Recently, it has been argued that by increasing the dielectric constant of materials used in OPVs, this strong interaction could be screened. In this work, we report the application of periodic density functional theory together with the coupled perturbed Kohn–Sham method to calculate the electronic contribution to the dielectric constant for fullerene C60 derivatives, a ubiquitous class of molecules in the field of OPVs. The results show good agreement with experimental data when available and also reveal an important undesirable outcome when manipulating the side chain to maximize the static dielectric constant: in all cases, the electronic contribution to the dielectric constant decreases as the side chain increases in size. This information should encourage both theoreticians and experimentalists to further investigate the relevance of contributions to the dielectric constant from slower processes like vibrations and dipolar reorientations for facilitating the charge separation, because electronically, enlarging the side chain of conventional fullerene derivatives only lowers the dielectric constant, and consequently, their electronic dielectric constant is upper bound by the one of C60.

1,2,4-Benzenetriol (BTO), sourced from the carbohydrate-derived platform chemical 5-hydroxylmethylfurfural (HMF), is an interesting starting point for the synthesis of various biobased aromatic products. However, BTO readily undergoes dimerization and other reactions under mild conditions, making analysis and isolation challenging. To both control and utilize the reactivity of BTO to produce biobased building blocks, its reactivity needs to be better understood. Here it was found that specific BTO aromatic C–H bonds are reactive toward deuterium exchange with D2O, which appears pronounced under acidic conditions at room temperature and can lead to the selective formation of BTO with an aromatic ring that contains one or two deuterium atoms, the first at the five and the second at the three position. By exposure to air, it was shown that BTO forms a 5,5′-linked BTO dimer [1,1′-biphenyl]-2,2′,4,4′,5,5′-hexaol (1) and subsequently a hydroxyquinone containing dimeric structure 2′,4,4′,5′-tetrahydroxy-[1,1′-biphenyl]-2,5-dione (2). Additionally, condensed dimer dibenzo[b,d]furan-2,3,7,8-tetraol (3) can be relatively easily accessed. The controlled formation of these symmetric and asymmetric multifunctional dimers illustrates diverse possibilities for BTO to be converted to valuable biobased aromatic compounds. Deuterium exchange was attributed to electrophilic aromatic substitution because this reactivity was found to be independent of oxygen and acid mediated. On the contrary, the dimerization was dependent on the presence of oxygen and thus likely involves radical intermediates. Thus this report overall displays different accessible reaction pathways for BTO that can be exploited for the production of BTO-derived compounds.