skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Unlocking enhanced thermal conductivity in polymer blends through active learning
Abstract Polymers play an integral role in various applications, from everyday use to advanced technologies. In the era of machine learning (ML), polymer informatics has become a vital field for efficiently designing and developing polymeric materials. However, the focus of polymer informatics has predominantly centered on single-component polymers, leaving the vast chemical space of polymer blends relatively unexplored. This study employs a high-throughput molecular dynamics (MD) simulation combined with active learning (AL) to uncover polymer blends with enhanced thermal conductivity (TC) compared to the constituent single-component polymers. Initially, the TC of about 600 amorphous single-component polymers and 200 amorphous polymer blends with varying blending ratios are determined through MD simulations. The optimal representation method for polymer blends is identified, which involves a weighted sum approach that extends existing polymer representation from single-component polymers to polymer blends. An AL framework, combining MD simulation and ML, is employed to explore the TC of approximately 550,000 unlabeled polymer blends. The AL framework proves highly effective in accelerating the discovery of high-performance polymer blends for thermal transport. Additionally, we delve into the relationship between TC, radius of gyration (Rg), and hydrogen bonding, highlighting the roles of inter- and intra-chain interactions in thermal transport in amorphous polymer blends. A significant positive association between TC andRgimprovement and an indirect contribution from H-bond interaction to TC enhancement are revealed through a log-linear model and an odds ratio calculation, emphasizing the impact of increasingRgand H-bond interactions on enhancing polymer blend TC.  more » « less
Award ID(s):
2332270
PAR ID:
10562601
Author(s) / Creator(s):
;
Publisher / Repository:
Nature
Date Published:
Journal Name:
npj Computational Materials
Volume:
10
Issue:
1
ISSN:
2057-3960
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Current research on ferroelectric polymers centers predominantly on poly(vinylidene fluoride) (PVDF)–based fluoropolymers because of their superior performance. However, they are considered “forever chemicals” with environmental concerns. We describe a family of rationally designed fluorine-free ferroelectric polymers, featuring a polyoxypropylene main chain and disulfonyl alkyl side chains with a C3 spacer: −SO2CH2CHRCH2SO2− (R = −H or −CH3). Both experimental and simulation results demonstrate that strong dipole-dipole interactions between neighboring disulfonyl groups induce ferroelectric ordering in the condensed state, which can be tailored by changing the R group: ferroelectric for R = −H or relaxor ferroelectric for R = −CH3. At low electric fields, the relaxor polymer exhibits electroactuation and electrocaloric performance comparable with those of state-of-the-art PVDF-based tetrapolymers. 
    more » « less
  2. Debenedetti, P (Ed.)
    Using angle-resolved X-ray photoelectron spectroscopy, sum-frequency generation vibrational spectroscopy, contact angle measurements, and molecular dynamics simulations, we verify that the glass transition temperature (Tg) of polymer glass is lower near the free surface. However, the experimentalTg-gradients showed a linear variation with depth (z) from the free surface, while the simulated equilibriumTg-gradients exhibited a double exponentialz-dependence. In typical simulations,Tgis determined based on the relaxation time of the system reaching a prescribed threshold value at equilibrium. Conversely, the experiments determinedTgby observing the unfreezing of molecular mobility during heating from a kinetically arrested, nonequilibrium glassy state. To investigate the impact of nonequilibrium effects on theTg-gradient, we reduced the thermal annealing time in simulations, allowing the system to fall out of equilibrium. We observe a decrease in the relaxation time and the emergence of a modifiedz-dependence consistent with a linearTg-gradient near the free surface. We further validate the impact of nonequilibrium effects by studying the dependence of theTgon the heating/cooling rate for polymer films of varying thickness (h). Our experimental results reveal significant variations in theTg-heating/cooling rate dependence withhbelow the bulkTg, which are also observed in simulation when the simulated system is not equilibrated. We explain our findings by the reduction in mass density within the inner region of the system under nonequilibrium conditions, as observed in simulation, and recent research indicating a decrease in the localTgof a polymer when placed next to a softer material. 
    more » « less
  3. Abstract Despite extensive research on piezoelectric polymers since the discovery of piezoelectric poly(vinylidene fluoride) (PVDF) in 1969, the fundamental physics of polymer piezoelectricity has remained elusive. Based on the classic principle of piezoelectricity, polymer piezoelectricity should originate from the polar crystalline phase. Surprisingly, the crystal contribution to the piezoelectric strain coefficientd31is determined to be less than 10%, primarily owing to the difficulty in changing the molecular bond lengths and bond angles. Instead, >85% contribution is from Poisson's ratio, which is closely related to the oriented amorphous fraction (OAF) in uniaxially stretched films of semicrystalline ferroelectric (FE) polymers. In this perspective, the semicrystalline structure–piezoelectric property relationship is revealed using PVDF‐based FE polymers as a model system. In melt‐processed FE polymers, the OAF is often present and links the crystalline lamellae to the isotropic amorphous fraction. Molecular dynamics simulations demonstrate that the electrostrictive conformation transformation of the OAF chains induces a polarization change upon the application of either a stress (the direct piezoelectric effect) or an electric field (the converse piezoelectric effect). Meanwhile, relaxor‐like secondary crystals in OAF (SCOAF), which are favored to grow in the extended‐chain crystal (ECC) structure, can further enhance the piezoelectricity. However, the ECC structure is difficult to achieve in PVDF homopolymers without high‐pressure crystallization. We have discovered that high‐power ultrasonication can effectively induce SCOAFin PVDF homopolymers to improve its piezoelectric performance. Finally, we envision that the electrostrictive OAF mechanism should also be applicable for other FE polymers such as odd‐numbered nylons and piezoelectric biopolymers. 
    more » « less
  4. The designability of orthogonal coiled coil (CC) dimers, which draw on well‐established design rules, plays a pivotal role in fueling the development of CCs as synthetically versatile assembly‐directing motifs for the fabrication of bionanomaterials. Here, we aim to expand the synthetic CC toolkit through establishing a “minimalistic” set of orthogonal, de novo CC peptides that comprise 3.5 heptads in length and a single buried Asn to prescribe dimer formation. The designed sequences display excellent partner fidelity, confirmed via circular dichroism (CD) spectroscopy and Ni‐NTA binding assays, and are corroborated in silico using molecular dynamics (MD) simulation. Detailed analysis of the MD conformational data highlights the importance of interhelical E@g‐N@ainteractions in coordinating an extensive 6‐residue hydrogen bonding network that “locks” the interchain Asn‐Asn′ contact in place. The enhanced stability imparted to the Asn‐Asn′ bond elicits an increase in thermal stability of CCs up to ~15°C and accounts for significant differences in stability within the collection of similarly designed orthogonal CC pairs. The presented work underlines the utility of MD simulation as a tool for constructing de novo, orthogonal CCs, and presents an alternative handle for modulating the stability of orthogonal CCs via tuning the number of interhelical E@g‐N@acontacts. Expansion of CC design rules is a key ingredient for guiding the design and assembly of more complex, intricate CC‐based architectures for tackling a variety of challenges within the fields of nanomedicine and bionanotechnology. 
    more » « less
  5. 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. 
    more » « less