Molecular dynamics simulation of isothermal crystallisation of polymer chains around single polymer lamella

The isothermal crystallisation of polyethylene (PE) chains around single PE lamella in vacuum is investigated by molecular dynamic simulation. The crystallisation process is analysed in terms of the orientational order parameters, principal moments of inertia for the simulated systems. The effects of charge interactions between the polymer chains and lamella are discussed. It is found that the crystallisation process for uncharged systems can be divided into three stages: (1) adsorption, (2) orientation and (3) arrangement. The single polymer lamella changes a little during the three stages. PE chains are arranged parallel to the chain direction of the stems in the crystalline state. When considering the effect of charge interactions between the polymer chains and lamella, a different crystallisation process appears. The single polymer lamella is affected by the charged polymer chains.     

Molecular dynamics simulation of melting and crystallization processes of polyethylene clusters confined in armchair single-walled carbon nanotubes

The confined interaction is important to understand the melting and crystallization of polymers within single-wall carbon tube (SWNT). However, it is difficult for us to observe this interaction. In the current work, the structures and behaviors of melting and crystallization for polyethylene (PE) clusters confined in armchair single-walled carbon nanotubes ((n,n)-SWNTs) are investigated and examined based on molecular dynamics (MD) simulations. The nonbonded energies, structures, Lindemman indices, radial density distributions, and diffusion coefficients are used to demonstrate the features of melting phase transition for PE clusters confined in (n,n)-SWNTs. The chain end-to-end distance (R _n) and chain end-to-end distribution are used to examine the flexibility of the PE chain confined in SWNT. The global orientational order parameter (P₂) is employed to reveal the order degree of whole PE polymer. The effect of polymerization degree on melting temperature and the influence of SWNT chirality on structure of PE cluster are examined and discussed. Results demonstrate that within the confined environment of SWNT, PE clusters adopt novel co-axial crystalline layer structure, in which parallel chains of each layer are approximately vertical to tube axis. The disordered-ordered transformation of PE chains in each layer is an important structural feature for crystallization of confined PE clusters. SWNTs have a considerable effect on the structures and stabilities of the confined PE clusters.     

Buckling analysis of defective cross-linked functionalized single- and double-walled carbon nanotubes with polyethylene chains using molecular dynamics simulations

Functionalized carbon nanotubes (CNTs) can be used for improving the mechanical properties and load transfer in nanocomposites. In this research, the buckling behavior of perfect and defective cross-linked functionalized CNTs with polyethylene (PE) chains is studied employing molecular dynamics (MD) simulations. Two different configurations with the consideration of vacancy defects, namely mapped and wrapped, are selected. According to the results, critical buckling force of cross-linked functionalized CNTs with PE chains increases as compared to pure CNTs, especially in the case of double-walled carbon nanotubes (DWCNTs). By contrast, it is demonstrated that critical strain of cross-linked functionalized CNTs decreases as compared to that of pristine CNTs. Also, it is observed that increasing the weight percentage leads to the higher increase and the decrease in critical buckling force and strain of cross-linked functionalized CNTs, respectively. Moreover, the presence of defect considerably reduces both critical buckling force and strain of cross-linked functionalized CNTs. Finally, it is shown that the critical buckling strain is more sensitive to the presence of defects as compared to critical buckling force.     

Adsorption of polypeptides on solid surfaces. I. Effect of chain stiffness

The general method of obtaining the partition function and thermodynamic characteristics of polymer chains near an adsorbing surface, simulated by random walks on a lattice, is developed. The method takes into account the effect of short-range interactions in polymer chains, in particular, the chain stiffness and secondary structure. The theory of adsorption of chains of different stiffness is developed, and the process of adsorption which occurs when the external conditions change is shown to be always a second-order phase transition. The critical adsorption energy decreases and the sharpness of transition grows when the chain stiffness increases. A simple model of a chain with “virtual” steps is proposed which simplifies the treatment; the results obtained are in good agreement with exact theories. A general scheme of analysis of adsorption of chains with a given secondary structure is set forth and the analogy between the stiffness of a noncooperative chain and the presence of helical segments in a polypeptide chain is discussed.     

Conformational changes in single carboxymethylcellulose chains on a highly oriented pyrolytic graphite surface under different salt conditions

Conformational changes in individual carboxymethylcellulose (CMC) chains deposited on a highly oriented pyrolytic graphite (HOPG) surface were investigated by atomic force microscopy (AFM). A small amount of CMC solution with various salt concentrations was deposited onto the HOPG surface. The CMC molecular chains adsorbed onto the HOPG surface were clearly visualized using tapping-mode AFM under ambient conditions, as compared with those on a hydrophilic mica surface. Each CMC chain was distinguishable at the molecular level based on the vertical profiles of the AFM images, and probably aligned along the HOPG crystal lattice. Higher NaCl concentrations brought about dramatic conformational changes from aligned single chains to globular aggregates via the molecular network structure only on the HOPG surface through electrostatic screening of the CM groups. Although CMC is a water-soluble hydrophilic polyelectrolyte, some interaction, possibly due to a CH-pi bonding between the glucopyranosic axial plane of CMC and the aromatic rings of HOPG, is considered to be effective and dominant for the unique molecular attachment. These phenomena would imply the potential use of HOPG as a substrate for not only molecular imaging, but also for nano-scale morphological control of cellulosic polymers and other structural polysaccharides.     

Monte carlo study of the percolation in two-dimensional polymer systems

The structure of a two-dimensional film formed by adsorbed polymer chains was studied by means of Monte Carlo simulations. The polymer chains were represented by linear sequences of lattice beads and positions of these beads were restricted to vertices of a two-dimensional square lattice. Two different Monte Carlo methods were employed to determine the properties of the model system. The first was the random sequential adsorption (RSA) and the second one was based on Monte Carlo simulations with a Verdier-Stockmayer sampling algorithm. The methodology concerning the determination of the percolation thresholds for an infinite chain system was discussed. The influence of the chain length on both thresholds was presented and discussed. It was shown that the RSA method gave considerably lower thresholds for longer chains. This behavior can be explained by a different pool of chain conformations used in the calculations in both methods under consideration.

Accurate general method for lattice approximation of three-dimensional structure of a chain molecule

An algorithm based on dynamic programming gives the lattice models having the minimal RMS deviations from the actual folds of protein (RNA, etc.) chains for a given lattice and a given orientation of the macromolecule relative to the lattice. The algorithm is applicable for 3-D lattices of any kind. The accuracy of the lattice approximation increases when the distance between neighbor chain links is not rigidly fixed. Special repulsive potentials facilitate generation of self-avoiding lattice chains. The results of model building show the efficiency and precisionof this proposed general method when compared with others. © 1995 Wiley-Liss, Inc.     

Modelling Studies of Crystalline PEEK

Abstract

The x-ray diffraction spectrum of the engineering thermoplastic poly (aryl-ether-ether-ketone), commonly known as PEEK, has been reported in the literature. The unit cell dimensions are consistent with two aryl units per unit cell, while the chemical repeat has three, suggesting that there is not a strict lateral register between chains. The polymer chain packing of crystalline PEEK has been investigated using lattice energy calculations and it is concluded that lateral chain ordering effects are rather unimportant in determining the crystal structure. The crystallographic equivalence of ether and ketone linkages arises due to essentially random chain stacking.     

The structure of adsorbed cyclic chains

In order to determine the structure of polymer films formed of cyclic chains (rings) we developed and studied a simple coarse-grained model. Our main goal was to check how the percolation and jamming thresholds in such a system were related to the thresholds obtained for linear flexible chains system, i.e., how the geometry of objects influenced both thresholds. All atomic details were suppressed and polymers were represented as a sequence of identical beads and the chains were embedded to a square lattice (a strictly 2D model). The system was athermal and the excluded volume was the only potential introduced. A random sequential adsorption algorithm was chosen to determine the properties of a polymer monolayer. It was shown that the percolation threshold of cyclic chains was considerably higher than those of linear flexible chains while the jamming thresholds for both chain architectures are very similar. The shape of adsorbed cyclic chains was found to be more prolate when compared to average single chain.     

Molecular-level insights into the surface-induced assembly of functional bacterial amyloid

Protein coating material is important in many technological fields. The interaction between carbon nanomaterial and protein is especially interesting since it makes the development of novel hybrid materials possible. Functional bacterial amyloid (FuBA) is promising as a coating material because of its desirable features, such as well-defined molecular structure, robustness against harsh conditions, and easily engineerable functionality. Here, we report the systematic assembly of the functional amyloid protein, CsgA, from Escherichia coli (E. coli) on graphite. We characterize the assemblies using scanning tunneling microscopy (STM) and show that CsgA forms assemblies according to systematic patterns, dictated by the graphite lattice. In addition, we show that graphite flakes induce the fibrillization of CsgA, in vitro, suggesting a surface-induced conformational change of CsgA facilitated by the graphite lattice. Using coarse-grained molecular dynamics simulations, we model the adhesion and lamellar formation of a CsgA-derived peptide and conclude that peptides are adsorbed both as monomers and smaller aggregates leading initially to unordered graphite-bound aggregates, which are followed by rearrangement into lamellar structures. Finally, we show that CsgA-derived peptides can be immobilized in very systematic assemblies and their molecular orientation can be tuned using a small chaperone-like molecule. Our findings have implications for the development of FuBA-based biosensors, catalysts, and other technologies requiring well-defined protein assemblies on graphite.     

Effect of chirality,length and diameter of carbon nanotubes on the adsorption of 20 amino acids: a molecular dynamics simulation study

We carried out molecular dynamics simulations to study the adsorption of all the 20 amino acids (AAs; aromatic, polar and non-polar) on the surface of chiral, zigzag and armchair single-walled carbon nanotubes. The adsorption was occurring in all systems. In the aromatic AAs, the π–π stacking and the semi-hydrogen bond formation cause a strong interaction with the carbon nanotubes (CNTs). We also investigated the chirality, length and diameter dependencies on adsorption energies. We found that all AAs have more tendency to adsorption on the chiral and zigzag CNTs over the armchair. The results show that increasing both the diameter and the length causes the enhancement of the adsorption energy. But, the effect of the length is more than of the diameter. For example, the adsorption energy of Trp on the surface of CNT (4,1), with 2 nm length, is 20.4 kcal/mol. When the length of CNT becomes twice, the adsorption energy increases by 24 ± 0.3%. But by doubling the diameter, the adsorption energy increased only by 9.8 ± 0.25%.     

Molecular dynamics simulations of the interaction of phospholipid bilayers with polycaprolactone

The molecular interaction between common polymer chains and the cell membrane is unknown. Molecular dynamics simulations offer an emerging tool to characterise the nature of the interaction between common degradable polymer chains used in biomedical applications, such as polycaprolactone, and model cell membranes. Herein we characterise with all-atomistic and coarse-grained molecular dynamics simulations the interaction between single polycaprolactone chains of varying chain lengths with a phospholipid membrane. We find that the length of the polymer chain greatly affects the nature of interaction with the membrane, as well as the membrane properties. Furthermore, we next utilise advanced sampling techniques in molecular dynamics to characterise the two-dimensional free energy surface for the interaction of varying polymer chain lengths (short, intermediate, and long) with model cell membranes. We find that the free energy minimum shifts from the membrane-water interface to the hydrophobic core of the phospholipid membrane as a function of chain length. Finally, we perform coarse-grained molecular dynamics simulations of slightly larger membranes with polymers of the same length and characterise the results as compared with all-atomistic molecular dynamics simulations. These results can be used to design polymer chain lengths and chemistries to optimise their interaction with cell membranes at the molecular level.     

Analysis of brush-particle interactions using self-consistent-field theory

Non-specific protein adsorption can be reduced by attaching polymer chains by one end to a sorbent surface. End-grafted polymer modified surfaces have also found application in size-based chromatographic bioseparations. To better understand how to tailor surfaces for these applications, a numerical SCF model has been used to calculate theoretical results for the polymer density distribution of interacting polymer chains around a solute particle positioned at a fixed distance from a surface. In addition, the excess energy required to move the particle into the polymer chains (interaction energy) is calculated using a statistical mechanical treatment of the lattice model. The effect of system variables such as particle size, chain length, surface density and Flory interaction parameters on density distributions and interaction energies is also studied. Calculations for the interaction of a solute particle with a surface covered by many polymer chains (a brush) show that the polymer segments will fill in behind the particle quite rapidly as it moves toward the surface. When there is no strong energetic attraction between the polymer and solute we predict that the interaction energy will be purely repulsive upon compression due to losses in conformational entropy of the polymer chains. Above a critical chain length, which depends upon particle size, a maximum in the force required to move the particle toward the surface is observed due to an engulfment of the particle as chains attempt to access the free volume behind the particle. If an attraction exists between the polymer and solute, such that a minimum in the interaction energy is seen, the optimum conditions for solute repulsion occur at the highest surface density attainable. Long chain length can lead to increased solute concentration within the polymer layer due to the fact that an increased number of favourable polymer–solute contacts are able to occur than with short chains at a similar entropic penalty.     

Effect of temperature on elastic properties of CNT-polyethylene nanocomposite and its interface using MD simulations

This paper investigates the effect of temperature on the elastic modulus of carbon nanotube-polyethylene (CNT-PE) nanocomposite and its interface using molecular dynamics (MD) simulations, by utilizing the second-generation polymer consistent force field (PCFF). Two CNTs—armchair and zigzag—were selected as reinforcing nano-fillers, and amorphous PE was used as the polymer matrix. For atomistic modelling of the nanocomposite, the commercially available code Materials Studio 8.0 was used and all other MD simulations were subsequently performed using the open source code Large-Scale Atomic/Molecular Massively Parallel Simulator (LAMMPS). To obtain the elastic modulus of the nanocomposite, stress-strain curves were drawn at different temperatures by performing uniaxial deformation tests on the nanocomposite material, whereas the curvatures of the interfacial interaction energy vs. strain curves were utilized to obtain Young’s modulus of the interface. In addition, the glass transition temperatures of the polymer matrix and nanocomposites were also evaluated using density-temperature curves. Based on the results, it is concluded that, irrespective of temperature condition, a nanocomposite reinforced with CNT of larger chirality (i.e., armchair) yields a higher value of Young’s modulus of the nanocomposite and its interface. It was also found that, at the phase transition (from a glassy to a rubbery state) temperature (i.e., glass transition temperature), Young’s moduli of the polymer matrix, nanocomposite, and its interface drop suddenly. The results obtained from MD simulations were verified with results obtained from continuum-based rule-of-mixtures.     

Modulation of the self-assembled structure of biomolecules: coarse grained molecular dynamics simulation

The mechanisms governing the self-assembled structure of biomolecules (single chain and bundle of chains) are studied with an AB copolymer model via the coarse grained molecular dynamics simulations. Non-local hydrophobic interaction is found to play a critical role in the pattern formation of the assembled structure of polymer chains. We show that the polymer structure could be controlled by adjusting the balance between local (short range) and non-local (long range) hydrophobic interaction which are influenced by various factors such as the sequences, chain length, stiffness, confinement, and the topology of polymers. In addition, the competition between the intrachain hydrophobic interaction and interchain hydrophobic interaction determines the structural transition of the chain bundles. This work may provide important insights into the fundamental physics in the structure control and the self-assembly of biomolecules for various practical applications.     

Chemical modification of wheat protein-based natural polymers: grafting and cross-linking reactions with poly(ethylene oxide) diglycidyl ether and ethyl diamine

Mobile poly(ethylene oxide) diglycidyl ether (PEODGE) segments were chemically grafted onto a soluble wheat protein (WP), and different network structures were formed via coupling reactions with ethyl diamine (EDA) in different PEODGE/EDA (PE) ratios. When the PE ratio was 1:1, linear PEs were the predominant segments grafted onto WP chains and the whole WP-PEODGE-EDA (WPE) system was still soluble with an increased molecular weight. Reducing the amount of EDA in the systems produced insoluble cross-linked WPE networks. The broad distribution of network structures and chain mobility resulted in a broad glass transition for the WPE materials. However, the glass transition started at lower temperatures, and the materials became flexible at room temperature. The PE segments were present in all rigid, intermediate, and mobile phases in WPE networks, while the proportion of mobile WP chains was increased as a result of the plasticization effect from the mobile PE segments. The mobility of the most mobile component lipid was also restricted to some extent when forming the cross-linked WPE networks. The study demonstrated that the formation of different network structures with PE segments could significantly improve the flexibility of WP materials, vary the solubility, and modify the mechanical performance of WP-based natural polymer materials.     

New ordered metastable phases between the gel and subgel phases in hydrated phospholipids

Formation of low-temperature ordered gel phases in several fully hydrated phosphatidylethanolamines (PEs) and phosphatidylcholines (PCs) with saturated chains as well as in dipalmitoylphosphatidylglycerol (DPPG) was observed by synchrotron x-ray diffraction, microcalorimetry, and densitometry. The diffraction patterns recorded during slow cooling show that the gel-phase chain reflection cooperatively splits into two reflections, signaling a transformation of the usual gel phase into a more ordered phase, with an orthorhombic chain packing (the Y-transition). This transition is associated with a small decrease (2-4 microl/g) or inflection of the partial specific volume. It is fully reversible with the temperature and displays in heating direction as a small (0.1-0.7 kcal/mol) endothermic event. We recorded a Y-transition in distearoyl PE, dipalmitoyl PE (DPPE), mono and dimethylated DPPE, distearoyl PC, dipalmitoyl PC, diC(15)PC, and DPPG. No such transition exists in dimyristoyl PE and dilauroyl PE where the gel L(beta) phase transforms directly into subgel L(c) phase, as well as in the unsaturated dielaidoyl PE. The PE and PC low-temperature phases denoted L(R1) and SGII, respectively, have different hydrocarbon chain packing. The SGII phase is with tilted chains, arranged in an orthorhombic lattice of two-nearest-neighbor type. Except for the PCs, it was also registered in ionized DPPG. In the L(R1) phase, the chains are perpendicular to the bilayer plane and arranged in an orthorhombic lattice of four-nearest-neighbor type. It was observed in PEs and in protonated DPPG. The L(R1) and SGII phases are metastable phases, which may only be formed by cooling the respective gel L(beta) and L(beta') phases, and not by heating the subgel L(c) phase. Whenever present, they appear to represent an indispensable intermediate step in the formation of the latter phase.     

Molecular dynamics simulation of the folding of single alkane chains with different lengths on single-walled carbon nanotubes and graphene

Chain folding is an important step during polymer crystallization. In order to study the effects of the surface on chain folding, molecular dynamics simulations of the folding of different alkane chains on three kinds of single-walled carbon nanotubes (SWCNTs) and graphene were performed. The folding behaviors of the single alkane chains on these surfaces were found to be different from their folding behaviors in vacuum. The end-to-end distances of the chains were calculated to explore the chain folding. An increasing tendency to fold into two or more stems with increasing alkane chain length was observed. This result indicates that the occurrence and the stability of chain folding are related to the surface curvature, the diameter of the SWCNT, and surface texture. In addition, the angle between the direction of the alkane chain segment and the direction of the surface texture was measured on different surfaces.     

The cell envelope of Thermoproteus tenax: three-dimensional structure of the surface layer and its role in shape maintenance

The sulphur-dependent archaebacterium Thermoproteus tenax has a cylindrical cell shape variable in length, but constant in diameter. Its whole surface is covered by a regular protein layer (S-layer). The lattice has p6 symmetry and a lattice constant of 32.8 nm. The three-dimensional reconstruction from a tilt series of isolated and negatively stained S-layer shows a complex mass distribution of the protein: a prominent, pillar-shaped protrusion is located at the 6-fold crystallographic axis with radiating arms connecting neighbouring hexamers in the vicinity of the 3-fold axis. The base vectors of the S-layer lattice have a preferred orientation with respect to the longitudinal axis of the cell. The layer can be seen as a helical structure consisting of a right-handed, two-stranded helix, with the individual chains running parallel. Supposing that new S-layer protein is inserted at lattice faults (wedge disclinations) near the poles, growing of the layer would then proceed by moving a disclination at the end of the helix. The constant shape of the cell, as well as the particular structure of the layer, strongly suggest that this S-layer has a shape-maintaining function.     

Atomistic simulations of interfacial mechanical characteristics of carbon nanotube/silicon nanocomposites

The interfacial mechanical properties of carbon nanotube (CNT)-reinforced silicon nanocomposites are investigated by using molecular dynamics simulation method. The hybrid potential that includes Tersoff_2 potential for Si–Si in the silicon matrix, AIREBO potential for C–C in the CNTs and the Lennard-Jones (LJ) potential for Si–C in the interface is used in the nanocomposite system. The effects of such parameters as the CNT chirality, the CNT diameter and the CNT embedded length, the defects (Vacancy defects, Stone–Wales defects), the size of model, the temperature, the bonding strength and the cut-off distance of the interfacial LJ potential of nanocomposites on the pull force and the average interfacial shear strength (ISS) are investigated and discussed. The results show that the toughness and the maximum tensile strength have been increased significantly by adding the CNTs into the Si matrix. Also by increasing the LJ bonding strength and the cut-off distance of the LJ potential, the pull force and the ISS are increased significantly. The CNT chirality, the CNT diameter and the CNT embedded length have a great influence on the pull force and the ISS, while the effects of temperature, the defects and the size of model are very slight.