Conformational behavior of polymers adsorbed on nanotubes

The importance of hydrophobic interactions in determining polymer adsorption and wrapping of carbon nanotubes is still under debate. In this work, we concentrate on the effect of short-ranged weakly attractive hydrophobic interactions between polymers and nanotubes (modeled as an infinitely long and smooth cylindrical surface), neglecting all other interactions apart for chain flexibility. Using coarse-grained Monte Carlo simulation of such simplified systems, we find that uniform adsorption and wrapping of the nanotube occur for all degrees of chain flexibility for tubes with sufficiently large outer radii. However, the adsorbed conformations depend on chain stiffness, ranging from randomly adsorbed conformations of the flexible chain to perfect helical or multihelical conformations (in the case of more concentrated solutions) of the rigid chains. Adsorption appears to occur in a sequential manner, wrapping the nanotube nearly one monomer at a time from the point of contact. Once adsorbed, the chains travel on the surface of the cylinder, retaining their helical conformations for the semiflexible and rigid chains. Our findings may provide additional insight to experimentally observed ordered polymer wrapping of carbon nanotubes.     

Polymer Nanocomposites Using Urchin‐Shaped Carbon Nanotube‐Silica Hybrids as Reinforcing Fillers

Summary: Carbon nanotubes (CNTs) have been grown on MCM‐41 supported Fe nanoparticles and the as‐prepared (no further purification) CNT‐silica hybrid was directly incorporated into nylon‐6 (PA6) by simple melt‐compounding. The urchin‐shaped CNT‐silica hybrid filler was observed to be homogeneously dispersed throughout the matrix by scanning electron and transmission electron microscopy. Compared with neat PA6, the tensile modulus and strength of the composite are greatly improved by about 110%, with incorporation of only 1 wt.‐% CNT‐silica filler.

SEM image and schematic representation showing polymer chains wrapping around the urchin‐shaped CNT‐silica hybrid filler.     

Morphology,rheological and crystallization behavior in non-covalently functionalized carbon nanotube reinforced poly(butylene succinate) nanocomposites with low percolation threshold

Multi-walled carbon nanotubes (CNTs) were non-covalently functionalized by surface wrapping of poly(sodium 4-styrenesulfonate) (PSS) with the aid of ultrasound. The functionalized CNTs were incorporated into poly(butylene succinate) (PBS) through solution coagulation to fabricate CNTs filled PBS nanocomposites. The morphologies of the PBS/CNT nanocomposites were studied by scanning electron microscope (SEM) and transmission electron microscope (TEM), and the effect of loading of functionalized CNT on the rheological behavior, electrical conductivity and mechanical properties of the nanocomposites was investigated systemically. SEM observation indicates that functionalized CNTs dispersed in PBS matrix without obvious aggregation and showed good interfacial adhesion with the PBS phase. TEM observation reveals that a CNT network was formed when the loading of CNTs increased from 0.1 to 0.3 wt%. Rheological investigation indicates the formation of a CNT network with a percolation threshold of only 0.3 wt%. Significant improvement in electrical conductivity occurred at CNT loading of 0.3 wt%, with the value of electrical conductivity increasing by six orders of magnitude compared to neat PBS. Differential scanning calorimetry indicates that the melt crystallization temperature of PBS was improved by ∼14 °C with addition of only 0.05 wt% functionalized CNTs. Tensile tests indicate that both the yield strength and Young's modulus of PBS were apparently reinforced by incorporation of functionalized CNTs, while the elongation at break was reduced gradually.     

Binding Mode Multiplicity and Multiscale Chirality in the Supramolecular Assembly of DNA and a π-Conjugated Polymer

Water-soluble π-conjugated polymers are increasingly considered for DNA biosensing. However, the conformational rearrangement, supramolecular organization and dynamics upon interaction with DNA have been overlooked, which prevents the rational design of such detection tools. To elucidate the binding of a cationic polythiophene (CPT) to DNA with atomistic resolution, we performed molecular simulations of their supramolecular assembly. Comparison of replicated simulations show a multiplicity of CPT binding geometries that contribute to the wrapping of CPT around DNA. The different binding geometries are stabilized by both electrostatic interactions between CPT lateral cations and DNA phosphodiesters and van der Waals interactions between the CPT backbone and the DNA grooves. Simulated circular dichroism (CD) spectra show that the induced CD signal stems from a conserved geometrical feature across the replicated simulations, i. e. the presence of segments of syn configurations between thiophene units along the CPT chain. At the macromolecular scale, we inspected the different shapes related to the CPT binding modes around the DNA through symmetry metrics. Altogether, molecular dynamics (MD) simulations, model Hamiltonian calculations of the CD spectra, and symmetry indices provide insights into the origin of induced chirality from the atomic to the macromolecular scale. Our multidisciplinary approach points out the hierarchical aspect of CPT chiral organization induced by DNA.     

Effects of size constraint on water filling process in nanotube

Molecular dynamics (MD) simulation and the potential of mean force (PMF) analysis are used to investigate the structural properties of water molecules near the end of nanotube for the whole process from the initial water filling up to the configuration stabilization inside the carbon nanotubes (CNTs). Numerical simulations showed that when a small-sized nanotube is immersed into the water bath, the size constraint will induce a prevailing orientation for the water molecule to diffuse into the tube and this effect can persist approximately 3.3 angstroms from the end of CNT. As the structure within the CNTs stabilizes, the ambient structural properties can indirectly reflect their corresponding properties inside the nanotube. Our results also showed that there exists a close correlation between the PMF analysis and the results of MD simulations, and the properties at nanometer scale are closely related to the size-constraint effect.     

Miscibility behavior and single chain properties in polymer blends: a bond fluctuation model study

Computer simulation studies on the miscibility behavior and single chain properties in binary polymer blends are reviewed. We consider blends of various architectures in order to identify important architectural parameters on a coarse grained level and study their qualitative consequences for the miscibility behavior. The phase diagram, the relation between the exchange chemical potential and the composition, and the intermolecular pair correlation functions for symmetric blends of linear chains, blends of cyclic polymers, blends with an asymmetry in cohesive energies, blends with different chain lengths, blends with distinct monomer shapes, and blends with a stiffness disparity between the components are discussed. For strictly symmetric blends the Flory‐Huggins theory becomes quantitatively correct in the long chain length limit, when the χ parameter is identified via the intermolecular pair correlation function. For small chain lengths composition fluctuations are important. They manifest themselves in 3D Ising behavior at the critical point and an upward parabolic curvature of the χ parameter from small‐angle neutron scattering close to the critical point. The ratio between the mean field estimate and the true critical temperature decreases like √χ/(ρb³) for long chain lengths. The chain conformations in the minority phase of a symmetric blend shrink as to reduce the number of energeticaly unfavorable interactions. Scaling arguments, detailed self‐consistent field calculations and Monte Carlo simulations of chains with up to 512 effective segments agree that the conformational changes decrease around the critical point like 1/√N. Other mechanisms for a composition dependence of the single chain conformations in asymmetric blends are discussed. If the constituents of the blends have non‐additive monomer shapes, one has a large positive chain‐length‐independent entropic contribution to the χ parameter. In this case the blend phase separates upon heating at a lower critical solution temperature. Upon increasing the chain length the critical temperature approaches a finite value from above. For blends with a stiffness disparity an entropic contribution of the χ parameter of the order 10^–3 is measured with high accuracy. Also the enthalpic contribution increases, because a back folding of the stiffer component is suppressed and the stiffer chains possess more intermolecular contacts. Two aspects of the single chain dynamics in blends are discussed: (a) The dynamics of short non‐entangled chains in a binary blend are studied via dynamic Monte Carlo simulations. There is hardly any coupling between the chain dynamics and the thermodynamic state of the mixture. Above the critical temperatures both the translational diffusion and the relaxation of the chain conformations are independent of the temperature. (b) Irreversible reactions of a small fraction of reactive polymers at a strongly segregated interface in a symmetric binary polymer blend are investigated. End‐functionalized homopolymers of different species react at the interface instantaneously and irreversibly to form diblock copolymers. The initial reaction rate for small reactant concentrations is time dependent and larger than expected from theory. At later times there is a depletion of the reactive chains at the interface and the reaction is determined by the flux of the chains to the interface. Pertinent off‐lattice simulations and analytical theories are briefly discussed.     

Structural and Dynamic Properties of Monoclonal Antibodies Immobilized on CNTs: A Computational Study

Due to the widespread application of carbon nanotube (CNT)‐based materials in nanomedicine, it is nowadays of paramount importance to unravel at the atomistic level of detail the structural properties of such bioconjugates in order to rationalize and predict the effect exerted by the graphitic framework on the bio‐active counterpart. In this paper, we report for the first time all‐atom explicit solvent molecular dynamics (MD) simulations investigating the structural and dynamic properties of a noncovalent bioconjugate in which the monoclonal Cetuximab antibody (Ctx) is adsorbed on a CNT surface. Upon selection of the three most representative adsorption modes as obtained by docking studies, force‐field MD and DFT simulations unambiguously showed that hydrophobic interactions mainly govern the adsorption of the protein on the graphitic surface. Two main adsorption poses have been predicted: a pose‐fab (p‐fab) and pose‐fc (p‐fc) (fab = fragment antigen binding region; fc = fragment crystallizable region), the former being favored with small‐diameter tubes (≤40 Å). In all the predicted poses, the secondary structure of Ctx is largely unaffected by the presence of the graphitic surface and, consistently with previous literature studies, our simulations reveal that positively charged amino acidic residues, such as Lys and Arg, predominantly contribute to the stabilization of the CNT?Ctx complex acting like surfactants. The predicted structural models are consistent with the experimental data, for which the immobilization of the antibody on CNTs does not disrupt the structural and recognition properties of the Ctx, consequently supporting the reliability of the used bioconjugation strategy for engineering stable and responsive hybrid nanomaterials for therapeutic applications. Moreover, a remarkable structural similarity of Ctx with antibodies of different isotypes suggests that in principle the CNT framework can interact in the same manner with all antibodies currently used in clinical applications.     

Solvent effects on electronic structures and chain conformations of alpha-oligothiophenes in polar and apolar solutions

Solvent effects on electronic structures and chain conformations of alpha-oligothiophenes nTs (n = 1 to 10) are investigated in solvents of n-hexane, 1,4-dioxane, carbon tetrachloride, chloroform, and water by using density functional theory (DFT) and molecular dynamics (MD) simulations. Both implicit and explicit solvent models are employed. The polarized continuum model (PCM) calculations and MD simulations demonstrate the weak solvent effects on the electronic structures of alpha-oligothiophenes. The lowest dipole-allowed vertical excitation energies of nTs, obtained from time-dependent DFT/PCM calculations at the B3LYP/6-31G(d) level, exhibit a red shift as the solvent polarity increases, in agreement with experiments. The studied solvents have little impact on the state order of the low-lying excited states provided that the nTs are kept in C2h or C2v symmetry. The MD simulations demonstrate that the chain conformations are distorted to some extent in polar and nonpolar solvents. A qualitative picture of the distribution of solvent molecules around the solvated nTs is drawn by means of radial and spatial distribution functions. The S...H-O and pi...H-O solute-solvent interactions are insignificant in aqueous solution.     

疏水蛋白吸附碳纳米管的分子动力学模拟

采用分子动力学方法模拟了水溶液中Ⅱ型疏水蛋白HFBI在单壁碳纳米管(SWNTs)表面的吸附过程, 考察了3种不同的HFBI初始取向, 并计算了结合自由能, 从累计240 ns的模拟轨迹中得到了不同的吸附结构. 结果表明, 当HFBI完全通过疏水面与SWNTs作用时, 其结合自由能最有利吸附, 且吸附最稳定. 另外, 由于HFBI含有4个二硫键, 因此吸附过程几乎不改变其二级结构.     

Using ONIOM calculations to investigate the abilities of simple and nitrogen,boron, sulfur-doped carbon nanotubes in sensing of carbon monoxide

In this work, geometries, stabilities, and electronic properties of the carbon monoxide (CO) molecule as an adsorbent in a simple carbon nanotube (CNT) and nitrogen (N), boron (B), sulfur (S)-doped CNTs (NCNT, BCNT, and SCNT) with parallel and perpendicular configurations are fully considered using ONIOM, natural bond orbital, and quantum theory of atom in molecule (QTAIM) calculations. The adsorption energies (E_ad) demonstrate that a CO molecule could be adsorbed on the surface of the simple CNT with parallel configuration and N-doped CNT with perpendicular configuration in an exothermic process. QTAIM calculations showed the close-shell (noncovalent) interactions between the CO molecule and CNT or N, B, S-doped CNTs. In addition, the energy gap (E_g) values between the highest occupied molecular orbital and the lowest unoccupied molecular orbital are calculated. In accordance with the results of energy gap, simple and N-doped CNTs could be used as CO sensors.     

Collapse transitions of a supersized neutral chain due to irreversibly adsorbed small colloidal particles

 We performed Monte Carlo simulations to study the destabilization processes of large neutral and flexible polymer chains due to irreversibly adsorbed colloidal particles attached to the chains like beads on a necklace. The particles are modeled as charged spherical units which interact with each other via repulsive electrostatic and attractive van der Waals (vdW) potentials. The usual Monte Carlo search procedure is extended and carefully checked to completely sample the chain conformational space and achieve dense conformations in the limit of both strong attractive and repulsive interaction potentials. Configurational properties, such as the radius of gyration, the end-to-end length, and the Kuhn length, are calculated as a function of the intensity of the vdW interactions and ionic strength values. It is observed that chains exhibit a new range of possible conformations compared to the classical random walk and self avoiding walk chains or polyelectrolytes. In the limit of low salt concentration, by gradually increasing vdW interactions, chains undergo a cascade of transitions from extended structures to dumbbells, from dumbbells to pearl necklaces, and from pearl necklaces to collapsed coils. Because of strong competition between the vdW and electrostatic forces, the distance along the chain between the interacting particles, and the sampling limitations, these transitions are found to sample metastable domains and to depend on the initial conformations. To gain insight into the spatial organization of the collapsed conformations, the pair correlation functions of both monomers and particles are calculated. It is shown that collapsed conformations which are the result of strong particle–particle interactions exhibit two distinct parts: a hard core mainly composed of particles and a surrounding polymeric shell composed of loops and tails. Possible effects of such a collapsed transition on the kinetics of flocculation of a mixture containing large flexible chains and small adsorbing colloidal particles are discussed. Received: 26 July 1999 Accepted in revised form: 9 November 1999     

Non-isothermal crystallization of ethylene-vinyl acetate copolymer containing a high weight fraction of graphene nanosheets and carbon nanotubes

The effect of the different geometrical dimensionality of two dimensional graphene nanosheets(2D GNSs) and one dimensional carbon nanotubes(1D CNTs) on the non-isothermal crystallization of an ethylene-vinyl acetate(EVA) copolymer at high loading(5 wt%) was studied.Transmission electron microscopy indicated a homogeneous dispersion of GNSs and CNTs in EVA obtained by a solution dispersion process.Fourier-transform infrared spectroscopy and differential scanning calorimetry measurements showed that 1D CNTs and 2D GNSs acted as effective nucleating agents,with a noticeably increased onset crystallization temperature of EVA.A high weight fraction of nano-fillers slowed the overall crystallization rate of composites.At the same crystallization temperature,the crystallization behavior of GNS/EVA composites was slowed compared to that of the CNT/EVA ones owing to larger nucleus barrier and activation energy of diffusion.Dynamic mechanical relaxation and rheology behavior of CNT/EVA and GNS/EVA composites demonstrated that the planar structure of the GNSs had an intensively negative effect on EVA chain mobility due to interactions between nanofillers and polymer chains,as well as spatial restriction.     

Single-chain compaction of long duplex DNA by cationic nanoparticles: modes of interaction and comparison with chromatin

The compaction of long duplex DNA by cationic nanoparticles (NP) used as a primary model of histone core particles has been investigated. We have systematically studied the effect of salt concentration, particle size, and particle charge by means of single-molecule observations-fluorescence microscopy (FM) and transmission electron microscopy (TEM)-and molecular dynamics (MD) simulations. We have found that the large-scale DNA compaction is progressive and proceeds through the formation of beads-on-a-string structures of various morphologies. The DNA adsorbed amount per particle depends weakly on NP concentration but increases significantly with an increase in particle size and is optimal at an intermediate salt concentration. Three different complexation mechanisms have been identified depending on the correlation between DNA and NPs in terms of geometry, chain rigidity, and electrostatic interactions: free DNA adsorption onto NP surface, DNA wrapping around NP, and NP collection on DNA chain.     

Polymer Structures and Glass Transition: A Molecular Dynamics Simulation Study

Full atomistic molecular dynamics (MD) simulations on five polymers with different chain backbone (C—C, Si—O, and C—O) and different side groups (—H, one —CH₃, and two —CH₃) are performed to study the effects of chain flexibility and side groups on the glass transition of polymers. Molecular dynamics simulations of NPT (constant pressure and constant temperature) dynamics are carried out to obtain specific volume as a function of temperature for polyethylene (PE), poly(propylene) (PP), polyisobutylene (PIB), poly(oxymethylene) (POM), and poly(dimethylsiloxane) (PDMS). The volumetric glass transition temperature has been determined as the temperature marking the discontinuity in slope of the plots of V–T simulation data. Various energy components at different temperatures of the polymers are investigated and their roles played in the glass transition process are analyzed. In order to understand the polymer chain conformations above and below the glass transition temperature, dihedral angle distributions of polymer chains at various temperatures are also studied.     

Enhanced sampling method in molecular simulations using genetic algorithm for biomolecular systems

We propose a molecular simulation method using genetic algorithm (GA) for biomolecular systems to obtain ensemble averages efficiently. In this method, we incorporate the genetic crossover, which is one of the operations of GA, to any simulation method such as conventional molecular dynamics (MD), Monte Carlo, and other simulation methods. The genetic crossover proposes candidate conformations by exchanging parts of conformations of a target molecule between a pair of conformations during the simulation. If the candidate conformations are accepted, the simulation resumes from the accepted ones. While conventional simulations are based on local update of conformations, the genetic crossover introduces global update of conformations. As an example of the present approach, we incorporated genetic crossover to MD simulations. We tested the validity of the method by calculating ensemble averages and the sampling efficiency by using two kinds of peptides, ALA3 and (AAQAA)₃. The results show that for ALA3 system, the distribution probabilities of backbone dihedral angles are in good agreement with those of the conventional MD and replica-exchange MD simulations. In the case of (AAQAA)₃ system, our method showed lower structural correlation of α-helix structures than the other two methods and more flexibility in the backbone ψ angles than the conventional MD simulation. These results suggest that our method gives more efficient conformational sampling than conventional simulation methods based on local update of conformations. © 2018 Wiley Periodicals, Inc.     

Effect of Size and Temperature on Water Dynamics inside Carbon Nano-Tubes Studied by Molecular Dynamics Simulation

Water transport inside carbon nano-tubes (CNTs) has attracted considerable attention due to its nano-fluidic properties, its importance in nonporous systems, and the wide range of applications in membrane desalination and biological medicine. Recent studies show an enhancement of water diffusion inside nano-channels depending on the size of the nano-confinement. However, the underlying mechanism of this enhancement is not well understood yet. In this study, we performed Molecular Dynamics (MD) simulations to study water flow inside CNT systems. The length of CNTs considered in this study is 20 nm, but their diameters vary from 1 to 10 nm. The simulations are conducted at temperatures ranging from 260 K to 320 K. We observe that water molecules are arranged into coaxial water tubular sheets. The number of these tubular sheets depends on the CNT size. Further analysis reveals that the diffusion of water molecules along the CNT axis deviates from the Arrhenius temperature dependence. The non-Arrhenius relationship results from a fragile liquid-like water component persisting at low temperatures with fragility higher than that of the bulk water.     

Physical association of polymers with nanotubes

We use a coarse-grained Monte Carlo model to further investigate the association of polymers with carbon nanotubes (CNTs). Previous studies have shown ordered helical wrapping conformations for a range of investigated parameters. Such adsorbed conformations allow the polymers to spiral up and down the surface of the nanotube, retaining their helical state. We analyze the helical pitch of such conformations, and relate it to nanotube radius and chain stiffness using a simple model. The model reveals that the helical pitch is approximately determined by the matching between the radius of curvature of the helix with the average bending angle of the polymer, determined by its persistence length. In addition, we simulate adsorption of block and triblock copolymers (BCPs) whose different blocks are differentiated by their degree of association with the nanotube (hydrophobic or polar). The hydrophobic blocks of the copolymers initially adsorb in both helical and random conformations of the hydrophobic block, depending on which part of the chain (center or ends) adsorbs first on the CNTs surface. In both configurations, however, the polar block extends away from the nanotube, forming loops and tails for triblock and diBCPs, respectively. Such configurations may improve the interfacial adhesion in polymer–CNTs composites. © 2008 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 46: 2711–2718, 2008     

Effect of dispersion and selective localization of carbon nanotubes on rheology and electrical conductivity of polyamide 6 (PA6), Polypropylene (PP), and PA6/PP nanocomposites

Nanocomposites were prepared by adding 1–3 vol % multiwalled carbon nanotubes (MWCNTs) to polyamide 6 (PA6), polypropylene (PP), and their co‐continuous blends of 60/40 and 50/50 volume compositions. Because of the good interaction and interfacial adhesion to the PA6, nanotubes were disentangled and distributed evenly through nanocomposites containing PA6. In contrast, lack of active interactions between the matrix and the CNTs resulted in poor tube dispersion in PP. These observations were then verified by studying the rheology and electrical conductivity of their respective nanocomposites. Absence of percolated CNT clusters and possible wrapping of the tubes by PA6 resulted in low electrical conductivity of PA6/CNT nanocomposites. On the other hand, despite the weak dispersion of the tubes, electrical conductivities of PP/CNT nanocomposites were much higher than all other counterparts. This could be the result of good three‐dimensional distribution of the agglomerated bundles and secondary aggregation of tubes in PP. Adding CNTs to blends of PA6/PP (60/40 and 50/50) resulted in almost full localization of carbon nanotubes in PA6, leading to their higher effective concentration. At the same CNT loadings, the blend nanocomposites had three to seven orders of magnitude higher electrical conductivity than pure PA6. © 2014 Wiley Periodicals, Inc. J. Polym. Sci., Part B: Polym. Phys. 2015 , 53, 368–378     

Phase separation in binary mixtures containing polymers: A quantitative comparison of single-chain-in-mean-field simulations and computer simulations of the corresponding multichain systems

We discuss a phenomenological, coarse-grained simulation scheme, single-chain-in-mean-field (SCMF) simulation, for investigating the kinetics of phase separation in dense polymer blends and mixtures of polymers and solvents. In the spirit of self-consistent-field calculations, we approximate the interacting multichain problem by that of a single chain in an external field, which, in turn, depends on the local densities of the components. To study the time evolution of the mixture, we perform an explicit Monte Carlo (MC) simulation of an ensemble of independent chains in the external field and periodically calculate the average densities and update the external field. Unlike dynamic self-consistent-field theory, these SCMF simulations do not assume that the chain conformations relax much more quickly than the density and incorporate the single-chain dynamics explicitly rather than via an Onsager coefficient. This allows us to study systems with large spatial inhomogeneities and dynamic asymmetries. To assess the accuracy and limitations of the simulation scheme, we compare the results of SCMF simulations using a discretized Edwards Hamiltonian with computer simulations of the corresponding multichain system for (1) the early stages of spinodal decomposition of a symmetric binary polymer blend in response to a quench from χN = 0.314 to χN = 5 (where χ is the Flory–Huggins parameter and N is the number of segments), for which the growth rate of composition fluctuations is compared with MC simulations of the bond fluctuation model and alternative dynamic self-consistent-field calculations, and (2) the evaporation of a solvent from a low-molecular-weight thin polymer film, for which a comparison is made with molecular dynamics (MD) simulations of a bead-necklace model with a monomeric solvent. In the latter case, the polymer conformations are extracted from MD simulations and modeled in the SCMF simulations by a discretized Edwards Hamiltonian augmented by a chain-bending potential. From the MD simulations of thin polymer films in equilibrium with its vapor, phase coexistence has been determined, and the second- and third-order virial coefficients in the SCMF simulations have been adjusted accordingly. Finally, MD simulations of bulk solutions of a polymer and a solvent over a range of compositions, as well as the pure solvent at various densities, have been performed to determine self-diffusion coefficients that enter the SCMF simulations in the form of density-dependent segmental mobilities. A comparison of the polymer and solvent profiles in a thin film as a function of time and the fraction of the solvent evaporating from a solvent-swollen film, as obtained from MD simulations and parameterized SCMF simulations, shows satisfactory agreement for this simple mapping procedure. © 2005 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 43: 934–958, 2005     

Molecular dynamics simulations of carbon nanotube/silicon interfacial thermal conductance

Using molecular dynamics simulations with Tersoff reactive many-body potential for Si-Si, Si-C, and C-C interactions, we have calculated the thermal conductance at the interfaces between carbon nanotube (CNT) and silicon at different applied pressures. The interfaces are formed by axially compressing and indenting capped or uncapped CNTs against 2 x 1 reconstructed Si surfaces. The results show an increase in the interfacial thermal conductance with applied pressure for interfaces with both capped and uncapped CNTs. At low applied pressure, the thermal conductance at interface with uncapped CNTs is found to be much higher than that at interface with capped CNTs. Our results demonstrate that the contact area or the number of bonds formed between the CNT and Si substrate is key to the interfacial thermal conductance, which can be increased by either applying pressure or by opening the CNT caps that usually form in the synthesis process. The temperature and size dependences of interfacial thermal conductance are also simulated. These findings have important technological implications for the application of vertically aligned CNTs as thermal interface materials.