1
|
Stroet M, Sanderson S, Sanzogni AV, Nada S, Lee T, Caron B, Mark AE, Burn PL. PyThinFilm: Automated Molecular Dynamics Simulation Protocols for the Generation of Thin Film Morphologies. J Chem Inf Model 2023; 63:2-8. [PMID: 36539938 DOI: 10.1021/acs.jcim.2c01334] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
The performance of organic optoelectronic devices, such as organic light-emitting diodes (OLEDs) and organic solar cells (OSCs), is intrinsically related to the molecular-scale morphology of the thin films from which they are composed. However, the experimental characterization of morphology at the molecular level is challenging due to the often amorphous or at best semicrystalline nature of these films. Classical molecular modeling techniques, such as molecular dynamics (MD) simulation, are increasingly used to understand the relationship between morphology and the properties of thin-film devices. PyThinFilm (github.com/ATB-UQ/PyThinFilm) is an open-source Python package which allows fully automated MD simulations of thin film growth to be performed using vacuum and/or solution deposition processes. PyThinFilm utilizes the GROMACS simulation package in combination with interaction parameters from the Automated Topology Builder (atb.uq.edu.au). Here, PyThinFilm is described along with an overview of applications in which PyThinFilm has been used to study the thin films of organic semiconductor materials typically used in OLEDs and OSCs.
Collapse
Affiliation(s)
- Martin Stroet
- Centre for Organic Photonics & Electronics, The University of Queensland, St Lucia Campus, Brisbane, Queensland4072, Australia.,Molecular Dynamics Group, School of Chemistry and Molecular Biosciences, The University of Queensland, St Lucia Campus, Brisbane, Queensland4072, Australia
| | - Stephen Sanderson
- The Australian Institute for Bioengineering and Nanotechnology, The University of Queensland, St Lucia Campus, Brisbane, Queensland4072, Australia.,College of Science and Engineering, James Cook University, Townsville, Queensland4811, Australia
| | - Audrey V Sanzogni
- Centre for Organic Photonics & Electronics, The University of Queensland, St Lucia Campus, Brisbane, Queensland4072, Australia.,Molecular Dynamics Group, School of Chemistry and Molecular Biosciences, The University of Queensland, St Lucia Campus, Brisbane, Queensland4072, Australia
| | - Sharif Nada
- Molecular Dynamics Group, School of Chemistry and Molecular Biosciences, The University of Queensland, St Lucia Campus, Brisbane, Queensland4072, Australia
| | - Thomas Lee
- Centre for Organic Photonics & Electronics, The University of Queensland, St Lucia Campus, Brisbane, Queensland4072, Australia.,Molecular Dynamics Group, School of Chemistry and Molecular Biosciences, The University of Queensland, St Lucia Campus, Brisbane, Queensland4072, Australia
| | - Bertrand Caron
- Molecular Dynamics Group, School of Chemistry and Molecular Biosciences, The University of Queensland, St Lucia Campus, Brisbane, Queensland4072, Australia
| | - Alan E Mark
- Molecular Dynamics Group, School of Chemistry and Molecular Biosciences, The University of Queensland, St Lucia Campus, Brisbane, Queensland4072, Australia
| | - Paul L Burn
- Centre for Organic Photonics & Electronics, The University of Queensland, St Lucia Campus, Brisbane, Queensland4072, Australia
| |
Collapse
|
2
|
Ronsin OJJ, Harting J. Formation of Crystalline Bulk Heterojunctions in Organic Solar Cells: Insights from Phase-Field Simulations. ACS APPLIED MATERIALS & INTERFACES 2022; 14:49785-49800. [PMID: 36282868 DOI: 10.1021/acsami.2c14319] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/16/2023]
Abstract
The performance of organic solar cells strongly depends on the bulk-heterojunction (BHJ) morphology of the photoactive layer. This BHJ forms during the drying of the wet-deposited solution, because of physical processes such as crystallization and/or liquid-liquid phase separation (LLPS). However, the process-structure relationship remains insufficiently understood. In this work, a recently developed, coupled phase-field-fluid mechanics framework is used to simulate the BHJ formation upon drying. For the first time, this allows to investigate the interplay between all the relevant physical processes (evaporation, crystal nucleation and growth, liquid demixing, composition-dependent kinetic properties), within a single coherent theoretical framework. Simulations for the model system P3HT-PCBM are presented. The comparison with previously reported in situ characterization of the drying structure is very convincing: The morphology formation pathways, crystallization kinetics, and final morphology are in line with experimental results. The final BHJ morphology is a subtle mixture of pure crystalline donor and acceptor phases, pure and mixed amorphous domains, which depends on the process parameters and material properties. The expected benefit of such an approach is to identify physical design rules for ink formulation and processing conditions to optimize the cell's performance. It could be applied to recent organic material systems in the future.
Collapse
Affiliation(s)
- Olivier J J Ronsin
- Helmholtz Institute Erlangen-Nürnberg for Renewable Energy, Forschungszentrum Jülich, Fürther Straße 248, 90429Nürnberg, Germany
- Department of Chemical and Biological Engineering, Friedrich-Alexander-Universität Erlangen-Nürnberg, Fürther Straße 248, 90429Nürnberg, Germany
| | - Jens Harting
- Helmholtz Institute Erlangen-Nürnberg for Renewable Energy, Forschungszentrum Jülich, Fürther Straße 248, 90429Nürnberg, Germany
- Department of Chemical and Biological Engineering, Friedrich-Alexander-Universität Erlangen-Nürnberg, Fürther Straße 248, 90429Nürnberg, Germany
- Department of Physics, Friedrich-Alexander-Universität Erlangen-Nürnberg, Fürther Straße 248, 90429Nürnberg, Germany
| |
Collapse
|
3
|
Ronsin OJJ, Harting J. Phase‐Field Simulations of the Morphology Formation in Evaporating Crystalline Multicomponent Films. ADVANCED THEORY AND SIMULATIONS 2022. [DOI: 10.1002/adts.202200286] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Affiliation(s)
- Olivier J. J. Ronsin
- Helmholtz Institute Erlangen‐Nürnberg for Renewable Energy Forschungszentrum Jülich Fürther Straße 248 90429 Nürnberg Germany
- Department of Chemical and Biological Engineering Friedrich‐Alexander‐Universität Erlangen‐Nürnberg Fürther Straße 248 90429 Nürnberg Germany
| | - Jens Harting
- Helmholtz Institute Erlangen‐Nürnberg for Renewable Energy Forschungszentrum Jülich Fürther Straße 248 90429 Nürnberg Germany
- Department of Chemical and Biological Engineering Friedrich‐Alexander‐Universität Erlangen‐Nürnberg Fürther Straße 248 90429 Nürnberg Germany
- Department of Physics Friedrich‐Alexander‐Universität Erlangen‐Nürnberg Fürther Straße 248 90429 Nürnberg Germany
| |
Collapse
|
4
|
Nguyen HTL, Huang DM. Systematic bottom-up molecular coarse-graining via force and torque matching using anisotropic particles. J Chem Phys 2022; 156:184118. [DOI: 10.1063/5.0085006] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We derive a systematic and general method for parametrizing coarse-grained molecular models consisting of anisotropic particles from fine-grained (e.g. all-atom) models for condensed-phase molecular dynamics simulations. The method, which we call anisotropic force-matching coarse-graining (AFM-CG), is based on rigorous statistical mechanical principles, enforcing consistency between the coarse-grained and fine-grained phase-space distributions to derive equations for the coarse-grained forces, torques, masses, and moments of inertia in terms of properties of a condensed-phase fine-grained system. We verify the accuracy and efficiency of the method by coarse-graining liquid-state systems of two different anisotropic organic molecules, benzene and perylene, and show that the parametrized coarse-grained models more accurately describe properties of these systems than previous anisotropic coarse-grained models parametrized using other methods that do not account for finite-temperature and many-body effects on the condensed-phase coarse-grained interactions. The AFM-CG method will be useful for developing accurate and efficient dynamical simulation models of condensed-phase systems of molecules consisting of large, rigid, anisotropic fragments, such as liquid crystals, organic semiconductors, and nucleic acids.
Collapse
|
5
|
Ronsin OJJ, Jang D, Egelhaaf HJ, Brabec CJ, Harting J. Phase-Field Simulation of Liquid-Vapor Equilibrium and Evaporation of Fluid Mixtures. ACS APPLIED MATERIALS & INTERFACES 2021; 13:55988-56003. [PMID: 34792348 DOI: 10.1021/acsami.1c12079] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
In solution processing of thin films, the material layer is deposited from a solution composed of several solutes and solvents. The final morphology and hence the properties of the film often depend on the time needed for the evaporation of the solvents. This is typically the case for organic photoactive or electronic layers. Therefore, it is important to be able to predict the evaporation kinetics of such mixtures. We propose here a new phase-field model for the simulation of evaporating fluid mixtures and simulate their evaporation kinetics. Similar to the Hertz-Knudsen theory, the local liquid-vapor (LV) equilibrium is assumed to be reached at the film surface and evaporation is driven by diffusion away from this gas layer. In the situation where the evaporation is purely driven by the LV equilibrium, the simulations match the behavior expected theoretically from the free energy: for evaporation of pure solvents, the evaporation rate is constant and proportional to the vapor pressure. For mixtures, the evaporation rate is in general strongly time-dependent because of the changing composition of the film. Nevertheless, for highly nonideal mixtures, such as poorly compatible fluids or polymer solutions, the evaporation rate becomes almost constant in the limit of low Biot numbers. The results of the simulation have been successfully compared to experiments on a polystyrene-toluene mixture. The model allows to take into account deformations of the liquid-vapor interface and, therefore, to simulate film roughness or dewetting.
Collapse
Affiliation(s)
- Olivier J J Ronsin
- Forschungszentrum Jülich GmbH, Helmholtz Institute Erlangen-Nürnberg for Renewable Energy (IEK-11), Fürther Straße 248, 90429 Nürnberg, Germany
- Department of Chemical and Biological Engineering, Friedrich-Alexander-Universität Erlangen-Nürnberg, Fürther Straße 248, 90429 Nürnberg, Germany
| | - DongJu Jang
- ZAE Bayern─Solar Factory of the Future, Energy Campus Nürnberg, Fürther Straße 250, 90429 Nürnberg, Germany
| | - Hans-Joachim Egelhaaf
- ZAE Bayern─Solar Factory of the Future, Energy Campus Nürnberg, Fürther Straße 250, 90429 Nürnberg, Germany
- Institute of Materials for Electronics and Energy Technology (i-MEET), Friedrich-Alexander-Universität Erlangen-Nürnberg, 91058 Erlangen, Germany
- Forschungszentrum Jülich GmbH, Helmholtz Institute Erlangen-Nürnberg for Renewable Energy (IEK-11), Immerwahrstrasse 2, 91058 Erlangen, Germany
| | - Christoph J Brabec
- Institute of Materials for Electronics and Energy Technology (i-MEET), Friedrich-Alexander-Universität Erlangen-Nürnberg, 91058 Erlangen, Germany
- Forschungszentrum Jülich GmbH, Helmholtz Institute Erlangen-Nürnberg for Renewable Energy (IEK-11), Immerwahrstrasse 2, 91058 Erlangen, Germany
| | - Jens Harting
- Forschungszentrum Jülich GmbH, Helmholtz Institute Erlangen-Nürnberg for Renewable Energy (IEK-11), Fürther Straße 248, 90429 Nürnberg, Germany
- Department of Chemical and Biological Engineering, Friedrich-Alexander-Universität Erlangen-Nürnberg, Fürther Straße 248, 90429 Nürnberg, Germany
- Department of Physics, Friedrich-Alexander-Universität Erlangen-Nürnberg, Fürther Straße 248, 90429 Nürnberg, Germany
| |
Collapse
|
6
|
Cohen AE, Jackson NE, de Pablo JJ. Anisotropic Coarse-Grained Model for Conjugated Polymers: Investigations into Solution Morphologies. Macromolecules 2021. [DOI: 10.1021/acs.macromol.1c00302] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Alexander E. Cohen
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, Illinois 60637, United States
| | - Nicholas E. Jackson
- Department of Chemistry, University of Illinois, Urbana, Illinois 61801, United States
| | - Juan J. de Pablo
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, Illinois 60637, United States
- Center for Molecular Engineering, Argonne National Laboratory, Lemont, Illinois 60439, United States
| |
Collapse
|
7
|
Marmolejo-Valencia AF, Mata-Pinzón Z, Amador-Bedolla C. Charge-transfer electronic states in organic solar cells: a TDDFT study. Phys Chem Chem Phys 2021; 23:16806-16815. [PMID: 34323261 DOI: 10.1039/d1cp00723h] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The prediction of new organic photovoltaic materials in organic solar cells (OSCs) must include a precise description of charge-transfer states because they are involved in electron-transfer processes such as charge separation and charge recombination which govern the device efficiency. Also, as the experimental performance of an optoelectronic device is measured for nonequilibrium nanostructures, computational approaches need models that can incorporate morphology effects. Usually, this aspect is treated by molecular dynamics simulation (MDS) methodologies; however, methodologies and formalisms to calculate the electron-transfer processes are still controversial and sometimes do not connect their information with the phase morphologies. In this work we propose a simple and fast characterization of electron-transfer processes to find the rate constants by analysing the distribution of vertical excitation energies of both local excitation (LE) and charge-transfer (CT) states using TD-DFT calculations in the donor-acceptor pair structures which were extracted from MDS. This proposal assumes that conformational changes are prevented and equilibria are not achieved while the electron-transfer events take effect, and thus the only pathway that connects the LE and CT states is their surface crossing point where an ideal distribution might exist. Different density functionals and dialectric models were tested. The results indicate a close relationship between the proposal and experimental data for electron-transfer events, suggesting the application of this method in the rational design of new photovoltaic materials.
Collapse
Affiliation(s)
- Andres F Marmolejo-Valencia
- Facultad de Química, Universidad Nacional Autónoma de México, Av. Universidad 3000, Coyoacán, CDMX 04510, Mexico.
| | | | | |
Collapse
|
8
|
Ronsin OJJ, Jang D, Egelhaaf HJ, Brabec CJ, Harting J. A phase-field model for the evaporation of thin film mixtures. Phys Chem Chem Phys 2020; 22:6638-6652. [PMID: 32159553 DOI: 10.1039/d0cp00214c] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The performance of solution-processed solar cells strongly depends on the geometrical structure and roughness of the photovoltaic layers formed during film drying. During the drying process, the interplay of crystallization and liquid-liquid demixing leads to structure formation on the nano- and microscale and to the final rough film. In order to better understand how the film structure can be improved by process engineering, we aim at theoretically investigating these systems by means of phase-field simulations. We introduce an evaporation model based on the Cahn-Hilliard equation for the evolution of the fluid concentrations coupled to the Allen-Cahn equation for the liquid-vapour phase transformation. We demonstrate its ability to match the experimentally measured drying kinetics and study the impact of the parameters of our model. Furthermore, the evaporation of solvent blends and solvent-vapour annealing are investigated. The dry film roughness emerges naturally from our set of equations, as illustrated through preliminary simulations of spinodal decomposition and film drying on structured substrates.
Collapse
Affiliation(s)
- Olivier J J Ronsin
- Helmholtz Institute Erlangen-Nürnberg for Renewable Energy, Forschungszentrum Jülich, Fürther Straße 248, 90429 Nürnberg, Germany.
| | - DongJu Jang
- ZAE Bayern-Solar Factory of the Future, Energy Campus Nürnberg, Fürther Straße 250, 90429 Nürnberg, Germany
| | - Hans-Joachim Egelhaaf
- ZAE Bayern-Solar Factory of the Future, Energy Campus Nürnberg, Fürther Straße 250, 90429 Nürnberg, Germany
| | - Christoph J Brabec
- Institute of Materials for Electronics and Energy Technology (i-MEET), Friedrich-Alexander-Universität Erlangen-Nürnberg, 91058 Erlangen, Germany and Helmholtz Institute Erlangen-Nürnberg for Renewable Energy, Forschungszentrum Jülich, Immerwahrstrasse 2, 91058 Erlangen, Germany
| | - Jens Harting
- Helmholtz Institute Erlangen-Nürnberg for Renewable Energy, Forschungszentrum Jülich, Fürther Straße 248, 90429 Nürnberg, Germany. and Department of Applied Physics, Eindhoven University of Technology, PO box 513, 5600MB Eindhoven, The Netherlands
| |
Collapse
|
9
|
Boehm BJ, Nguyen HTL, Huang DM. The interplay of interfaces, supramolecular assembly, and electronics in organic semiconductors. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2019; 31:423001. [PMID: 31212263 DOI: 10.1088/1361-648x/ab2ac2] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
Organic semiconductors, which include a diverse range of carbon-based small molecules and polymers with interesting optoelectronic properties, offer many advantages over conventional inorganic semiconductors such as silicon and are growing in importance in electronic applications. Although these materials are now the basis of a lucrative industry in electronic displays, many promising applications such as photovoltaics remain largely untapped. One major impediment to more rapid development and widespread adoption of organic semiconductor technologies is that device performance is not easily predicted from the chemical structure of the constituent molecules. Fundamentally, this is because organic semiconductor molecules, unlike inorganic materials, interact by weak non-covalent forces, resulting in significant structural disorder that can strongly impact electronic properties. Nevertheless, directional forces between generally anisotropic organic-semiconductor molecules, combined with translational symmetry breaking at interfaces, can be exploited to control supramolecular order and consequent electronic properties in these materials. This review surveys recent advances in understanding of supramolecular assembly at organic-semiconductor interfaces and its impact on device properties in a number of applications, including transistors, light-emitting diodes, and photovoltaics. Recent progress and challenges in computer simulations of supramolecular assembly and orientational anisotropy at these interfaces is also addressed.
Collapse
Affiliation(s)
- Belinda J Boehm
- Department of Chemistry, School of Physical Sciences, The University of Adelaide, SA 5005, Australia
| | | | | |
Collapse
|
10
|
Meng X, Zhang L, Xie Y, Hu X, Xing Z, Huang Z, Liu C, Tan L, Zhou W, Sun Y, Ma W, Chen Y. A General Approach for Lab-to-Manufacturing Translation on Flexible Organic Solar Cells. ADVANCED MATERIALS (DEERFIELD BEACH, FLA.) 2019; 31:e1903649. [PMID: 31423693 DOI: 10.1002/adma.201903649] [Citation(s) in RCA: 43] [Impact Index Per Article: 8.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/09/2019] [Revised: 07/31/2019] [Indexed: 06/10/2023]
Abstract
The blossoming of organic solar cells (OSCs) has triggered enormous commercial applications, due to their high-efficiency, light weight, and flexibility. However, the lab-to-manufacturing translation of the praisable performance from lab-scale devices to industrial-scale modules is still the Achilles' heel of OSCs. In fact, it is urgent to explore the mechanism of morphological evolution in the bulk heterojunction (BHJ) with different coating/printing methods. Here, a general approach to upscale flexible organic photovoltaics to module scale without obvious efficiency loss is demonstrated. The shear impulse during the coating/printing process is first applied to control the morphology evolution of the BHJ layer for both fullerene and nonfullerene acceptor systems. A quantitative transformation factor of shear impulse between slot-die printing and spin-coating is detected. Compelling results of morphological evolution, molecular stacking, and coarse-grained molecular simulation verify the validity of the impulse translation. Accordingly, the efficiency of flexible devices via slot-die printing achieves 9.10% for PTB7-Th:PC71 BM and 9.77% for PBDB-T:ITIC based on 1.04 cm2 . Furthermore, 15 cm2 flexible modules with effective efficiency up to 7.58% (PTB7-Th:PC71 BM) and 8.90% (PBDB-T:ITIC) are demonstrated with satisfying mechanical flexibility and operating stability. More importantly, this work outlines the shear impulse translation for organic printing electronics.
Collapse
Affiliation(s)
- Xiangchuan Meng
- College of Chemistry, Nanchang University, 999 Xuefu Avenue, Nanchang, 330031, China
- Institute of Polymers and Energy Chemistry, Nanchang University, 999 Xuefu Avenue, Nanchang, 330031, China
| | - Lin Zhang
- State Key Laboratory for Mechanical Behavior of Materials, Xi'an Jiaotong University, 28 Xianning West Road, Xi'an, 710049, China
| | - Yuanpeng Xie
- School of Chemistry, Beihang University, 37 Xueyuan Road, Beijing, 100191, China
| | - Xiaotian Hu
- College of Chemistry, Nanchang University, 999 Xuefu Avenue, Nanchang, 330031, China
- Institute of Polymers and Energy Chemistry, Nanchang University, 999 Xuefu Avenue, Nanchang, 330031, China
| | - Zhi Xing
- College of Chemistry, Nanchang University, 999 Xuefu Avenue, Nanchang, 330031, China
| | - Zengqi Huang
- College of Chemistry, Nanchang University, 999 Xuefu Avenue, Nanchang, 330031, China
| | - Cong Liu
- College of Chemistry, Nanchang University, 999 Xuefu Avenue, Nanchang, 330031, China
| | - Licheng Tan
- College of Chemistry, Nanchang University, 999 Xuefu Avenue, Nanchang, 330031, China
- Institute of Polymers and Energy Chemistry, Nanchang University, 999 Xuefu Avenue, Nanchang, 330031, China
| | - Weihua Zhou
- College of Chemistry, Nanchang University, 999 Xuefu Avenue, Nanchang, 330031, China
- Institute of Polymers and Energy Chemistry, Nanchang University, 999 Xuefu Avenue, Nanchang, 330031, China
| | - Yanming Sun
- School of Chemistry, Beihang University, 37 Xueyuan Road, Beijing, 100191, China
| | - Wei Ma
- State Key Laboratory for Mechanical Behavior of Materials, Xi'an Jiaotong University, 28 Xianning West Road, Xi'an, 710049, China
| | - Yiwang Chen
- College of Chemistry, Nanchang University, 999 Xuefu Avenue, Nanchang, 330031, China
- Institute of Polymers and Energy Chemistry, Nanchang University, 999 Xuefu Avenue, Nanchang, 330031, China
- Institute of Advanced Scientific Research (iASR), Jiangxi Normal University, 99 Ziyang Avenue, Nanchang, 330022, China
| |
Collapse
|
11
|
Greco C, Melnyk A, Kremer K, Andrienko D, Daoulas KC. Generic Model for Lamellar Self-Assembly in Conjugated Polymers: Linking Mesoscopic Morphology and Charge Transport in P3HT. Macromolecules 2019; 52:968-981. [PMID: 30792553 PMCID: PMC6376450 DOI: 10.1021/acs.macromol.8b01863] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/29/2018] [Revised: 12/17/2018] [Indexed: 01/27/2023]
Abstract
We develop a generic coarse-grained model of soluble conjugated polymers, capable of describing their self-assembly into a lamellar mesophase. Polymer chains are described by a hindered-rotation model, where interaction centers represent entire repeat units, including side chains. We introduce soft anisotropic nonbonded interactions to mimic the potential of mean force between atomistic repeat units. The functional form of this potential reflects the symmetry of the molecular order in a lamellar mesophase. The model can generate both nematic and lamellar (sanidic smectic) molecular arrangements. We parametrize this model for a soluble conjugated polymer poly(3-hexylthiophene) (P3HT) and demonstrate that the simulated lamellar mesophase matches morphologies of low molecular weight P3HT, experimentally observed at elevated temperatures. A qualitative charge-transport model allows us to link local chain conformations and mesoscale order to charge transport. In particular, it shows how coarsening of lamellar domains and chain extension increase the charge carrier mobility. By modeling large systems and long chains, we can capture transport between lamellar layers, which is due to rare, but thermodynamically allowed, backbone bridges between neighboring layers.
Collapse
Affiliation(s)
- Cristina Greco
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - Anton Melnyk
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - Kurt Kremer
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - Denis Andrienko
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - Kostas Ch. Daoulas
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| |
Collapse
|