Nanoscale finite element models for vibrations of single-walled carbon nanotubes: atomistic versus continuum

By the atomistic and continuum finite element models, the free vibration behavior of single-walled carbon nanotubes （SWCNTs） is studied. In the atomistic finite element model, the bonds and atoms are modeled by the beam and point mass elements, respectively. The molecular mechanics is linked to structural mechanics to determine the elastic properties of the mentioned beam elements. In the continuum finite element approach, by neglecting the discrete nature of the atomic structure of the nanotubes, they are modeled with shell elements. By both models, the natural frequencies of SWCNTs are computed, and the effects of the geometrical parameters, the atomic structure, and the boundary conditions are investigated. The accuracy of the utilized methods is verified in comparison with molecular dynamic simulations. The molecular structural model leads to more reliable results, especially for lower aspect ratios. The present analysis provides valuable information about application of continuum models in the investigation of the mechanical behaviors of nanotubes.     

Atomistic simulations of elastic deformation and dislocation nucleation during nanoindentation

Nanoindentation experiments have shown that microstructural inhomogeneities across the surface of gold thin films lead to position-dependent nanoindentation behavior [Phys. Rev. B (2002), to be submitted]. The rationale for such behavior was based on the availability of dislocation sources at the grain boundary for initiating plasticity. In order to verify or refute this theory, a computational approach has been pursued. Here, a simulation study of the initial stages of indentation using the embedded atom method (EAM) is presented. First, the principles of the EAM are given, and a comparison is made between atomistic simulations and continuum models for elastic deformation. Then, the mechanism of dislocation nucleation in single crystalline gold is analyzed, and the effects of elastic anisotropy are considered. Finally, a systematic study of the indentation response in the proximity of a high angle, high sigma (low symmetry) grain boundary is presented; indentation behavior is simulated for varying indenter positions relative to the boundary. The results indicate that high angle grain boundaries are a ready source of dislocations in indentation-induced deformation.     

Non-equilibrium grain boundary structure and inelastic deformation using atomistic simulations

Grain boundary influence on material properties becomes increasingly significant as grain size is reduced to the nanoscale. Nanostructured materials produced by severe plastic deformation techniques often contain a higher percentage of high-angle grain boundaries in a non-equilibrium or energetically metastable state. Differences in the mechanical behavior and observed deformation mechanisms are common due to deviations in grain boundary structure. Fundamental interfacial attributes such as atomic mobility and energy are affected due to a higher non-equilibrium state, which in turn affects deformation response. In this research, atomistic simulations employing a biased Monte Carlo method are used to approximate representative non-equilibrium bicrystalline grain boundaries based on an embedded atom method potential, leveraging the concept of excess free volume. An advantage of this approach is that non-equilibrium boundaries can be instantiated without the need of simulating numerous defect/grain boundary interactions. Differences in grain boundary structure and deformation response are investigated as a function of non-equilibrium state using Molecular Dynamics. A detailed comparison between copper and aluminum bicrystals is provided with regard to boundary strength, observed deformation mechanisms, and stress-assisted free volume evolution during both tensile and shear simulations.     

Multiscale modeling of crack initiation and propagation at the nanoscale

Fracture occurs on multiple interacting length scales; atoms separate on the atomic scale while plasticity develops on the microscale. A dynamic multiscale approach (CADD: coupled atomistics and discrete dislocations) is employed to investigate an edge-cracked specimen of single-crystal nickel, Ni, (brittle failure) and aluminum, Al, (ductile failure) subjected to mode-I loading. The dynamic model couples continuum finite elements to a fully atomistic region, with key advantages such as the ability to accommodate discrete dislocations in the continuum region and an algorithm for automatically detecting dislocations as they move from the atomistic region to the continuum region and then correctly “converting” the atomistic dislocations into discrete dislocations, or vice-versa. An ad hoc computational technique is also applied to dissipate localized waves formed during crack advance in the atomistic zone, whereby an embedded damping zone at the atomistic/continuum interface effectively eliminates the spurious reflection of high-frequency phonons, while allowing low-frequency phonons to pass into the continuum region.The simulations accurately capture the essential physics of the crack propagation in a Ni specimen at different temperatures, including the formation of nano-voids and the sudden acceleration of the crack tip to a velocity close to the material Rayleigh wave speed. The nanoscale brittle fracture happens through the crack growth in the form of nano-void nucleation, growth and coalescence ahead of the crack tip, and as such resembles fracture at the microscale. When the crack tip behaves in a ductile manner, the crack does not advance rapidly after the pre-opening process but is blunted by dislocation generation from its tip. The effect of temperature on crack speed is found to be perceptible in both ductile and brittle specimens.     

Simulation and study of the modulus of elasticity of nanocrystalline materials

ＩｎｔｒｏｄｕｃｔｉｏｎＩｎｔｈｅ１９８０’ｓ,ｎａｎｏｃｒｙｓｔａｌｉｎｅｍａｔｅｒｉａｌｓ（ＮＣＭ）ｈａｖｅｂｅｅｎｐｒｏｐｏｓｅｄｂｙｐｒｏｆｅｓｏｒＨ．ＧｌｅｉｔｅｒｉｎＧｅｒｍａｎ．Ｔｈｅａｃａｄｅｍｉｃｗｏｒｌｄｗａｓｉｎｔｅｒｅｓｔｅｄ...     

Molecular/cluster statistical thermodynamics methods to simulate quasi-static deformations at finite temperature

The rapid evolution of nanotechnology appeals for the understanding of global response of nanoscale systems based on atomic interactions, hence necessitates novel, sophisticated, and physically based approaches to bridge the gaps between various length and time scales. In this paper, we propose a group of statistical thermodynamics methods for the simulations of nanoscale systems under quasi-static loading at finite temperature, that is, molecular statistical thermodynamics (MST) method, cluster statistical thermodynamics (CST) method, and the hybrid molecular/cluster statistical thermodynamics (HMCST) method. These methods, by treating atoms as oscillators and particles simultaneously, as well as clusters, comprise different spatial and temporal scales in a unified framework. One appealing feature of these methods is their “seamlessness” or consistency in the same underlying atomistic model in all regions consisting of atoms and clusters, and hence can avoid the ghost force in the simulation. On the other hand, compared with conventional MD simulations, their high computational efficiency appears very attractive, as manifested by the simulations of uniaxial compression and nanoindenation.     

Evolution of nanoscale defects to planar cracks in a brittle solid

Fracture of a solid is a highly multiscale process that associates atomic scale bond breaking with macroscopic crack propagation, and the process can be dramatically influenced by the presence of defects in materials. In a nanomaterial, defect formation energy decreases with the reduction of material size, and therefore, the role of defects in crack formation and subsequent crack growth in such materials may not be understood from the classical laws of fracture mechanism. In this study, we investigated the crack formation process of a defective (with missing atoms) nanostructured material (NaCl) using a series of molecular dynamics (MD) simulations. It was demonstrated that simple defects in the form of several missing atoms in the material could develop into a planar crack. Subsequently, MD simulations on failures of nanosized NaCl with pre-defined planar atomistic cracks of two different lengths under prescribed tensile displacement loads were performed. These failure loads were then applied on the equivalent continuum models, separately, to evaluate the associated fracture toughness values using the finite element analysis. For small cracks, the fracture toughness thus obtained is cracksize dependent and the corresponding critical energy release rate is significantly smaller than Griffith’s theoretical value. Explanation for this discrepancy between LEFM and the atomistic model was attempted.     

A generalized self-consistent polycrystal model for the yield strength of nanocrystalline materials

Inspired by recent molecular dynamic simulations of nanocrystalline solids, a generalized self-consistent polycrystal model is proposed to study the transition of yield strength of polycrystalline metals as the grain size decreases from the traditional coarse grain to the nanometer scale. These atomic simulations revealed that a significant portion of atoms resides in the grain boundaries and the plastic flow of the grain-boundary region is responsible for the unique characteristics displayed by such materials. The proposed model takes each oriented grain and its immediate grain boundary to form a pair, which in turn is embedded in the infinite effective medium with a property representing the orientational average of all these pairs. We make use of the linear comparison composite to determine the nonlinear behavior of the nanocrystalline polycrystal through the concept of secant moduli. To this end an auxiliary problem of Christensen and Lo (J. Mech. Phys. Solids 27 (1979) 315) superimposed on the eigenstrain field of Luo and Weng (Mech. Mater. 6 (1987) 347) is first considered, and then the nonlinear elastoplastic polycrystal problem is addressed. The plastic flow of each grain is calculated from its crystallographic slips, but the plastic behavior of the grain-boundary phase is modeled as that of an amorphous material. The calculated yield stress for Cu is found to follow the classic Hall-Petch relation initially, but as the gain size decreases it begins to depart from it. The yield strength eventually attains a maximum at a critical grain size and then the Hall-Petch slope turns negative in the nano-range. It is also found that, when the Hall-Petch relation is observed, the plastic behavior of the polycrystal is governed by crystallographic slips in the grains, but when the slope is negative it is governed by the grain boundaries. During the transition both grains and grain boundaries contribute competitively.     

Hyper-QC: An accelerated finite-temperature quasicontinuum method using hyperdynamics

The quasicontinuum (QC) method is a spatial multiscale method that extends the length scales accessible to fully atomistic simulations (like molecular dynamics (MD)) by several orders of magnitude. While the recent development of the so-called “hot-QC method” enables dynamic simulations at finite temperature, the times accessible to these simulations remain limited to the sub-microsecond time scale due to the small time step required for stability of the numerical integration. To address this limitation, we develop a novel finite-temperature QC method that can treat much longer time scales by coupling the hot-QC method with hyperdynamics—a method for accelerating time in MD simulations. We refer to the new approach as “hyper-QC”. As in the original hyperdynamics method, hyper-QC is targeted at dynamical systems that exhibit a separation of time scales between short atomic vibration periods and long waiting times in metastable states. Acceleration is achieved by modifying the hot-QC potential energy to reduce the energy barriers between metastable states in a manner that ensures that the characteristic dynamics of the system are preserved. First, the high accuracy of hot-QC in reproducing rare event kinetics is demonstrated. Then, the hyper-QC methodology is validated by comparing hyper-QC results with those of hot-QC and full MD for a one-dimensional chain of atoms interacting via a Lennard–Jones potential.     

Traction–separation laws and stick–slip shear phenomenon of interfaces between cellulose nanocrystals

Cellulose nanocrystals (CNCs) are one of nature's most abundant structural material building blocks and possess outstanding mechanical properties including a tensile modulus comparable to Kevlar. It remains challenging to upscale these properties in CNC neat films and nanocomposites due to the difficulty of characterizing interfacial bonding between CNCs that governs stress transfer under deformation. Here we present new analyses based on atomistic simulations of shear and tensile failure of the interfaces between Iβ CNCs, providing new insight into factors governing the mechanical behavior of hierarchical nanocellulose materials. We compare the two most relevant crystal interfaces and find that hydrogen bonded surfaces have greater tensile strength compared to the surfaces governed by weaker interactions. On the contrary, shearing simulations reveal that friction between the atomic interfaces depends not only on surface energy but also the energy landscape along the shear direction. While being a weaker interface, the intersheet plane exhibits greater energy barriers to shear. The molecular roughness of this interface, characterized by a greater energy barrier, exhibits stick–slip deformation behavior as opposed to a more continuous sliding and rebonding mechanism observed for the interfaces with hydrogen bonds. Analytical models to describe the energy landscapes are developed using energy scaling relations for van der Waals surfaces in combination with a modification of the Prandtl–Tomlinson model for atomic friction. Our simulations pave the way for tailoring hierarchical CNC materials by taking a similar approach to techniques employed for describing metals, where mechanical properties can be tuned through a deeper understanding of grain boundary physics and nanoscale interfaces.     

Deformation gradients for continuum mechanical analysis of atomistic simulations

We present an expression developed for calculating an atomic-scale deformation gradient within atomistic simulations. This expression is used to analyze the deformation fields for a one-dimensional atomic chain, a biaxially stretched thin film containing a surface ledge, and a FCC metal subject to indentation loading from a nanometer-scale indenter. The analyses presented show that the metric established here is consistent with the continuum mechanical concept of deformation gradient (which is known to have a zero curl for compatible deformations) in most instances. However, our metric does yield non-zero values of curl for atoms near loaded geometric inhomogeneities, such as those that form the ledges themselves and those beneath or adjacent to the indentation contact region. Also, we present expressions for higher order gradients of the deformation field and discuss the requirements for their calculation. These expressions are necessary for linking atomistic simulation results with advanced continuum mechanics theories such as strain gradient plasticity, thereby enabling fundamental, atomic-scale information to contribute to the formulation and parameterization of such theories.     

Multiscale diffusion method for simulations of long-time defect evolution with application to dislocation climb

In many problems of interest to materials scientists and engineers, the evolution of crystalline extended defects (dislocations, cracks, grain boundaries, interfaces, voids, precipitates) is controlled by the flow of point defects (interstitial/substitutional atoms and/or vacancies) through the crystal into the extended defect. Precise modeling of this behavior requires fully atomistic methods in and around the extended defect, but the flow of point defects entering the defect region can be treated by coarse-grained methods. Here, a multiscale algorithm is presented to provide this coupling. Specifically, direct accelerated molecular dynamics (AMD) of extended defect evolution is coupled to a diffusing point defect concentration field that captures the long spatial and temporal scales of point defect motion in the presence of the internal stress fields generated by the evolving defect. The algorithm is applied to study vacancy absorption into an edge dislocation in aluminum where vacancy accumulation in the core leads to nucleation of a double-jog that then operates as a sink for additional vacancies; this corresponds to the initial stages of dislocation climb modeled with explicit atomistic resolution. The method is general and so can be applied to many other problems associated with nucleation, growth, and reaction due to accumulation of point defects in crystalline materials.     

Dislocation nucleation from bicrystal interfaces and grain boundary ledges: Relationship to nanocrystalline deformation

Molecular dynamics simulations are used to evaluate the primary interface dislocation sources and to estimate both the free enthalpy of activation and the critical emission stress associated with the interfacial dislocation emission mechanism. Simulations are performed on copper to study tensile failure of a planar Σ5 {2 1 0} 53.1° interface and an interface with the same misorientation that contains a ledge. Simulations reveal that grain boundary ledges are more favorable as dislocation sources than planar regions of the interface and that their role is not limited to that of simple dislocation donors. The parameters extracted from the simulations are utilized in a two-phase composite mesoscopic model for nanocrystalline deformation that includes the effects of both dislocation emission and dislocation absorption mechanisms. A self-consistent approach based on the Eshelby solution for grains as ellipsoidal inclusions is augmented by introduction of stress concentration in the constitutive law of the matrix phase to account for more realistic grain boundary effects. Model simulations suggest that stress concentration is required in the standard continuum theory to activate the coupled grain boundary dislocation emission and absorption mechanisms when activation energy of the dislocation source is determined from atomistic calculation on grain boundaries without consideration of impurities or other extrinsic defects.     

Atomistic simulation study on key factors dominating dislocation nucleation from a crack tip in two FCC materials: Cu and Al

Dislocation nucleations from crack tips in FCC copper and aluminum are studied using atomistic simulations. It is shown that the critical load for dislocation nucleation predicted by Rice’s model (Rice, 1992) based on the Peierls concept of dislocation can either be under- or over-estimated in reference to the simulation results. Such discrepancies have not been fully resolved by existing improved nucleation models, due to the complicated atomic environments at crack tips. Based on our simulation results, it is proposed that such discrepancies can be reconciled by the competition of two coupling processes at a crack tip: the tension-shear coupling, which facilitates the dislocation nucleation, and the nucleation-debonding coupling, which retards the dislocation nucleation. In addition, the two couplings are applied to explain the paradoxical observation: easy dislocation nucleation at a blunted crack tip. The present work provides a detailed picture to justify future improvements on Rice’s model for dislocation nucleation and to accurately predict intrinsic brittle to ductile transition for crystalline materials.     

Estimation of anisotropic elastic properties of nanocomposites using atomistic-continuum interphase model

We have revised classical micromechanics by accounting for the effect of interface to predict the effective anisotropic elastic properties of heterogeneous materials containing nano-inhomogeneities. In contrast to sharp interface between the matrix and inhomogeneity, we introduce the concept of interphase to account for the interfacial-stress effect at the nano-scale. The interphase’s constitutive properties are derived from atomistic simulations and then incorporated in a micromechanics-based interphase model to compute effective properties of nanocomposites. This scale transition approach bridges the gap between discrete atomic level interactions and continuum mechanics. An advantage of this approach is that it combines atomistic with continuum models that consider inhomogeneity and interphase morphology. It thereby enables us to account simultaneously for both the shape and the anisotropy of a nano-inhomogeneity and interphase at the continuum level when we compute material’s overall properties. In so doing, it frees us from making any assumptions about the interface characteristics between matrix and the nano-inhomogeneity.     

Influence of anisotropic elasticity on the mechanical properties of fivefold twinned nanowires

Previous atomistic simulations and experiments have shown an increased Young's modulus and yield strength of fivefold twinned (FT) face-centered cubic metal nanowires (NWs) when compared to single crystalline (SC) NWs of the same orientation. Here we report the results of atomistic simulations of SC and FT Ag, Al, Au, Cu and Ni NWs with diameters between 2 and 50 nm under tension and compression. The simulations show that the differences in Young's modulus between SC and FT NWs are correlated with the elastic anisotropy of the metal, with Al showing a decreased Young's modulus. We develop a simple analytical model based on disclination theory and constraint anisotropic elasticity to explain the trend in the difference of Young's modulus between SC and FT NWs. Taking into account the role of surface stresses and the elastic properties of twin boundaries allows to account for the observed size effect in Young's modulus. The model furthermore explains the different relative yield strengths in tension and compression as well as the material and loading dependent failure mechanisms in FTNWs.     

Continuum Mechanics Modeling and Simulation of Carbon Nanotubes

The understanding of the mechanics of atomistic systems greatly benefits from continuum mechanics. One appealing approach aims at deductively constructing continuum theories starting from models of the interatomic interactions. This viewpoint has become extremely popular with the quasicontinuum method. The application of these ideas to carbon nanotubes presents a peculiarity with respect to usual crystalline materials: their structure relies on a two-dimensional curved lattice. This renders the cornerstone of crystal elasticity, the Cauchy–Born rule, insufficient to describe the effect of curvature. We discuss the application of a theory which corrects this deficiency to the mechanics of carbon nanotubes (CNTs). We review recent developments of this theory, which include the study of the convergence characteristics of the proposed continuum models to the parent atomistic models, as well as large scale simulations based on this theory. The latter have unveiled the complex nonlinear elastic response of thick multiwalled carbon nanotubes (MWCNTs), with an anomalous elastic regime following an almost absent harmonic range.