skip to main content

Title: Assessment of molecular dynamics simulations for amorphous poly(3-hexylthiophene) using neutron and X-ray scattering experiments
The molecular morphology and dynamics of conjugated polymers in the bulk solid state play a significant role in determining macroscopic charge transport properties. To understand this relationship, molecular dynamics (MD) simulations and quantum mechanical calculations are used to evaluate local electronic properties. In this work, we investigate the importance of system and simulation parameters, such as force fields and equilibration methods, when simulating amorphous poly(3-hexylthiophene) (P3HT), a model semiconducting polymer. An assessment of MD simulations for five different published P3HT force fields is made by comparing results to experimental wide-angle X-ray scattering (WAXS) and to a broad range of quasi-elastic neutron scattering (QENS) data. Moreover, an in silico analysis of force field parameters reveals that atomic partial charges and torsion potentials along the backbone and side chains have the greatest impact on structure and dynamics related to charge transport mechanisms in P3HT.
Authors:
; ; ; ; ; ; ; ;
Award ID(s):
1708317
Publication Date:
NSF-PAR ID:
10143430
Journal Name:
Soft Matter
Volume:
15
Issue:
25
Page Range or eLocation-ID:
5067 to 5083
ISSN:
1744-683X
Sponsoring Org:
National Science Foundation
More Like this
  1. Evaluating new, promising organic molecules to make next-generation organic optoelectronic devices necessitates the evaluation of charge carrier transport performance through the semi-conducting medium. In this work, we utilize quantum chemical calculations (QCC) and kinetic Monte Carlo (KMC) simulations to predict the zero-field hole mobilities of ∼100 morphologies of the benchmark polymer poly(3-hexylthiophene), with varying simulation volume, structural order, and chain-length polydispersity. Morphologies with monodisperse chains were generated previously using an optimized molecular dynamics force-field and represent a spectrum of nanostructured order. We discover that a combined consideration of backbone clustering and system-wide disorder arising from side-chain conformations are correlated with hole mobility. Furthermore, we show that strongly interconnected thiophene backbones are required for efficient charge transport. This definitively shows the role “tie-chains” play in enabling mobile charges in P3HT. By marrying QCC and KMC over multiple length- and time-scales, we demonstrate that it is now possible to routinely probe the relationship between molecular nanostructure and device performance.
  2. The field-effect electron mobility of aqueous solution-processed indium gallium oxide (IGO) thin-film transistors (TFTs) is significantly enhanced by polyvinyl alcohol (PVA) addition to the precursor solution, a >70-fold increase to 7.9 cm2/Vs. To understand the origin of this remarkable phenomenon, microstructure, electronic structure, and charge transport of IGO:PVA film are investigated by a battery of experimental and theoretical techniques, including In K-edge and Ga K-edge extended X-ray absorption fine structure (EXAFS); resonant soft X-ray scattering (R-SoXS); ultraviolet photoelectron spectroscopy (UPS); Fourier transform-infrared (FT-IR) spectroscopy; time-of-flight secondary-ion mass spectrometry (ToF-SIMS); composition-/processing-dependent TFT properties; high-resolution solid-state1H,71Ga, and115In NMR spectroscopy; and discrete Fourier transform (DFT) analysis with ab initio molecular dynamics (MD) liquid-quench simulations. The71Ga{1H} rotational-echo double-resonance (REDOR) NMR and other data indicate that PVA achieves optimal H doping with a Ga···H distance of ∼3.4 Å and conversion from six- to four-coordinate Ga, which together suppress deep trap defect localization. This reduces metal-oxide polyhedral distortion, thereby increasing the electron mobility. Hydroxyl polymer doping thus offers a pathway for efficient H doping in green solvent-processed metal oxide films and the promise of high-performance, ultra-stable metal oxide semiconductor electronics with simple binary compositions.

  3. Molybdenum disulfide (MoS 2 ) is a layered material with outstanding electrical and optical properties. Numerous studies evaluate the performance in sensors, catalysts, batteries, and composites that can benefit from guidance by simulations in all-atom resolution. However, molecular simulations remain difficult due to lack of reliable models. We introduce an interpretable force field for MoS 2 with record performance that reproduces structural, interfacial, and mechanical properties in 0.1% to 5% agreement with experiments. The model overcomes structural instability, deviations in interfacial and mechanical properties by several 100%, and empirical fitting protocols in earlier models. It is compatible with several force fields for molecular dynamics simulation, including the interface force field (IFF), CVFF, DREIDING, PCFF, COMPASS, CHARMM, AMBER, and OPLS-AA. The parameters capture polar covalent bonding, X-ray structure, cleavage energy, infrared spectra, bending stability, bulk modulus, Young's modulus, and contact angles with polar and nonpolar solvents. We utilized the models to uncover the binding mechanism of peptides to the MoS 2 basal plane. The binding strength of several 7mer and 8mer peptides scales linearly with surface contact and replacement of surface-bound water molecules, and is tunable in a wide range from −86 to −6 kcal mol −1 . The binding selectivitymore »is multifactorial, including major contributions by van-der-Waals coordination and charge matching of certain side groups, orientation of hydrophilic side chains towards water, and conformation flexibility. We explain the relative attraction and role of the 20 amino acids using computational and experimental data. The force field can be used to screen and interpret the assembly of MoS 2 -based nanomaterials and electrolyte interfaces up to a billion atoms with high accuracy, including multiscale simulations from the quantum scale to the microscale.« less
  4. Intracellular transport of chloride by members of the CLC transporter family involves a coupled exchange between a Cl− anion and a proton (H+), which makes the transport function dependent on ambient pH. Transport activity peaks at pH 4.5 and stalls at neutral pH. However, a structure of the WT protein at acidic pH is not available, making it difficult to assess the global conformational rearrangements that support a pH-dependent gating mechanism. To enable modeling of the CLC-ec1 dimer at acidic pH, we have applied molecular dynamics simulations (MD) featuring a new force field modification scheme—termed an Equilibrium constant pH approach (ECpH). The ECpH method utilizes linear interpolation between the force field parameters of protonated and deprotonated states of titratable residues to achieve a representation of pH-dependence in a narrow range of physiological pH values. Simulations of the CLC-ec1 dimer at neutral and acidic pH comparing ECpH-MD to canonical MD, in which the pH-dependent protonation is represented by a binary scheme, substantiates the better agreement of the conformational changes and the final model with experimental data from NMR, cross-link and AFM studies, and reveals structural elements that support the gate-opening at pH 4.5, including the key glutamates Gluin and Gluex.
  5. Peptide nucleic acids (PNAs) are a promising group of synthetic analogues of DNA and RNA that offer several distinct advantages over the naturally occurring nucleic acids for applications in biosensing, drug delivery, and nanoelectronics. Because of its structural differences from DNA/RNA, methods to analyze and assess the structure, conformations, and dynamics are needed. In this work, we develop synergistic techniques for the study of the PNA conformation. We use CuQ2, a Cu(II) complex with 8-hydroxyquinoline (HQ), as an alternative base pair and as a spin label in electron paramagnetic resonance (EPR) distance methods. We use molecular dynamics (MD) simulations with newly developed force field parameters for the spin labels to interpret the distance constraints determined by EPR. We complement these methods by UV–vis and circular dichroism measurements and assess the efficacy of the Cu(II) label on a PNA duplex whose backbone is based on aminoethylglycine and a duplex with a hydroxymethyl backbone modification. We show that the Cu(II) label functions efficiently within the standard PNA and the hydroxymethyl-modified PNA and that the MD parameters may be used to accurately reproduce our EPR findings. Through the combination of EPR and MD, we gain new insights into the PNA structure and conformationsmore »as well as into the mechanism of orientational selectivity in Cu(II) EPR at X-band. These results present for the first time a rigid Cu(II) spin label used for EPR distance measurements in PNA and the accompanying MD force fields for the spin label. Our studies also reveal that the spin labels have a low impact on the structure of the PNA duplexes. The combined MD and EPR approach represents an important new tool for the characterization of the PNA duplex structure and provides valuable information to aid in the rational application of PNA at large.« less