1. Introduction
One of the long-standing problems in materials research has been the development and improvement of materials capable of withstanding ballistic/blast impact. These materials, often referred to as blast and penetration resistant materials (BPRMs), require the ability to mitigate damage and dissipate energy and momentum. These materials are designed to stop fragments and/or projectiles while preventing behind armor debris or spall. Functionally graded ceramic particle reinforced metal-matrix composites (MMCs) are an emerging class of materials that show significant promise for applications in protective structures due to their improved damage mitigation response in addition to their superior strength, stiffness and wear resistance. The improved damage mitigation is attributed to the compositional/structural gradients introduced through the variation of the concentration of the reinforcing ceramic particles (SiC, Al2O3, etc.) in the matrix [Mortensen & Suresh, 1995; Suresh & Mortensen, 1997]. These functionally graded materials show promise to reduce the need for layered structures that are susceptible to delamination and the formation of axial cracks initiated at the interface of the ceramic and the more ductile backing material. The mechanisms responsible for plastic deformation and failure of these metal-matrix composites, however, are complex and are affected by multiple factors, such as the distribution and location/size of the reinforcing particles, grain size in the metal-matrix, the structural characteristics of the interface between the particles and matrix grains, and loading conditions [Li & Ramesh, 1998, Li et al., 2007]. The design and optimization of BPRM nanocomposite materials can be significantly accelerated by improvements in the understanding of the deformation and failure mechanisms under conditions of ballistic impact. An optimized design of these materials, therefore, requires a fundamental understanding of the links between the length scales (particle size/distribution, layer thickness) and properties (structure, chemistry, etc.) of these composites that enable enabling a combination of high strength and good ductility, and improved failure resistance in impact environments.
While there has been significant progress toward understanding the behavior of these materials, there are still major knowledge gaps to be filled with respect to the identification and characterization of the elementary processes responsible for plastic deformation and failure. These knowledge gaps are attributed to the inability to fully integrate theory, models, and experiments because of disparate time scales and the heterogeneous nature of these materials. The time scales of the processes in these experiments can range from picoseconds to micro-seconds and the length scales can range from nanometers to micrometers [Remington et al., 2004]. These small time scales of the processes make it difficult to identify and characterize the elementary processes responsible for plastic deformation and failure using experiments alone. The molecular dynamics (MD) simulation technique has the ability to provide the atomic level structural information on the relevant micromechanisms at the time-scales of the spall experiments and may be instrumental in physical interpretation of experimental observations [Dongare et al., 2010a]. The applicability of MD simulations as an analysis/predictive tool, however, relies on the accuracy of the interatomic potentials to model the interactions between atoms in the system. The complex heterogeneity arises from the multicomponent interactions that results in complex bonding characteristics at the metal-ceramic interfaces. This heterogeniety limits the applicability of the MD technique due to inavailability of interatomic potentials capable of providing an accurate description of the structural and energy characteristics of metal-ceramic interfaces at the atomic scale. A new class of interatomic potentials referred as the angular dependent Embedded Atom Method (A-EAM) [Dongare et al., 2009a; Dongare et al., 2009b] has been recently developed that extends the ability of MD simulations to MMCs by combining the potentials conventionally used for metals with those conventionally used for ceramics in one functional form.
This chapter provides an overview of the current availability of interatomic potentials to investigate deformation and failure behavior in metal-matrix composites at the atomic scales using MD simulations in Section 2. An in-depth discussion on the formulation of interatomic potentials for complex metal-ceramic interactions is provided in Section 3 and an example parameterization is provided in Section 4. The applicabilty the A-EAM is demonstrated by investigating the strengthening and failure behavior of Si-particle reinforced Al-matrix nanocrystalline composites at high strain rates in Section 5.
2. Interatomic potentials for Molecular Dynamics simulations
The MD simulations study the evolution of the system as a function of time using classical atomistic models defined as interatomic potentials [Allen & Tildesley (1987)]. The interatomic potentials,
The interatomic potentials typically are analytical functions that define the energy of an atom based on the coordinates of its neighboring atoms. The parameters of the analytical functions are chosen to describe properties of the material system as observed experimentally (crystal structure, cohesive energy, density, elastic constants, phase transformations, etc.) or as predicted using
The applicability of MD simulations to investigate deformation and failure micromechanisms in ceramic particle reinforced MMCs requires development of new computationally efficient interatomic potenials that are able to accurately reproduce the structural and energy characteristics of the metal matrix, the ceramic phase, as well as the metal-ceramic interface. An easier solution would be to design new alloy potentials by combining the well-established and thoroughly tested potentials developed for pure components within a unified approach. A number of studies have focused on the relationships between various potentials to justify this unified approach. For example, a potential for the Pt-C system was developed with an analytical form that reduces to the bond-order Brenner potential for C and an EAM-like potential for Pt [Albe et al., 2002]. Similar connections between the EAM formalisms and the bond-order scheme of the Tersoff potential [Brenner, 1989] and between the SW and MEAM potentials [Thijsse, 2002] have been discussed. These studies suggest that despite the apparent distinctions between the functional forms and underlying physical arguments used in the description of interatomic bonding in metallic and covalent systems, a unified approach is feasible.
This chapter discusses the idea of developing a unified interatomic potential for systems with mixed type of bonding is extended to combine the EAM potential for metals with the SW potential commonly used in simulations of Si [Dongare et al., 2009a] and the Tersoff potential for Si and C [Dongare et al., 2009b]. The combined Angular-dependent EAM (A-EAM) potential incorporates a description of the angular dependence of the interatomic interactions into the framework of the EAM potential to make it compatible with the SW/Tersoff potentials. The A-EAM potential retains all the properties of the pure components as predicted by the original SW/Tersoff and EAM potentials and therefore eliminates the need for extensive testing for pure systems. In addition, the scope of the potential parameterization is limited to cross-interaction between the components. An example of the parameterization of the cross-interaction in the A-EAM potential is given for the Al-Si system here to study the micromechanisms related to deformation and failure of functionally graded Si reinforced Al-matrix composites (A359 alloy). Experimentally, the as-cast Al-Si (A359) alloy consists of a dendritic microstructure containing fine Si particles (wt % ~ 9 % Si) [Li et al., 1998, 2007]. The EAM potential for Al in the form suggested in Ref. [Mishin et al., 1999] is chosen as it provides a good description of the unstable and stable stacking fault energies and the surface and grain boundary energies in Al. In addition, the advantage of this potential is in the simplicity of the functional form and the availability of parameterizations for many metals, allowing for an easy implementation for a broad range of alloy systems. For pure Si, the original parameterization of SW potential [Stillinger & Weber, 1985] is used as it has been found to provide a good description of the crystalline, liquid, and amorphous phases of Si. The choice of the EAM and SW potentials renders the A-EAM potential to be computationally efficient and also provides an accurate description of the pure systems. The functional formulation of the A-EAM potential and the parameterization for the Al/Si cross-interactions are described in Section 3 below.
3. The Angular-dependent Embedded Atom Method (A-EAM)
A unified alloy potential based on a reformulation of the EAM and SW potentials in a compatible functional form [Dongare et al, 2009a] is discussed in this section. A reformulation of the conventional EAM potential into a form that includes three-body terms in the expression for the total electron density function is presented first, followed by the description of an approach for incorporation of the angular dependence compatible with the SW potential.
3.1 Reformulation of the EAM potential
The energy of an atom in the EAM formulation [Daw & Baskes, 1984; Foiles, 1985; Johnson, 1988], the energy of an atom is expressed as
where,
As discussed in [Dongare et al, 2009a] and [Dongare et al, 2009b], the ability to combine the EAM potentials with the SW or the TF potentials requires the addition of an angular dependence to the electron density. This chapter provides a more efficient formulation of the equations in [Dongare et al, 2009a]. The electron density term no longer requires a square-root dependence. Instead, the electron density for an atom in the A-EAM formulation is written as
For modeling pure metals, the parameter,
3.2 Reformulation of the SW potential
The SW formation for the energy of an [Stillinger & Weber, 1985] is given by
The two-body (
Here,
Here, the terms
Here, the parameter
3.3 Angular-dependent EAM compatible with SW potential
The A-EAM formulation [Dongare et al, 2009a] that combines both the EAM and SW formulations is given by:
Here the embedding energy,
The parameter,
Here, the form of the two-body term
4. A-EAM Potential for Si-particle reinforced Al-matrix composites
The A-EAM formulation given by equation (9) enables the modeling of metal-covalent interactions in alloys wherein the energy of an Al atom due to a neighboring Si atom is modified by only the pair term in the electron density fuction and a pair energy term for the Al-Si cross-interaction. Similarly, the energy of a Si atom due to a neighboring Al atom is modified by only through the angular term in the electron density function for the Si atom and the pair energy term for the Al-Si cross-interaction. The capability of the A-EAM potential is demonstrated using an example parameterization of the Al/Si system below.
The pair energy term and the partial electron density contribution for the Al-Si cross-interactions is chosen as:
where
Since there are no stable AlSi alloys (in the solid phase) observed experimentally, the parameters for the electron density function and the pair energy function for Al-Si cross-interactions, given by Eqs. (10)-(11), are selected based on the results of DFT calculations performed for several representative Al-Si bulk alloy structures. Since the number of the energy and structural parameters evaluated in DFT calculations exceeds the number of the fitting parameters, the fitting procedure is not aimed at reproducing the exact values of the material properties predicted in the DFT calculations. Rather, the results of the DFT calculations and the experimental data on the enthalpy of mixing of liquid Al-Si solutions [Kanibolotsky et al., 2002] are considered together as a target in the optimization of the overall agreement. Although several sets of parameters were obtained that provided a good agreement (less than 10 % deviation) with DFT results, the prediction of the enthalpy of mixing of the liquid alloy was the true test to study material behavior in extreme environments (high temperatures, pressures, and strain rates). The final parameter set chosen was found to provide an overall satisfactory description of the DFT results in addition to excellent prediction of the mixing behavior of the liquid alloy.
The plots of the embedding energy, the pair energy, and the partial electron density for Si (red lines) and Al (green lines) are shown in Fig. 1.(a), (b), and (c) respectively. The plots of the the pair energy functions and the partial electron density functions and are shown by the dashed blue lines in Fig. 1.(b), and (c), respectively. Thus, according to Eq. (2) and the plots in Fig. 1., the presence of Al atoms in the neighborhood of a Si atom causes a reduction in the partial electron density at the location of a Si atom and a reduced pair energy and therefore results in the weakening of the strength of the Si-Si interactions. This weakening
increases the tolerance of the local configuration in Si to deviations from the perfect diamond lattice structure in the presence of Al atoms. This observation is consistent with the results of recent DFT calculations which show a weakening of the Si-Si bonds at Al/Si interface [Nakayama et al., 2006]. The prediction of the experimental dependence of the enthalpy of mixing for the liquid Al-Si alloy at 1870 K by the A-EAM potential is shown in Fig. 1.(d) in comparison to that obtained using experiments [Kanibolotsky et al., 2002]. The plot shows excellent agreement of the predication using the A-EAM potential that shows a minimum of -5.36 kJ/mol at a concentration of 49 at. % Si with experimental data exhibiting a minimum of -5.4 kJ/mol at a composition of 50 at. % Si. The accurate prediction of the mixing behavior across the composition range and the overall prediction of structural characteristics and energetics from DFT calculations enables the A-EAM to predict the relevant intermixing/melting behavior at Al/Si interfaces in extreme environments such as shock loading and radiation damage. More details of the fitting procedure and a quantitative analysis of the performance of the A-EAM potential in comparison to DFT calculations is provided elsewhere [Dongare et al., 2011].
5. Deformation and failure behavior of MMCs at high strain rates
MD simulations with the newly developed A-EAM Al/Si potential are used to investigate the strengthening and failure behavior of functionally graded Al/Si composites. Experimentally, the as-cast Al-Si (A359) alloy consists of a dendritic microstructure containing fine Si particles (wt % ~ 9 % Si) [Li et al., 2007]. While it is beyond the capabilities of MD simulations to model length scales that incorporate the gradation in composition as observed experimentally, the response of the material at local regions with different compositions can be investigated individually using MD simulations. Nanocrystalline Al/Si systems with varying volume fractions of Si grains are used as the first approximations of the local microstructure in the functionally graded Al/Si alloy. The initial Al/Si nanocrystalline systems with an average grain size of 6 nm and containing 43 grains are created using the Voronoi construction method [Derlet & Van Swygenhoven 2003]. The distribution of the Si grains is chosen at random with random grain orientations depending on the desired composition of the nanocomposites. To study the effect of volume fraction, five Al/Si nanocrystalline systems are created with volume fractions of ~ 5 %, ~ 16 %, ~ 28 %, ~ 35 %, and ~40% Si. The volume fraction is defined as the percentage of Si grains in the nanocrystalline Al/Si system. To facilitate easy comparison, the nanocrystalline composites of different compositions are created with the same grain morphologies and the same grain orientations. Periodic boundary conditions are used in all the three directions. The as-created systems are first relaxed to have zero pressure, and then equilibrated at 300 K for 100 ps. The elements of the atomic-level stress tensor are calculated as
where α and β label the Cartesian components, Ω0 is the atomic volume, Fij is the force on atom i due to atom j, Mi is the mass of atom i, and vi is the velocity of atom i. The strain (ε) calculated here is the engineering strain and the mean stress (
where
The common neighbor analysis (CNA) [Honneycutt & Andersen, 1987 in addition to coordination number (CN) is used to identify local deformations in the FCC and diamond-cubic lattice. The initial nanocrystalline Al/Si systems with an average grain size of 6 nm are shown in Fig. 2. with the atoms colored according to the CNA/CN values. The contour for the Al atoms colored according to the CNA/CN characterization is as follows: the green atoms represent bulk fcc stacking, the red colored atoms represent local hexagonal close-packed order (stacking faults), the yellow atoms represent a coordination greater than 12, the light blue atoms represent a coordination of 12 other than fcc, and the dark blue atoms represent a coordination less than 12. The contour for the Si atoms colored according to the CNA/CN characterization is as follows: the purple atomsrepresent the tetrahedral bonding in the diamond cubic lattice, the white atomsrepresent coordination greater than 4, and the dark blue atoms represent a coordination of less than 4. The time step for all of the MD simulation runs was chosen to be 4 fs. The temperature was allowed to evolve during the deformation process.
5.1. Strengthening behavior of the Al/Si nanocrystalline composite
At the macro scale, the strengths of materials are typically analyzed using the commonly used phenomenological yield criteria such as von Mises and Tresca [Hertzberg, 1996[29]. While, these criteria have been found to be appropriate to study deformation behavior in metals under quasi-static loading conditions, recent MD simulations suggest that several effects result in asymmetries in the tensile and compressive strengths and render these criteria to be inappropriate under dynamic loading conditions [Dongare et al., 2010a[30] ; Dongare et al., 2010b. As a result, the. MD simulations of uniaxial stress loading were carried out to compute the flow stress of the composite at high strain rates. Deformation simulations were carried out at constant strain rate by deforming the sample in the loading direction while maintaining zero stress conditions in the lateral directions. The scaling parameter in the loading directions is chosen so as to achieve a strain rate of 109 s-1. For each nanocrystalline system, three deformation simulations are carried out independently for loading in the X, Y, and Z directions to eliminate the effect of loading direction on the orientation of the Si grains in the composite.
The stress-strain curves obtained for the composite with a volume fraction of ~ 35 % Si are plotted in Fig. 3.(a) and (b) for loading in tension and compression, respectively. The curves are initially linear and lie on top of each other up to the yield point, after which they start to deviate from elastic behavior. The flow stress (σf) is defined as the peak value of the stress in the stress-strain curve [Dongare et al., 2010b; Dongare et al., 2010c]. It can be seen from Fig. 3. that the nanocrystalline composite is stronger in compression than in tension. The flow stress for the nanocrystalline Al/Si composite for various compositions are tabulated in Table 1 for conditions of uniaxial tensile stress and uniaxial compressive stress loading at a constant strain rate of 109 s-1.
The snapshots for the nanocrystalline Al system at 10% strain are shown in Fig. 4.(a) for conditions of tensile loading and in Fig. 4.(b) for conditions of compressive loading at a strain of 10 %. Intermediate snapshots of composite containing ~ 35 % Si grains are shown in Fig. 4.(c) for conditions of tensile loading and in Fig. 4.(d) for conditions of compressive loading at a strain of 10 % with the atoms colored according to the CNA/CN characterization. The snapshots do not show significant dislocation activity in the nanocrystalline Al and the Al/Si composite due to the small grain size (6 nm) of the nanocrystalline systems. While no dislocations are observed in the Si grains, it can be seen that the presence of the reinforcing phase does not significantly alter the dislocation density of the Al grains. The higher strengths in the composite are therefore attributed to reduced grain boundary rotation/sliding at the Al/Si grain boundaries. The fraction of Al/Si grain boundaries increases with increasing fractions of Si grains and therefore results in higher strength of the nanocrystalline composite.
The average values of the flow stress in tension and compression are plotted in Fig. 5.(a) as a function of Si composition. While there are some deviations, the plot shows a general trend that the flow stress values increase with increasing Si fraction. The deviations in the trend are attributed to the small system size and the small number of grains. The asymmetry in the strength of the nanocomposites in tension and compression is plotted in Fig. 5.(b) for various Si compositions. The strength asymmetry for the pure nanocrystalline Al system is predicted to be very small (~ 1.6 %) and addition of Si results in an increase in the asymmetry as the Si composition in the nanocomposites increases. The higher asymmetry in the strength values at higher volume fractions results from the limited grain boundary rotation/sliding behavior of the Al/Si grain boundaries in compression as compared to tension as observed for pure nanocrystalline systems at higher strain rates [Dongare et al., 2010b; Dongare et al., 2010c]. These results are consistent with results from MD simulations using the MEAM potential which suggest that the inclusion of Si particles reduces the relative motion of the Al grains by suppressing both sliding/shearing along the grain boundaries and dislocation nucleation from the boundaries that results in enhanced strengths of the Al–Si nanocomposites [Ward et al., 2006].
5.2. Failure behavior of the Al/Si nanocrystalline composite
The failure behavior of nanocomposites formed by adding Si particles to Al nanocrystalline materials has been investigated using molecular dynamics simulations with the MEAM potential [Ward et al., 2006]. These results suggest that failure initiates by the nucleation of a void at the interface between the Al/Si grains. Experimental studies of deformation and failure of the A359 alloy at high strain rates also suggest that the fracture propagates along the interdendritic network of Si [Li et al., 2007]. A valid test of the A-EAM potential, therefore, is to study the micromechanisms related to ductile failure in nanocrystalline Al/Si composites consisting of an embedded Si nanocrystal in a nanocrystalline Al matrix. Deformation is carried out under conditions of uniaxial expansion (εx = εy = 0, and εz = ε) at a constant strain rate of 109 s-1 to relate to experimental dynamic loading conditions of dynamic failure [Dongare et al., 2009c]. Three deformation simulations are carried out independently for loading in the X, Y, and Z directions to eliminate the effect of loading direction on the orientation of the Si grains in the composite.
The stresses averaged over the entire system in the X, Y, and Z direction and the are plotted in Fig. 6.(a) as a function of engineering strain (ε) for uniaxial tensile strain loading in the Z direction. The evolution of mean stress (σm) as a function of engineering strain (ε) in the loading direction is plotted in Fig. 6.(b) as obtained in the three simulation runs during uniaxial tensile strain loading independently in the X, Y, and Z directions. The stress-strain curves in Fig. 6 show three stages of deformation. The first stage corresponds to elastic deformation of the nanocrystalline system leading to a linear evolution of the curves for σx, σy, and σz up to the yield point. Continued deformation results in deviation from linearity due to GB sliding/rotation processes that are accommodated by the nucleation of dislocations [Dongare et al., 2009c]. Continued deformation results in an increase in the stress to a peak value after which, nucleation of voids results in the relaxation of stresses as indicated by a sharp drop in the stresses in the X, Y, and Z directions as shown in Fig. 6.(a).
The stress state transitions to an almost hydrostatic stress condition
Continued deformation results in the stress state to transition back to that corresponding to uniaxial strain.
Intermediate snapshots of a section of the system are illustrated in Fig. 7. with the atoms colored according to CAN/CN contour. The snapshot in Fig. 7.(a) corresponds to a time just before the peak in the stress value (ε ~ 9.0 %). A single void is observed to nucleate at a grain boundary interface between the embedded Si nanocrystal and the neighboring Al grains in Fig. 7.(b). Also, there is relatively little dislocation activity at the GB interface between the Al/Si grains. This indicates that void nucleation cannot be attributed to a dislocation pile-up process, but rather results from the mechanical separation/sliding at the Al/Si grain boundary interface. The creation of the void is accommodated by the shearing of the nearby atoms leading to the formation of a disordered shell of atoms around the void as shown in Fig. 7.(c). Continued deformation increases the size of this shell of disordered atoms around the voids and the void shape changes to near spherical as shown in Fig. 7.(b)-(d). Although a few dislocations are observed to propagate into the grains surrounding the void, most of the material surrounding the void consists of disordered atoms and void growth is observed to occur by the shearing of the disordered regions. The creation of disordered atoms surrounding the void increases until a minima in stresses is reached in Fig. 6. This minima in stresses corresponds to the snapshot shown in Fig. 7.(e). A transition in the stress state back to that corresponding to uniaxial strain conditions results in a slower growth rate of the void that is accompanied by the recrystallization of the surrounding disordered regions. The recrystallization of the disordered region around the void is attributed to the increased temperature of the system due to nucleation and growth of the void [Dongare et al., 2009c]. The results from these MD simulations thus agree with those predicted by the MEAM potential [Ward et al., 2006] that the failure of the nanocomposite is initiated by the nucleation of the void at the Al/Si interface.
The value of the maximum mean tensile stress in Fig. 6.(b) can be compared to the peak spall pressure during spallation experiments. The plots suggest that the spall strength changes slightly for the nanocrystalline Al/Si composite with a change in loading direction. Intermediate snapshots of a section of the system at a strain of ε = 11 % during loading under conditions of uniaxial tensile strain in the X direction, the Y direction, and the Z direction are shown in Fig. 8.(a), (b), and (c), respectively at a constant strain rate of 109 s-1. The snapshots suggest that failure (nucleation of void) initiates along the Al/Si grain boundary interface. Thus, the embedding of a Si nanocrystal results in the initiation of failure at the Al/Si interface irrespective of the loading conditions. These results suggest that the failure resistance behavior of these metal-matrix composites is dependent on the characteristics of the metal-ceramic interfaces in the nanocomposites. More work is needed to be done to obtain the links between the length scales (particle size/distribution, layer thickness) and properties (structure, chemistry, etc.) of these composites for enabling high strength behavior and improved failure resistance so as to allow for an optimized microstructural design of the MMnCs for use in extreme environments.
6. Conclusion
A new angular-dependent EAM interatomic potential is developed that allows for the study of metal-ceramic interfaces at the atomic scale. The A-EAM potential is obtained by combining the EAM potential for Al with the SW potential for Si. The A-EAM potential retains the original functional forms and parameters and therefore all the properties of the pure components as predicted by the original SW and EAM potentials. The formulation of the A-EAM potential limits the parameterization to only the cross-interaction between the Al-Si components. An example parameterization is demonstrated for the Al/Si system based on the on the results of DFT calculations performed for artificial bulk Al/Si alloys and the compositional dependence of the enthalpy of mixing in liquid Al-Si alloys at 1870 K. Large scale molecular dynamics simulations using the A-EAM potential are carried out to study the strengthening behavior of the ceramic (Si) particle reinforced nanocrystalline metal (Al) matrix composites for varying volume fractions of the ceramic phase. The results suggest that the strength of the nanocrystalline composite increases with an increase in the volume fraction of the ceramic phase. The nanocrystalline composite is observed to be stronger in compression loading as compared to tensile loading and this strength asymmetry is observed to increase with increasing volume fraction of the ceramic phase. The microscopic mechanisms related dynamic failure predicted using MD simulations suggest that failure initiates at the metal-ceramic interface. These results suggest that the performance of the nanocomposites in extreme environments can be improved by a tailored design of a metal-ceramic interfaces in the nanocomposite. Future work will therefore focus on a fundamental understanding of the effect of the structure, size and distribution of the reinforcing ceramic phase that will enable high strengths of the nanocomposite in extreme environments.
Acknowledgments
This work has been done during AMD’s tenure as a National Research Council (NRC) Research Associate funded by the US Army Research Office.
References
- 1.
Albe K. Nordlund K. Averback R. S. 2002 Modeling the metal-semiconductor interaction: Analytical bond-order potential for platinum-carbon . ,65 19 195124 1098-0121 - 2.
Allen M. P. Tildesley D. J. 1987 Computer simulation of liquids , Clarendon Press,0-19855-645-4 - 3.
Baskes M. I. 1992 Modified embedded-atom potentials for cubic materials and impurities . ;46 5 2727 2742 ,0163-1829 - 4.
Brenner D. W. 1989 Relationship between the embedded-atom method and Tersoff potentials . ,63 9 1022 1022 ,0031-9007 - 5.
Brenner D. W. 1990 Empirical potential for hydrocarbons for use in simulating the chemical vapor deposition of diamond films. ,42 15 9458 9471 ,0163-1829 - 6.
Daw M. S. Baskes M. I. 1984 Embedded-atom method: Derivation and application to impurities, surfaces, and other defects in metals . . Rev. B,29 12 6443 6453 ,1098-0121 - 7.
Derlet P. M. Van Swygenhoven H. 2003 Atomic positional disorder in fcc metal nanocrystalline grain boundaries . ,67 1 014202 1098-0121 - 8.
Ding K. Andersen H. 1986 Molecular-dynamics simulation of amorphous germanium . Phys. Rev. B,34 10 6987 6991 ,0163-1829 - 9.
Dongare A. M. Neurock M. Zhigilei L. V. 2009 Angular-dependent embedded atom method potential for atomistic simulations of metal-covalent systems . . 80,18 184106 1098-0121 - 10.
Dongare A. M. Zhigilei L. V. Rajendran A. M. La Mattina B. 2009 Atomic scale modeling of multi-component interactions in nanocomposites, ,40 6 461 467 ,1359-8368 - 11.
Dongare A. M. Rajendran A. M. La Mattina B. Zikry M. A. Brenner D. W. 2009 Atomic scale simulations of ductile failure micromechanisms in nanocrystalline Cu at high strain rates . ,80 10 104108 1098-0121 - 12.
Dongare A. M. Rajendran A. M. La Mattina B. Zikry M. A. Brenner D. W. 2010 Atomic scale study of plastic-yield criterion in nanocrystalline metals using molecular dynamics simulations. ,41A 2 523 531 ,1073-5623 - 13.
Dongare A. M. Rajendran A. M. La Mattina B. Zikry M. A. Brenner D. W. 2010 Tension-Compression Asymmetry in Nanocrystalline Cu: High Strain Rate vs Quasi-Static Deformation. ,49 2 260 265 ,0927-0256 - 14.
Dongare A. M. Rajendran A. M. La Mattina B. Zikry M. A. Brenner D. W. 2010 Atomic scale studies of spall behavior in nanocrystalline Cu . ,108 11 113518 0021-8979 - 15.
Dongare A. M. Irving D. L. Rajendran A. M. La Mattina B. Zikry M. A. Brenner D. W. 2011 Physical Rev. B, In review. - 16.
Foiles S. M. 1985 Application of the embedded-atom method to liquid transition metals . . Rev. B,32 32 3409 3415 ,0163-1829 - 17.
Germann T. C. Kadau K. 2008 Trillion-atom molecular dynamics becomes a reality . ,19 9 1315 1319 ,0129-1831 - 18.
Hertzberg R. W. 1996 Deformation and fracture mechanics of engineering materials , 4th ed., Wiley,0-47101-214-9 York (NY), USA - 19.
Honneycutt D. J. Andersen H. C. 1987 Molecular dynamics study of melting and freezing of small Lennard-Jones clusters. ,91 19 4950 4963 ,0022-3654 - 20.
Johnson R. A. 1988 Analytic nearest-neighbor model for fcc metals . . Rev. B,37 37 3924 3931 ,0163-1829 - 21.
Kanibolotsky D. S. Bieloborodova O. A. Kotova N. V. Lisnyak V. V. 2002 Thermodynamic properties of liquid Al-Si and Al-Cu alloys . ,70 3 975 983 ,1418-2874 - 22.
Li Y. Ramesh K. T. 1998 Influence of particle volume fraction, shape, and aspect ratio on the behavior of particle-reinforced metal-matrix composites at high rates of strain . ,46 6 5633 5646 ,1359-6454 - 23.
Li Y. Ramesh K. T. Chin E. S. C. 2007 Plastic Deformation and failure in A359 aluminum and an A359-SiCp MMC under quasistatic and high-strain-rate tension .,41 1 27 40 ,0021-9983 - 24.
Mishin Y. Farkas D. Mehl M. J. Papaconstantopoulos D. A. 1999 Interatomic potentials for monoatomic metals from experimental data and ab initio calculations . ,59 5 3393 3407 ,1098-0121 - 25.
Mortensen A. Suresh S. 1995 Functionally graded metals and metal-ceramic composites: Part 1 Processing . ,40 6 239 265 ,0950-6608 - 26.
Nakayama T. Itaya S. Murayama D. 2006 Nano-scale view of atom intermixing at metal/semiconductor interfaces . ,38 216 219 , Maui, HI, Nov 27-Dec 02, 2005,1742-6588 - 27.
Remington B. A. Bazan G. Belak J. Bringa E. Caturla M. Colvin J. D. Edwards M. J. Glendinning S. G. Ivanov D. S. Kad B. Kalantar D. H. Kumar M. Lasinski B. F. Lorenz K. T. Mc Naney J. M. Meyerhofer D. D. Meyers M. A. Pollaine S. M. Rowley D. Schneider M. Stölken J. S. Wark J. S. Weber S. V. Wolfer W. G. Yaakobi B. Zhigilei L. V. 2004 Materials science under extreme conditions of pressure and strain rate.35 9 2587 2607 ,1073-5623 - 28.
Stillinger F. H. Weber T. A. 1985 Computer simulation of local order in condensed phases of silicon .31 5262 ,0163-1829 - 29.
Suresh S. Mortensen A. 1997 Functional graded metals and metal-ceramic composites: Part 2 Thermomechanical Behavior. ,42 3 85 116 ,0950-6608 - 30.
Tersoff J. 1988 New empirical approach for the structure and energy of covalent systems . ,37 12 6991 7000 ,0163-1829 - 31.
Tersoff J. 1988 Empirical interatomic potential for carbon, with applications to amorphous carbon . ,61 25 2879 2882 ,0031-9007 - 32.
Thijsse B. J. 2005 Silicon potentials under (ion) attack: Towards a new MEAM model . ;228 Sp. Iss. SI,198 211 ,0016-8583 X - 33.
Thijsse B. J. 2002 Relationship between the modified embedded-atom method and Stillinger-Weber potentials in calculating the structure of silicon . ,65 19 195207 1098-0121 - 34.
Ward D. K. Curtin W. A. Qi Y. 2006 Mechanical behavior of aluminum-silicon nanocomposites: A molecular dynamics study . ,54 17 4441 4451 ,1359-6454