High temperature rise dominated cracking mechanisms in ultra-ductile and tough titanium alloy

Extensive use of titanium alloys is partly hindered by a lack of ductility, strain hardening, and fracture toughness. Recently, several β-metastable titanium alloys were designed to simultaneously activate both transformation-induced plasticity and twinning-induced plasticity effects, resulting in significant improvements to their strain hardening capacity and resistance to plastic localization. Here, we report an ultra-large fracture resistance in a Ti-12Mo alloy (wt.%), that results from a high resistance to damage nucleation, with an unexpected fracture phenomenology under quasi-static loading. Necking develops at a large uniform true strain of 0.3 while fracture initiates at a true fracture strain of 1.0 by intense through-thickness shear within a thin localized shear band. Transmission electron microscopy reveals that dynamic recrystallization occurs in this band, while local partial melting is observed on the fracture surface. Shear band temperatures of 1250–2450 °C are estimated by the fusible coating method. The reported high ductility combined to the unconventional fracture process opens alternative avenues toward Ti alloys toughening. Specific titanium alloys combine transformation-induced plasticity and twinning-induced plasticity for improved work hardening. Here, the authors show that these alloys also have an ultra-large fracture resistance and an unexpected fracture mechanism via dynamic recrystallization and local melting in a deformation band.

R ecently, a class of β-metastable Ti alloys exhibiting exceptional levels of strain hardening and uniform elongation resulting from the synergetic activation of TRIP (TRansformation-Induced Plasticity) and TWIP (TWinning-Induced Plasticity) effects has attracted major attention [1][2][3][4][5][6][7][8] . However, the post-necking and fracture properties of such alloys together with the underlying failure mechanisms have never been addressed in the literature, even though toughness and fracture resistance are critical properties for numerous structural applications.
The mechanical performances of structural materials can be defined by several properties. The yield stress, σ y , sets the resistance to the onset of plastic deformation. The uniform elongation, ϵ u , (expressed hereafter as true strain) quantifies the resistance to plastic localization and is directly related to the strain hardening capacity, with σ u being the corresponding stress and W n ¼ R ϵ u 0 σdϵ the work per unit volume dissipated until necking, sometimes called "toughness". The true fracture strain, ϵ f , characterizes the resistance to damage, with σ f being the true fracture stress and W f ¼ R ϵ f 0 σdϵ the work per unit volume dissipated until final fracture in the fracture zone. The stress intensity factor, K IC , the critical J integral, J IC , and/or corresponding critical crack tip opening displacement, δ c , characterize the resistance to crack initiation in a pre-cracked material, known as fracture toughness.
The development of alloys over the last decades has been often driven by the quest of increased strength levels σ y and/or σ u (by playing on hardening mechanisms) and, more recently, of enhanced strength-ductility balance σ u À ϵ u (by boosting the strain hardening capacity). However, microstructure engineering studies to tackle damage and fracture resistance are more rarely reported.
Depending on the stress state and on the mechanical properties and microstructure, plastic localization or damage accumulation is predominant and dictates the onset of fracture of the material 9 . Plastic localization usually dominates in materials with low strain hardening like pure metals, as cp-Ti 10 . In uniaxial tension, macroscopic plastic localization occurs first through the necking process when the Considère criterion is reached, i.e., when the cross section reduction is not compensated anymore by the strain hardening. The initially diffuse localization can then continue through a more localized shear banding process until the material fails by pure plastic flow without any occurrence of damage as in pure metals with no inclusions, or until damage initiates in the localized band as in pure metals containing non-metallic inclusions 11 . In materials with almost inexistent work hardening like nanocrystalline materials or bulk metallic glasses, shear banding is often activated very soon after the onset of plasticity, resulting in a catastrophic failure 12,13 . On the other hand, damage accumulation usually dominates failure (or occurs soon after the onset of necking) in the most recent advanced metallic alloys (e.g., high strength Ti alloys like β-CEZ, TRIP and TWIP steels, Q&P steels, …), where numerous interfaces are integrated within the microstructure to increase strength, but often at the expense of ductility [14][15][16][17][18] . The introduction of a large density of interfaces is in principle detrimental as it provides more opportunities for damage nucleation before plastic localization, resulting from the strain incompatibility at these interfaces. The accumulated porosity can partly counteract the strain-hardening capacity and accelerate plastic localization, or even directly lead to failure if void coalescence occurs before any plastic instability 9 . However, Wang et al. 19 recently postulated that the activation of the TWIP effect (mechanically induced nanotwins) could improve the fracture resistance of TRIP-TWIP steels by acting as a barrier against the growth of nanovoids nucleated at the interface of martensite transformed from austenite (TRIP effect). Therefore, a simultaneous activation of TRIP and TWIP effects in Ti alloys could possibly have the same positive influence on the resistance to damage. This is certainly essential with respect to numerous applications requiring high fracture toughness, not only for critical structural components but also in some forming processes limited by damage and cracking issues.
Here, we report the large damage resistance, as quantified by the work of fracture W f and confirmed by an estimation of the fracture toughness, of the generic Ti-12Mo (wt.%) alloy despite the huge density of interfaces originating from the TRIP and TWIP effects. In other words, this alloy not only involves a high resistance to plastic localization in terms of high ϵ u and high ϵ u σ u indicators as proved in earlier recent reports 1,2 , but also an exceptional high resistance to fracture that overtops classical titanium alloys. This is a unique combination of performances allowing both high formability, high crash resistance and high resistance to cracking. In addition, a fundamental study of the underlying failure mechanisms demonstrates large temperature rise during the quasi-static deformation of the investigated alloy with multiple local modifications of the microstructure. The resulting damage and cracking mechanisms are also found to be totally unusual for a polycrystalline alloy exhibiting a high workhardening rate.

Results
Tensile properties. Figure 1a presents the true stress-true strain tensile curves of the Ti-12Mo (wt.%) alloy (see in "Methods" for material information) compared with other classical Ti alloys, i.e., TA6V, Ti5553, and Ti LCB alloys 20 . While the yield strength of the Ti-12Mo alloy is lower, the true uniform strain is by far the largest since the macroscopic strain localization has been postponed owing to a high strain-hardening capacity 1,2 . Figure 1a also illustrates the magnitude of the properties at fracture. The true fracture stress (σ f ¼ F f =A f ) and true fracture strain (ϵ f ¼ ln A 0 =A f ) were computed from the last force value F f recorded before fracture, from the fracture area A f and from the initial section area of the tensile specimen A 0 . The fracture area was estimated from the smallest width at fracture, measured from the reconstituted fractured parts, and from the thickness of the neck, averaged over ten measurements along the width of the polished fractured parts. The tensile properties between the neck and the fracture point are extrapolated as a dashed line according to a power law relationship. The true fracture strain of the Ti-12Mo alloy is twice larger than the one of other reference Ti alloys with a comparable fracture stress. The association of such large levels of strength and ductility leads to exceptionally high levels of work of fracture W f , representing the total amount of energy to be brought in the fracture zone to activate ductile failure, approximated here by 1 2 ðσ y þ σ u Þϵ u þ 1 2 ðσ u þ σ f Þðϵ f À ϵ u Þ in order to allow the comparison with results from the literature. The work of fracture of Ti-12Mo, 865 MJ.m −3 , significantly overtops the one of other classical Ti alloys, i.e., 360 MJ.m −3 for Ti-5553, 440 MJ.m −3 for Ti-LCB, and 620 MJ.m −3 for TA6V. The properties of the present Ti-12Mo alloy are also compared with other high performance metallic materials in Fig. 1b in terms of combination of W f and ϵ u . This graph combines the two indicators that dictate the resistance to failure of ductile metals: W f for damage controlled fracture and ϵ u for plastic localization controlled fracture. The figure shows that, more than largely overpassing classical titanium alloys, the present Ti-12Mo alloy exhibits a combination of W f and ϵ u comparable to the currently best materials like specific TRIP steels 21,22 , high entropy alloys and TWIP steels [23][24][25] . The high ductility and failure resistance of Ti-12Mo, combined with the well-known interesting properties of Ti alloys (low density, high resistance to corrosion, biocompatibility), increases the range of potential applications of Ti alloys.
Moreover, a high level of work of fracture is a good indicator of anticipated large fracture toughness. The fracture toughness J IC is indeed defined as the work per unit area spent in the fracture process zone (FPZ). The work of separation in the FPZ is given by where u is the displacement at the boundaries of the FPZ, σ the opening stress and X 0 the thickness of the FPZ 26 . The fracture toughness is thus directly related to the work of fracture, J IC = X 0 W f . However, a limitation to this reasoning is the change of stress triaxiality between a uniaxial tensile test and the propagation of a crack. The work of fracture measured in uniaxial tension will be higher than the work of fracture evaluated at the tip of the crack, and the stress concentration at the tip of a crack could also trigger a change of fracture mechanism. Consequently, a better estimation of the fracture toughness is given below by the critical crack tip opening displacement δ c .
Fracture toughness. The critical crack tip opening displacement, δ c , was estimated from the propagation of pre-cracks in three double edge notched tensile specimens (DENT) under near plane stress conditions as δ c ¼ 146 ± 29 μm. Based on a Shih factor d n in plane stress equal to 0.4 27 , the fracture toughness of Ti-12Mo can be estimated as Even though the fracture toughness derived from δ c is only a rough estimation, the obtained value for Ti-12Mo is larger than the fracture toughness of classical Ti alloys such as TA6V(K IC $ 40À100 MPa:m 1 2 ) 28,29 .
Damage resistance. The origin of the exceptional ductility of the Ti-12Mo alloy is associated with its surprisingly large resistance to damage nucleation. While no cavities could be observed through the thickness of fractured specimens by X-ray computed tomography (see Supplementary Fig. 1), SEM micrographs of the fractured specimens revealed the presence of some cavities, 1-3 μm thick and 2-6 μm long, aligned as a micro-crack after the void coalescence step underneath the fracture surface (see Fig. 2a-d). The few cavities found inside the fractured specimen stand in line within localized deformation bands of same thickness (1-3 μm). The two shear offsets indicated by the red arrows in Fig. 2d show that these bands result from an intense shearing process. The nucleation of these cavities occurs just before fracture, as no cavity was detected in a tensile specimen deformed to a true strain of 1.03, and stopped just before fracture (see Supplementary  Fig. 2).
Transmission electron microscopy (TEM) analysis. A TEM specimen was machined perpendicularly to the band by focused ion beam (FIB) milling (as illustrated by the white rectangle in Fig. 2d) in order to investigate the microstructure changes occurring within these bands. The bright-field TEM micrograph of Fig. 2e as well as the selected area electron diffraction (SAED) patterns of Fig. 2f, g illustrate that the band presents a different microstructure from the surrounding. Out of this band, very small nanoblocks of about 30 nm can be observed, in agreement with the continuous diffraction rings shown in the corresponding SAED. Within the band, both TEM contrast and SAED show the presence of well-defined small equiaxed grains (118 ± 24 nm). Furthermore, the SAED patterns of Fig. 2f, g do not present exactly the same diffraction spots. Indeed, while no other spots than the ones corresponding to the β phase can be found in Fig. 2g, extra spots corresponding to the athermal ω (ω ath ) are present in Fig. 2f. The presence of nanoprecipitates of ω ath has also been confirmed by high-resolution TEM as shown in Fig. 2h. This means that nanoprecipitates of ω ath initially present in the β phase of quenched Ti-12Mo (wt.%) 1 were dissolved by the large level of plastic deformation, but appeared again in the band leading to fracture.
The evolution of the microstructure during the tensile deformation is thus schematically represented in Fig. 3. While the initial grains of Ti-12Mo are quite large (40 ± 9 μm) (see Fig. 3a), the high density of mechanically induced twins and martensite dynamically reduces the size of apparent grains 2 . Figure 3b presents the microstructure at a true strain of 0.05, where the size of the β untwinned blocks is already reduced to about 4 μm. TEM analyses carried out for small levels of deformation showed the formation of submicron "blocks" inside the twins as a result of the activation of secondary martensite transformation and twinning mechanisms 1,2,6 . At necking, nanoblocks about 50-100 nm in size can be observed as illustrated in Fig. 3c. The size of these nanoblocks is further decreased for a larger level of deformation near the fracture surface (where the deformation level is close to the true fracture strain (see Fig. 3d). A few localized shear bands eventually form Ti-LCB  and consist of larger equiaxed grains containing again ω ath nanoprecipitates, as already presented in Fig. 2.
Fractography. Figure 4 provides a multiscale characterization of the fracture surface, with an overall view in Fig. 4a. It mostly consists of flat dimples such as shown in Fig. 4b, with a diameter of 3.1 ± 0.1 μm and an average height of the walls of the order of half of the thickness of the localized shear bands (1.5 ± 0.3 μm). The observation of the walls of the dimples at high magnification in Fig. 4c reveals a granular appearance except at the top of the wall which exhibits a smooth, viscous-like appearance, resulting from melting and resolidification. On the other hand, Fig. 4d shows that a stretched zone is present at the edge of the fracture surface, resulting from the shearing of the specimen at 45°from the tensile axis and from the thickness direction (z axis). This zone is present along the entire width of the specimen (x direction). The length of this stretched surface, or shear offset, was estimated to be 27 ± 8 μm. Figure 4e shows that a specific fracture pattern can be observed at the transition between the stretched surface and the area covered with the dimples. This pattern is very similar to a so-called Taylor meniscus instability pattern 30,31 , which consists of several interconnected ridges caused by a fingerlike perturbation growth starting at the edges of the specimen. Such pattern is often found in metallic glasses [32][33][34] , polymers 35 , and liquid interfaces 36 . It was also hardly reported in the case of nanocrystalline/amorphous composite microstructure as the Taylor meniscus instability is caused by the penetration into a viscous layer of a less viscous material like air 31 , but it was never reported for crystalline metals. Consequently, the presence of such pattern on the fracture surface of the investigated alloy implies the development of a locally highly viscous state. A local crystalline-amorphous transition could induce the formation of Taylor patterns, however no amorphous transition was observed by high-resolution TEM underneath the fracture surface, as presented in Supplementary Fig. 3. A similar vein pattern was observed on the fracture surface of a TRIP β metastable Ti alloy 37 , mistakenly identified as shallow dimples, and on the fracture surface of a nanostructured Ti alloy failing by shear banding 38 , even though it was not for the same conditions of testing or the same level of deformation. It can therefore be anticipated that the original fracture mechanism of Ti-12Mo (wt.%) is not restricted to this alloy but is shared by more βÀ metastable Ti alloys. It can also be argued that the presence of Taylor patterns on the fracture surface finds its origin in a strong localization of the deformation in a shear band and a local thermal softening. A local remelting is also indicated by the presence of droplets on the Taylor patterns, pointed out by white arrows in Fig. 4e. More micrographs of the Taylor meniscus patterns are provided in Supplementary Fig. 4.
Fusible coating method. The fusible coating method originally proposed by Lewandowski et al. 39 was used to determine the temperature rise related to the fracture process. The frontier of the melting zone of a deposited coating provides a good indication of the temperature attained locally. A simple one-dimensional Fig. 2 Damage in localized shear bands. a After testing, a fractured tensile specimen was polished to a third of its thickness. b Secondary electrons SEM micrograph shows that there is almost no cavity in the fractured tensile specimen. c The only cavities observed are aligned in a micro-crack. d Backscattered electrons contrast shows that the cavities are formed within localized shear bands located underneath the fracture surface. The two red arrows indicate the shear offsets resulting from the shearing process. The white rectangle indicates the location where the TEM foil was extracted by FIB. e TEM bright-field micrograph revealing the microstructure inside the band (equiaxed grains) and outside the band (nanoblocks). The limits of the band are highlighted by dotted lines. f-g Selected area electron diffraction (SAED) patterns taken inside the band and outside the band, respectively. While ω spots are detected in f together with the β spots, the ω spots are not visible in g. h High-resolution TEM micrograph revealing the presence of ω phase inside the band (f). The HRTEM micrograph was taken along the <111> zone axis of the β phase. The measured distances of 0.42 and 0.28 nm within the precipitate correspond to the interplanar spacings of the (110) and (100) planes of the ω phase. differential equation for heat conduction (see Supplementary Equations (1a) and (1b)) was then used to estimate the distribution and time variation of temperature inside the fracture zone. Figure 5a shows that a 50-nm-thick Cu coating deposited on the tensile specimen prior to straining melted near the fracture surface. Two zones can be distinguished: the non-melted region stretched by the deformation of the tensile specimen (consequently divided into several "islands"); and the melted region that can be distinguished owing to its smoother appearance. The coating melted over a total width l of 1.1 ± 0.3 μm, i.e. along the tensile axis, as indicated in Fig. 5a. A correction must be applied to obtain the smallest distance x between the melted region and the fracture surface: x ¼ cos π 4 À Á l ¼ 0:8 ± 0:2μm. The melted region can itself be divided into two zones as observed in Fig. 5a: a first zone h ¼ 0:7 ± 0:1μm, where the specimen was deformed the most and which could correspond to the location of the shear band; and a second zone where the Cu coating presents a smaller wettability with the Ti alloy substrate, as the wettability of the melted coating decreases with the decrease of temperature.

Discussion
The outstanding ductility of generic β-metastable Ti-12Mo alloy (wt.%) is demonstrated both in terms of resistance to plastic localization and resistance to damage accumulation and cracking. As already developed in earlier works 1,2 , a large resistance to plastic localization was observed due to the high strain hardening caused by the simultaneous TRIP and TWIP effects, resulting in a large true uniform strain. When Considere necking condition is reached at a true strain of 0.3, the alloy can be considered as nanostructured. The creation of a high density of interfaces can indeed accelerate the process of nanostructuration as highlighted by Zafari et al. for TRIP Ti-5553 alloy, where a nanostructure is obtained after a true strain of 3, while an ultra-fined grain microstructure is obtained only after a true strain of 121 for Ti-20Mo (wt.%) alloy deforming by slip 40 . However, the work hardening in Ti-12Mo alloy (wt.%) at a true strain of 0.3 is still large enough to prevent catastrophic failure and a diffuse neck forms. Intensification of the plastic localization at a true strain of about 1.0 brings the formation of a few thin shear bands, in which the final fracture occurs. Even though the large work hardening of Ti-12Mo alloy has postponed this intense localization of the deformation to a large level of strain, it is not enough to ensure a large ductility at fracture as resistance to cracking in metals is controlled by the weakest resistance amongst plastic localization and resistance to damage. In a TRIP-TWIP steel, Wang et al. 19 proposed that the presence of nanotwins increases the damage resistance of the alloy due to the confinement of the nucleated nanovoids. In the present work, the resistance to damage is even more efficient as the damage nucleation itself is only activated at very large levels of deformation. The interfaces between the martensite, the β phase, and the twin boundaries do not constitute weak spots. Otherwise, multiple cavities would have nucleated earlier and more homogeneously. This extraordinary resistance to damage nucleation results in a fracture mechanism dominated by the plastic localization rather than the nucleation, growth and coalescence of cavities. Associated to a large work hardening, it leads to a large true fracture strain ϵ f , a large work of fracture W f and a very large estimated fracture toughness J IC .
The unusual flatness of the fracture surface can be explained by the development of damage inside a localized deformation band shearing through the entire section of the tensile specimen. The nucleated cavities are all aligned within the bands, and the final thickness of the dimples corresponds to half of the localized band thickness, suggesting that the entire damage mechanism, from nucleation to coalescence, occurs within these bands. Observation of Taylor instability patterns together with the smooth cap of the dimple walls and the presence of droplets on the fracture surface constitute concomitant evidences of a large increase of the temperature up to melting, thus of a viscous liquid state where fracture occurs. Moreover, a local large increase of temperature could trigger dynamic recrystallization and explain the presence of small equiaxed grains and nanoprecipitates of ω ath inside the shear band (see Fig. 2e-f).
A physically-based model has been developed in order to relate the measured width of the melted coating region with the temperature rise associated to plastic dissipation. As a quasi-static deformation rate was used during the deformation, no significant temperature rise could occur before the localization of the deformation within the shear bands. The heat rate, _ ω [kJ:m À3 :s À1 ], can therefore be estimated from the mechanical work dissipated during shearing at 45°from the tensile axis through the shear band: where β is the fraction of mechanical work dissipated into heat (taken here as 0.9); τ f ¼ σ f 2 (=640 MPa) is the true shear stress at fracture computed from the macroscopic fracture stress measurement; and _ γ ¼ v s 2h is the true shear strain rate, with v s being the shearing velocity and h is the half-thickness of the shear band (taken as 0.5 μm, the corrected width of the first zone as indicated in Fig. 5a). This heat rate _ ω is produced during the shearing time δt, which is the ratio between the shear offset measured on the fracture surface (27 μm) and the shearing speed v s . The only unknown parameter of the model that describes the temperature rise is therefore the shearing speed v s . This parameter is identified to bring the best correspondence with the width of the melted coating zone. Two boundaries can however be set for the shearing speed. The minimum shearing speed ν s is set by the applied straining rate v applied : v s > v applied cos π=4 $ 2:10 À5 m:s À1 . The maximum shearing speed on the other hand can be estimated as 10% of the maximum speed of sound in the material 41   temperature needed for dynamic recrystallization ð0:4T TiÀMo m Þ. The estimated range of shearing speed needed to melt the Cu coating over a distance of 0.6-1.0 μm is 55-300 m.s −1 , which corresponds to a maximum temperature reached at the fracture surface in the range of 1250-2450°C. This temperature range fits with the melting temperature of Ti-12Mo alloy, explaining the melted features and the Taylor meniscus pattern observed on the fracture surface. Moreover, the estimated temperature distribution also correlates quite well with the dynamic recrystallization occurring within the localized shear bands, as a temperature over 0.4 T m is predicted to be reached over a distance of 2 μm from the center of the shear band. The important temperature rise observed in the Ti-12Mo alloy is linked to its high shear cracking velocity and to its small thermal diffusivity: 200 cm 2 .s −1 for Ti-12Mo at room temperature, compared with 400 cm 2 .s −1 for 316 L stainless steel, 1300 cm 2 .s −1 for 1018 steel, or 6000 cm 2 .s −1 for 7xxx Al alloys. A small thermal diffusivity is also measured in other Ti alloys: 220 cm 2 .s −1 for Ti-6-6-2 or 290 cm 2 .s −1 for TA6V. However, in these alloys, the smaller damage resistance induces a damage-dominated ductile failure and prevents the unusual fracture sequence observed in the Ti-12Mo alloy. The same fracture mechanism involving local melting and dynamic recrystallization could nevertheless happen in many other highly ductile Ti alloys including the recently designed TRIP-TWIP Ti alloys [3][4][5]7,43 as implied by the similar fracture surface of a TRIP β-metastable Ti alloy 37 , opening a path for fracture properties improvements.
In conclusion, the Ti-12Mo alloy shows an outstanding capacity of dissipating mechanical work during deformation and fracture by combining the damage resistance and failure mechanisms of very ductile materials with the large strainhardening capacity of TRIP-TWIP alloys providing strength at large strains. The large ductility is due to an extraordinary high resistance to damage nucleation, bringing the full use of the large strain hardening capacity until the localization of the deformation at 45°from the tensile axis into few shear bands. A large temperature rise up to 1250-2450°C is estimated by the fusible coating method in the center of the shear band during its fracture, causing a local melting of the alloy at the fracture surface and dynamic recrystallization inside the shear band demonstrated by the re-precipitation of (ω ath ) in equiaxed β grains.

Methods
Alloy processing. Ingot of Ti-12 wt.% Mo was processed from commercially pure titanium and molybdenum by the self-consumable melting technique. The chemical composition after casting, measured by Inductively Coupled Plasma (ICP), gives Ti-11.8Mo-0.02Fe-0.01Cu-0.0115 O (in wt.%). A 10 mm thick slice of the ingot was homogenized at 900°C during 15 min, followed by water quenching. It was then cold rolled down to 1.1 mm, corresponding to a reduction level of 89%. A recrystallization treatment was carried out at 900°C for 15 min in high-purity Ar atmosphere, followed by water quenching to retain the fully β microstructure at room temperature. It is worth noting that nano-sized athermal ω ath phase is also present in the microstructure, as shown in previous work 1 .
Material characterization. Uniaxial tension was carried out at 1 mm.min −1 on 4 flat tensile specimens with a calibrated gauge length of 27 mm and gage section of 1.1 × 6.1 mm², and, at the same displacement rate, on 2 double edge notched tensile (DENT) specimens of dimensions 0.9 × 20 × 70 mm³, with 7 mm long and 300 μm thick edge notches on each side. Pre-cracks with an opening of <140 μm were produced at the tip of the notches of one DENT specimen with a razor blade and pre-cracks with an opening of less than 1 μm were produced by fatigue on the second DENT specimen. δ c was estimated from the difference between the opening of a blunted pre-crack, taken at 45°from the crack, and the opening of a new crack, for a limited new crack length of 400 μm, following Pardoen and Delannay 44 . The opening of the blunted crack and of the new crack were measured five times through the thickness by polishing two cracked DENT specimens. SEM was used to characterize the fracture surfaces and the melted Cu coating regions. A TEM sample was lifted from the localized band using a dual-beam focused ion beam (FIB) instrument, and was characterized using TEM operating at 200 kV. X-ray computed tomography was carried on four cylinders of 300 μm in diameter, and minimum 2 mm long cut from the fracture surface along the tensile direction, with a voxel size of 0.3 μm. The micrographs were analyzed by ImageJ Fiji 45 : the initial images were transformed into binary images by setting a threshold chosen to enhance the contrast between the matrix and the crack coming from the fracture surface, and the ImageJ algorithm "despeckle" was used to remove the noise. The height of the walls of the dimples was measured by white light interferometry, with a depth resolution of 89 nm. The density was measured by Archimedes' principle in water. The specific heat capacity was measured by differential scanning calorimetry under Ar atmosphere. The thermal diffusivity α was measured by the laser flash technique, in agreement with the ASTM E-1461 standard.
Cu coating. Tensile specimens were polished down to 1 μm diamond paste and OPS suspension. A 50-nm-thick coating Cu was deposited by Physical Vapor Deposition at a rate of 1 Å/s.

Data availability
The datasets generated during and/or analysed during the current study are available from the corresponding author on reasonable request.