Tack of Entangled Molecular Bottlebrushes.

Polymer softness (i.e., low elastic modulus) is a known requirement for good tack in adhesives. We assess the tack performance of α-olefin molecular bottlebrushes having an elastic modulus from 4 to 30 times lower than a linear polyolefin (polypropylene). Monotonic increases of the tack parameters are observed as the bottlebrush side chain length (Nsc) increases and the modulus decreases. All-atom molecular dynamics simulations reveal that the monomeric bonding energy increases with Nsc due to a strong van der Waals interaction between the side chains and the aluminum sheet, which overcomes the energy penalty imposed by side-chain bending.

Solution properties of spherical gold nanoparticles with grafted DNA chains from simulation and theory.

There has been a rapidly growing interest in the use of functionalized Au nanoparticles (NPs) as platforms in multiple applications in medicine and manufacturing. The sensing and targeting characteristics of these NPs, and the realization of precisely organized structures in manufacturing applications using such NPs, depend on the control of their surface functionalization. NP functionalization typically takes the form of polymer grafted layers, and a detailed knowledge of the chemical and structural properties of these layers is required to molecularly engineer the particle characteristics for specific applications. However, the prediction and experimental determination of these properties to enable the rational engineering of these particles is a persistent problem in the development of this class of materials. To address this situation, molecular dynamic simulations were performed based on a previously established coarse-grained single-stranded DNA (ssDNA) model to determine basic solution properties of model ssDNA-grafted NP-layers under a wide range of conditions. In particular, we emphasize the calculation of the hydrodynamic radius for ssDNA-grafted Au NPs as a function of structural parameters such as ssDNA length, NP core size, and surface coverage. We also numerically estimate the radius of gyration and the intrinsic viscosity of these NPs, which in combination with hydrodynamic radius estimates, provide valuable information about the fluctuating structure of the grafted polymer layers. We may then understand the origin of the commonly reported variation in effective NP "size" by different measurement methods, and then exploit this information in connection to material design and characterization in connection with the ever-growing number of applications utilizing polymer-grafted NPs.

The interfacial zone in thin polymer films and around nanoparticles in polymer nanocomposites.

We perform coarse-grained simulations of model unentangled polymer materials to quantify the range over which interfaces alter the structure and dynamics in the vicinity of the interface. We study the interfacial zone around nanoparticles (NPs) in model polymer-NP composites with variable NP diameter, as well as the interfacial zone at the solid substrate and free surface of thin supported polymer films. These interfaces alter both the segmental packing and mobility in an interfacial zone. Variable NP size allows us to gain insight into the effect of boundary curvature, where the film is the limit of zero curvature. We find that the scale for perturbations of the density is relatively small and decreases on cooling for all cases. In other words, the interfaces become more sharply defined on cooling, as naively expected. In contrast, the interfacial mobility scale ξ for both NPs and supported films increases on cooling and is on the order of a few nanometers, regardless of the polymer-interfacial interaction strength. Additionally, the dynamical interfacial scale of the film substrate is consistent with a limiting value for polymer-NP composites as the NP size grows. These findings are based on a simple quantitative model to describe the distance dependence of relaxation that should be applicable to many interfacial polymer materials.

Influence of knot complexity on glass-formation in low molecular mass ring polymer melts.

We perform molecular dynamics simulations on a coarse-grained polymer melt to study the dynamics of glass-formation in ring polymer melts of variable knot complexity. After generating melts of non-concatenated polymeric rings having a range of minimum crossing number values, mc, we compute the coherent intermediate scattering function, the segmental α-relaxation time, fragility, and the glass transition temperature as a function of mc. Variation of knot complexity is found to have a pronounced effect on the dynamics of polymer melts since both molecular rigidity and packing are altered, primary physical factors governing glass-formation in polymeric materials.

Communication: A comparison between the solution properties of knotted ring and star polymers.

We computationally investigate the good solvent solution properties of knotted ring and star polymers by combining molecular dynamics (MD) simulation and path-integral calculations. We consider knotted rings having a minimal crossing number m c in the range, 0 ≤ m c ≤ 9, and star polymers having a range of f star arms, 2 ≤ f ≤ 20, attached to a common core monomer particle. After generating configurational ensembles of these polymers by MD, we use the path-integration program ZENO to calculate basic configurational properties, i.e., radius of gyration, hydrodynamic radius, intrinsic viscosity, as well as fluctuations in these properties. Our simulations indicate that the configurational properties of knotted rings and star polymers in solution show a similar decrease with increasing m c and f. Moreover, fluctuations in these properties also decrease with increasing topological complexity. Our findings should be helpful in polymer characterization and more generally for understanding the role of polymer topology in polymer material properties.

Structure and Dynamics of a Graphene Melt.

We explore the structural and dynamic properties of bulk materials composed of graphene nanosheets using coarse-grained molecular dynamics simulations. Remarkably, our results show clear evidence that bulk graphene materials exhibit a fluid-like behavior similar to linear polymer melts at elevated temperatures and that these materials transform into a glassy-like "foam" state at temperatures below the glass-transition temperature ( Tg) of these materials. Distinct from an isolated graphene sheet, which exhibits a relatively flat shape with fluctuations, we find that graphene sheets in a melt state structurally adopt more "crumpled" configurations and correspondingly smaller sizes, as normally found for ordinary polymers in the melt. Upon approaching the glass transition, these two-dimensional polymeric materials exhibit a dramatic slowing down of their dynamics that is likewise similar to ordinary linear polymer glass-forming liquids. Bulk graphene materials in their glassy foam state have an exceptionally large free-volume and high thermal stability due to their high Tg (≈ 1600 K) as compared to conventional polymer materials. Our findings show that graphene melts have interesting lubricating and "plastic" flow properties at elevated temperatures, and suggest that graphene foams are highly promising as high surface filtration materials and fire suppression additives for improving the thermal conductivities and mechanical reinforcement of polymer materials.

Fiber Network Formation in Semi-Flexible Polymer Solutions: An Exploratory Computational Study.

The formation of gels through the bundling of semi-flexible polymer chains into fiber networks is ubiquitous in diverse manufactured and natural materials, and, accordingly, we perform exploratory molecular dynamics simulations of a coarse-grained model of semi-flexible polymers in a solution with attractive lateral interchain interactions to understand essential features of this type of gel formation. After showing that our model gives rise to fibrous gels resembling real gels of this kind, we investigate how the extent of fiber bundling influences the "melting" temperature, T m , and the emergent rigidification of model bundled fibers having a fixed number of chains, N, within them. Based on our preliminary observations, we suggest the fiber size is kinetically selected by a reduced thermodynamic driving force and a slowing of the dynamics within the fibers associated with their progressive rigidification with the inclusion of an increasing number of chains in the bundle.

Molecular rigidity and enthalpy-entropy compensation in DNA melting.

Enthalpy-entropy compensation (EEC) is observed in diverse molecular binding processes of importance to living systems and manufacturing applications, but this widely occurring phenomenon is not sufficiently understood from a molecular physics standpoint. To gain insight into this fundamental problem, we focus on the melting of double-stranded DNA (dsDNA) since measurements exhibiting EEC are extensive for nucleic acid complexes and existing coarse-grained models of DNA allow us to explore the influence of changes in molecular parameters on the energetic parameters by using molecular dynamics simulations. Previous experimental and computational studies have indicated a correlation between EEC and changes in molecular rigidity in certain binding-unbinding processes, and, correspondingly, we estimate measures of DNA molecular rigidity under a wide range of conditions, along with resultant changes in the enthalpy and entropy of binding. In particular, we consider variations in dsDNA rigidity that arise from changes of intrinsic molecular rigidity such as varying the associative interaction strength between the DNA bases, the length of the DNA chains, and the bending stiffness of the individual DNA chains. We also consider extrinsic changes of molecular rigidity arising from the addition of polymer additives and geometrical confinement of DNA between parallel plates. All our computations confirm EEC and indicate that this phenomenon is indeed highly correlated with changes in molecular rigidity. However, two distinct patterns relating to how DNA rigidity influences the entropy of association emerge from our analysis. Increasing the intrinsic DNA rigidity increases the entropy of binding, but increases in molecular rigidity from external constraints decreases the entropy of binding. EEC arises in numerous synthetic and biological binding processes and we suggest that changes in molecular rigidity might provide a common origin of this ubiquitous phenomenon in the mutual binding and unbinding of complex molecules.

Knot Energy, Complexity, and Mobility of Knotted Polymers.

The Coulomb energy E C is defined by the energy required to charge a conductive object and scales inversely to the self-capacity C, a basic measure of object size and shape. It is known that C is minimized for a sphere for all objects having the same volume, and that C increases as the symmetry of an object is reduced at fixed volume. Mathematically similar energy functionals have been related to the average knot crossing number 〈m〉, a natural measure of knot complexity and, correspondingly, we find E C to be directly related to 〈m〉 of knotted DNA. To establish this relation, we employ molecular dynamics simulations to generate knotted polymeric configurations having different length and stiffness, and minimum knot crossing number values m for a wide class of knot types relevant to the real DNA. We then compute E C for all these knotted polymers using the program ZENO and find that the average Coulomb energy 〈E C〉 is directly proportional to 〈m〉. Finally, we calculate estimates of the ratio of the hydrodynamic radius, radius of gyration, and the intrinsic viscosity of semi-flexible knotted polymers in comparison to the linear polymeric chains since these ratios should be useful in characterizing knotted polymers experimentally.

Universal interrelation between measures of particle and polymer size.

The characterization of many objects involves the determination of a basic set of particle size measures derived mainly from scattering and transport property measurements. For polymers, these basic properties include the radius of gyration Rg, hydrodynamic radius Rh, intrinsic viscosity [η], and sedimentation coefficient S, and for conductive particles, the electric polarizability tensor αE and self-capacity C. It is often found that hydrodynamic measurements of size deviate from each other and from geometric estimates of particle size when the particle or polymer shape is complex, a phenomenon that greatly complicates both nanoparticle and polymer characterizations. The present work explores a general quantitative relation between αE, C, and Rg for nanoparticles and polymers of general shape and the corresponding properties η, Rh, and Rg using a hydrodynamic-electrostatic property interrelation.

Giant Surface Conductivity Enhancement in a Carbon Nanotube Composite by Ultraviolet Light Exposure.

Carbon nanotube composites are lightweight, multifunctional materials with readily adjustable mechanical and electrical properties-relevant to the aerospace, automotive, and sporting goods industries as high-performance structural materials. Here, we combine well-established and newly developed characterization techniques to demonstrate that ultraviolet (UV) light exposure provides a controllable means to enhance the electrical conductivity of the surface of a commercial carbon nanotube-epoxy composite by over 5 orders of magnitude. Our observations, combined with theory and simulations, reveal that the increase in conductivity is due to the formation of a concentrated layer of nanotubes on the composite surface. Our model implies that contacts between nanotube-rich microdomains dominate the conductivity of this layer at low UV dose, while tube-tube transport dominates at high UV dose. Further, we use this model to predictably pattern conductive traces with a UV laser, providing a facile approach for direct integration of lightweight conductors on nanocomposite surfaces.

Hydrodynamic Radius Fluctuations in Model DNA-Grafted Nanoparticles.

We utilize molecular dynamics simulations (MD) and the path-integration program ZENO to quantify hydrodynamic radius (Rh) fluctuations of spherical symmetric gold nanoparticles (NPs) decorated with single-stranded DNA chains (ssDNA). These results are relevant to understanding fluctuation-induced interactions among these NPs and macromolecules such as proteins. In particular, we explore the effect of varying the ssDNA-grafted NPs structural parameters, such as the chain length (L), chain persistence length (lp), NP core size (R), and the number of chains (N) attached to the nanoparticle core. We determine Rh fluctuations by calculating its standard deviation (σRh ) of an ensemble of ssDNA-grafted NPs configurations generated by MD. For the parameter space explored in this manuscript, σRh shows a peak value as a function of N, the amplitude of which depends on L, lp and R, while the broadness depends on R.

Intrinsic conductivity of carbon nanotubes and graphene sheets having a realistic geometry.

The addition of carbon nanotubes (CNTs) and graphene sheets (GSs) into polymeric materials can greatly enhance the conductivity and alter the electromagnetic response of the resulting nanocomposite material. The extent of these property modifications strongly depends on the structural parameters describing the CNTs and GSs, such as their shape and size, as well as their degree of particle dispersion within the polymeric matrix. To model these property modifications in the dilute particle regime, we determine the leading transport virial coefficients describing the conductivity of CNT and GS composites using a combination of molecular dynamics, path-integral, and finite-element calculations. This approach allows for the treatment of the general situation in which the ratio between the conductivity of the nanoparticles and the polymer matrix is arbitrary so that insulating, semi-conductive, and conductive particles can be treated within a unified framework. We first generate ensembles of CNTs and GSs in the form of self-avoiding worm-like cylinders and perfectly flat and random sheet polymeric structures by using molecular dynamics simulation to model the geometrical shapes of these complex-shaped carbonaceous nanoparticles. We then use path-integral and finite element methods to calculate the electric and magnetic polarizability tensors (αE, αM) of the CNT and GS nanoparticles. These properties determine the conductivity virial coefficient σ in the conductive and insulating particle limits, which are required to estimate σ in the general case in which the conductivity contrast Δ between the nanoparticle and the polymer matrix is arbitrary. Finally, we propose approximate relationships for αE and αM that should be useful in materials design and characterization applications.

Dimensional reduction of duplex DNA under confinement to nanofluidic slits.

There has been much interest in the dimensional properties of double-stranded DNA (dsDNA) confined to nanoscale environments as a problem of fundamental importance in both biological and technological fields. This has led to a series of measurements by fluorescence microscopy of single dsDNA molecules under confinement to nanofluidic slits. Despite the efforts expended on such experiments and the corresponding theory and simulations of confined polymers, a consistent description of changes of the radius of gyration of dsDNA under strong confinement has not yet emerged. Here, we perform molecular dynamics (MD) simulations to identify relevant factors that might account for this inconsistency. Our simulations indicate a significant amplification of excluded volume interactions under confinement at the nanoscale due to the reduction of the effective dimensionality of the system. Thus, any factor influencing the excluded volume interaction of dsDNA, such as ionic strength, solution chemistry, and even fluorescent labels, can greatly influence the dsDNA size under strong confinement. These factors, which are normally less important in bulk solutions of dsDNA at moderate ionic strengths because of the relative weakness of the excluded volume interaction, must therefore be under tight control to achieve reproducible measurements of dsDNA under conditions of dimensional reduction. By simulating semi-flexible polymers over a range of parameter values relevant to the experimental systems and exploiting past theoretical treatments of the dimensional variation of swelling exponents and prefactors, we have developed a novel predictive relationship for the in-plane radius of gyration of long semi-flexible polymers under slit-like confinement. Importantly, these analytic expressions allow us to estimate the properties of dsDNA for the experimentally and biologically relevant range of contour lengths that is not currently accessible by state-of-the-art MD simulations.

Confronting the complexity of CNT materials.

The morphology of commercially available carbon nanotube materials is often much more complex than the term "carbon nanotube" (CNT) would imply. Commercial CNT materials are typically composed of roughly spherical CNT domains having a highly ramified internal structure and a size on the order of microns. Clearly, such structures cannot reasonably be modeled as "rods". To address this problem, we first perform molecular dynamics simulations (MD) to generate structures similar to those measured experimentally, based on the presumptions that CNT domains are composed of worm-like cylinders having observed persistence lengths and that these CNTs are confined to spherical domains having the observed average domain size. This simple model generates structures remarkably similar to those observed experimentally. We then consider numerical path-integral computations to calculate the self-capacitance C and intrinsic conductivity [σ]∞ of these CNT rich domains. This information is then incorporated in a generalized effective medium theory to estimate the conductivity of bulk composite materials composed of these complex-shaped "particles". We term these CNT structures "tumbleweeds", given their evident morphological similarity to this naturally occurring growth form. Based on this model, we find that the conductivity percolation threshold of the tumbleweeds can be quite low, despite their quasi-spherical average shape. We also examine the structure factor S(q) of the CNT-rich domains as function of the number N of CNTs within them, to aid in the structural characterization of CNT nanocomposites. The structure factor S(q) of our model tumbleweed is found to resemble that of hyperbranched, star and dendrimer polymers, and also domain structures observed in polyelectrolytes. Commercial CNT materials at high loading should then have physical features in common with suspension of "soft" colloidal particles by virtue of their deformability and roughly spherical shape.

High-speed, high-purity separation of gold nanoparticle-DNA origami constructs using centrifugation.

DNA origami is a powerful platform for assembling gold nanoparticle constructs, an important class of nanostructure with numerous applications. Such constructs are assembled by the association of complementary DNA oligomers. These association reactions have yields of <100%, requiring the development of methods to purify the desired product. We study the performance of centrifugation as a separation approach by combining optical and hydrodynamic measurements and computations. We demonstrate that bench-top microcentrifugation is a simple and efficient method of separating the reaction products, readily achieving purities of >90%. The gold nanoparticles play a number of critical roles in our system, functioning not only as integral components of the purified products, but also as hydrodynamic separators and optical indicators of the reaction products during the purification process. We find that separation resolution is ultimately limited by the polydispersity in the mass of the gold nanoparticles and by structural distortions of DNA origami induced by the gold nanoparticles. Our study establishes a methodology for determining the design rules for nanomanufacturing DNA origami-nanoparticle constructs.

Internal structure of nanoparticle dimers linked by DNA.

We construct nanoparticle dimers linked by DNA. These dimers are basic units in a possible multiscale, hierarchical assembly and serve as a model system to understand DNA-mediated interactions, especially in the nontrivial regime when the nanoparticle and DNA are comparable in their sizes. We examine the structure of nanoparticle dimers in detail by a combination of scattering experiments and molecular simulations. We find that, for a given DNA length, the interparticle separation within the dimer is controlled primarily by the number of linking DNA. We summarize our findings in a simple model that captures the interplay of the number of DNA bridges, their length, the particle's curvature, and the excluded volume effects. We demonstrate the applicability of the model to our results, without any free parameters. As a consequence, the increase of dimer separation with increasing temperature can be understood as a result of changing the number of connecting DNA.