Morphological and nanomechanical changes in tungsten in high heat flux conditions

Morphological and nanomechanical alteration of tungsten in extreme environments, like those in edge localized modes in nuclear fusion environments, up to 46.3 GWm−2 heat fluxes were experimentally simulated using electrothermal plasma. Surface and subsurface damage to the tungsten is seen mainly in the form of pore formation, cracks, and resolidified melt instabilities. Mirco voids, rosette-type microfeatures, core-shell structure, particle enrichment, and submicron channels all manifest in the damaged subsurface. The formation of voids in the subsurface was determined to originate from the ductile fracture of hot tungsten by plastic flow but not developed to cracking. The voids were preferentially settled in grain boundaries, interfaces. The directionality of elongated voids and grains is biased to the heat flow vector or plasma pathway, which is the likely consequence of the thermally driven grain growth and sliding in the high-temperature conditions. The presence of a border between the transient layer and heat-affected zone is observed and attributed to plasma shock and thermal spallation of fractural tungsten at high temperature. Plasma peening-like hardening effects in tungsten were observed in the range of 22.7–46.3 GWm−2 but least in the case of the lowest heat flux, 12.5 GWm−2.


INTRODUCTION
Tungsten is a leading material for the divertor of tokamak-style fusion reactors due to its high melting point (~3422°C), good thermal conductivity (~175 Wm −1°C−1 at room temperature and 120 Wm −1°C−1 at 800°C), reasonably good mechanical strength (E Y~4 50 Gpa, nanohardness~8 GPa), high sputtering threshold energy (E th~2 00 eV for deuterium), resistance to hydride formation, and low tritium retention 1 . Despite these advantages, there are significant concerns as to the survivability of tungsten in a fusion environment. Damage to tungsten by hydrogen and helium ions, thermal shock from plasma exposure, and high heat effects are all important to understand for predicting the materials performance in nuclear fusion reactor environments [2][3][4] .
In particular, an extremely high heat flux environment known as an edge localized mode (ELM) 5 deliver 1-10 GWm −2 flux over submillisecond time frames for 10-100 sec, and is cause for concern for long term reactor operation. This amount of heat flux is as much as ten million times the intensity of solar radiation (~1 kWm −2 ) on the Earth's surface. A laboratory-scale method for physically simulating such high heat flux is by using a capillary discharging device. This device functions through discharging high voltage across a capillary, generates a pulse of electrothermal plasma with density of 10 23 -10 27 m −3 , temperature of 1-5 eV, particle velocities of~6 kms −1 , pressure of~300 MPa 6 . Previous studies with capillary discharging devices have reached heat flux ranges from 10 to~140 GWm −2 making them excellent tools for simulating ELM conditions 7 .
Depending upon the heat flux intensity, macroscopic melting can be observed in tungsten plates in fusion reactors 8 . Melting severely damages and loosens the material and can cause tungsten ejection and mass loss from the reactor wall surface 9,10 . Moreover, resolidifed melts deform the plasma contact surface thereby change the operational environment permanently. At the microscopic level, noticeable morphological, phasic, and structural changes emerge under ELM like conditions [11][12][13][14] . Resolidification of tungsten exposed to ELM condition has been studied 14 , but details of morphological changes in surface are still not fully understood. Generally, the resolidification process of plasma exposed molten metals has been widely studied in the field of welding 15 . This is the basis for understanding transient melts of tungsten. A high heat flux induced temperature gradient of welded joints is stratified roughly into three zones: the fusion zone, partially melted zone, and heat-affected zone. The fusion zone experiences full melting, whereas the partially melted zone undergoes partial phase change where solid and liquid material coexist briefly. As liquid phase smears into grains, intragranular liquid film separates grains, thereby localizing the tensile stress, which triggers cracking after cooling 16 . The heat-affected zone experiences high-temperature conditions but does not undergo melting. Despite the large body of research on welds, the mechanism of phase change related damage in plasma exposed tungsten is not fully understood in the nuclear fusion context. As fusion energy relies on the deuterium-tritium reaction, which produces helium particles, many studies have been conducted on helium effects in tungsten, particularly the formation of bubbles and nanostructures 17,18 . Nanostructures form in tungsten where helium energy is higher than 20 eV and temperature are in the 1273-2273°C range. Swelling is associated with the coalescence of helium into bubbles, which in turn responsible for embrittlement, and reduction in thermal conductivity. Therefore, mechanical properties such as hardness before and after helium implantation with high radiation damage has been extensively studied by others 19,20 . In the case of ELM conditions, helium plasma irradiation on tungsten showed that the nanostructures were affected by the high heat load 21 . Collective damage of ELMs and tungsten fuzz effects all lead to erosion 22 . However, limited numbers of studies have been done on the mechanical properties of tungsten in both high heat flux conditions and with helium particle implantation. Also, tungsten behavior in the extreme intensity of high heat fluxes of ELM condition is not well established.
The goal of this study is to understand the damaging mechanisms of material interactions at the tungsten surface in extremely high heat flux, plasma discharge environments, at low contact angle as shown in Fig. 1a. To investigate these mechanisms, morphological changes of the surface and subsurface of samples in Fig. 1b were observed by using focused ion beam (FIB) and scanning electron microscopes (SEM) and behaviors were interpreted based on stopping and range of interaction in matter (SRIM), and nanoindentation measurements. Figure 2 displays morphologies of the surface and subsurface of tungsten as a function of heat fluxes. The observed cracks in undamaged surface are likely introduced from the manufacturing process or during mechanical cutting. In the undamaged subsurface, grains of polycrystalline tungsten were revealed by polishing. The 12.5 GWm −2 damaged surface has multiple numbers of micron sized pores, located in the resolidified melts. A crosssectional view of 12.5 GWm −2 heat flux exposed sample has cluster of pores with a size of few microns to tens of microns. Some of the resolidified melts were near these pores. There is a recognizable 10-μm thick layer between the porous and the bulk side of the tungsten. Next, a panoramic view of the 46.3 Gm −2 damaged surface has tungsten matrix (bright) and copper features (dark) with the size of 50-100 μm implanted from the electrode of HELIOS. Cracks across the surface, micro-size pores, and resolidified melts are irregularly distributed. A cluster of pores on the surface were likely formed by prompt boiling during the moment of plasma interaction. This would imply that temperature was beyond the boiling point of 5930°C and was heterogeneously distributed from the non-uniform nature of the delivered heat flux and/or a fluctuating molten surface. In addition, wrinkles and flow of melts due to the high speed plasma flow can be seen in the resolidified surfaces. These characteristics show a fingerprint for plasma flow direction. This suggests a relatively slow cooling rate, and a viscosity high enough to retain fluidity during cooling.

Morphology
The cross section of damaged tungsten shows three layers; namely, a resolidified zone, an intermediate layer, and the bulk tungsten. Several tens of μm thick resolidified melts lie on the top of the damaged surface. It is hard to tell the uniformity of thickness due to the fluidic melt behavior. About a 10-μm thick transient layer was observed below the resolidified melts and is uniformly distributed over the subsurface. A cross-sectional view of bulk tungsten implies sporadic macroscopic cracks may be introduced during the plasma exposure or from the manual cutting process. A closer look at this 46.   layer, which is indicated by the red arrow. The portion that was not covered by the tungsten and copper in EDS maps is a carbon layer for FIB treatment. A more detailed SEM image is shown in Fig. 4 to clarify the transient layer of Fig. 3. The transient layer not only appears in Fig. 3 of the 46.3 GWm −2 exposure, but also in Fig. 5, meaning that this damage is not region specific. A transient layer is divided into an upper 5 μm or greater thickness region (layer 1), a lower 5-10 μm region (layer 2), and base tungsten (layer 3). Layer 1 contains a couple of μm-size voids, dispersed nano-tomicro-size particulates, irregular patterns of resolidified melts, intergranular nanochannels, and dendrite microfeatures of grains. Scattered nanospheres are copper enrichments according to EDS results. Copper enrichment is mostly located at the interface of two different grains throughout the entire matrix. In line with copper positioning, some of the copper atoms were condensed at the corner of voids. The shape of resolidified melts was like a river delta by making local grains as islands. The feature of these spotted phases is sometimes more like a rosette shape rather than a thinner dendrite. The nanochannels link previously liquified material in the dendrites to the long strand of resolidified melts. It is observable that the copper nanoparticles mostly dwell in that stream of resolidified melts, thus giving an indication for how far and fine melts penetrated in the matrix. Layer 2 is with voids at a triple junction. A vertical strand of sausage-like voids with copper particulates extended across the transient layer alongside the 1-2μm-sized cavity all describing the bigger size voids was preferred. Figure 5c shows an egg-shaped tungsten core copper shell microfeature about the size of 2 μm. The border of layer 2 and layer 3 has a bubbly line across the entire sample that clearly divides the transient layer and base tungsten. Layer 3 showed equiaxial grains and no visible voids within. Undamaged tungsten has a grain size of 2.28 ± 0.43 μm while 46.3 GWm −2 damaged sample showed 4.37 ± 0.45 μm, about a 90% increase in size. The average grain size of the transient layer between layer 1 and layer 2 is 0.85 ± 0.08 μm as can seen in Fig. 5b. The grain size has decreased by about 67.2, 80.5% compared to layer 3 and undamaged tungsten, respectively.

Void directionality
The angular distribution of voids in layer 2 of Figs. 4 and 5a and presented in Fig. 6. For the voids in Fig. 4, the results are fitted to Gaussian distribution with a broad range from 17°to 93°but with two different peaks. One sits at an average of 38.19°with a standard deviation of 9.14°and the other curve at an average of 77.09°with a standard deviation of 8.24°. A rough cutoff angle between two Gaussian distributions is at 59°, which divides 47 counts into a lower angle group for 27 voids whereas 20 voids grouped in a higher angle. In the case of Fig. 5a, 133 voids out of 220 voids were scattered with orientation from 4°to 169°c orresponding 60% chance of having directionality. The Gaussian pattern closely fitted to the angles from 4°to 93°where the average angle is 57.48°with a standard deviation of 19.2°. It is also observable that the multiple voids sitting at higher angles (>93°) still counted meaningfully. The average of two curves in Fig. 4 is 54.75°that closely matches with the average angle 57.48°in Fig. 5a. Table 1 listed measured nanoindentation and Young's modulus information. Figure 7 shows the decreasing pattern of nanohardness and fluctuating Young's modulus along the distance. Nanohardness decreased 38.8% for the 12.5 GWm −2 damaged sample and 24.0% for 46.3 GWm −2 damaged one from the transient layer to 1265 μm in depth. Likewise, Young's modulus decreased to 6.8, 3.3% for 12.5 GWm −2 , 46.3 GWm −2 damaged sample from the transient layer to 1265 μm, respectively. Most of the nanomechanical properties of samples and positions are showing a hardening effect relative to the undamaged sample. However, the 12.5 GWm −2 exposure's damaged spot at 1265 μm has nanoindentation 12.3% lower than the value of undamaged sample.

SRIM
Ion distribution of 5 keV helium and copper plasma in molten tungsten with a 1°incident angle is shown in Fig. 8. Figure 8a shows the ion range of helium and its farthest traveled length to      be 18.1 nm and~50 nm, respectively. Figure 8b shows the ion range of copper and its farthest penetration length to be 2.3 nm, 9 nm, respectively. This indicates that high-Z tungsten would easily hamper the particles. The level of radiation damage from bombarding plasma here would far exceed the melting limit of the lattice. Asymmetry of the ion distribution with respect to the normal vector of surface is due to the low angle bombardment. Since ions travel three order of magnitude less than the thickness of transient layer (~50 μm), initial kinematics cannot allow particles past the 50 μm depth.

DISCUSSION
A descriptive and inclusive damage process of tungsten by plasma and high heat fluxes in extreme ELMs (46.3 GWm −2 ) is illustrated in Fig. 9. The first interaction of tungsten with plasma is melting. A fluidic behavior of tungsten melts is governed by Raleigh Taylor and Kelvin Helmholtz instabilities such that waves of melts wobble on the viscous surface 23 . Raleigh Taylor instability is for normal directional forces acting on the liquid surface while Kelvin Helmholtz instability is for shear caused by high velocity plasma to the surface. These melt instabilities are undesirable in real operation since they alter the surface properties of tungsten significantly unless effectively healed. On top of that, the boiling of tungsten and copper is possible due to the high-temperature condition thus leading to surface instability. The formation of isolated copper regions is due to the immiscibility of copper and tungsten in any temperatures of compositional ranges based on the Cu-W phase diagram. Based on the Hume-Rothery rule, the ionic radius of copper is 1.57 Å while tungsten is 2.07 Å so that 22% of dimensional difference hampers the alloying. Also, the bcc structure of the tungsten and fcc structure of copper are unfavorable to form an intermetallic compound. Furthermore, 36.6 kJmol −1 of high positive enthalpy of mixing between tungsten makes insolubility of two metals 24 . Macroscopic cracking on the surface is due to the rapid solidification of melts as the cooling time is slower than the time of liquid relaxation thus tensile stress remained 25 . Unrelaxed structures remain and follow Griffith theory-based crack propagation at low temperatures where tungsten is brittle.
The subsurface damage is much more complicated as stratified interactions take place over the depth dimension before and after the damage process. Layer 1 can be considered as similar to the welding fusion zone and the melt pool hence having lowest viscosity compared to other layers during the damage process. As soon as the body cools down, layer 1 experiences accompanying nucleation and growth of new grains governed by shear stress, time during solidification, and cooling rate accounts for the formation of rosette and dendrite microfeatures 26 . The microstructural evolution at near melting point and quenching was previously investigated by MD simulation 27 . For example, partial     11 . With the temperature near at melting point and sweeping plasma, the transient layer possibly undergoes ductile fracture by plastic flow. First, voids nucleate at inclusions sites, such as impurities, and interfaces in tungsten, and continue to grow by plasticity. The nucleation of voids at 3073°C with shear stress obeys a simple relation of gas free void nucleation 28 : where r c is a critical radius, γ is the surface energy and σ is applied stress. As voids size exceeds than the critical radius, nucleation turns into growth. A continuously applied shear then triggers voids to coarsening and leads to crack formation. However, unlike the surface damage, subsurface does not show any major cracks. Despite the void size seen in the SEM images, the linkages between voids that potentially become a seed for the cracking are yet to found 29 . Copper-localized sites can act as void nucleation sites since they provide an interface as shown in Fig. 3. However, the heterogeneous void nucleation and growth were preferable since the temperature of subsurface dropped from the melting points with a relatively long duration that gives a chance for voids to migrate to nucleation sites. Because of the short stopping range of helium and copper in tungsten, their transport should have been assisted by fluidic transportation in nanochannels as copper halted in the middle of the transient layer. This indicates a gradual change of viscosity and permeability along the temperature gradient or depth.
Layer 2 is an intermediate layer with a higher solid fraction in the semi-solid slurry with relatively high viscosity. It is inappropriate to set an exact border between layer 1 and layer 2 because of the vagueness in the boundary. The deformation of dendrite microfeatures readily occurs by showing thixotropic behavior in the mush zone 30 . A mush zone is a feature found in castings of molten metal characterized by a formation of dendrite network in higher solid fraction at the dendrite coherence point. Figure 5a shows a similar but slightly wider transient layer compared to Fig. 4, which indicates the inequality of plasma damage from region to region. Melts from layer 1 were smeared into preexisting grain boundaries and grain interior 17 . Likewise, the columnar grains were preferentially grown to the melted zone in the welded tungsten 31 . Temperature gradient influenced grains in layer 2 were tilted but not layer 1 because of the higher viscosity. For example, the viscosity of tungsten temperature from 3000°C to the melting point is 10-5 mPas 32 . Melt instability also easily erases the patterns of directionality information whereas in layer 2 the heat source direction remains, meaning that the directionality of grains in a certain layer may belong to the flexibility or robustness of grains against the external thermal energy. Yet the duration of exposure was not enough to alter the grains to fully align with the heat source direction but the angle of voids sway at 57°. Also, the skewed pathway of plasma in the tungsten body as shown in Fig. 8 is responsible for the heat source direction. The angle varies due to multiple factors such as the variance of plasma intensity, deviation in the duration of interaction, and the morpho-chemical environment. Like layer 1, the origin of micro to submicron size voids in layer 2 is due to the shear and thermal stress-induced void formation but with smaller size, likely because the plasmainduced plastic flow is not as intense as the layer 1. It is unclear if helium was able to travel down to the voids in layer 2 but helium or possibly voids would prefer to settle down the elongated layers and interfaces 33 . Grain boundary sliding at an arbitrary boundary can occur due to the local stress concentration in layer 2. The triple junction of the grain boundary is the place where stress builds up by grain boundary sliding induced dislocation emission under high-stress loading 34 . The grain boundary gets weak as temperature goes up due to easier vacancy migration and a decrease of strength of the material. Figure 5c shows immiscible copper and tungsten can form an egg-like core-shell microfeatures in the liquid state according to the Marangoni effect 35 . As the minor liquid phase always forms the core of the egg, the enrichment of copper in a specific region is necessary. The tendency does harmonize with the copper chunk formation in Figs. 2 and 3. The formation of the core-shell structure is based on the surface energy difference between copper and tungsten phases. Yet the boiling temperature of copper (~2562°C) is lower than the melting point of tungsten, there will be at least 900°C temperature difference from the outer shell to the inner core. Although the tungsten droplet formation is less desirable than tungsten absorbing into the matrix, the matrix is no more mobile as it cools down to 2562°C. Inside the shell, the condition for the volumetric growth speed of the core was faster than the resolidification of copper. The Marangoni movement of droplet velocity in steady state is the following: where r is the radius of the droplet, μ W and μ Cu are the viscosities of the tungsten 32 and copper shell, respectively. σ is the interfacial energy, and x is the distance. The interfacial energy between two liquid phases must deal with the formation of the core microstructure. The interfacial energy between two liquid phases in the copper and tungsten binary system is simply proportional to the square of the concentration difference.
N 0 is Avogadro's number, N * is the number of front atoms per unit area of the interface, z 0 is the coordination number, z * is the number of cross bonds per front atom, L CuW is the interaction energy between copper and tungsten, and (C L2 Cu À C L1 Cu ) corresponds to the miscibility gap in the liquid phase. The driving force is the interfacial energy between the copper and tungsten liquid phases that is evaluated with the metastable miscibility gap in the copper tungsten binary system. The interfacial energy increases considerably with decreasing temperature because of the large temperature dependence of (C L2 Cu À C L1 Cu ). This suggests that if there is a temperature gradient in the droplet, an interfacial tension gradient should also exist. Microchannels are observed in Fig. 4 as a branch of liquid tungsten smeared into subjected grains. Again, copper was transported by the flow of liquid tungsten that is commensurate with particles delivered by microchannel analog to microfluidics of liquid metal 36 . A combination of tungsten and copper was suggested to use as a thermomechanically improved material. Evidently, nanohardness of tungsten in transient layer is greater than layer 3 and undamaged side as shown in Fig. 10 as hardening effect attributed to accumulation of residual stress, higher void density and clusters of copper nanoparticles.
A border between layer 2 and layer 3 distinguished the transient layer and heat-affected zone. This is shock and thermal induced spallation when laser driven shock wave can introduce a spallation in tungsten 37 . Shock-induced spallation is possible since HELIOS exerts pressures higher than 100 MPa to the tungsten surface during the plasma exposure. It is also known that yield strength of tungsten weakens at high-temperature ranges and will drop to 100 MPa at near melting point as 2300°C showed 140 Mpa 38 . Increasing strain and strain rate induces microfractures in tungsten 39 . Therefore, this border is the place where fractures happen from the exterior tensile stress by a shock wave near melting temperature. That is why SEM images showed a clear distinction between partially melted layer 2 and solid layer 3 in morphology. Thermal spallation might be another reason that gives a border because of the strong temperature and stress gradient along the depth axis 40 . A rapid temperature rises and drop afford to generate thermal stress as different parts of the tungsten expand and shrink. This effect gets severe as the temperature of tungsten drops to ductile to brittle transition temperature acting as like a fragile glass. Fusion relevant transient heat loads on tungsten would lead crack formation by thermal stress 41 . The thermal shock brings tensile and compressive stress along the way it propagates then spallation type voids were nucleated and enlarged after the shock.
The heat-affected zone experienced a high temperature but below the melting point. There are no more copper particles that exist in layer 3, so heat and pressure driven damage is the only factor. The recrystallization temperature of tungsten is around at the range of 1300-1500°C according to the supplier. Thus, the tungsten grain growth was affordable by comparing the size of the grains in undamaged subsurface of Fig. 2 and 46.3 GWm −2 damaged subsurface of Fig. 4 possibly accompanying recovery and recrystallization behind. Figure 10 shows a relationship between grain size and nanohardness for each sample but not all samples agree on classical Hall Petch relations: where σ Y is the yield stress, σ 0 is a constant for starting stress causing dislocation movement in tungsten, k Y is the strengthening coefficient, D is the average grain diameter. Normally yield stress has proportional relations to nanohardness. By comparing layers 1 and 2 to layer 3, nanohardness increased as grain size decreased. However, layer 3 and the heat-affected zone show reverse Hall Petch type relations, which points out that the grain size is not the only factor affecting hardness, rather residual stress such as dislocation density was affected. Although grain size of pristine tungsten is about 100 times larger than the critical grain size (~20 nm) of normally suggested reverse Hall Petch relations in metals 34 , nanohardness has increased probably because of the accumulated residual stress. As the plasma-induced shock wave propagates to the tungsten sample, residual stress has increased by the formation of dislocations. Simply, the residual stress generation and high-temperature annealing compete together counteractively along the depth dimension. Due to the overwhelming plasma shock-induced hardening over annealing effects, mechanical strength is higher in damaged samples than the undamaged sample. This is similar to the plasma peening process often used to enhance the hardness 42 . Compressive stress and hardening occur as long as the pressure plastically deform the structure 43 . A higher number of voids in layer 1 and layer 2 impede the propagation of dislocation causing hardening 28 . The dislocation density in the metal increased by laser shock for the large portion of the body 44 .
The intensity of hardness slightly decreases along the distance is due to the heating effects by the shock wave 45 and diffused over larger volume 46 . The hardening effects turns into softening as hightemperature annealing influenced more than shock peening on 12.5 GWm −2 damaged sample at 1265 μm meaning that plasma shock decays faster than the mitigation of shock damage along the depth dimension. Higher the heat flux, the more aggressive shock stimulates the more damage.
In conclusion, tungsten was exposed to high heat flux up to 46.3 GWm −2 in order to investigate ELM-relevant damage in the near surface. The morphological surface damages showed frozen melt instability with high porosity and macro-cracks that originated from the fast-cooled brittle tungsten. Below the surface, resolidified melts (layer 1), partially melted zone (layer 2), and heat-affected zone (layer 3) were found.
• Layer 1 contained rosette-type microfeatures that formed from shear stress, time during solidification, and cooling rate under melt instability. Nanochannels of liquid tungsten transported copper particles of short stopping range to deeper depth. Voids in this layer were formed by plasma and thermal stress assisted nucleation and growth under plasticity but did not evolve to the point of fracture.
• Layer 2 has columnar grains and void directionality aligned to the heat source with an average angle of 57°, determined by the robustness of semi-solid mainly viscosity. Void localization in grain boundaries, stress-induced sliding, and higher mobility by vacancy migration all weaken the material strength in high-temperature conditions. The hardening effect in layer 1 and layer 2 is from to small grain size, copper enrichment, higher void numbers, and residual stress. • The sharp border between layer 2 and layer 3 showed shock and thermal spallation took place on the fractured tungsten during the damage process.

•
In layer 3, grain size has increased compared to the undamaged material through by high-temperature grain growth. Plasma peening-like hardening overwhelmed thermal annealing for all depth dimensions except for the 12.5 GWm −2 damaged tungsten at the far edge region.

Sample preparations and plasma exposure experiments
The tungsten used in this work were supplied by Midwest Tungsten Service Inc. and machined by the supplier into hollow cylinders, or sleeves, with a 4 mm inner diameter, 6.9 mm outer diameter, and 110 mm in length. The products were manufactured by powder sintering under high temperature and pressure. The sample density is 99.6% theoretical density and compositional purity is 99.95% as shown in Table 2. The sleeve geometry is a requirement for samples installed in the electrothermal plasma generator used for this work, the heat exposure limit, ion outflow system (HELIOS). A schematic of this device is shown in Fig. 1a. Prior to discharge, a capacitor bank is charged to high energy. When discharged, high voltage jumps from the device cathode to anode and this current arc delivers energy to the walls and generates a plume of plasma from the dissociated wall materials. During the plasma exposure experiments, a copper electrode was used and discharges were initiated at room temperature. The device was evacuated and backfilled with helium gas to generate high energy helium. To achieve ELM like conditions, a single 1-2 µs pulse with peak voltages up to 5 kV and maximum currents of 30 kA were applied. Because the exposed surface area is small, high heat loads can be easily achieved. Four separate heat fluxes were investigated: 12.5 GWm −2 , 14.8 GWm −2 , 22.7 GWm −2 , and 46.3 GWm −2 . The plasma profiles as a function of time were calculated based on the electrothermal conditions in the previous study 7 . The samples were carefully unloaded after cooling at room temperature and sectioned into smaller pieces for analysis as shown in Fig. 1b. FIB techniques were employed for undamaged and 46.3 GWm −2 damaged samples in order to excavate the cross-sectional surface. The grains of undamaged tungsten in Fig. 2  Void directionality Forty-seven voids in Fig. 4 and 220 voids in Fig. 5a were manually identified and their angles respect to the bottom bubbly border were analyzed, see inset image of Fig. 6. The aspect ratio of voids needed to be at least higher than 2 so that the directionality is recognizable. Tolerance of measurement was 1-2°the slight tilt of the SEM images was accounted for, Figs. 4 and 5a. Nanoindentation Nanoindentation was performed on the polished surface in order to quantify the nanomechanical properties by using a Hysitron TI-900 TriboIndenter. The tip is a Berkovich-type indenter with a size of 200 nm. The samples were prepared in solid epoxy resin as shown in Fig. 1b so that the subjected surfaced is mechanically stable during the indentation. For the undamaged samples, 16 measurements were taken from a wide sampling area to get an average nanomechanical value. For 12.5, 22.7 and 46.3 GWm −2 damaged samples, indentations were conducted linearly from 0, 25, 115, 215, 815, and 1265 μm along the depth direction as shown in Fig. 7. The 0 μm site is simply a transient layer. These positions were verified by an optical microscope. For the 0 μm region, indentation was conducted four times to get 16 indentations with 10 μm spacing. For other bulk side positions, 16 indentations with 10 μm spacing were conducted, and obtained vales were averaged. Reduced Young's modulus was calculated: where E * is reduced Young's modulus, E is the modulus of a specimen, E' is the modulus of the indenter, and ν is Poisson's ratio of tungsten 0.28.

SRIM
Stopping and range of interaction in matter (SRIM) was used to understand copper and helium particle collisions with the tungsten matrix, 1 × 10 5 50 keV copper and helium ions were bombarded on the tungsten sample with low impact angle, 1°, relative to the impact surface. Displacement energy of tungsten was 68 eV and the collision mechanism follows the Kinchin-Pease model. The phase of tungsten would not matter for the simulation nonetheless the density of liquid tungsten is 17.6 gcm −3 .

DATA AVAILABILITY
The data that support the findings of this study are available from the corresponding author upon reasonable request.