Toughening of thermoset polymers by rigid crystalline particles

The toughening of an aromatic amine-cured diglycidyl ether of bisphenol-A epoxy with particles of crystalline polymers was studied. The crystalline polymers were poly (butylene terephthalate), nylon 6, and poly(vinylidene fluoride). Nylon 6 and poly(vinylidene fluoride) were found to toughen the epoxy about as well as did an equivalent amount of CTBN rubber. Poly(butylene terephthalate) was found to toughen the epoxy about twice as well as did the rubber. The toughness of poly(butylene terephthalate)-epoxy blends was independent of particle size for sizes in the range of tens of micrometres, but the toughness of the nylon 6-epoxy blends decreased with increasing particle size for sizes smaller than about 40 μm. There was no loss of either Young's modulus or yield strength of the epoxy with the inclusion of either nylon 6 or poly(butylene terephthalate) and less loss of these with the inclusion of poly(vinylidene fluoride) than with the inclusion of rubber. Toughness seems to have arisen from a combination of mechanisms. The poly(butylene terephthalate)-epoxy blends alone seem to have gained toughness from phase-transformation toughening. Crack path alteration and the formation of steps and welts and secondary crack bridging seem to have accounted for an especially large part of the fracture energy of the poly(vinylidene fluoride)-epoxy blends. Secondary crack nucleation contributed to the toughness of the nylon 6-epoxy blends.     

Tight-binding molecular dynamics for materials simulations

Summary Tight-binding molecular dynamics has recently emerged as a useful method for atomistic simulation of the structural, dynamical and electronic properties of realistic materials. The method incorporates quantum-mechanical calculations into molecular dynamics through an empirical tight-binding Hamiltonian and bridges the gap between ab initio molecular dynamics and simulations using empirical classical potentials. In this paper, we review the accuracy, efficiency, and predictive power of the method and discuss some opportunities and challenges for future development.     

Million-atom molecular dynamics simulations of magnetic iron

The problem of large-scale molecular dynamics simulations of iron has recently attracted attention in connection with the need to understand the microscopic picture of radiation damage in ferritic steels. In this paper we review the development of a new interatomic potential for magnetic iron, and describe the first large-scale atomistic simulations performed using the new method. We investigate the structure and thermally activated mobility of self-interstitial atom clusters and show that the spatial distribution of magnetic moments around a cluster is well correlated with the distribution of hydrostatic pressure, highlighting the significant part played by magneto-elasticity in the treatment of radiation damage. We show that self-interstitial atom clusters exhibit a transition from relatively immobile configurations containing 〈1 1 0〉-like groups of atoms to 〈1 1 1〉-like configurations occurring at a critical cluster size N_c ∼ 5 atoms. We discuss implications of this finding for the treatment of cascade damage effects, and the possibility of observing new low-temperature resistivity recovery stages in neutron-irradiated α-iron.     

Computational limits of classical molecular dynamics simulations

The system sizes and time scales accessible by classical molecular dynamics techniques on current-generation parallel supercomputers are briefly discussed. The implications for simulation of glasses and glass-forming materials now and in the near future are highlighted.     

A multisurface constitutive model for highly cross-linked polymers with yield data obtained from molecular dynamics simulations

Constitutive properties for highly cross-linked glassy polymers are currently determined by molecular dynamics (MD) simulations. This avoids the need for ad-hoc experimentation. Constitutive data in functional form, such as yield surfaces, still require identification and correspondence to an existing function. In addition, loss of information occurs with fitting procedures. The present alternative consists in directly defining piecewise-linear yield functions from a set of points obtained by MD simulations. To prevent the algorithmic issues of multisurface plasticity, we propose an alternative to active-set strategies by simultaneously including all yield functions regardless of being active. We smooth the complementarity conditions using the Chen–Mangasarian function. In addition, extrapolation is proposed for slowly-evolving quantities such as the effective plastic strain while fully implicit integration is adopted for rapidly-evolving constitutive quantities. Since polymers exhibit finite-strain behavior, we propose a semi-implicit integration algorithm which allows a small number of steps to be used up to very large strains. Experimentally-observed effects herein considered are: thermal effects on strain (i.e. thermal expansion), Young’s modulus dependence on temperature and the effects of strain rate and temperature on the yield stress. A prototype model is first studied to assess the performance of the integration algorithm, followed by a experimental validation and a fully-featured, thermally-coupled 2D example.     

Vibrational study of a molecular device using molecular dynamics simulations

A Potential scenario for the implementation of molecular electronic systems is introduced by applying digital signal processing techniques to results from classical molecular dynamics simulations of a molecular system interconnected by nanosize gold clusters. Under this new scenario, signals can be introduced, processed, and read through interactions with the internal vibrational modes of the small molecular unit. We use modulation operations intrinsically inherent to any molecular system as a concept proof. As an example of this type of analysis, we focus on the individual oscillations between C-H and C-C bonds and cluster-cluster displacements.     

Coarse-grained molecular dynamics simulations of ionic polymer networks

The stress-strain behavior of cross-linked polymeric networks was investigated using molecular dynamics simulations with a coarse-grained representation of the repeating units. The network structure was formed by dynamically cross-linking the reactants placed between two rigid layers comprised of particles of the same type. We studied two types of networks which differ only by one containing ionic pairs that amount to 7% of the total number of bonds present. The stress-strain curves were obtained after imposing deformation in tensile and shear modes to the networks and measuring their stress response. Under both forms of deformations there was improvement in the level of stress that the material could bear. Moreover, the time dependent behavior of the improvement in mechanical properties signified a self-healing mechanism.     

Large-scale molecular dynamics simulations of hyperthermal cluster impact

Using multimillion-atom classical molecular dynamics simulations, we have studied the impact dynamics of solid and liquid spherical copper clusters (10–30 nm radius) with a solid surface, at velocities ranging from 100 m/s to 2 km/s. The resulting shock, jetting, and fragmentation processes are analyzed, demonstrating three distinct mechanisms for fragmentation. At early times, shock-induced ejection and hydrodynamic jetting produce fragments in the normal and tangential directions, respectively, while sublimation (evaporation) from the shock-heated solid (liquid) surface produces an isotropic fragment flux at both early and late times.     

Large-scale molecular dynamics simulations of fracture and deformation

Summary We have discussed the prospects of applying massively parallel molecular dynamics simulation to investigate brittle versus ductile fracture behaviors and dislocation intersection. This idea is illustrated by simulating dislocation emission from a three-dimensional crack. Unprecedentedly, the dislocation loops emitted from the crack fronts have been observed. It is found that dislocation-emission modes, jogging or blunting, are very sensitive to boundary conditions and interatomic potentials. These 3D phenomena can be effectively visualized and analyzed by a new technique, namely, plotting only those atoms within the certain ranges of local potential energies.     

Modified morphology of graphene sheets by Argon-atom bombardment: molecular dynamics simulations

By a molecular dynamics method, we simulated the process of Argon-atom bombardment on a graphene sheet with 2720 carbon atoms. The results show that, the damage of the bombardment on the graphene sheet depends not only on the incident energy but also on the particle flux density of Argon atoms. To compare and analyze the effect of the incident energy and the particle flux density in the Argon-atom bombardment, we defined the impact factor on graphene sheet by calculating the broken-hole area. The results indicate that, there is an exponential accumulated-damage for the impact of both the incident energy and the particle flux density and there is a critical incident energy ranging from 20-30 eV/atom in Argon-atom bombardment. Different configurations, such as sieve-like and circle-like graphene can be formed by controlling of different particle flux density as the incident energy is more than the critical value. Our results supply a feasible method on fabrication of porous graphene-based materials for gas-storages and molecular sieves, and it also helps to understand the damage mechanism of graphene-based electronic devices under high particle radiation.     

Study of protein adsorption on octacalcium phosphate surfaces by molecular dynamics simulations

This study numerically studies absorption of human serum albumin (HSA) and basic protein lysozyme (LSZ) on crystallographic planes of octacalcium phosphate (OCP), an essential bioactive calcium phosphate. The molecular simulations include constructing atomic structure of OCP crystallographic planes and representative segments of HSA and LSZ with three different initiate orientations respect to OCP planes. The simulation reveals the dynamic process of the protein absorption. The absorption behavior of proteins is quantified by the interaction energy between proteins and OCP planes and the strain energy of proteins in absorption. The results show that absorption interaction energy of basic LSZ is higher than that of acidic HSA, which indicates that LSZ is more favorable to adsorb onto OCP surface than HSA. The interaction energies change with the OCP crystallographic planes, the trend of changes for both proteins are similar, that is OCP (001) > OCP (111) > OCP (110) > OCP (100), which is corrected with surface energy variation of crystallographic planes. The strain energy strongly depends on the orientations of the proteins before absorption, but weakly depends on crystallographic planes. The simulation results provide useful significant information for predicting/designing interface between bioceramic materials and organic tissues as well as for understanding the mechanism of the osteoinductivity at an atomic level.     

Simulating scratch behavior of polymers with mesoscopic molecular dynamics

Part replacement and repair is needed in structures with moving parts because of scratchability and wear. In spite of some accumulation of experimental evidence, scratch resistance is still not well understood. We have applied molecular dynamics to study scratch resistance of amorphous polymeric materials through computer simulations. As a first approach, a coarse grain model was created for high density polyethylene at the mesoscale. We have also extended the traditional approach and used real units rather than reduced units (to our knowledge, for the first time), which enable an improved quantification of simulation results. The obtained results include analysis of penetration depth, residual depth and recovery percentage related to indenter force and size. Our results show there is a clear effect from these parameters on the tribological properties. We also discuss a "crooked smile" effect on the scratched surface and the reasons for its appearance.     

Million atom molecular dynamics simulations of materials on parallel computers

Recent advances in computing tecnology — parallel computer architectures, portable software and development of robust O(N) algorithms — have revolutionized the field of computer simulation. Using the space-time multiresolution molecular dynamics algorithms it is possible to carry out multimillion atom simulations of materials in different ranges of density, temperature and uniaxial strain.     

Fracture properties of metals and alloys from molecular dynamics simulations

Molecular dynamics simulations of fracture have been performed on the metals Al and Nb, and the intermetallic alloys RuAl, Nb₃Al and NiAl. The forces and energies were modelled with embedded atom method potentials. The increasing external stress was applied using displacements of the outer boundaries of the array, calculated by anisotropic elasticity theory, until the pre-existing cracks propagated or dislocation nucleation occurred. The resulting critical stress intensity factor was calculated at various orientations and temperatures, and the results compared with theory. Observations of slip systems are reported, as well as values for surface energies and “unstable stacking” energies.     

A molecular dynamics approach in simulations of fluid-solid particles flow

In the paper a discrete system of particles carried by fluid is considered in a planar motion. The volumetric density of particles is assumed to be small enough such that they can be treated within the framework of a molecular dynamics model. The fluid is then considered as a carrier of particles. The Landau-Lifshitz concept of turbulence is used to describe the fluctuating part of fluid velocity. This approach is applied to simulate different regimes (laminar and turbulent) and various states of particle motion (moving bed, heterogeneous flow, and homogeneous flow) using only two parameters, which have to be determined experimentally. These two parameters, found for a particular pipe and for a particular velocity from a simple experiment, then can be used for other pipe diameters and different velocities. The computer simulations performed for the flow of particles in pipes at different flow velocities and different pipe diameters agree favorably with experimental observations of the type of flow and critical velocities identifying transitions from one type to another. Received: 8 January 1999     

Thickness and chirality effects on tensile behavior of few-layer graphene by molecular dynamics simulations

The mechanical response of few-layer graphene (FLG), consisting of 2-7 atomic planes and bulk graphite is investigated by means of molecular dynamics simulations. By performing uniaxial tension tests at room temperature, the effects of number of atomic planes and chirality angle on the stress-strain response and deformation behavior of FLG were studied using the Tersoff potential. It was observed that by increasing of the FLG number of layers, the increase of bonding strength between neighboring layers reduce the elastic modulus and ultimate strength. It was found that, while the chirality angle of FLG showed a significant effect on the elastic modulus and ultimate tensile strength of two and three graphene layers, it turns to be less significant when the numbers of layers are more than four. Finally, by plotting the deformation behavior, it was concluded that FLGs present brittle performance.     

Generation of single-walled carbon nanotubes from alcohol and generation mechanism by molecular dynamics simulations

Recent advances in high-purity and high-yield catalytic chemical vapor deposition (CVD) generation of single-walled carbon nanotubes (SWNTs) from alcohol are comprehensively presented and discussed on the basis of results obtained from both experimental and numerical investigations. We have uniquely adopted alcohol as a carbon feedstock, and this has resulted in high-quality, low-temperature synthesis of SWNTs. This technique can produce SWNTs even at a very low temperature of 550 degrees C, which is about 300 degrees C lower than the conventional CVD methods in which methane or acetylene is typically used. We demonstrate the excellence of the proposed alcohol catalytic CVD method for high-yield production of SWNTs when Fe-Co on USY-zeolite powder was used as a catalyst. At optimum CVD conditions, a SWNT yield of more than 40 wt % was achieved over the weight of the catalytic powder within the reaction time of 120 min. In addition to the advantages for mass production, this method is also suitable for the direct synthesis of high-quality SWNTs on Si and quartz substrates when combined with the newly developed liquid-based "dip-coat" technique to mount catalytic metals on the surface of substrates. This method allows easy and costless loading of catalytic metals without the need for any support or underlayer materials that were usually required in previous studies for the generation of a sufficient quantity of SWNTs on an Si surface. Finally, the result of molecular dynamics simulation for the SWNT growth process is presented to obtain a fundamental insight into the initial growth mechanism on the catalytic particles.     

Possible phase transformation toughening of thermoset polymers by poly(butylene terephthalate)

Mechanisms were explored by which particles of poly(butylene terephthalate) (PBT) are able to toughen a brittle epoxy. The epoxy studied was an aromatic amine-cured diglycidyl ether of bisphenol-A, which was toughened at about twice the rate with particles of poly(butylene terephthalate) as with particles of nylon 6, poly(vinylidene fluoride), or CTBN rubber. Many of the mechanisms of toughening are visible on the fracture surface of the PBT-epoxy blend, but a mechanism suggested to account for perhaps half of the increased toughness with PBT, phase transformation toughening, is not. The two types of experiment performed to detect phase transformation toughening were: (1) measurements of the rubber cavitation zone in PBT-CTBN rubber-epoxy ternary blends, which would detect an expansion of the PBT particles during fracture if it occurred, and (2) measurements of the fracture energy in PBT-epoxy blends in which the various mechanisms of toughening were selectively suppressed. Both types of experiment indicated the occurrence of phase transformation toughening in these PBT-epoxy blends.