Abstract:The conductive polymer complex poly(3,4-ethylenedioxythiophene) polystyrene sulfonate (PEDOT:PSS) is integral to many next-generation polymeric devices. All-atom molecular dynamics (MD) is increasingly employed to probe film microstructure, but the commonly used Generalized AMBER Force Field (GAFF) parameters remain unvalidated for this system. Here, we develop all-atom force fields for both undoped and highly doped PEDOT compatible with GAFF. Molecular geometries, vibrational energies, torsional profiles, and conformational energetics calculated with density functional theory comprise the training data. The optimized force fields improve upon GAFF for all the training data, including notable improvements in torsional barriers, and capture the finite size of charge carriers in doped PEDOT. Although GAFF predicts similar behavior between undoped and doped polymers, our models identify meaningful differences in ground-state geometry, vibrational spectra, and torsional barriers between the two. Furthermore, vacuum MD simulations show that the neutral polymer is more flexible than the highly doped polymer, primarily due to the weaker inter-monomer dihedral potentials. The foregoing trends are directly attributable to the aromatic-to-quinoid transition PEDOT experiences upon doping. This coupling between doping and molecular flexibility highlights the importance of charge–geometry interactions in atomistic modeling and should be considered while modeling conjugated polymers beyond PEDOT.The Supporting Information is available free of charge at <a class="ext-link" href="/doi/10.1021/acs.macromol.1c00860?goto=supporting-info">https://pubs.acs.org/doi/10.1021/acs.macromol.1c00860</a>.Force field parameters (<a class="ext-link" href="/doi/suppl/10.1021/acs.macromol.1c00860/suppl_file/ma1c00860_si_001.zip">ZIP</a>),Objective function details for parameterization, assessment of vibrational mode overlap for large molecules, comparison of neutral and doped PEDOT vibrational spectra, energy-matching results for neutral PEDOT<sub>6</sub>, and effect of bonded and nonbonded parameters on flexibility of doped PEDOT (<a class="ext-link" href="/doi/suppl/10.1021/acs.macromol.1c00860/suppl_file/ma1c00860_si_002.pdf">PDF</a>)This article has not yet been cited by other publications.

Data-Driven Parametrization of All-atom force fields for Organic Semiconductors

PhyNEO: A Neural-Network-Enhanced Physics-Driven Force Field Development Workflow for Bulk Organic Molecule and Polymer Simulations.

A Scalable Graph Neural Network Method For Developing An Accurate Force Field Of Large Flexible Organic Molecules

Construction of atomistic force field model for perfluorosulfonic acid polymer

Force Field Development for Organic Molecules: Modifying Dihedral and 1-N Pair Interaction Parameters.

Data-Driven Parametrization of Molecular Mechanics Force Fields for Expansive Chemical Space Coverage

An Extendible, Graph-Neural-Network-Based Approach for Accurate Force Field Development of Large Flexible Organic Molecules

Parameterization of General Organic Polymers within the Open Force Field Framework

Parametrization of Non-Bonded Force Field Terms for Metal-Organic Frameworks Using Machine Learning Approach

De Novo Calculation of the Charge Carrier Mobility in Amorphous Small Molecule Organic Semiconductors

AUTOMATED FORCE FIELD PARAMETERIZATION FOR NON-POLARIZABLE AND POLARIZABLE ATOMIC MODELS BASED ON AB INITIO TARGET DATA.

Strategies for the Development of Conjugated Polymer Molecular Dynamics Force Fields Validated with Neutron and X-ray Scattering

Automated Force Field Developer and Optimizer Platform: Torsion Reparameterization

System-specific parameter optimization for non-polarizable and polarizable force fields

Atomistic Modeling of PEDOT:PSS Complexes II: Force Field Parameterization

Generating High-Precision Force Fields for Molecular Dynamics Simulations to Study Chemical Reaction Mechanisms using Molecular Configuration Transformer

Development and Comprehensive Benchmark of a High Quality AMBER-Consistent Small Molecule Force Field with Broad Chemical Space Coverage for Molecular Modeling and Free Energy Calculation

Broadening access to small-molecule parameterization with the force field toolkit

Open Force Field BespokeFit: Automating Bespoke Torsion Parametrization at Scale

Parameterization of prototype organic small molecules suitable for OPVs and molecular dynamics simulations: the BTT and BPT cases

MACE-OFF23: Transferable Machine Learning Force Fields for Organic Molecules