Elastic Energy Storage in Al–Al₄C₃ Composites: Effects of Dislocation Character and Interfacial Graphite Formation

1. Introduction

Aluminum-based metal matrix composites (MMCs) have attracted significant attention in the aerospace, automotive, and defense industries due to their excellent combination of properties, including high specific strength, high specific stiffness, good thermal conductivity, wear resistance, and outstanding corrosion resistance compared to monolithic materials. These characteristics, together with their high strength-to-density ratio, good ductility, and manufacturability, make them ideal candidates for critical structural applications [1]. Among fabrication techniques, high-energy mechanical milling stands out as an effective method for producing MMCs through powder metallurgy. This solid-state process promotes a homogeneous distribution of reinforcements through repeated welding and fracturing cycles of particles, resulting in nanocrystalline microstructures that enhance mechanical strength [2]. Optimization of milling parameters, such as rotational speed, makes it possible to tune strength and ductility, obtaining ultrahigh-strength or ultraductile nanocomposites from the same composition [3]. Given these materials exposure to dynamic loads—such as impacts in aerospace structures or armor—it is essential to evaluate properties like fracture toughness and energy-absorption capacity [4]. Recent studies have shown that nanoreinforcements not only improve strength and stiffness, but also simultaneously increase energy absorption and the ability to dissipate vibrations—critical aspects for preventing brittle failures in service [4]. The reinforcement type and content are determining factors. For example, a 40% volume fraction of SiCp/2024Al exhibits ductile behavior, whereas a 50% volume fraction yields predominantly brittle behavior, evidencing a critical content for the transition between these regimes [5]. Likewise, the incorporation of nanometric reinforcements such as CNTs, graphene, or fullerenes improves tribological behavior and wear resistance, but excessive additions can cause agglomeration and loss of mechanical properties [6]. The fabrication technique also directly affects performance. In carbon-fiber–reinforced aluminum (CF/Al) composites, indirect extrusion has shown significant increases in compressive strength, fracture toughness, and fracture work compared with direct extrusion, due to better reinforcement distribution and intermediate alloy layers that facilitate energy absorption [7]. In hybrid systems—such as Al–Mg composites reinforced with ammonium polyphosphate (APP) via spark plasma sintering (SPS)—a 5 wt% APP content has been identified as optimal for compressive strength (432 MPa, +222% relative to the unreinforced matrix), whereas higher concentrations increase brittleness [8]. These results show that microstructural engineering, control of reinforcement content, and selection of the processing technique are critical variables for optimizing toughness and energy absorption in Al-based MMCs. The elastic energy density per unit volume stored (or elastic energy density) in a material (Ee​), mainly within the elastic stress fields generated by dislocations during plastic deformation, is given by:

$$E_e=A·G·b²·\rho ·ln(Re/r₀)$$

(1)

In this expression, G is the shear modulus (26 GPa for Al), b is the magnitude of the Burgers vector (which quantifies the atomic displacement associated with the dislocation and is 0.286 nm for aluminum), and ρ is the dislocation density, defined as the total dislocation line length per unit volume of the crystal. The terms Re (external cut-off radius) and r_o​ (internal or core cut-off radius) define the limits of the elastic field of the dislocation, where r_o​ is usually approximated as b, while Re​ can be associated with a characteristic distance such as the grain diameter in nanocrystalline materials or the average dislocation spacing [9]. The constant A varies depending on the type of dislocation: A=1/(4π) for a screw dislocation and A=1/(4π(1−ν)) for an edge dislocation [10]. This stored energy constitutes a fundamental parameter in dislocation physics and microstructural restoration processes, as its release acts as the driving force for phenomena such as recovery and recrystallization [11, 12]. The practical applications of elastic energy stored in the stress fields generated by dislocations span various areas of materials science. First, the quantification and characterization of deformed microstructures are achieved through the dislocation density (ρ), which allows the estimation of stored energy and the evaluation of microstructural evolution; in this context, X-ray diffraction methods such as complete line profile analysis (CMWP) are employed [13]. Likewise, in modeling strain hardening, the stored energy explains the increase in strength in metals and plays a central role in models such as that of Kocks–Mecking [14]. Similarly, in the design of nanocrystalline materials, the balance between dislocation density and stored energy decisively influences their mechanical properties [15]. The strain energy density (or energy per unit volume) that a material can store up to the yield point (also known as resilience) corresponds to the area under the stress–strain curve in the linear elastic region [16]. For proportional uniaxial loading within the proportional limit (just below the yield point), it can be expressed as:

$$E_s={\displaystyle \frac{1}{2}}·{\sigma }_p·{\epsilon }_p$$

(2)

In this expression, σ_p y ε_p represent the stress and the engineering strain within the proportional region of the stress–strain curve. The Es (J·m⁻³) quantity, quantifies the material’s ability to absorb and release energy without permanent damage; it is useful for comparing materials in scenarios involving impact, vibration, or transient loads.The mechanical performance of nanocrystalline materials is strongly conditioned by the nature of their interfaces. In particular, the contact zones between metallic and ceramic phases concentrate microstructural heterogeneities that generate local deformations, which in turn have a direct impact on the overall response of the material [17].To characterize this behavior, the Convolutional Multiple Whole Profile (CMWP) method has been established as a reliable approach, allowing the determination of key microstructural parameters such as dislocation density, dislocation character (edge, screw, or mixed), and the crystallite size distribution [13]. This procedure is based on a mathematical fitting of the full diffraction profile, using models that relate peak broadening to the presence of dislocations and grain size dispersion. Among the derived parameters, the factor q is fundamental for identifying the nature of dislocations by comparing experimental values with those theoretically calculated through contrast factors [18]. The theoretical value of parameter q is later compared with the experimental value obtained using the CMWP program. When the experimental q closely matches the value estimated for edge dislocations, these are classified accordingly. Conversely, if the experimental q approaches the characteristic value for screw dislocations, it is assigned that character, as described by T. Ungár and collaborators. Finally, if the experimental q value lies within an intermediate range between the edge and screw limits, the dislocations are considered to be of mixed character. Furthermore, the spatial arrangement of dislocations can be evaluated through the dimensionless parameter:

$$M=Re\sqrt{\rho }$$

(3)

which describes the degree of correlation among them: values of M > 1 indicate dislocations distributed in a random manner, whereas values of M < 1 reflect the formation of correlated configurations with short-range strain fields [19]. The effective external cut-off radius of dislocations (Rₑ) is a fundamental parameter in modeling diffraction profiles broadened by dislocations. In the Krivoglaz–Wilkens theory, Re​ is interpreted as the maximum distance over which the stress field associated with a dislocation retains physical significance before being compensated by other dislocations within the crystal [20]. This radius does not represent a real geometric limit, but rather an effective value that integrates the collective effect of long-range dislocation interactions. In practice, Re​ is determined by fitting X-ray diffraction profiles using CMWP. Its magnitude is directly related to the elastic energy stored in the lattice and to the material’s capacity to accumulate internal strains [13]. High Re​ values suggest that dislocation stress fields extend over longer distances, reflecting a less correlated arrangement of defects; in contrast, low values indicate that the fields are rapidly screened due to more ordered configurations with strong spatial correlation. The parameter q corresponds to the dislocation character (edge or screw). The theoretical q value for dislocations is taken from [19] and compared with the experimental q value obtained from XRD analyses. According to the methodology described by T. Ungár et al., the reference values are q = 1.31 for pure screw dislocations and q = 0.33 for pure edge. In this work the microstructural analysis focused on quantifying the parameters governing the elastic energy stored in the aluminum matrix. Using XRD–CMWP analysis, key dislocation-related parameters (ρ, q, and Re) were determined. Furthermore, the strain field distribution at the Al–Al4C3 nanorod interface, obtained from HRTEM-based strain mapping, was correlated with these microstructural parameters, thereby establishing a direct link between local interfacial deformation and the specific dislocation arrangements that control the material’s elastic energy storage capacity.

2.2. Materials and Methods

The reinforcement powder of Al₄C₃ was obtained through a reaction between Al and C, following the methodology described in a previous study [20]. The reinforcing powder (R) is composed of approximately 51 wt.% Al₄C₃ particles with an average size of 13 nm, about 3 wt.% of C phase with an average size of 20 nm, and the remainder consisting of Al. The Al-based composites were produced by mixing Al powder (99.5% purity, Sigma-Aldrich) with 1 and 2 wt.% of R. Each Al–R mixture was subjected to mechanical milling in a high-energy Simoloyer mill for 8 h, using argon as the milling atmosphere and approximately 4 ml of methanol as a process control agent. The device and milling media were made of hardened steel, establishing a ball-to-powder weight ratio of 50:1. Consolidated samples with a diameter of 6 mm and a length of 12 mm were obtained by uniaxial pressing of the powder mixtures for 2 min at approximately 1250 MPa. Subsequently, the samples were sintered at 550 °C for different time intervals: 2 h, 4 h, and 6 h. Table 1 describes the nomenclature, composition, and sintering times of the composites; the first number corresponds to the wt.% of the R powder mixture and the second to the sintering time. The composites were characterized by X-ray diffraction (XRD) and transmission electron microscopy (TEM). The diffraction profiles were recorded using a Philips X’pert powder diffractometer with a Cu cathode (λ = 0.15406 nm), a step size of 0.02°, and a counting time of 5 s per step. The diffraction profile analysis was applied to determine the crystallite size distribution and dislocation density of the nanocomposites using the CMWP fitting procedure. TEM characterization was performed using a JEOL JEM-2200FS transmission electron microscope equipped with a 200 kV field emission gun (FEG) and an energy-dispersive spectrometer (EDS). TEM samples were prepared using the focused ion beam (FIB) technique (JEM9320FIB). Compression tests were carried out using a 5-ton universal testing machine at room temperature and a strain rate of 1×10⁻³ s⁻¹ on cylindrical specimens of 6 mm diameter and 12 mm height, in accordance with ASTM E9-09 [21].

3. Results

3.1. Elastic Energy Density Analysis

Table 2 presents the values obtained for the samples containing 1, 2, and 3 wt.% of R, derived from the X-ray diffraction patterns using the CMWP technique: the average dislocation density (ρ), the experimental parameter (q), the average dislocation character (M), and the effective outer cut-off radius of dislocations (Rₑ). The table also includes the results of the elastic energy density (Es), obtained using Eq. (1).

Figure 1 presents the calculated Ee, derived from XRD data, alongside the dislocation density (ρ) and effective outer cut-off radius (Re), illustrating their variation with composite composition and sintering time.

3.2. Strain Energy Density Analysis

The graphs in Figure 2 show the behavior of the specimens under compression conditions at different reinforcement contents and sintering temperatures.

Table 3 shows the proportional stress (σ_p), proportional train (ε_p) and strain energy density (Es) obtained from the compression tests.

Figure 3a shows the strain energy density (Es) calculated using Eq. (2) and the compression test data as a function of composition and for different sintering times. For comparison purposes, the graph also includes the elastic energy density (Ee) values prior to the compression test, obtained from the XRD analyses. Graphic of Figure 3b, shows the energy absorption efficiency (EAE) as a function of composition and sintering time. Where EAE is expressed as:

$$EAE\left(\mathrm{\%}\right)={\displaystyle \frac{E_s}{E_e+E_s}}\mathrm{x}100$$

(4)

The graph also displays the q values, and it can be seen that, for most samples, the EAE exhibits a clear correspondence with the q values. As a first approximation, a mathematical model is proposed that considers only the proportional increment of Ee to estimate the value of Es. The calculated stored energy density (Ecs) is expressed as:

Ecs =Ee

(5)

Where α = 3.2159. Figure 3a also incorporates the curve corresponding to Ecs. The graphs show that Ecs curve does not exhibit a clear correlation with Es at several points; that is, Es does not increase proportionally with the initial Ee values. Notable discrepancies are observed in certain intervals; for example, at the points corresponding to samples Al-12, Al-22, Al-32, and Al-16, the model fails to adequately fit the experimental data.

Figure 4a shows the points corresponding to the difference ΔEsc = (Es - Esc) and compares them with the behavior of the q parameter values for each sample, under each processing condition.

3.3. Mathematical Model for Estimating Es

Figure 5 shows the comparison between the experimental ΔEsc curve and the fitted curve as a function of the parameter q.

This behavior suggests a parabolic-type model that explicitly incorporates the initial condition through the stored energy (Ee), and also accounts for the effect of the character of the pre-existing dislocations (parameterized by (q)) as a determining factor in the energy absorption capacity. The proposed model is:

$$Esq=\alpha Ee+{\beta }_1+{\beta }_{2}q+{\beta }_3{q}^2$$

(6)

Where α = 3.2159, β₁ = 3.2906, β₂ = –11.784, and β₃ = 7.6255. Figure 6 compares the Es curve with the Esqc curve calculated using Eq. (6). The graph shows a close correlation between both (R² ≈ 0.93). This result supports the feasibility of predicting the strain energy density (Es) in compression tests based on the initial stored energy, estimated through microstructural analysis by XRD.

A high-resolution transmission electron microscopy (HRTEM) image corresponding to the interface of an Al₄C₃ nanorod from the Al-34 sample is shown in Figure 7a. The inserted digital diffraction pattern, obtained through the fast Fourier transform (FFT) of region A, reveals the presence of two distinct zones: one corresponding to the Al phase, oriented along the [0 1 1] zone axis, and another associated with the Al₄C₃ phase, oriented along the [1 0 0] zone axis. The image also shows the region corresponding to the Al₄C₃ and Al phases, outlined by a dashed line. For the geometric phase analysis (GPA), two primary reciprocal lattice vectors of the Al₄C₃ phase, g₁ = 0 0 6 and g₂ = 0 1 6, were selected based on the digital diffraction pattern. The distribution of the strain field (εₓᵧ), influenced by the presence of dislocations at the interface, is shown in Figure 7b, highlighting the local strain variations induced by the lattice mismatch between both phases. Figure 8 shows an HRTEM image of the Al-36 sample, in which an Al₄C₃ nanorod can be observed. The image reveals the presence of a graphite phase located on both sides of the nanofiber, while the Al phase is observed in a more distant region. In the HRTEM image in Figure 9a, a region near the Al₄C₃ nanorod can be observed. The inserted digital diffraction pattern, obtained from the FFT of region A, reveals two distinct zones: one corresponding to the Al₄C₃ phase,, oriented along the [1 1 0] zone axis, and the other corresponding to the graphite phase, oriented along the [1 2 0] zone axis. For the GPA, two primary reciprocal lattice vectors of the Al₄C₃ phase were selected, g₁ = 0 0 3 and g₂ = 1 1 0, determined from the digital diffraction pattern. The strain field distribution (εₓᵧ), due to the presence of dislocations at the interface, is presented in Figure 9b.

A further insight into the local deformation mechanisms at the graphite/Al₄C₃ interface, the strain fields obtained from the geometric phase analysis are examined in detail. Figure 10a shows an amplified view of the (εxy) strain field extracted from the region highlighted in Figure 9, together with the corresponding (εyy) strain distribution (see Fig 9b) and in Figure 9c, the raw phase image, in which the presence of intrinsic stacking faults is clearly revealed.

TEM with Z-contrast imaging was employed to further characterize the composite microstructure, including the spatial distribution of the reinforcing phase and the graphite. Figure 11 presents a representative micrograph of the Al-32 sample, highlighting the arrangement of the Al₄C₃ nanorods (Figure 11a) and the graphite phase surrounding them, as revealed by the Z-contrast technique (Figure 11b). In the HAADF–STEM image, the bright halo surrounding the Al₄C₃ nanorod is consistent with the presence of graphite. As reported for graphitized carbon black, carbon layers locally oriented parallel to the electron beam increase the effective mass–thickness and therefore appear brighter [34].

3.4. Results Discussion

The elastic energy density Ee​ is mainly attributed to the microstructural parameters ρ, G, q, Re​, and b contained in Eq. (1), of which ρ, q, and Re​ were experimentally determined by XRD. In the case of the analyzed samples, it was observed that the dislocation density exerts a significant effect on the value of Ee​. The high dislocation density observed is attributed to the mechanical milling process, in which numerous linear defects are generated as a consequence of the dispersion of the reinforcement within the aluminum matrix. Additionally, the mismatch in physical properties between the matrix and the reinforcement, particularly in the coefficients of thermal expansion and elastic modulus, generates residual stresses during the sintering and cooling stages. These stresses are mainly relaxed through the nucleation and multiplication of dislocations at the matrix/reinforcement interface. As a result, a higher reinforcement content increases the number of active interfaces and, consequently, the dislocation density in the material. Figure 1 shows that larger fractions of reinforcement are associated with higher ρ values and, consequently, with increases in Ee​, in accordance with Eq. (1). For example, samples with 3 wt.% reach values between ~2 and 3 MJ/m³, which are approximately one order of magnitude higher than those corresponding to compositions with 1 wt.% and 2 wt.% (~ 0.24 and 0.55 MJ/m³, respectively). The lowest Ee​ value was verified for the Al-16 sample, sintered for 6 hours. This result is consistent, since the sample exhibits a very low dislocation density (~1.8 × 10¹⁴ m⁻²), and although the Re​ value is quite high (~102 nm), its contribution is marginal due to the logarithmic dependence of Ee​ on Re (Note: in the graph, Re is plotted as 30 nm for visual clarity). Conversely, the maximum Ee​ value corresponds to the Al-34 sample, resulting from its high dislocation density (~21.3 × 10¹⁴ m⁻²). Although Re is low (~6.5 nm), its effect on Ee​ is negligible compared to the dominant contribution of ρ. In fact, a clear inverse relationship is observed between Re​ and the dislocation density (ρ), where low Re values are systematically associated with high ρ values (see Figure 1). This correlation is particularly notable in the Al-32, Al-34, and Al-36 samples. Furthermore, a progressive reduction of the Re parameter is observed in all samples as the reinforcement fraction increases, for each of the sintering times evaluated (2, 4, and 6 h). Another factor influencing Ee​ is the parameter q, which quantifies the character of the dislocations present in the microstructure (edge/screw ratio). The results of the q parameter are presented in Table 2. The table shows that samples with q values between 1.29 and 1.58 are associated with a predominantly screw character, while the only sample exhibiting a pure edge character (q=0.33) is the Al-34 sample. The remaining samples, with q values ranging from 0.89 to 1.11, are classified as mixed character. The character of the dislocations is incorporated through the factor A in Eq. (1) and is presented for each q value in Table 2. When the dislocation configuration is dominated by screw-character segments, the factor A takes a value of approximately A≈0.079, which leads to a more pronounced reduction in Ee​ compared with microstructures in which edge dislocations predominate, for which A≈0.118. Under configurations where edge and screw dislocations coexist (A≈0.099), the value of Ee ​ falls within an intermediate range between the cases dominated by either pure character. The dislocation density (ρ) is a fundamental parameter that quantifies the degree of distortion in the crystal lattice. An increase in ρ leads to a proportional rise in the total dislocation line length per unit volume, which directly translates into an increase in the elastic energy density (Ee​) in the material. On the other hand, the effective range (Re​) of the stress fields associated with dislocations, determines the scale of long-range lattice distortion. A reduction in Re spatially limits the influence of these elastic fields, resulting in a decrease in Ee. However, when both the dislocation density and the dislocation character (q) are approximately similar, the effect of Re​ on Ee​ becomes practically negligible. This is evidenced by the Al-12 and Al-22 samples, where Re=20.02 nm and Re=14.5 nm, respectively, yet the corresponding Ee values remain very close (0.311 and 0.464 MJ m-3). Additionally, q parameter quantifies the influence of the dislocation character on the elastic energy storage capacity of the material. For instance, edge dislocations exhibit a higher ability to store elastic energy than screw dislocations, as their stress fields are more intense and localized [23]. This greater intensity promotes stronger interactions with microstructural obstacles (solute atoms, precipitates, grain boundaries), which increases local distortion and contributes to hardening and energy storage. In contrast, screw dislocations generate a more distributed pure shear field and can rearrange via cross-slip, so their contribution to energy storage is lower [24]. Overall, microstructures with a higher fraction of edge dislocations tend to exhibit greater accumulated elastic energy, whereas a predominance of screw dislocations reduces it [25]. The graphs in Figure 2 show the behavior of the specimens under compression. In the samples with 1 wt.% and 2 wt.% Al₄C₃, a progressive increase in stress is observed as the load on the specimen increases, without a clearly defined fracture point. This indicates a greater capacity for plastic deformation prior to failure. In contrast, the samples with 3 wt.% reinforcement exhibit a different behavior: the stress–strain curves display a well-defined fracture point, evidencing a sudden loss of load-bearing capacity. This result suggests that a higher reinforcement content promotes microstructural mechanisms that limit the ductility of the material. Table 3 presents the values of strain energy density, Es​, calculated using Eq. 2 and employing the proportional yield stress (σ_p) and proportional strain (ε_p) from the stress–strain curve. Figure 3 summarizes these results and, for comparison, includes the Ee​ values. At this stage, the deformation is essentially elastic, compression displaces the atoms from their equilibrium positions, storing potential energy that allows full recovery upon load removal. According to the results shown in Figure 3a, the energy Es​ is influenced by the initial elastic energy conditions (Ee​). That is, samples with a high reinforcement content (3 wt.% R) and a high Ee​ value (average of 2.5 MJ·m−3) also exhibit high Es​ values (average of 8.7 MJ·m−3). The results show that although the energy stored during compression, Es​, is related to the initial Ee​, its increase is not proportional to the initial Ee​ values. For example, although samples Al-22 and Al-12 showed similar Ee​ values (~0.48 MJ·m−3), sample Al-22 reached a significantly higher Es​, with 5.36 MJ·m−3 versus 3.7 MJ·m−3 for sample Al-12. Another case is sample Al-16 which, despite having a very low Ee​ value (0.24 MJ·m−3), achieved a considerably high Es, close to 5.32 MJ·m−3. To facilitate understanding of this behavior, the energy absorption efficiency (EAE) attained in the elastic region during the compression test was calculated for each sample. Figure 3b shows the EAE values for each sample under different processing conditions and also compares them with the q parameter values. The graph shows that the samples with the highest EAE correspond to Al-22 and Al-16, with approximate values of 92.0% and 95.7%, respectively. In contrast, samples Al-32, Al-24, Al-34, and Al-26 exhibit the lowest EAE values, around 70 - 77%. Therefore, directly applying a proportional model between Es​ and Ee​ (Eq. 5) does not adequately reproduce the experimentally observed behavior. The behavior of the calculated curve (Esc) according to Eq. 5 is shown in Figure 3a. This non-proportionality indicates the suitability of examining the microstructural factors that govern energy storage in the elastic regime for each sample in particular (e.g., dislocation character and density q, ρ and the effective outer cut-off radius Re, etc.). Figure 4a presents the difference curve between the experimentally measured stored energy (Es) and the calculated value (Esc​) from Eq. 5, defined as ΔEsc=(Es−Esc) for each sample. In Figure 4b, q ranges are included to characterize the microstructural state according to the predominant dislocation type: mainly screw dislocations (q≈1.3), mainly edge dislocations (q≈0.33), or a balanced mixture of both (q≈0.8). Comparing both plots, it is noted that samples Al-32, Al-14, Al-24, Al-34 and Al-26 exhibit the lowest absolute ΔEsc​ values, between 0.27 and 0.48 MJ·m−3. These same samples also exhibit low q values close to 1 (0.8–1.3), as shown in Figure 4b, which correspond to the aforementioned samples. Such q values are associated with a microstructural state dominated by dislocations of mixed character. Conversely, samples Al-12, Al-22, Al-16, and Al-36 exhibit a substantially elevated ΔEs, with values between ~1.9 and 4.45 MJ·m−3, together with markedly higher q values in the 1.51–1.56 range. This q range is associated with a microstructural condition characterized by a greater presence of screw dislocations. The plot also reveals a high correlation between the ΔEsc and q curves, suggesting that the observed difference may be governed by q parameter. Figure 5 shows the ΔEs curve and its fit as a function of q, where the data exhibit a clear parabolic trend. Figure 6 shows the Es​ curve and its comparison with the stored energy as a function of q (Esq), calculated using Eq. 6. The plot shows a close correlation between Es and Esq, with R²≈0.95. According to the results obtained, the parameter q in Eq. 1 (which quantifies the influence of dislocation character) is a determining factor for predicting the behavior observed in the composites. Although this parameter reflects the different energy-storage capacities between edge dislocations, which exhibit more intense and localized stress fields (greater energy-storage capacity), and screw dislocations, whose pure shear stress fields are more diffuse and contribute less to energy storage, the results show that, in this particular system, other microstructural mechanisms related to q predominate, in accordance with the proposed model (Eq. 6). The response may be related to the presence of dislocation dipoles and stacking faults. In FCC metals, dislocation motion (particularly that of screw dislocations undergoing cross-slip) produces jogs and generates dipoles formed by edge dislocations of opposite sign. These dipoles act as short-range elastic obstacles that reduce dislocation mobility and contribute to strain hardening [26]. A decrease in Re with increasing reinforcement indicates a stronger tendency for dislocation dipole formation. Low Re​ values imply rapid cancellation of strain fields, which limits dislocation motion and multiplication; consequently, the density of mobile dislocations decreases, resulting in reduced stored elastic energy [27]. This stored energy plays a critical role in redistributing and dissipating stresses during deformation; therefore, its reduction leads to strain localization and shortens the plastic stage prior to fracture. As a result, the total energy absorbed before failure decreases, leading to reduced material toughness. This behavior is observed in the 3 wt. % samples subjected to compression (see Figure 2c). In particular, the microstructural characterization of the Al-34 alloy shows a predominance of edge dislocations (see Table 2). Since dipoles are mainly formed by pure edge dislocations [28], it is reasonable to expect that this sample contains a high density of dislocation dipoles. This condition limits the ability of the material to store elastic energy, which is reflected in the EAE value of the Al-34 sample (~74%) the lowest among the specimens reinforced with 3 wt.% (see Figure 3b). Figure 7b, shows the HRTEM image of the Al-34 sample where is observable the strain field distribution (ε_xy), which arises primarily from dislocation dipoles distributed around an Al₄C₃ nanorod. On the other hand, the Al-36 sample exhibits an improvement in EAE of ~82.6%. According to the XRD analysis, the microstructure of this sample is associated with a higher dislocation density (~2.7×10¹⁵ m⁻²), with a predominantly screw character (see Table 2). According with the theory, higher fraction of screw dislocations reduces the formation of dislocation dipoles, which would otherwise hinder dislocation multiplication during plastic deformation. According to dislocation theory, a higher fraction of screw dislocations is generally associated with a reduced formation of dislocation dipoles. However, the observed decrease in Re with increasing reinforcement content indicates a stronger tendency for dipole formation. This sample exhibits the lowest Re value (~3.58 nm) and would therefore be expected to show enhanced dipole formation, which in turn hinders dislocation multiplication during plastic deformation. However, this same sample presents the highest dislocation density, suggesting that additional mechanisms are involved in the generation and accumulation of defects. A distinctive feature of the Al-36 sample is the presence of graphite at the Al–Al₄C₃ interface. HRTEM analysis confirms this, revealing graphite located at the interface of the Al₄C₃ nanorods (see Figure 8). The formation of this graphite phase can be attributed to diffusion processes associated with the prolonged sintering times. Therefore, the high dislocation density at the interface can reasonably be attributed primarily to the thermal mismatch between the graphite and the metallic matrix that develops during cooling.

The HRTEM image in Figure 9a shows a region near an Al₄C₃ nanorod; the phase present at the interface corresponds to graphite, whose observed planes are associated with the basal (002) plane. In turn, Figure 9b shows the distribution of the strain field (ε_xy), associated with the presence of well-defined dislocation dipoles, as well as dipoles that exhibit reduced strain fields. Figure 10a shows the distribution of strain fields (ε_xy) in a specific region of imagen of Figure 9b. In the image, arrows indicate the presence of strain concentrations that are barely perceptible in the (ε_yy) map shown in Figure 10b. The presence of these local strain concentrations, which are only faintly discernible in the ε_xy map, is consistent with the existence of extended dislocations, possibly associated with the dissociation of screw segments into partial dislocations separated by stacking faults. In the raw phase image (see Figure 10c), regions consistent with stacking faults are identified, and these correlate with the presence of such dislocations. These configurations may correspond to partial dislocations of mixed character (edge + screw), in which the local elastic field combines a dominant rotational shear component with weaker normal contributions. This behavior has been reported both experimentally and through atomistic simulations [27,29]. The presence of screw dislocations in systems containing graphite or in the vicinity of interfaces has been previously reported in the literature [30]. The dissociation of these dislocations into mixed-character partials, observed in the graphite phase, may account for both the high dislocation density and the reduction in the system’s energy, through the formation of stacking faults that relax the core field and redistribute the shear strain over a wider glide region [31,32]. This condition accounts for the improved energy storage efficiency and the increased toughness of the material in samples that exhibit a predominantly screw-type dislocation character. Previous studies have documented the addition of elements such as lithium to aluminum alloys with the aim of reducing the stacking fault energy (SFE) and promoting the emission of partial dislocations. This mechanism facilitates the nucleation of stacking faults and twinning planes, which subsequently activate additional plastic deformation modes, thereby enhancing ductility [33]. In Figure 11b, the graphite phase surrounding the Al₄C₃ nanorods is clearly revealed by the Z-contrast technique. This microstructural configuration, characterized by well-dispersed Al₄C₃ nanorods embedded in a graphite-containing matrix, is expected to play a key role in the mechanical response of the composite, particularly in terms of energy absorption and damage tolerance.

4. Conclusions

Compression tests were conducted to determine the strain energy density, Es in the Al/Al₄C₃ composites. These results were compared with the elastic energy density Ee​, calculated from the parameters obtained by XRD. The development of a mathematical model based on the dislocation character made it possible to establish that the difference observed between Es and the proportional increase in Ee associated with the elastic deformation of the Al matrix, depends on the distribution of stored energy through microstructural defects generated during the fabrication processes. These defects are linked to the dislocation character and are reflected in the parameter q. The samples with a dislocation density dominated by screw segments exhibited a more efficient redistribution of energy during processing, attributable to the dissociation into partial dislocations and the formation of stacking faults improving the EAE. In contrast, samples containing predominantly edge or mixed-character dislocations presented a higher proportion of dislocation dipoles, which generate highly localized strain fields reducing the EAE. The increased presence of screw dislocations in these samples is associated with the occurrence of graphite at the interface of the Al₄C₃ nanorods, a phase that develops as a result of prolonged sintering times. The formation of dislocations at this interface, as well as at the Al4C3/Al matrix interface, arises from the thermal mismatch between the graphite phase, Al₄C₃, and the aluminum matrix.