1
|
Bhat V, Callaway CP, Risko C. Computational Approaches for Organic Semiconductors: From Chemical and Physical Understanding to Predicting New Materials. Chem Rev 2023. [PMID: 37141497 DOI: 10.1021/acs.chemrev.2c00704] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/06/2023]
Abstract
While a complete understanding of organic semiconductor (OSC) design principles remains elusive, computational methods─ranging from techniques based in classical and quantum mechanics to more recent data-enabled models─can complement experimental observations and provide deep physicochemical insights into OSC structure-processing-property relationships, offering new capabilities for in silico OSC discovery and design. In this Review, we trace the evolution of these computational methods and their application to OSCs, beginning with early quantum-chemical methods to investigate resonance in benzene and building to recent machine-learning (ML) techniques and their application to ever more sophisticated OSC scientific and engineering challenges. Along the way, we highlight the limitations of the methods and how sophisticated physical and mathematical frameworks have been created to overcome those limitations. We illustrate applications of these methods to a range of specific challenges in OSCs derived from π-conjugated polymers and molecules, including predicting charge-carrier transport, modeling chain conformations and bulk morphology, estimating thermomechanical properties, and describing phonons and thermal transport, to name a few. Through these examples, we demonstrate how advances in computational methods accelerate the deployment of OSCsin wide-ranging technologies, such as organic photovoltaics (OPVs), organic light-emitting diodes (OLEDs), organic thermoelectrics, organic batteries, and organic (bio)sensors. We conclude by providing an outlook for the future development of computational techniques to discover and assess the properties of high-performing OSCs with greater accuracy.
Collapse
Affiliation(s)
- Vinayak Bhat
- Department of Chemistry & Center for Applied Energy Research, University of Kentucky, Lexington, Kentucky 40506-0055, United States
| | - Connor P Callaway
- Department of Chemistry & Center for Applied Energy Research, University of Kentucky, Lexington, Kentucky 40506-0055, United States
| | - Chad Risko
- Department of Chemistry & Center for Applied Energy Research, University of Kentucky, Lexington, Kentucky 40506-0055, United States
| |
Collapse
|
2
|
Borzdun N, Glova A, Larin S, Lyulin S. Influence of Asphaltene Modification on Structure of P3HT/Asphaltene Blends: Molecular Dynamics Simulations. NANOMATERIALS (BASEL, SWITZERLAND) 2022; 12:2867. [PMID: 36014732 PMCID: PMC9413297 DOI: 10.3390/nano12162867] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 07/20/2022] [Revised: 08/08/2022] [Accepted: 08/18/2022] [Indexed: 06/15/2023]
Abstract
Further development and commercialization of bulk heterojunction (BHJ) solar cells require the search for novel low-cost materials. The present study addresses the relations between the asphaltenes' chemical structure and the morphology of the poly(3-hexylthiohene) (P3HT)/asphaltene blends as potential materials for the design of BHJ solar cells. By means of all-atom molecular dynamics simulations, the formation of heterophase morphology is observed for the P3HT-based blends with carboxyl-containing asphaltenes, as well as the aggregation of the asphaltenes into highly ordered stacks. Although the π-π interactions between the polyaromatic cores of the asphaltenes in solutions are sufficient for the molecules to aggregate into ordered stacks, in a blend with a conjugated polymer, additional stabilizing factors are required, such as hydrogen bonding between carboxyl groups. It is found that the asphaltenes' aliphatic side groups may improve significantly the miscibility between the polymer and the asphaltenes, thereby preventing the formation of heterophase morphology. The results also demonstrate that the carboxyl-containing asphaltenes/P3HT ratio should be at least 1:1, as a decrease in concentration of the asphaltenes leads to the folding of the polymer chains, lower ordering in the polymer phase and the destruction of the interpenetrating 3D structure formed by P3HT and the asphaltene phases. Overall, the results of the present study for the first time reveal the aggregation behavior of the asphaltenes of varying chemical structures in P3HT, as well the influence of their presence and concentration on the polymer phase structure and blend morphology, paving the way for future development of BHJ solar cells based on the conjugated polymer/asphaltene blends.
Collapse
Affiliation(s)
- Natalia Borzdun
- Institute of Macromolecular Compounds, Russian Academy of Sciences, Bolshoi pr. 31 (V.O.), 199004 St. Petersburg, Russia
| | - Artyom Glova
- Institute of Macromolecular Compounds, Russian Academy of Sciences, Bolshoi pr. 31 (V.O.), 199004 St. Petersburg, Russia
| | - Sergey Larin
- Institute of Macromolecular Compounds, Russian Academy of Sciences, Bolshoi pr. 31 (V.O.), 199004 St. Petersburg, Russia
- Faculty of Physics, St. Petersburg State University, Ulyanovskaya str. 1–3, Peterhof, 198504 St. Petersburg, Russia
| | - Sergey Lyulin
- Institute of Macromolecular Compounds, Russian Academy of Sciences, Bolshoi pr. 31 (V.O.), 199004 St. Petersburg, Russia
| |
Collapse
|
3
|
Wu ECK, Salamat CZ, Tolbert SH, Schwartz BJ. Molecular Dynamics Study of the Thermodynamics of Integer Charge Transfer vs Charge-Transfer Complex Formation in Doped Conjugated Polymers. ACS APPLIED MATERIALS & INTERFACES 2022; 14:26988-27001. [PMID: 35657331 DOI: 10.1021/acsami.2c06449] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
Molecular dopants such as 2,3,5,6-tetrafluoro-7,7,8,8-tetracyanoquinodimethane (F4TCNQ) can interact with conjugated polymers such as poly(3-hexylthiophene-2,5-diyl) (P3HT) in two different ways: they can undergo integer charge transfer (ICT) or they can form a partial-charge-transfer complex (CTC). Both are seen experimentally, but the CTC has been challenging to characterize, making it difficult to answer questions such as the following. Which polymorph is more stable? Do they have similar barriers for formation? Is there a thermodynamic route to convert one to the other? Here, we study the structure and the thermodynamics of bulk F4TCNQ-doped P3HT with all-atom molecular dynamics simulations, using thermodynamic integration to calculate the relative free energies. We find that the ICT and CTC polymorphs have similar thermodynamic stabilities. The barrier to create the ICT polymorph, however, is lower than that to make the CTC polymorph, because the ICT polymorph has a small critical nucleus, but the critical nucleus for the CTC polymorph is larger than what we can simulate. Moreover, simulated thermal annealing shows that the activation barrier for converting the CTC polymorph to the ICT polymorph is relatively modest. Overall, the simulations explain both the observed structures and the thermodynamics of F4TCNQ-doped P3HT and offer guidelines for targeting the production of a desired polymorph for different applications.
Collapse
Affiliation(s)
- Eric Chih-Kuan Wu
- Department of Chemistry and Biochemistry University of California, Los Angeles Los Angeles, California 90095-1569, United States
| | - Charlene Z Salamat
- Department of Chemistry and Biochemistry University of California, Los Angeles Los Angeles, California 90095-1569, United States
| | - Sarah H Tolbert
- Departments of Chemistry and Biochemistry and Materials Science and Engineering University of California, Los Angeles Los Angeles, California 90095-1569, United States
| | - Benjamin J Schwartz
- Department of Chemistry and Biochemistry University of California, Los Angeles Los Angeles, California 90095-1569, United States
| |
Collapse
|
4
|
Casalegno M, Famulari A, Meille SV. Modeling of Poly(3-hexylthiophene) and Its Oligomer’s Structure and Thermal Behavior with Different Force Fields: Insights into the Phase Transitions of Semiconducting Polymers. Macromolecules 2022. [DOI: 10.1021/acs.macromol.2c00131] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Mosè Casalegno
- Dipartimento di Chimica, Materiali e Ingegneria Chimica “G. Natta”, Politecnico di Milano, via Mancinelli 7, I-20131 Milano (MI), Italy
| | - Antonino Famulari
- Dipartimento di Chimica, Materiali e Ingegneria Chimica “G. Natta”, Politecnico di Milano, via Mancinelli 7, I-20131 Milano (MI), Italy
| | - Stefano Valdo Meille
- Dipartimento di Chimica, Materiali e Ingegneria Chimica “G. Natta”, Politecnico di Milano, via Mancinelli 7, I-20131 Milano (MI), Italy
| |
Collapse
|
5
|
Temperature-Dependent Conformation Behavior of Isolated Poly(3-hexylthiopene) Chains. Polymers (Basel) 2022; 14:polym14030550. [PMID: 35160539 PMCID: PMC8840214 DOI: 10.3390/polym14030550] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/03/2021] [Revised: 01/18/2022] [Accepted: 01/25/2022] [Indexed: 02/01/2023] Open
Abstract
We use atomistic as well as coarse-grained molecular dynamics simulations to study the conformation of a single poly(3-hexylthiopene) chain as a function of temperature. We find that mainly bundle and toroid structures appear with bundles becoming more abundant for decreasing temperatures. We compare an atomistic and a Martini-based coarse-grained model which we find in very good agreement. We further illustrate how the temperature dependence of P3HT can be connected to that of simple Lennard–Jones model polymers in a vacuum. Upon adding solvent (THF) we observe the occurrence of a prominent swelling of the molecular size at a temperature of about 220 K. This swelling is in close agreement with the interpretation of recent spectroscopic experiments which allows us to explain the experimental observations by an increased frequency of bundle structures.
Collapse
|
6
|
Atomistic simulation of poly (lactic acid) of different regioregularity. POLYMER 2021. [DOI: 10.1016/j.polymer.2021.123577] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
|
7
|
Tsourtou FD, Peristeras LD, Apostolov R, Mavrantzas VG. Molecular Dynamics Simulation of Amorphous Poly(3-hexylthiophene). Macromolecules 2020. [DOI: 10.1021/acs.macromol.0c00454] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Flora D. Tsourtou
- Department of Chemical Engineering, University of Patras & FORTH-ICE/HT, GR 26504 Patras, Greece
| | - Loukas D. Peristeras
- Molecular Thermodynamics and Modelling of Materials Laboratory, Institute of Nanoscience and Nanotechnology, National Center for Scientific Research “Demokritos”, GR-15310 Agia Paraskevi
Attikis, Athens, Greece
| | | | - Vlasis G. Mavrantzas
- Department of Chemical Engineering, University of Patras & FORTH-ICE/HT, GR 26504 Patras, Greece
- Department of Mechanical and Process Engineering, Particle Technology Laboratory, ETH Zürich, CH-8092 Zürich, Switzerland
| |
Collapse
|
8
|
Glova AD, Melnikova SD, Mercurieva AA, Larin SV, Lyulin SV. Grafting-Induced Structural Ordering of Lactide Chains. Polymers (Basel) 2019; 11:polym11122056. [PMID: 31835722 PMCID: PMC6961058 DOI: 10.3390/polym11122056] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/07/2019] [Revised: 12/06/2019] [Accepted: 12/09/2019] [Indexed: 01/05/2023] Open
Abstract
The structure of a grafted layer of lactide chains in the “dry brush” regime immersed in a melt of chemically similar polymer was examined while varying graft lengths. To this end, microsecond atomistic molecular dynamics simulations were performed. Almost no influence of graft length on the fraction of the grafted chains backfolded to the grafting surface was found. However, a structural ordering was unexpectedly observed in the system when the length of the grafted lactide chains was close to approximately 10 Kuhn segments. This ordering of the grafts is characterized by the formation of helical fragments whose structure is in good agreement with the experimental data for the α crystal of the lactide chains. Both the backfolding and the structural ordering may be viewed as the initial stage of the crystallization of the layer of grafted lactide chains. In contrast to the known behavior for conventional polymer brushes in the “dry brush” regime, the structure of the grafted lactide chains can be either amorphous or ordered, depending on the graft length N and the grafting density σ when their product Nσ is fixed.
Collapse
Affiliation(s)
- Artyom D. Glova
- Institute of Macromolecular Compounds, Russian Academy of Sciences, Bolshoj pr. 31 (V.O.), 199004 St. Petersburg, Russia; (A.D.G.); (A.A.M.); (S.V.L.)
| | - Sofya D. Melnikova
- Institute of Physics, Nanotechnology and Telecommunications, Peter the Great St. Petersburg Polytechnic University, Polytechnicheskaya st. 29, 195251 St. Petersburg, Russia;
| | - Anna A. Mercurieva
- Institute of Macromolecular Compounds, Russian Academy of Sciences, Bolshoj pr. 31 (V.O.), 199004 St. Petersburg, Russia; (A.D.G.); (A.A.M.); (S.V.L.)
| | - Sergey V. Larin
- Institute of Macromolecular Compounds, Russian Academy of Sciences, Bolshoj pr. 31 (V.O.), 199004 St. Petersburg, Russia; (A.D.G.); (A.A.M.); (S.V.L.)
| | - Sergey V. Lyulin
- Institute of Macromolecular Compounds, Russian Academy of Sciences, Bolshoj pr. 31 (V.O.), 199004 St. Petersburg, Russia; (A.D.G.); (A.A.M.); (S.V.L.)
- Correspondence: ; Tel.: +7-812-323-0216
| |
Collapse
|
9
|
Miller ED, Jones ML, Henry MM, Chery P, Miller K, Jankowski E. Optimization and Validation of Efficient Models for Predicting Polythiophene Self-Assembly. Polymers (Basel) 2018; 10:E1305. [PMID: 30961230 PMCID: PMC6401914 DOI: 10.3390/polym10121305] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/02/2018] [Revised: 11/20/2018] [Accepted: 11/21/2018] [Indexed: 01/27/2023] Open
Abstract
We develop an optimized force-field for poly(3-hexylthiophene) (P3HT) and demonstrate its utility for predicting thermodynamic self-assembly. In particular, we consider short oligomer chains, model electrostatics and solvent implicitly, and coarsely model solvent evaporation. We quantify the performance of our model to determine what the optimal system sizes are for exploring self-assembly at combinations of state variables. We perform molecular dynamics simulations to predict the self-assembly of P3HT at ∼350 combinations of temperature and solvent quality. Our structural calculations predict that the highest degrees of order are obtained with good solvents just below the melting temperature. We find our model produces the most accurate structural predictions to date, as measured by agreement with grazing incident X-ray scattering experiments.
Collapse
Affiliation(s)
- Evan D Miller
- Micron School of Materials Science and Engineering, Boise State University, Boise, ID 83705, USA.
| | - Matthew L Jones
- Micron School of Materials Science and Engineering, Boise State University, Boise, ID 83705, USA.
| | - Michael M Henry
- Micron School of Materials Science and Engineering, Boise State University, Boise, ID 83705, USA.
| | - Paul Chery
- Physics, Macalester College, St. Paul, MN 55105, USA.
| | - Kyle Miller
- Physics, University of Puget Sound, Tacoma, WA 98416, USA.
| | - Eric Jankowski
- Micron School of Materials Science and Engineering, Boise State University, Boise, ID 83705, USA.
| |
Collapse
|
10
|
Tsourtou FD, Peroukidis SD, Peristeras LD, Mavrantzas VG. Monte Carlo Algorithm Based on Internal Bridging Moves for the Atomistic Simulation of Thiophene Oligomers and Polymers. Macromolecules 2018. [DOI: 10.1021/acs.macromol.8b01344] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Flora D. Tsourtou
- Department of Chemical Engineering, University of Patras and FORTH-ICE/HT, GR-26504 Patras, Greece
| | - Stavros D. Peroukidis
- Department of Chemical Engineering, University of Patras and FORTH-ICE/HT, GR-26504 Patras, Greece
- School of Science/Technology, Natural Sciences, Hellenic Open University, GR-26335 Patras, Greece
| | - Loukas D. Peristeras
- Institute of Nanoscience and Nanotechnology, Molecular Thermodynamics and Modelling of Materials Laboratory, National Center for Scientific Research “Demokritos”, GR-15310 Agia Paraskevi Attikis, Greece
| | - Vlasis G. Mavrantzas
- Department of Chemical Engineering, University of Patras and FORTH-ICE/HT, GR-26504 Patras, Greece
- Particle Technology Laboratory, Department of Mechanical and Process Engineering, ETH Zürich, CH-8092 Zürich, Switzerland
| |
Collapse
|
11
|
Glova AD, Falkovich SG, Dmitrienko DI, Lyulin AV, Larin SV, Nazarychev VM, Karttunen M, Lyulin SV. Scale-Dependent Miscibility of Polylactide and Polyhydroxybutyrate: Molecular Dynamics Simulations. Macromolecules 2017. [DOI: 10.1021/acs.macromol.7b01640] [Citation(s) in RCA: 35] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Affiliation(s)
- Artyom D. Glova
- Institute
of Macromolecular Compounds, Russian Academy of Sciences, Bolshoj
pr. V.O., 31, 199004 St. Petersburg, Russia
| | - Stanislav G. Falkovich
- Institute
of Macromolecular Compounds, Russian Academy of Sciences, Bolshoj
pr. V.O., 31, 199004 St. Petersburg, Russia
| | - Daniil I. Dmitrienko
- Faculty
of Physics, Saint-Petersburg University, Ulyanovskaya str. 1, Petrodvorets, 198504 St. Petersburg, Russia
| | - Alexey V. Lyulin
- Theory
of Polymers and Soft Matter Group, Technische Universiteit Eindhoven, PO Box 513, 5600 MB Eindhoven, The Netherlands
| | - Sergey V. Larin
- Institute
of Macromolecular Compounds, Russian Academy of Sciences, Bolshoj
pr. V.O., 31, 199004 St. Petersburg, Russia
| | - Victor M. Nazarychev
- Institute
of Macromolecular Compounds, Russian Academy of Sciences, Bolshoj
pr. V.O., 31, 199004 St. Petersburg, Russia
| | - Mikko Karttunen
- Institute
of Macromolecular Compounds, Russian Academy of Sciences, Bolshoj
pr. V.O., 31, 199004 St. Petersburg, Russia
- Department
of Chemistry and Department of Applied Mathematics, Western University, 1151 Richmond St., London, Ontario, Canada N6A 5B7
| | - Sergey V. Lyulin
- Institute
of Macromolecular Compounds, Russian Academy of Sciences, Bolshoj
pr. V.O., 31, 199004 St. Petersburg, Russia
- Faculty
of Physics, Saint-Petersburg University, Ulyanovskaya str. 1, Petrodvorets, 198504 St. Petersburg, Russia
| |
Collapse
|