Subscriber access provided by HACETTEPE UNIVERSITESI KUTUPHANESI
Article
Thermal Conductance in Crosslinked Polymers: Effects of Non-Bonding Interactions Vahid Rashidi, Eleanor J. Coyle, Katherine Sebeck, John Kieffer, and Kevin P. Pipe J. Phys. Chem. B, Just Accepted Manuscript • DOI: 10.1021/acs.jpcb.7b01377 • Publication Date (Web): 31 Mar 2017 Downloaded from http://pubs.acs.org on April 3, 2017
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
The Journal of Physical Chemistry B is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 33
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
The Journal of Physical Chemistry
Thermal Conductance in Crosslinked Polymers: Effects of Non-Bonding Interactions Vahid Rashidi,† Eleanor J. Coyle,‡ Katherine Sebeck,‡ John Kieffer,‡ and Kevin P. Pipe∗,†,¶ Department of Mechanical Engineering, University of Michigan, Ann Arbor, MI, 48109, Department of Materials Science and Engineering, University of Michigan, Ann Arbor, MI, 48109, and Department of Electrical Engineering and Computer Science, University of Michigan, Ann Arbor, MI, 48109 E-mail:
[email protected] Abstract Weak inter-chain interactions have been considered to be a bottleneck for heat transfer in polymers, while covalent bonds are believed to give a high thermal conductivity to polymer chains. For this reason, crosslinkers have been explored as a means to enhance polymer thermal conductivity; however, results have been inconsistent. Some studies show an enhancement in the thermal conductivity for polymers upon crosslinking, while others show the opposite trend. In this work we study the mechanisms of heat transfer in crosslinked polymers in order to understand the reasons for these discrepancies, in particular examining the relative contributions of covalent (referred to here as “bonding”) and non-bonding (e.g., van der Waals and electrostatic) interactions. ∗
To whom correspondence should be addressed Department of Mechanical Engineering, University of Michigan, Ann Arbor, MI, 48109 ‡ Department of Materials Science and Engineering, University of Michigan, Ann Arbor, MI, 48109 ¶ Department of Electrical Engineering and Computer Science, University of Michigan, Ann Arbor, MI, 48109 †
1
ACS Paragon Plus Environment
The Journal of Physical Chemistry
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Our results indicate crosslinkers enhance thermal conductivity primarily when they are short in length and thereby bring polymer chains closer to each other; this increases inter-chain heat transfer by enhancing non-bonding interactions between the chains (non-bonding interactions being highly dependent on inter-chain distance). This suggests that enhanced non-bonding interactions (rather than thermal pathways through crosslinker covalent bonds) are the primary transport mechanism by which thermal conductivity is increased. We further illustrate this by showing that energy from THz acoustic waves travels significantly faster in polymers when non-bonding interactions are included versus when only covalent interactions are present. These results help to explain prior studies that measure differing trends in thermal conductivity for polymers upon crosslinking with various species.
1. Introduction Polymers have been extensively studied in recent decades for many applications, such as organic field effect transistors and other electronic devices, 1,2 thermoelectric devices, 3–7 photovoltaic devices, 8 thermal devices, 9 lithography and patterning, 10 drug delivery, 11 photonic devices, 12 and plastics products. 13,14 In many of these applications, the rate of thermal transport through the polymer plays a significant role in device reliability and performance. 15,16 For example, enhancing the thermal conductivity of polymers (which is on the order of 0.1 W/mK for most amorphous polymers 17,18 ) could improve their applications in solar cells 16 and the plastics industry, 19 while for a better thermoelectric figure of merit a polymer with a lower thermal conductivity is desired. 20 Previous computational 17,19,21–24 and experimental 25–34 studies have largely shown that the thermal conductivity of aligned linear polymer chains in the direction of chain alignment is significantly higher than that of bulk amorphous polymers. This has been attributed to a relatively large thermal conductivity for covalent bonds. 17,29,35 Meanwhile, the low thermal conductivities of amorphous bulk polymers are generally believed to be due to a bottleneck associated with weak inter-chain forces such as 2
ACS Paragon Plus Environment
Page 2 of 33
Page 3 of 33
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
The Journal of Physical Chemistry
van der Waals (vdw) interactions, 18 which have been calculated to be an order of magnitude less efficient at transferring heat than covalent bonds. 36 Mechanisms to enhance thermal transport in polymers have traditionally been based on the inclusion of thermally conductive fillers 37–39 and more recently on enhanced inter-chain interactions. 18,40–43 Crosslinking of polymer chains through covalent bonds is a common synthesis approach for polymers and would seem to be a very compelling technique to enhance their thermal conductivity. 18,41,44–47 However, computational and experimental studies of thermal conductivity in crosslinked polymers have shown significant discrepancies. For example, experimental studies have shown 50% enhancement for polyethylene samples 41 at high crosslinking densities and approximately 30% enhancement for polystyrene at 20% crosslinking density, 46 while computational studies have predicted a threefold enhancement in the bulk thermal conductivity of polyethylene upon crosslinking 18 and almost no change in the thermal conductivity of polystyrene at 20% crosslinking density. Other studies have measured or predicted only a very small enhancement 18 or even a reduction in thermal conductivity upon crosslinking. 48–50 For example, experimental results by Yu et al. 49 show a nearly 30% reduction in the thermal conductivity of polyethylene upon crosslinking. Numerical simulations by Ni et al. 48 also suggest that crosslinking 10% of the carbon atoms in polyethylene chains results in a 44.2% reduction in bulk thermal conductivity. While both Yu et al. 49 and Ni et al. 48 predict a monotonic decrease in thermal conductivity with increasing crosslinking density for polyethylene, Kikugawa et al. 18 predict an opposite trend. Given the many applications of crosslinked polymers 11,12,51–60 and the importance of thermal conductivity to many of these applications, it is crucially important to build a comprehensive understanding of the effects of crosslinking on thermal transport. Thermal conduction in electrically insulating crosslinked polymers occurs through four mechanisms, as shown in Figure 1. These mechanisms may be closely coupled. For example, the length of the crosslinking agent is expected to strongly influence the spacing between linked polymer chains. This spacing in turn greatly influences heat flow via non-bonding in-
3
ACS Paragon Plus Environment
The Journal of Physical Chemistry
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment
Page 4 of 33
Page 5 of 33
The Journal of Physical Chemistry
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment
The Journal of Physical Chemistry
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 6 of 33
fixed in order to control the chain endpoints. The system was then evolved as an NVT ensemble for 1.5 ns at 300 K using the Langevin 81 thermostat followed by a 1.5 ns relaxation at 300K using the Nose-Hoover 82,83 thermostat to reach thermal equilibrium. Sample temperature and potential energy profiles during thermal equilibrium are shown in the Supporting Information. A temperature gradient was then established in the system under NVE conditions by pumping "heat" into one chain and removing it from the other chain at the same rate. This heat flux was generated by rescaling the velocities of atoms in each chain to stimulate or suppress thermal motion. The system was allowed to evolve for 3 ns to reach equilibrium. Thermal data was then collected from the system. For example, thermal conductance (TC) between two chains was calculated using Fourier's law of conduction: 63,64,67,84,85
TC =
q ∆T
(1)
where q is the thermal energy transferred between the two chains and ∆T is the difference between the average temperatures of the chains. This technique has been used for intertube thermal transport in one dimensional systems such as CNTs. 84 Sample temperature profiles are shown in the Supporting Information. All relaxation steps used free, non-periodic boundary conditions. After relaxation, boundaries were fixed for thermal data collection and trajectory analysis. The OPLS force field was used to model the interactions of PMMA chains. 86–88 The parameters are reported in the Supporting Information. This force field has been successfully used in other studies for heat transfer modeling. 89–94 The force fields used for PVA and PE are described in their relevant sections, and their parameters are reported in the Supporting Information. We have used a variety of force fields to ensure the results are not force field dependent. For non-bonding interactions, a cutoff distance longer than the size of the simulation box was used in order to account for all interactions between the atoms in the system and to avoid any drift in the total energy of the system. This inclusion of all interactions in our simulations is crucial for drawing clear conclusions regarding the 6
ACS Paragon Plus Environment
Page 7 of 33
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
The Journal of Physical Chemistry
effects of these interactions. Inter-chain distances were measured by calculating the average distance between the centers of mass in each chain using the coordinates of the atoms in each chain.
3. Results and Discussion 3.1. Effects of Crosslinking Species and Non-bonding Interactions on Inter-chain Thermal Transport To study the effects of crosslinking on thermal transport in polymers and to understand the differing results reported in the literature, we simulated systems of two parallel PMMA polymer chains consisting of 20 monomers in each chain (an icosamer), terminated with hydrogens (similar to Figure 2). We first varied the crosslinker species in PMMA chains with 100% crosslinking to assess the maximum impact of the crosslinking agent on interchain thermal conductance. Although 100% crosslinking density for PMMA may not be practically achievable, we include it in our study as it gives us an understanding of the theoretical maximum by which crosslinkers may enhance heat transfer. In the following sections we also study lower crosslinking densities that are achievable in practice. Note that such high crosslinking densities are often found in epoxies and resins. 95 The species used (shown in Figure 3) span a range of shapes, masses, bond strengths, lengths, and vibrational characteristics, all of which potentially play a role in thermal transport. Although most of the crosslinkers are real chemical structures, some such as "light benzene-1,4-diyl" (a benzene-1,4-diyl crosslinker in which the carbon masses are reduced to half of the true carbon mass) do not have a real-world chemical representation. In order to isolate the effect of crosslinker mass, we performed a parametric study in which all of the other parameters (chemical structure, charge, force field, etc.) were kept constant and only the mass of the crosslinker was varied. We performed a similar simulation to elucidate the effect of bond strength by representing a carbon atom that has only two (rather than four) covalent bonds 7
ACS Paragon Plus Environment
The Journal of Physical Chemistry
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment
Page 8 of 33
Page 9 of 33
The Journal of Physical Chemistry
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment
The Journal of Physical Chemistry
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
for systems in which non-bonding interactions are excluded (Figure 4(d)). Although the polymer chains get closer in this case, because there are no non-bonding interactions in the system, inter-chain thermal transport does not benefit from this reduction in the interchain distance. We note that the reduction in inter-chain distance due to the absence of the repulsive portion of the vdW interactions is more significant for systems crosslinked with benzene-1,4-diyl, since these repulsive forces are strong in benzene-1,4-diyl and the additional degrees of freedom associated with its ring structure enable a more significant collapse upon their removal. For systems in which all non-bonding interactions are included (squares, Figure 4(b)), chains crosslinked with CH2 (carbene) have the highest enhancement for inter-chain thermal conductance. This is because the crosslinkers bring the chains close to each other and hence enhance the non-bonding inter-chain interactions between the chains. The interchain distance for this crosslinker is smaller than the radius of effectiveness for non-bonding interactions (∼8-10Å or 2.5 times the Lennard Jones σ parameter; widely used in molecular simulations as a distance beyond which vdW interactions are negligible 61 ) and hence these interactions significantly enhance heat transfer between the chains. Moreover, there is an overlap between the vibrational density of states of the crosslinker and that of the main polymer chain, which further contributes to inter-chain heat transfer. This is discussed in more detail in the Supporting Information. When all non-bonding interactions are removed from the system, the only path for heat transfer between the chains is through the crosslinker covalent bonds. In this case we observe a large drop in inter-chain thermal conductance (triangles in Figure 4(d)). For example, a factor of ∼5 reduction in thermal conductance is observed for the system crosslinked with CH2 . This drastic drop in thermal conductance suggests that the majority of inter-chain thermal transport in crosslinked polymers occurs through enhanced inter-chain non-bonding interactions rather than thermal conduction along crosslinkers. Quantifying the number of bonding and non-bonding interactions offers another perspective on this effect, and is
10
ACS Paragon Plus Environment
Page 10 of 33
Page 11 of 33
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
The Journal of Physical Chemistry
discussed in Section 3.4. Above we have shown that crosslinking species that bring the chains closer have a higher impact on thermal conductance between the chains. It is interesting to consider whether inter-chain distance also has an influence on the portion of thermal transport between the chains that occurs through covalent bonds. To study whether the presence or absence of non-bonding interactions influences heat transfer through crosslinker covalent bonds, we calculated the thermal conductance between polymer chains crosslinked with carbene as a function of inter-chain distance in the absence of non-bonding interactions. Here, changing the inter-chain distance is meant to approximate the physical effect that non-bonding interactions may have on the crosslinked chains. These results are shown in Figure 5. For these simulations we manually moved the chains with respect to each other in order to change the inter-chain distance. This will result in some residual stress in the system. The range of inter-chain distances in Figure 5 is chosen to cover the range of inter-chain distances in Figure 4(b). In the absence of non-bonding interactions, we did not see any clear relation or enhancement in thermal conductance between the chains for different inter-chain distances. Thus, based on Figures 4(d) and 5, we conclude that the thermal conductance between the chains in the absence of non-bonding interactions is not a function of inter-chain distance. This indicates that heat conduction through crosslinker covalent bonds is not greatly perturbed by the presence or absence of non-bonding interactions, and thus can be considered approximately to be independent. Therefore, the large drop in thermal conductance observed in Figure 4(a) when non-bonding interactions are removed is expected to be due primarily to removal of heat transfer through non-bonding interactions rather than an alteration of heat conduction through crosslinker covalent bonds. Together, the results of Figures 4 and 5 suggest that the primary thermal function of a crosslinker is to bring chains closer together so that non-bonding interactions are increased. Eiermann 36 calculated that thermal conductance through vdW interactions is 10 times weaker than thermal transport through covalent bonds. While this applies for a single
11
ACS Paragon Plus Environment
The Journal of Physical Chemistry
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
vdW bond in comparison with a single covalent bond, our results indicate that inter-chain heat transfer through vdW interactions in larger systems is greater than that through covalent interactions due to the much larger number of vdW interactions. In addition, short crosslinker length enhances non-bonding interactions and hence increases their contribution to heat transfer. For example, for each PMMA monomer that forms a crosslink with another PMMA monomer, there is only one covalent bond between them to conduct heat. However, since each PMMA monomer has 15 atoms, there are 225 (15 × 15) heat transfer pathways through non-bonding interactions between the two crosslinked monomers. Although each non-bonding interaction transfers less heat than a covalent bond, the significantly larger number of them has a greater effect.
Figure 5: Change in thermal conductance (TC) between two PMMA polymer chains crosslinked with carbene crosslinkers in the absence of non-bonding interactions, with respect to the distance between the two chains. The importance of non-bonding interactions revealed here could explain results such as that of Kikugawa et al. 18 in which the thermal conductivity of simulated polystyrene was not calculated to increase even for high crosslinking densities. In this case the benzene side group in polystyrene may prevent the inter-chain distance from decreasing even in the presence of a large density of crosslinkers. Likewise, Yu et al. 49 reported a reduction in thermal conductivity for polyethylene sam12
ACS Paragon Plus Environment
Page 12 of 33
Page 13 of 33
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
The Journal of Physical Chemistry
ples crosslinked by dicumyl peroxide. The length of this crosslinking agent may lead to increased inter-chain distance and hence a reduction in the contribution of non-bonding interactions to thermal transport. This is supported by the reduced density they measure.
3.2. Acoustic Wave Propagation Velocity Analysis In Figure 4(b) we showed that the thermal conductance between two polymer chains crosslinked with CH2 (carbene) crosslinkers is approximately 5 times greater when non-bonding interactions are included. The long-range nature of such non-bonding interactions (versus the short-range nature of covalent bonds) has important implications for thermal transport. In this section, we show that acoustic waves (which carry heat) travel faster when non-bonding interactions are included versus when only covalent bonds are present. Figure 6 shows a crosslinked polymer system for simple acoustic wave propagation analysis, with input wave location and "probe" locations at which we detect the energy of the wave after it propagates from its initial point. We first reduce the temperature of the structure to 0K to remove all molecular motion. We then create a 1 THz oscillation input in the transverse direction with respect to the polymer backbone. The 1 THz oscillation is chosen to represent a wave that carries heat in the materials at representative temperatures. 63,96 The method used to detect the wave at the probe locations is discussed in detail in the Supplementary Information. In Table 1 we show the normalized wave propagation velocities for a 1 THz wave. The results show that a wave at 1 THz frequency travels approximately 9 times faster from one chain to another chain when non-bonding interactions are included. This wave propagation analysis agrees with our previous results in which we showed that thermal conductance between two chains is 5 times higher when non-bonding interactions are included versus when only covalent bonds conduct heat between the chains. This analysis
13
ACS Paragon Plus Environment
The Journal of Physical Chemistry
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment
Page 14 of 33
Page 15 of 33
The Journal of Physical Chemistry
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment
The Journal of Physical Chemistry
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Above 65% crosslinking density, the system with all interactions included shows a nonlinear increase in thermal conductance. This is due to the fact that the PMMA polymer chains simulated here have an intrinsically syndiotactic conformation, and when they are 100% crosslinked, the monomer side groups in the two chains all face each other and inter-chain distance is reduced drastically. As further discussed in the Supporting Information, we also find that inter-chain distance for the results shown in Figure 7 are inversely proportional with the inter-chain thermal conductance. This means that reducing crosslinking density results in a larger inter-chain distance which in turn weakens the non-bonding interactions between the chains.
3.4. Inter-chain Thermal Transport through Non-bonding Interactions In the previous section we have shown that non-bonding interactions have a more significant effect on thermal transport in crosslinked polymers than covalent bonds when the crosslinking agent reduces the inter-chain distance. Additionally, in Figure 7 we show that the thermal conductance between polymer chains scales approximately linearly with crosslinking density whether the non-bonding interactions are included or not. This means that the difference between the red and blue curves in Figure 7 also scales linearly with crosslinking density. The values obtained by subtracting the two curves in Figure 7 are proportional to the portion of thermal transport through non-bonding interactions. Therefore, even in polymers without crosslinkers, there should be a linear relationship between the number of interactions between polymer chains (which we quantify below) and the inter-chain thermal conductance. Here we simulate relatively long (120 monomers in each chain, terminated with hydrogens) pairs of different polymer chains (PMMA, PVA, and PE) to study the effect of inter-chain distance and entanglement on inter-chain thermal transport for polymers without crosslinkers. We used nanoHUB Polymer Modeler 78 to generate the polymer chains geometry and the force field provided to study PE chains. The force-field developed by 16
ACS Paragon Plus Environment
Page 16 of 33
Page 17 of 33
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
The Journal of Physical Chemistry
Muller 97 was used to model the PVA chains; however, since the lengths of the bonds in their model were fixed, we used bond stiffness values from the OPLS force field 86 to allow the bonds to vibrate. The average bond length using the OPLS bond stiffness values is the same as the original bond length set in Muller’s force field 97 (post-simulation bond length distributions are given in the Supplementary Information). The diversity of force fields 78,88,97 that we have used here reinforces the generality of the results. Polymers were relaxed, and the thermal transport between the chains was measured similar to previous simulations using NEMD. A time step of 0.5 fs was used for all polymer simulations in this section. The cutoff values for the non-bonding interactions were set to be larger than the simulation box in order to account for all of them. After each simulation we analyzed the atomic trajectories to find the degree of inter-chain interaction between the monomers in the polymer chains. We define this degree of interaction as the number of monomers from one chain that are within a specific distance from a monomer in another chain. This distance is approximately equal to 2.5 times the average value of σ parameter in Lennard Jones potential for carbon atoms (∼8Å) in each polymer. This 2.5 σ factor is widely used in molecular simulations as the distance beyond which vdW interactions are negligible. 61 Figure 8 shows the method we have used to determine the number of interactions per monomer in our simulations. We then calculate the average number of interactions per monomer for all monomers in a chain and plot the inter-chain thermal conductance with respect to this average. When calculating the number of interactions for a given monomer, we measure the distance between the center of mass of the monomer and the centers of mass of monomers in the other chain. If the distance is smaller than the specified distance cut-off then we consider the monomers as interacting. To control the level of entanglement of polymer chains, we let them relax for different periods of time ranging from a few femtoseconds to several thousands of femtoseconds, before we fix the terminating hydrogen atoms. The results for inter-chain conductance as a function of average number of interacting monomers are shown in Figure 9. The inter-chain thermal conductance exhibits a linear relation with the average number
17
ACS Paragon Plus Environment
The Journal of Physical Chemistry
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment
Page 18 of 33
Page 19 of 33
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
The Journal of Physical Chemistry
of monomers that are interacting. Each pair of polymer chains has a unique slope that could be used to inform a new theoretical model for the thermal conductivity of bulk samples or to improve existing models. 47,98 The relation observed here is in agreement with a previous study on the effect of chain length on thermal transport in polymers. 99 For different polymers, thermal transport between the chains scales linearly with the average number of interacting monomers and hence scales with inter-chain distance. This is also in agreement with the linear regime observed in Figure 7. Also, the slopes of the lines in Figure 9 are strongly related to the strength and number of non-bonding interactions between the chains. For example, the PMMA monomer has 15 atoms, and hence if two PMMA monomers interact with each other there will be 225 (15 × 15) interactions that can transfer heat. This number for two interacting PE monomers is only 9 (3 × 3) and 49 (7 × 7) for two PVA monomers.
3.5. Effects of Different Inter-chain Interactions on Thermal Transport along Polymer Chains High thermal conductivity along polymer chains is usually associated with robust thermal transport through covalent bonds in the polymer. 17,29,35 In the previous section we showed the significant contribution of non-bonding interactions on inter-chain thermal transport in crosslinked polymers. Single polymer chains could also be considered as several monomers that are sequentially linked together through covalent bonds. Thus, one may expect a significant contribution from non-bonding interactions on thermal transport along the chain. To investigate this possibility, we ran two sets of simulations using the NEMD method detailed in Section 2. In one of them we removed the non-bonding interactions from the PMMA chain and then measured the thermal conductance along the chain. In the other we kept all of the non-bonding interactions included and measured the thermal conductance. These results are reported in Table 2. We found a factor of ten increase in thermal transport along the PMMA polymer chain when non-bonding interactions are included versus when they are not. Similar to the previous section, we observed that non-bonding interactions 19
ACS Paragon Plus Environment
The Journal of Physical Chemistry
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
significantly enhance the intra-chain thermal transport. These results suggest that thermal transport along a single polymer chain could be enhanced by improving the non-bonding interactions along the chain through introduction of polar side groups or hydrogen bonding within the chain. Table 2: Effects of non-bonding interactions on intra-chain thermal conductance in polymers. With non-bonding interactions 0.288 (nW/K) Without non-bonding interactions 0.029 (nW/K)
4. Conclusions Our investigation of heat transfer in polymers through non-bonding and bonding interactions suggests that non-bonding interactions are very important to overall thermal transport. Despite previous beliefs that strong covalent bonds increase heat transfer in polymers, studies have shown both enhancement 18 and reduction 49 in the thermal conductivities of polymers upon crosslinking. We showed that conduction along covalent crosslinkers between polymer chains is not the primary mechanism for heat transfer enhancements in crosslinked polymers. If the crosslinker is short enough to bring the chains close to each other, the large number of non-bonding interactions between the chains is enhanced, leading to a large contribution to heat transfer. This is reasonable since the number of crosslinkers is small compared to the number of non-bonding interactions. Thus, enhancing non-bonding interactions could result in a large improvement in thermal conductivity. For example, recent measurements of polymers with inter-chain hydrogen bonds have shown an order of magnitude enhancement in thermal conductivity. 40 We further illustrated the importance of non-bonding interactions by showing that acoustic waves that carry heat travel significantly faster when non-bonding interactions are included compared to when only covalent bonds are included, because of the long-range nature of non-bonding interactions. The effects of crosslinking species chemical structure on inter20
ACS Paragon Plus Environment
Page 20 of 33
Page 21 of 33
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
The Journal of Physical Chemistry
chain thermal transport were also studied. Our results suggest that in order to synthesize crosslinked polymers with high thermal conductivity, one should choose a short crosslinker and a polymer that is capable of forming strong non-bonding interactions between its chains.
Supporting Information 1) Sample temperature and potential energy profiles 2) Sample force fields used 3) Bond length distribution 4) Inter-chain distance for various crosslinking densities 5) Bulk density and thermal conductivity predictions for force field validation 6) Effects of crosslinker species vibrational density of states of crosslinker agent on interchain thermal transport 7) Method of acoustic wave detection at the probe location
Acknowledgement This work used the Extreme Science and Engineering Discovery Environment (XSEDE), which is supported by National Science Foundation grant number OCI-1053575.
References 1. Klauk, H. Organic Thin-Film Transistors. Chem. Soc. Rev. 2010, 39, 2643–2666. 2. Hu, Z.; Tian, M.; Nysten, B.; Jonas, A. M. Regular Arrays of Highly Ordered Ferroelectric Polymer Nanostructures for Non-Volatile Low-Voltage Memories. Nat. Mater. 2009, 8, 62–67. 21
ACS Paragon Plus Environment
The Journal of Physical Chemistry
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
3. Culebras, M.; Gómez, C. M.; Cantarero, A. Review on Polymers for Thermoelectric Applications. Materials 2014, 7, 6701–6732. 4. Shi, W.; Chen, J.; Xi, J.; Wang, D.; Shuai, Z. Search for Organic Thermoelectric Materials with High Mobility: The Case of 2,7-Dialkyl[1]benzothieno[3,2-b][1]benzothiophene Derivatives. Chemistry of Materials 2014, 26, 2669–2677. 5. Kim, G.-H.; Shao, L.; Zhang, K.; Pipe, K. P. Engineered Doping of Organic Semiconductors for Enhanced Thermoelectric Efficiency. Nat. Mater. 2013, 12, 719–723. 6. See, K. C.; Feser, J. P.; Chen, C. E.; Majumdar, A.; Urban, J. J.; Segalman, R. A. Water-Processable Polymer-Nanocrystal Hybrids for Thermoelectrics. Nano Lett. 2010, 10, 4664–4667, PMID: 20923178. 7. Bubnova, O.; Khan, Z. U.; Malti, A.; Braun, S.; Fahlman, M.; Berggren, M.; Crispin, X. Optimization of the Thermoelectric Figure of Merit in the Conducting Polymer Poly(3,4ethylenedioxythiophene). Nat. Mater. 2011, 10, 429–433. 8. Dayal, S.; Kopidakis, N.; Olson, D. C.; Ginley, D. S.; Rumbles, G. Photovoltaic Devices with a Low Band Gap Polymer and CdSe Nanostructures Exceeding 3% Efficiency. Nano Lett. 2010, 10, 239–242, PMID: 20000623. 9. Pal, S.; Puri, I. K. Thermal Rectification in a Polymer-Functionalized Single-Wall Carbon Nanotube. Nanotechnology 2014, 25, 345401. 10. Tang, C.; Lennon, E. M.; Fredrickson, G. H.; Kramer, E. J.; Hawker, C. J. Evolution of Block Copolymer Lithography to Highly Ordered Square Arrays. Science 2008, 322, 429–432. 11. Ryu, J.-H.; Chacko, R. T.; Jiwpanich, S.; Bickerton, S.; Babu, R. P.; Thayumanavan, S. Self-Cross-Linked Polymer Nanogels: A Versatile Nanoscopic Drug Delivery Platform. J. Am. Chem. Soc. 2010, 132, 17227–17235, PMID: 21077674. 22
ACS Paragon Plus Environment
Page 22 of 33
Page 23 of 33
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
The Journal of Physical Chemistry
12. He, G. S.; Kim, K.-S.; Yuan, L.; Cheng, N.; Prasad, P. N. Two-Photon Pumped Partially Cross-Linked Polymer Laser. Appl. Phys. Lett. 1997, 71, 1619–1621. 13. Seymour, R. B. Plastics. Ind. Eng. Chem. 1966, 58, 61–73. 14. Hashemi, S.; Williams, J. A Fracture Toughness Study on Low Density and Linear Low Density Polyethylenes. Polymer 1986, 27, 384–392. 15. Weathers, A.; Khan, Z. U.; Brooke, R.; Evans, D.; Pettes, M. T.; Andreasen, J. W.; Crispin, X.; Shi, L. Significant Electronic Thermal Transport in the Conducting Polymer Poly(3,4-ethylenedioxythiophene). Adv. Mater. 2015, 27, 2101–2106. 16. Guo, Z.; Lee, D.; Liu, Y.; Sun, F.; Sliwinski, A.; Gao, H.; Burns, P. C.; Huang, L.; Luo, T. Tuning the Thermal Conductivity of Solar Cell Polymers through Side Chain Engineering. Phys. Chem. Chem. Phys. 2014, 16, 7764–7771. 17. Henry, A.; Chen, G. High Thermal Conductivity of Single Polyethylene Chains Using Molecular Dynamics Simulations. Phys. Rev. Lett. 2008, 101, 235502. 18. Kikugawa, G.; Desai, T. G.; Keblinski, P.; Ohara, T. Effect of Crosslink Formation on Heat Conduction in Amorphous Polymers. J. Appl. Phys. 2013, 114, 034302. 19. Zhang, T.; Wu, X.; Luo, T. Polymer Nanofibers with Outstanding Thermal Conductivity and Thermal Stability: Fundamental Linkage between Molecular Characteristics and Macroscopic Thermal Properties. J. Phys. Chem. C 2014, 118, 21148–21159. 20. Yu, C.; Kim, Y. S.; Kim, D.; Grunlan, J. C. Thermoelectric Behavior of SegregatedNetwork Polymer Nanocomposites. Nano Lett. 2008, 8, 4428–4432, PMID: 18959450. 21. Freeman, J. J.; Morgan, G. J.; Cullen, C. A. Thermal Conductivity of a Single Polymer Chain. Phys. Rev. B 1987, 35, 7627–7635. 22. Choy, C. L.; Wong, S. P.; Young, K. Model Calculation of the Thermal Conductivity of Polymer Crystals. J. Polym. Sci. Pol. Phys. 1985, 23, 1495–1504. 23
ACS Paragon Plus Environment
The Journal of Physical Chemistry
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
23. Zhang, T.; Luo, T. Morphology-Influenced Thermal Conductivity of Polyethylene Single Chains and Crystalline Fibers. J. Appl. Phys. 2012, 112, 094304. 24. Zhang, T.; Luo, T. High-Contrast, Reversible Thermal Conductivity Regulation Utilizing the Phase Transition of Polyethylene Nanofibers. ACS Nano 2013, 7, 7592–7600, PMID: 23944835. 25. Rojo, M. M.; Martin, J.; Grauby, S.; Borca-Tasciuc, T.; Dilhaire, S.; Martin-Gonzalez, M. Decrease in Thermal Conductivity in Polymeric P3HT Nanowires by Size-Reduction Induced by Crystal Orientation: New Approaches towards Thermal Transport Engineering of Organic Materials. Nanoscale 2014, 6, 7858–7865. 26. Singh, V.; Bougher, T. L.; Weathers, A.; Cai, Y.; Bi, K.; Pettes, M. T.; McMenamin, S. A.; Lv, W.; Resler, D. P.; Gattuso, T. R. et al. High Thermal Conductivity of Chain-Oriented Amorphous Polythiophene. Nat. Nano. 2014, 9, 384–390. 27. Choy, C. L.; Wong, Y. W.; Yang, G. W.; Kanamoto, T. Elastic Modulus and Thermal Conductivity of Ultradrawn Polyethylene. J. Polym. Sci. Pol. Phys. 1999, 37, 3359– 3367. 28. Wang, Z.; Carter, J. A.; Lagutchev, A.; Koh, Y. K.; Seong, N.-H.; Cahill, D. G.; Dlott, D. D. Ultrafast Flash Thermal Conductance of Molecular Chains. Science 2007, 317, 787–790. 29. Shen, S.; Henry, A.; Tong, J.; Zheng, R.; Chen, G. Polyethylene Nanofibres with Very High Thermal Conductivities. Nat. Nano. 2010, 5, 251–255. 30. Cao, B.-Y.; Li, Y.-W.; Kong, J.; Chen, H.; Xu, Y.; Yung, K.-L.; Cai, A. High Thermal Conductivity of Polyethylene Nanowire Arrays Fabricated by an Improved Nanoporous Template Wetting Technique. Polymer 2011, 52, 1711–1715.
24
ACS Paragon Plus Environment
Page 24 of 33
Page 25 of 33
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
The Journal of Physical Chemistry
31. Fujishiro, H.; Ikebe, M.; Kashima, T.; Yamanaka, A. Drawing Effect on Thermal Properties of High-Strength Polyethylene Fibers. Jpn. J. Appl. Phys. 1998, 37, 1994–1995. 32. Wang, X.; Ho, V.; Segalman, R. A.; Cahill, D. G. Thermal Conductivity of High-Modulus Polymer Fibers. Macromolecules 2013, 46, 4937–4943. 33. Canetta, C.; Guo, S.; Narayanaswamy, A. Measuring Thermal Conductivity of Polystyrene Nanowires Using the Dual-Cantilever Technique. Rev. Sci. Instrum. 2014, 85, 104901. 34. Zhong, Z.; Wingert, M. C.; Strzalka, J.; Wang, H.-H.; Sun, T.; Wang, J.; Chen, R.; Jiang, Z. Structure-Induced Enhancement of Thermal Conductivities in Electrospun Polymer Nanofibers. Nanoscale 2014, 6, 8283–8291. 35. Asegun, H. Thermal Transport in Polymers. Ann. Rev. Heat Transfer 2014, 17, 485–520. 36. Eiermann, K. Modellmäßige Deutung der Wärmeleitfähigkeit von Hochpolymeren. Colloid. Polym. Sci. (Kolloid-Zeitschrift und Zeitschrift für Polymere) 1964, 198, 5–16. 37. Liao, Q.; Liu, Z.; Liu, W.; Deng, C.; Yang, N. Extremely High Thermal Conductivity of Aligned Carbon Nanotube-Polyethylene Composites. Sci. Rep. 2015, 5, 16543. 38. Marconnet, A. M.; Yamamoto, N.; Panzer, M. A.; Wardle, B. L.; Goodson, K. E. Thermal Conduction in Aligned Carbon Nanotube–Polymer Nanocomposites with High Packing Density. ACS Nano 2011, 5, 4818–4825, PMID: 21598962. 39. Xiao, Y.; Wang, W.; Lin, T.; Chen, X.; Zhang, Y.; Yang, J.; Wang, Y.; Zhou, Z. Largely Enhanced Thermal Conductivity and High Dielectric Constant of Poly(Vinylidene Fluoride)/Boron Nitride Composites Achieved by Adding a Few Carbon Nanotubes. J. Phys. Chem. C 2016, 120, 6344–6355. 40. Kim, G.-H.; Lee, D.; Shanker, A.; Shao, L.; Kwon, M. S.; Gidley, D.; Kim, J.; Pipe, K. P.
25
ACS Paragon Plus Environment
The Journal of Physical Chemistry
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
High Thermal Conductivity in Amorphous Polymer Blends by Engineered Interchain Interactions. Nat. Mater. 2015, 14, 295–300. 41. Tonpheng, B.; Yu, J.; Andersson, O. Effects of Cross-Links, Pressure and Temperature on the Thermal Properties and Glass Transition Behaviour of Polybutadiene. Phys. Chem. Chem. Phys. 2011, 13, 15047–15054. 42. Zhang, L.; Chen, T.; Ban, H.; Liu, L. Hydrogen Bonding-Assisted Thermal Conduction in [Small Beta]-Sheet Crystals of Spider Silk Protein. Nanoscale 2014, 6, 7786–7791. 43. Xie, X.; Li, D.; Tsai, T.-H.; Liu, J.; Braun, P. V.; Cahill, D. G. Thermal Conductivity, Heat Capacity, and Elastic Constants of Water-Soluble Polymers and Polymer Blends. Macromolecules 2016, 49, 972–978. 44. Tomlinson, J. N.; Kline, D. E.; Sauer, J. A. Effect of Nuclear Radiation on the Thermal Conductivity of Polyethylene. Polym. Eng. Sci. 1965, 5, 44–48. 45. Yamamoto, O. Thermal Conductivity of Cross-Linked Polymers. Polym. J. 1971, 2, 509–517. 46. Yamamoto, O.; Kambe, H. Thermal Conductivity of Cross-Linked Polymers. A Comparison between Measured and Calculated Thermal Conductivities. Polym. J. 1971, 2, 623–628. 47. Kanari, K.; Ozawa, T. Thermal Conductivity of Epoxy Resins Cured with Aliphatic Amines. Polym. J. 1973, 4, 372–378. 48. Ni, B.; Watanabe, T.; Phillpot, S. R. Thermal Transport in Polyethylene and at Polyethylene–Diamond Interfaces Investigated Using Molecular Dynamics Simulation. J. Phys-Condens. Mat. 2009, 21, 084219. 49. Yu, S.; Park, C.; Hong, S. M.; Koo, C. M. Thermal Conduction Behaviors of Chemically Cross-Linked High-Density Polyethylenes. Thermochim. Acta 2014, 583, 67–71. 26
ACS Paragon Plus Environment
Page 26 of 33
Page 27 of 33
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
The Journal of Physical Chemistry
50. Knappe, W.; Yamamoto, O. Effects of Crosslinking and Chain Degradation on the Thermal Conductivity of Polymers. Colloid. Polym. Sci. (Kolloid-Zeitschrift und Zeitschrift für Polymere) 1970, 240, 775–783. 51. Perkin, K. K.; Turner, J. L.; Wooley, K. L.; Mann, S. Fabrication of Hybrid Nanocapsules by Calcium Phosphate Mineralization of Shell Cross-Linked Polymer Micelles and Nanocages. Nano Lett. 2005, 5, 1457–1461, PMID: 16178257. 52. Srinivasan, K.; Pohl, C.; Avdalovic, N. Cross-Linked Polymer Coatings for Capillary Electrophoresis and Application to Analysis of Basic Proteins, Acidic Proteins, and Inorganic Ions. Anal. Chem. 1997, 69, 2798–2805. 53. Ohara, T. J.; Rajagopalan, R.; Heller, A. Glucose Electrodes Based on CrossLinked Bis(2,2’-Bipyridine)Chloroosmium(+/2+) Complexed Poly(1-Vinylimidazole) Films. Anal. Chem. 1993, 65, 3512–3517. 54. Gosselin, M. A.; Guo, W.; Lee, R. J. Efficient Gene Transfer Using Reversibly CrossLinked Low Molecular Weight Polyethylenimine. Bioconjugate Chem. 2001, 12, 989–994, PMID: 11716690. 55. Peterson, A. M.; Jensen, R. E.; Palmese, G. R. Reversibly Cross-Linked Polymer Gels as Healing Agents for Epoxy-Amine Thermosets. ACS Appl. Mater. Inter. 2009, 1, 992–995, PMID: 20355883. 56. Scott, T. F.; Schneider, A. D.; Cook, W. D.; Bowman, C. N. Photoinduced Plasticity in Cross-Linked Polymers. Science 2005, 308, 1615–1617. 57. Trollsås, M.; Orrenius, C.; Sahlén, F.; Gedde, U. W.; Norin, T.; Hult, A.; Hermann, D.; Rudquist, P.; Komitov, L.; Lagerwall, S. T. et al. Preparation of a Novel Cross-Linked Polymer for Second-Order Nonlinear Optics. J. Am. Chem. Soc. 1996, 118, 8542–8548.
27
ACS Paragon Plus Environment
The Journal of Physical Chemistry
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
58. Gregg, B. A.; Heller, A. Cross-Linked Redox Gels Containing Glucose Oxidase for Amperometric Biosensor Applications. Anal. Chem. 1990, 62, 258–263. 59. Roberts, M. E.; Queraltó, N.; Mannsfeld, S. C. B.; Reinecke, B. N.; Knoll, W.; Bao, Z. Cross-Linked Polymer Gate Dielectric Films for Low-Voltage Organic Transistors. Chem. Mater. 2009, 21, 2292–2299. 60. Ge, H.; Wu, W.; Li, Z.; Jung, G.-Y.; Olynick, D.; Chen, Y.; Liddle, J. A.; Wang, S.Y.; Williams, R. S. Cross-Linked Polymer Replica of a Nanoimprint Mold at 30 nm Half-pitch. Nano Lett. 2005, 5, 179–182, PMID: 15792435. 61. Toxvaerd, S.; Dyre, J. C. Communication: Shifted Forces in Molecular Dynamics. J. Chem. Phys. 2011, 134 . 62. Liu, J.; Wang, X.; Li, D.; Coates, N. E.; Segalman, R. A.; Cahill, D. G. Thermal Conductivity and Elastic Constants of PEDOT:PSS with High Electrical Conductivity. Macromolecules 2015, 48, 585–591. 63. Rashidi, V.; Pipe, K. P. Contributions of Strain Relaxation and Interface Modes to Thermal Transport in Superlattices. Comp. Mater. Sci. 2015, 107, 151–156. 64. Schelling, P. K.; Phillpot, S. R.; Keblinski, P. Comparison of Atomic-Level Simulation Methods for Computing Thermal Conductivity. Phys. Rev. B 2002, 65, 144306. 65. Luo, T.; Lloyd, J. R. Non-Equilibrium Molecular Dynamics Study of Thermal Energy Transport in Au-SAM-Au Junctions. Int. J. Heat Mass Transfer 2010, 53, 1–11. 66. Zhang, H.; Fonseca, A. F.; Cho, K. Tailoring Thermal Transport Property of Graphene through Oxygen Functionalization. J. Phys. Chem. C 2014, 118, 1436–1442. 67. Liu, J.; Yang, R. Length-Dependent Thermal Conductivity of Single Extended Polymer Chains. Phys. Rev. B 2012, 86, 104307.
28
ACS Paragon Plus Environment
Page 28 of 33
Page 29 of 33
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
The Journal of Physical Chemistry
68. Varshney, V.; Patnaik, S. S.; Roy, A. K.; Froudakis, G.; Farmer, B. L. Modeling of Thermal Transport in Pillared-Graphene Architectures. ACS Nano 2010, 4, 1153–1161, PMID: 20112924. 69. Bagri, A.; Kim, S.-P.; Ruoff, R. S.; Shenoy, V. B. Thermal Transport across Twin Grain Boundaries in Polycrystalline Graphene from Nonequilibrium Molecular Dynamics Simulations. Nano Lett. 2011, 11, 3917–3921, PMID: 21863804. 70. Chen, X.; Munjiza, A.; Zhang, K.; Wen, D. Molecular Dynamics Simulation of Heat Transfer from a Gold Nanoparticle to a Water Pool. J. Phys. Chem. C 2014, 118, 1285–1293. 71. Varshney, V.; Patnaik, S. S.; Roy, A. K.; Farmer, B. L. Modeling of Thermal Conductance at Transverse CNT-CNT Interfaces. J. Phys. Chem. C 2010, 114, 16223–16228. 72. Ren, C.; Zhang, W.; Xu, Z.; Zhu, Z.; Huai, P. Thermal Conductivity of Single-Walled Carbon Nanotubes under Axial Stress. J. Phys. Chem. C 2010, 114, 5786–5791. 73. Pal, S.; Balasubramanian, G.; Puri, I. K. Modifying Thermal Transport in Electrically Conducting Polymers: Effects of Stretching and Combining Polymer Chains. J. Chem. Phys. 2012, 136, 044901. 74. Zhao, J.; Jiang, J.-W.; Wei, N.; Zhang, Y.; Rabczuk, T. Thermal Conductivity Dependence on Chain Length in Amorphous Polymers. J. Appl. Phys. 2013, 113, 184304. 75. Plimpton, S. Fast Parallel Algorithms for Short-Range Molecular Dynamics. J. Comput. Phys. 1995, 117, 1–19. 76. Hanwell, M.; Curtis, D.; Lonie, D.; Vandermeersch, T.; Zurek, E.; Hutchison, G. Avogadro: an Advanced Semantic Chemical Editor, Visualization, and Analysis Platform. J. Cheminform. 2012, 4 .
29
ACS Paragon Plus Environment
The Journal of Physical Chemistry
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
77. Humphrey, W.; Dalke, A.; Schulten, K. VMD: Visual Molecular Dynamics. J. Mol. Graphycs 1996, 14, 33–38. 78. Haley, B. P.; Wilson, N.; Li, C.; Arguelles, A.; Jaramillo, E.; Strachan, A. Polymer Modeler. https://nanohub.org/resources/polymod 2017, (DOI: 10.4231/D3JD4PQ7T). 79. Persistence of Vision Pty. Ltd. (2004). Persistence of Vision (TM) Raytracer. Persistence of Vision Pty. Ltd., Williamstown, Victoria, Australia. http://www.povray.org/ 80. Waters, M. J.; Sebeck, K.; Coyle, E. J.; Bregman, A. Kirke. Kieffer Group, Ann Arbor, MI. 2012. 81. Schneider, T.; Stoll, E. Molecular-Dynamics Study of a Three-Dimensional OneComponent Model for Distortive Phase Transitions. Phys. Rev. B 1978, 17, 1302–1322. 82. Nosé, S. A Molecular Dynamics Method for Simulations in the Canonical Ensemble. Mol. Phys. 1984, 52, 255–268. 83. Hoover, W. G. Canonical Dynamics: Equilibrium Phase-Space Distributions. Phys. Rev. A 1985, 31, 1695–1697. 84. Xu, Z.; Buehler, M. J. Nanoengineering Heat Transfer Performance at Carbon Nanotube Interfaces. ACS Nano 2009, 3, 2767–2775, PMID: 19702296. 85. Yang, N.; Zhang, G.; ; Li, B. Ultralow Thermal Conductivity of Isotope-Doped Silicon Nanowires. Nano Lett. 2008, 8, 276–280, PMID: 18095735. 86. Jorgensen, W. L.; Maxwell, D. S.; ; Tirado-Rives, J. Development and Testing of the OPLS All-Atom Force Field on Conformational Energetics and Properties of Organic Liquids. J. Am. Chem. Soc. 1996, 118, 11225–11236. 87. Kahn, K.; Bruice, T. C. Parameterization of OPLS–AA Force Field for the Conformational Analysis of Macrocyclic Polyketides. J. Comput. Chem. 2002, 23, 977–996. 30
ACS Paragon Plus Environment
Page 30 of 33
Page 31 of 33
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
The Journal of Physical Chemistry
88. Chen, C.; Maranas, J. K.; García-Sakai, V. Local Dynamics of Syndiotactic Poly(Methyl Methacrylate) Using Molecular Dynamics Simulation. Macromolecules 2006, 39, 9630– 9640. 89. Zhang, T.; Luo, T. Role of Chain Morphology and Stiffness in Thermal Conductivity of Amorphous Polymers. J. Phys. Chem. B 2016, 120, 803–812, PMID: 26751002. 90. Kikugawa, G.; Ohara, T.; Kawaguchi, T.; Torigoe, E.; Hagiwara, Y.; Matsumoto, Y. A Molecular Dynamics Study on Heat Transfer Characteristics at the Interfaces of Alkanethiolate Self-Assembled Monolayer and Organic Solvent. J. Chem. Phys. 2009, 130, 074706. 91. Chiavazzo, E.; Asinari, P. Enhancing Surface Heat Transfer by Carbon Nanofins: Towards an Alternative to Nanofluids? Nanoscale Res. Lett. 2011, 6, 1–13. 92. Konatham, D.; Papavassiliou, D.; Striolo, A. Thermal Boundary Resistance at the Graphene-Graphene Interface Estimated by Molecular Dynamics Simulations. Chem. Phys. Lett. 2012, 527, 47–50. 93. Acharya, H.; Mozdzierz, N. J.; Keblinski, P.; Garde, S. How Chemistry, Nanoscale Roughness, and the Direction of Heat Flow Affect Thermal Conductance of Solid–Water Interfaces. Ind. Eng Chem. Res. 2012, 51, 1767–1773. 94. Gao, W.; Kong, L.; Hodgson, P. Atomic Interaction of Functionalized Carbon NanotubeBased Nanofluids with a Heating Surface and its Effect on Heat Transfer. Int. J. Heat Mass Tran. 2012, 55, 5007–5015. 95. Chen, X.; Dam, M. A.; Ono, K.; Mal, A.; Shen, H.; Nutt, S. R.; Sheran, K.; Wudl, F. A Thermally Re-Mendable Cross-Linked Polymeric Material. Science 2002, 295, 1698– 1702.
31
ACS Paragon Plus Environment
The Journal of Physical Chemistry
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
96. Robbins, A. B.; Minnich, A. J. Crystalline Polymers with Exceptionally Low Thermal Conductivity Studied Using Molecular Dynamics. Appl. Phys. Lett. 2015, 107, 201908. 97. Muller-Plathe, F.; van Gunsteren, W. F. Solvation of Poly(Vinyl Alcohol) in Water, Ethanol and an Equimolar Water-Ethanol Mixture: Structure and Dynamics Studied by Molecular Dynamics Simulation. Polymer 1997, 38, 2259–2268. 98. Zhong, C.; Yang, Q.; Wang, W. Correlation and Prediction of the Thermal Conductivity of Amorphous Polymers. Fluid Phase Equilib. 2001, 181, 195–202. 99. Lohe, P. Wärmeleitfähigkeit Hochpolymerer Schmelzen. Colloid. Polym. Sci. (KolloidZeitschrift und Zeitschrift für Polymere) 1965, 204, 7–10.
32
ACS Paragon Plus Environment
Page 32 of 33
Page 33 of 33
The Journal of Physical Chemistry
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment