Thermal conductivity of aligned polymers with kinks

This paper theoretically demonstrates that thermal conductivity in aligned polymers with kinks exhibits superdiffusive scaling at long lengths and non-monotonic behavior at shorter scales due to a transition from ballistic transport to Anderson localization, providing a microscopic explanation consistent with experimental and simulation data.

The Gist

Imagine a long, flexible rope made of tiny beads connected by strong springs. In the world of materials science, this rope represents a polymer chain (like the stuff in plastic or DNA). When you line up millions of these ropes perfectly straight, they become incredible at conducting heat—sometimes even better than metals. This is because heat travels through them as vibrations (called "phonons") bouncing along the strong springs connecting the beads.

However, real ropes aren't perfect. They have little twists and turns called "kinks." This paper investigates what happens to the heat flow when these kinks are scattered randomly along the rope.

Here is the story of their findings, broken down into simple concepts:

The Three Stages of a Heat Journey

The researchers found that how well heat travels depends entirely on how long the rope is. The behavior changes in three distinct stages, like a runner facing different types of terrain:

1. The Short Run: The "Highway" Effect
When the rope is very short, almost all the heat vibrations can zoom through without hitting any obstacles. It's like driving on a perfectly clear, straight highway. Because the vibrations travel so fast and unimpeded (a state called "ballistic transport"), the ability to conduct heat actually increases as the rope gets slightly longer. More road means more traffic can flow.

2. The Middle Run: The "Traffic Jam"
As the rope gets longer, the random kinks start to cause problems. Imagine driving down a road where the pavement suddenly shifts left or right at random spots. The heat vibrations start to bounce around, get confused, and eventually get stuck in one spot. In physics, this is called "Anderson localization."
Instead of flowing forward, the heat gets trapped. Because of this, the rope's ability to conduct heat drops dramatically (by about four times) as it gets longer. It's like a traffic jam where the cars (heat) can't move forward no matter how long the road is.

3. The Long Run: The "Super-Runner"
If you make the rope incredibly long, something surprising happens again. The heat finds a way to sneak through. The very slow, long-wavelength vibrations (like a giant, slow wave rolling through the whole rope) are less affected by the small kinks. They manage to bypass the traffic jams.
At this extreme length, the heat flow starts to increase again, but it follows a specific, slow mathematical rule (scaling with the cube root of the length). It's not a superhighway anymore, but it's a steady, super-efficient path that only the "super-runners" (long waves) can use.

The "Fence" Analogy

To understand this, the authors used a model they call a "fence." Imagine a fence where the posts are the atoms in the polymer.

• The Kinks: Sometimes, the fence posts are twisted, forcing the fence line to bend.
• The Constraint: The fence is built inside a narrow corridor. The posts can wiggle a little, but they can't wander too far off the path.
• The Result: The researchers found that if the fence wiggles too much (large sideways movement) or the twists are too sharp, the heat flow gets crushed. But if the fence stays relatively straight, the heat can eventually find its way through, even with the twists.

Why This Matters (According to the Paper)

The paper explains why some experiments show heat flowing better in longer chains, while others show it getting worse. The answer is: it depends on the length of the chain.

• Short chains: Heat flows better as they get longer.
• Medium chains: Heat flow gets worse because the kinks trap the energy.
• Very long chains: Heat flow improves again, but only for specific types of vibrations.

The authors also note that the "kinks" (twists in the molecular chain) are the main culprit. If you can control how much these chains twist or how much they wiggle sideways, you can control how much heat they carry. This helps explain why aligning polymer fibers (making them straighter) makes them conduct heat so much better.

In summary: Heat traveling through a twisted polymer chain is like a traveler navigating a road with random detours. At first, the road is clear. Then, the detours cause a massive jam. But if the road is long enough, the traveler finds a secret, slow path that gets them to the destination anyway.

Technical Summary

1. Problem Statement

Aligned polymer molecules exhibit exceptionally high thermal conductivity along their alignment direction due to strong covalent bonding, potentially rivaling metals. However, this property is highly sensitive to molecular conformation. Specifically, gauche kinks (structural defects where the polymer backbone rotates) can suppress heat transport by orders of magnitude.

While experiments and molecular dynamics (MD) simulations suggest that thermal conductivity ($\kappa$) in 1D systems often follows a superdiffusive scaling law ($\kappa \propto L^\beta$ with $0 < \beta < 1$), the microscopic mechanisms governing this behavior in the presence of kinks remain unclear. Existing theories struggle to explain:

• Why some experiments show a decrease in conductivity at very long chain lengths (contradicting simple superdiffusive models).
• The non-monotonic behavior of conductivity observed in certain polymers like Poly(methyl acrylate) (PMA).
• The dominance of defect scattering over anharmonic effects at room temperature.

2. Methodology

The authors developed a theoretical framework to investigate phonon transport in kinked polymers using the following approach:

• The "Fence" Model: They utilized a 2D fence model (previously introduced in Ref. 60) extended to include gauche kinks.
  • The model represents an alkane chain where transverse displacements are restricted to a width $2W$.
  • Kinks are modeled as local rotations of the molecular axis by angles $\pm \phi_1$ at randomly selected sites.
  • Interactions are harmonic, defined by force constants $A$ (nearest neighbor) and $B$ (next-nearest neighbor), derived from Density Functional Theory (DFT) calculations of alkane chains.
• Simulation Parameters:
  • Kink angles tested: $\phi_1 = \pi/4$ and $\pi/3$.
  • Inter-kink distances ($l$) and transverse widths ($W$) were varied (e.g., $l=W=a$ and $l=W=2a$, where $a$ is the lattice period).
  • Configurations violating transverse constraints were discarded to simulate aligned systems.
• Computational Technique:
  • Transfer Matrix Method: Used to calculate phonon transmission coefficients ($T$) and reflection coefficients ($R$) for longitudinal (LA) and transverse (TA) acoustic phonon branches. This method was chosen for its computational efficiency and superior convergence at low frequencies compared to Green's function methods.
  • Landauer Formalism: Thermal conductivity was calculated using the generalized Landauer formula, integrating transmission over the phonon bandwidth. Anharmonic interactions were neglected, assuming low temperatures where defect scattering dominates.
• Defect Modeling: To interpret results, the authors modeled kinks as independent, non-intersecting defects with a specific scattering profile, approximating transmission through a chain of length $L$ using the formula $T(\omega) \approx 1/\cosh^2(nLR(\omega)/2)$.

3. Key Contributions

• Identification of Non-Monotonic Behavior: The study reveals that thermal conductivity in kinked polymers does not follow a single power law. Instead, it exhibits a non-monotonic dependence on chain length, characterized by three distinct regimes.
• Mechanism Elucidation: The paper provides a microscopic explanation for the competition between ballistic transport and Anderson localization. It demonstrates how the interplay between these mechanisms, mediated by the scattering of longitudinal and transverse phonons, dictates the scaling behavior.
• Universal Scaling at Long Lengths: The authors confirm that at sufficiently long lengths, heat transport becomes superdiffusive with a universal scaling of $\kappa \propto L^{1/3}$. This regime is governed by the scattering of longitudinal phonons into transverse modes by defects.
• Role of Transverse Constraints: The study highlights that the transverse displacement limit ($W$) and the kink angle ($\phi_1$) are critical parameters. Increasing these parameters strongly suppresses thermal conductivity, offering a theoretical basis for controlling heat transport via molecular alignment.

4. Key Results

The length dependence of thermal conductivity ($\kappa$) was found to follow three distinct regimes:

1. Short Lengths (Ballistic Regime):

   • $\kappa$ increases with length ($L$).
   • Mechanism: Most phonons propagate ballistically because the chain is shorter than the mean free path. The number of conducting modes increases with length.
   • Scaling: Consistent with experimental observations of $\kappa \propto L^{0.38}$ to $L^{0.46}$.

2. Intermediate Lengths (Localization Regime):

   • $\kappa$ decreases with length (by a factor of ~4 in simulations).
   • Mechanism: Anderson localization dominates. Scattering by kinks suppresses the transmission of the majority of phonon modes (particularly transverse modes).
   • Observation: This explains the decrease in conductivity observed in PMA at long chain lengths, which contradicts simple ballistic or purely superdiffusive models.

3. Long Lengths (Superdiffusive Regime):

   • $\kappa$ increases again with length.
   • Mechanism: Transport is dominated by a narrow window of low-frequency, long-wavelength phonons that are weakly scattered.
   • Scaling: The system enters the universal superdiffusive regime where $\kappa \propto L^{1/3}$. This is governed by longitudinal-to-transverse phonon scattering.

Additional Findings:

• Transverse Width Sensitivity: Chains with larger transverse widths ($W=2a$) showed lower overall conductivity and a less pronounced initial increase compared to narrower chains ($W=a$).
• Defect Density: The density of kinks and the specific geometry of the "defect" (modeled as a four-kink bridge) significantly alter the reflection coefficients and the transition points between regimes.

5. Significance

• Resolving Experimental Discrepancies: The paper reconciles conflicting experimental data (e.g., increasing vs. decreasing conductivity with length) by showing that the observed behavior depends on the specific length scale relative to the localization length and the dominance of specific phonon modes.
• Material Design: The findings suggest that thermal management in polymer-based materials can be tuned by controlling molecular alignment and the density of gauche kinks. This is crucial for applications in thermal management, heat exchangers, and energy conversion.
• Theoretical Advancement: By establishing the $L^{1/3}$ scaling as the asymptotic limit for kinked polymers, the work provides a rigorous theoretical foundation for 1D heat transport that accounts for both structural disorder and phonon mode coupling.
• Biological Implications: The authors note that similar mechanisms (kink formation and misalignment) could explain thermal conductivity changes in biological membranes (lipid phase transitions), suggesting new avenues for investigation using 2D infrared spectroscopy.

In conclusion, this work demonstrates that thermal transport in aligned polymers is a complex, length-dependent phenomenon where Anderson localization initially suppresses conductivity before superdiffusive transport of long-wavelength phonons takes over at macroscopic scales.