## 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, Al_{2}O_{3}, 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, *N* atoms as a function of the positions ^{12} atoms (~ 10 micron size cubic samples) using hundreds of thousands of processors on one of the world-largest supercomputers [Germann & Kadau, 2008], the duration of most of the simulations in the area of materials research do not exceed tens of nanoseconds. The length-scales and the time scales of the MD simulations is largely defined by the computational efficiency of the interatomic potentials used to define the interatomic interactions.

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 *ab initio* simulations where experimental data is not available. Metals and metallic alloys have been studied using the embedded atom method (EAM) [Daw & Baskes, 1984; Foiles, 1985; Johnson, 1988] potentials, whereas covalently bonded systems have been studied using the Stillinger – Weber (SW) potential [Stillinger & Weber, 1985; Ding & Anderson, 1986] for Si and Ge, Tersoff potential [Tersoff, 1988a, 1988b] for Si and C, and the Brenner potential [Brenner, 1990] for hydrocarbon systems. Each of these potentials has a unique functional form and underlying physical concept that allows them to model distinguishing characteristics of the material system. For example, the EAM potential uses an embedding energy term that simulates the effect of an atom embedded in an electron density in metals. The SW potential uses an explicit angular term that stabilizes the diamond cubic lattice for Si and Ge. The Tersoff potential uses a bond order interaction (which includes an angular dependence) that defines the strength of a bond based on the surrounding environment for Si, Ge, and C. These unique functional forms and concepts limit the transferability of these potentials from metallic systems to covalent systems and vice versa and hence limits the atomistic studies to pure systems or alloys with similar bonding characteristics. The modified embedded atom method (MEAM) potential [Baskes, 1992] is an attractive potential to model mixed metal-covalent systems as it includes parameterization for many cubic metals, as well as Si, Ge, C, H, N, and O. However, the formulation uses a many-body angular screening function as a cutoff mechanism that makes the potential computationally expensive. In addition, the parameterization of the mixed cross interactions in mixed systems also requires an adjustment of the parameters of the pure components [Thijsse, 2005] that will require testing the performance of the potentials for the pure systems.

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 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

In the EAM potentials [Daw & Baskes, 1984; Foiles, 1985; Johnson, 1988], the energy of an atom is expressed as

where *r*_{ij} is the distance between atoms *i* and *j*, *i*, and the summation is over all atoms interacting with atom *i*. The electron density

where *j* at the location of atom *i*. Since only interatomic distances *r*_{ij} are needed to calculate the energy and forces in the system, the EAM calculations are nearly as simple and computationally efficient as the ones with pair potentials. The lack of explicit three-body terms, however, makes the conventional EAM inappropriate for covalently bonded materials.

To make the connections to SW or Tersoff potentials and to allow for the introduction of the angular dependence of the interatomic interactions, the linear sum of partial electron density contributions in Eq. (2) can be expressed through the sum of products of partial electron densities,

The sum on the right hand side of the above equation includes two-body terms with identical pairs of atoms (*j = k*) and three-body terms (*j ≠ k*) that can be separated from each other. The three-body terms can be written in the form of a sum over unique triplets of atoms (*i, j, k*) [Thijsse, 2002]:

where in the first (two-body) term under the square root, the summation is over all atoms interacting with atom *i*, and in the second (three-body) term under the square root, the summation is over all pairs of atoms *j* and *k* that form unique triplets with atom *i*. This formulation includes an explicit dependence on triplets of neighboring atoms and, as shown in Section 3.2, allows for incorporation of the angular dependence of the interatomic interactions in a form compatible with SW potential. Alternatively, the three-body terms can be written in a form of the “bond order” dependence:

This formulation can be used to design an angular-dependent EAM potential compatible with the Tersoff potential [Dongare et al., 2009b].

### 3.2 Reformulation of the SW potential

The energy of an atom in a system described by SW potential [Stillinger & Weber, 1985] is defined as

The potential consists of a two-body (*U*_{2}) and a three-body (*U*_{3}) term, with the summation of the three-body terms being taken over all atom pairs *j* and *k* that form unique triplets with atom *i*. The two-body term has a Lennard-Jones form terminated at a distance *r*_{c} by a cut-off function:

The three-body term is defined as

were *θ*_{jik} is an angle between vectors *i* and directed to atoms *j* and *k*. The parameters *A*, *B*, *p*, *q*, *λ*, ε, *σ*, *γ*, and *r*_{c} are adjustable parameters that are chosen to reproduce the properties of crystalline, liquid, and amorphous phases, as well as surface structures for Si. The sum of the three-body terms in Eq. (8) can be rewritten in a form of embedding energy,

where the electron density and partial electron density contributions are defined as

where *i* and *j* of different type [Dongare et al., 2009a]. This definition provides more flexibility in fitting the potential to the experimental data and the results of *ab initio* calculations, as discussed in Section 4, where an example of the parameterization procedure for cross-interactions is given for Al-Si system.

### 3.3 Angular-dependent EAM compatible with SW potential

The reformulation of the three-body term of the SW potential in the functional form of the EAM embedding function, Eq. (9) with the inclusion of the first (two-body) term under the square root in Eq. (4) can be used to wite a combined electron density as

Two parameters, *δ*_{i} and *c*_{i}, are added to ensure that the combined electron density given by Eq. (12) reduces to the conventional EAM potential for pure metals and to the original SW potential for pure Si. For metals, *δ*_{i} *= c*_{i} *= 0* excludes the angular dependence and includes the radial contributions to electron density, thus, yielding the original EAM potential. For Si, *δ*_{i} *= 1* and *c*_{i} *= 2* transform the electron density function into the three-body function of the SW potential, Eq. (10), that retains the angular dependence. The combined potential that reduces to the conventional SW and EAM potentials for pure components can be then formulated:

Here, the form of the two-body term of the SW potential, _{ij}(r_{ij}) = U_{2}(r_{ij}). As will be discussed in the next sections, the parameterization of the A-EAM potential for Al-Si interactions does not require an angular dependence of the electron density function for metal atoms. As a result, the computational efficiency of the evaluation of the A-EAM potential can be improved by rewriting the total electron density function in a form that does not involve the summation over triplets of neighboring atoms around a metal atom:

Similarly to Eq. (13), the parameters are set to δ_{i} = c_{i} = 0 for metals (the triplet contribution to electron density is eliminated) and δ_{i} = 1 and c_{i} = 2 for Si (the radial contribution to the electron density is eliminated). The additional parameter n_{i} is equal to 1 for metals and 1/2 for Si. While in this work we use the formulation of the A-EAM potential given by Eq. (13a), the alternative formulation of Eq. (13) is useful when an angular dependence needs to be included in the electron density of metal atoms in order to provide an adequate representation of the interatomic interactions in metal-Si alloys.

The functional form and parameters of the embedding energy functions,

## 4. A-EAM Potential for Si-particle reinforced Al-matrix composites

To test the ability of the combined A-EAM potential to reproduce the properties of systems with mixed metallic-covalent bonding, in this section we provide an example of parameterization of the potential fo Al-Si systems. The formulation of the A-EAM potential, given by Eq. (13a), implies that the presence of a Si atom contributes to the electron density of a neighboring metal atom only through the two-body term of the electron density function that does not have the angular dependence. On the other hand, the presence of a metal atom contributes to the electron density of a neighboring Si atom through the three-body (triplet interaction) term of the electron density function.

For Al-Si cross-interactions, the pair energy term is defined as:

where A_{C}, B_{C}, α_{C}, β_{C}, σ_{C}, R_{M}, P_{C}, and Q_{C} are fitting parameters. It should be noted that the partial electron density contributions

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. (14)-(15), 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, F_{ij} is the force on atom i due to atom j, M_{i} is the mass of atom i, and v_{i} 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 10^{9} 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 10^{9} 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 10^{9} 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 10^{9} 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.