The directed, head-to-tail self-assembly of microtubule filaments may be generalized in the context of Janus colloidal rods. Specifically, their assembly at the tens of micron-length scale involves a careful balance between long-range electrostatic repulsion and short-range attractive forces. Here we show that the addition of counterion salts increases the rate of directed assembly by screening the electrostatic forces and enhancing the effectiveness of short-range interactions at the microtubule ends.
The aim of this study was to alter polymerization chemistry to improve network homogeneity in free-radical crosslinked systems. It was hypothesized that a reduction in heterogeneity of the network would lead to improved mechanical performance. Experiments and simulations were carried out to investigate the connection between polymerization chemistry, network structure and mechanical properties. Experiments were conducted on two different monomer systems - the first is a single monomer system, urethane dimethacrylate (UDMA), and the second is a two-monomer system consisting of bisphenol A glycidyl dimethacrylate (BisGMA) and triethylene glycol dimethacrylate (TEGDMA) in a ratio of 70/30 BisGMA/TEGDMA by weight. The methacrylate systems were crosslinked using traditional radical polymeriza- tion (TRP) with azobisisobutyronitrile (AIBN) or benzoyl peroxide (BPO) as an initiator; TRP systems were used as the control. The monomers were also cross-linked using activator regenerated by electron transfer atom transfer radical polymerization (ARGET ATRP) as a type of controlled radical polymerization (CRP). FTIR and DSC were used to monitor reac- tion kinetics of the systems. The networks were analyzed using NMR, DSC, X-ray diffraction (XRD), atomic force microscopy (AFM), and small angle X-ray scattering (SAXS). These techniques were employed in an attempt to quantify differences between the traditional and controlled radical polymerizations. While a quantitative methodology for characterizing net- work morphology was not established, SAXS and AFM have shown some promising initial results. Additionally, differences in mechanical behavior were observed between traditional and controlled radical polymerized thermosets in the BisGMA/TEGDMA system but not in the UDMA materials; this finding may be the result of network ductility variations between the two materials. Coarse-grained molecular dynamics simulations employing a novel model of the CRP reaction were carried out for the UDMA system, with parameters calibrated based on fully atomistic simulations of the UDMA monomer in the liquid state. Detailed metrics based on network graph theoretical approaches were implemented to quantify the bond network topology resulting from simulations. For a broad range of polymerization parameters, no discernible differences were seen between TRP and CRP UDMA simulations at equal conversions, although clear differences exist as a function of conversion. Both findings are consistent with experiments. Despite a number of shortcomings, these models have demonstrated the potential of molecular simulations for studying network topology in these systems.
We calculated the force-extension curves for a flexible polyelectrolyte chain with varying charge separations by performing Monte Carlo simulations of a 5000 bead chain using a screened Coulomb interaction. At all charge separations, the force-extension curves exhibit a Pincus-like scaling regime at intermediate forces and a logarithmic regime at large forces. As the charge separation increases, the Pincus regime shifts to a larger range of forces and the logarithmic regime starts are larger forces. We also found that force-extension curve for the corresponding neutral chain has a logarithmic regime. Decreasing the diameter of bead in the neutral chain simulations removed the logarithmic regime, and the force-extension curve tends to the freely jointed chain limit. This result shows that only excluded volume is required for the high force logarithmic regime to occur.
We perform molecular dynamics simulations of a coarse-grained model of ionomer melts in an applied oscillating electric field. The frequency-dependent conductivity and susceptibility are calculated directly from the current density and polarization density, respectively. At high frequencies, we find a peak in the real part of the conductivity due to plasma oscillations of the ions. At lower frequencies, the dynamic response of the ionomers depends on the ionic aggregate morphology in the system, which consists of either percolated or isolated aggregates. We show that the dynamic response of the model ionomers to the applied oscillating field can be understood by comparison with relevant time scales in the systems, obtained from independent calculations.
A coarse-grained (CG) model is developed for the thermoresponsive polymer poly(N-isopropylacrylamide) (PNIPAM), using a hybrid top-down and bottom-up approach. Nonbonded parameters are fit to experimental thermodynamic data following the procedures of the SDK (Shinoda, DeVane, and Klein) CG force field, with minor adjustments to provide better agreement with radial distribution functions from atomistic simulations. Bonded parameters are fit to probability distributions from atomistic simulations using multi-centered Gaussian-based potentials. The temperature-dependent potentials derived for the PNIPAM CG model in this work properly capture the coil-globule transition of PNIPAM single chains and yield a chain-length dependence consistent with atomistic simulations.
The performance and reliability of many mechanical and electrical components depend on the integrity of po lymer - to - solid interfaces . Such interfaces are found in adhesively bonded joints, encapsulated or underfilled electronic modules, protective coatings, and laminates. The work described herein was aimed at improving Sandia's finite element - based capability to predict interfacial crack growth by 1) using a high fidelity nonlinear viscoelastic material model for the adhesive in fracture simulations, and 2) developing and implementing a novel cohesive zone fracture model that generates a mode - mixity dependent toughness as a natural consequence of its formulation (i.e., generates the observed increase in interfacial toughness wi th increasing crack - tip interfacial shear). Furthermore, molecular dynamics simulations were used to study fundamental material/interfa cial physics so as to develop a fuller understanding of the connection between molecular structure and failure . Also reported are test results that quantify how joint strength and interfacial toughness vary with temperature.
For more than a decade now, biomolecular systems have served as an inspiration for the development of synthetic nanomaterials and systems that are capable of reproducing many of unique and emergent behaviors of living systems. One intriguing element of such systems may be found in a specialized class of proteins known as biomolecular motors that are capable of performing useful work across multiple length scales through the efficient conversion of chemical energy. Microtubule (MT) filaments may be considered within this context as their dynamic assembly and disassembly dissipate energy, and perform work within the cell. MTs are one of three cytoskeletal filaments in eukaryotic cells, and play critical roles in a range of cellular processes including mitosis and vesicular trafficking. Based on their function, physical attributes, and unique dynamics, MTs also serve as a powerful archetype of a supramolecular filament that underlies and drives multiscale emergent behaviors. In this review, we briefly summarize recent efforts to generate hybrid and composite nanomaterials using MTs as biomolecular scaffolds, as well as computational and synthetic approaches to develop synthetic one-dimensional nanostructures that display the enviable attributes of the natural filaments.
We present atomistic simulations of a single PNIPAM-alkyl copolymer surfactant in aqueous solution at temperatures below and above the LCST of PNIPAM. We compare properties of the surfactant with pure PNIPAM oligomers of similar lengths, such as the radius of gyration and solvent accessible surface area, to determine the differences in their structures and transition behavior. We also explore changes in polymer-polymer and polymer-water interactions, including hydrogen bond formation. The expected behavior is observed in the pure PNIPAM oligomers, where the backbone folds onto itself above the LCST in order to shield the hydrophobic groups from water. The surfactant, on the other hand, does not show much conformational change as a function of temperature, but instead folds to bring the hydrophobic alkyl tail and PNIPAM headgroup together at all temperatures. The atomic detail available from these simulations offers important insight into understanding how the transition behavior is changed in PNIPAM-based systems.
Designing acid- and ion-containing polymers for optimal proton, ion, or water transport would benefit profoundly from predictive models or theories that relate polymer structures with ionomer morphologies. Recently, atomistic molecular dynamics (MD) simulations were performed to study the morphologies of precise poly(ethylene-co-acrylic acid) copolymer and ionomer melts. Here, we present the first direct comparisons between scattering profiles, I(q), calculated from these atomistic MD simulations and experimental X-ray data for 11 materials. This set of precise polymers has spacers of exactly 9, 15, or 21 carbons between acid groups and has been partially neutralized with Li, Na, Cs, or Zn. In these polymers, the simulations at 120 °C reveal ionic aggregates with a range of morphologies, from compact, isolated aggregates (type 1) to branched, stringy aggregates (type 2) to branched, stringy aggregates that percolate through the simulation box (type 3). Excellent agreement is found between the simulated and experimental scattering peak positions across all polymer types and aggregate morphologies. The shape of the amorphous halo in the simulated I(q) profile is in excellent agreement with experimental I(q). The modified hard-sphere scattering model fits both the simulation and experimental I(q) data for type 1 aggregate morphologies, and the aggregate sizes and separations are in agreement. Given the stringy structure in types 2 and 3, we develop a scattering model based on cylindrical aggregates. Both the spherical and cylindrical scattering models fit I(q) data from the polymers with type 2 and 3 aggregates equally well, and the extracted aggregate radii and inter- and intra-aggregate spacings are in agreement between simulation and experiment. Furthermore, these dimensions are consistent with real-space analyses of the atomistic MD simulations. By combining simulations and experiments, the ionomer scattering peak can be associated with the average distance between branches of type 2 or 3 aggregates. This direct comparison of X-ray scattering data to the atomistic MD simulations is a substantive step toward providing a comprehensive, predictive model for ionomer morphology, gives substantial support for this atomistic MD model, and provides new credibility to the presence of stringy, branched, and percolated ionic aggregates in precise ionomer melts.
Cao, Zhen; Stevens, Mark J.; Carrillo, Jan M.; Dobrynin, Andrey V.
We use a combination of the molecular dynamics simulations and scaling analysis to study interactions between gel-like nanoparticles and substrates covered with rectangular shape posts. Our simulations have shown that nanoparticles in contact with substrate undergo a first-order transition between the Cassie-Baxter and Wenzel states, which depends on nanoparticle shear modulus, the strength of nanoparticle-substrate interactions, height of the substrate posts, and nanoparticle size, Rp. There is a range of system parameters where these two states coexist such that the average indentation δ produced by substrate posts changes with nanoparticle shear modulus, Gp. We have developed a scaling model that describes deformation of nanoparticle in contact with patterned substrate. In the framework of this model, the effect of the patterned substrate can be taken into account by introducing an effective work of adhesion, Weff, which describes the first-order transition between Wenzel and Cassie-Baxter states. There are two different shape deformation regimes for nanoparticles with shear modulus Gp and surface tension γp. The shape of small nanoparticles with size Rp < γp3/2Gp-1Weff-1/2 is controlled by capillary forces, while deformation of large nanoparticles, Rp > γp3/2Gp-1Weff-1/2, is determined by nanoparticle elastic and contact free energies. The model predictions are in good agreement with simulation results.