Abstract
Porous materials have peculiar characteristics that are relevant for inertial confinement fusion (ICF). Among them, chemically produced foams are proved to be able to smooth the laser inhomogeneities and to increase the coupling of the laser with the target. Foams realized with other elements and techniques may prove useful as well for ICF applications. In this work, we explore the potential of a novel class of porous materials for ICF, namely, carbon nanofoams produced with the pulsed laser deposition (PLD) technique, by means of hydrodynamic numerical simulations. By comparison with a simulation of solid-density carbon, PLD nanofoams show a higher pressure at the shock front, which could make them potential good candidates as ablators for a capsule for direct-drive fusion.
1. Introduction
Porous materials, or foams, are the subject of an intense research activity in the context of inertial confinement fusion (ICF) and, more generally, of laser-matter interaction. Foams of low-Z elements have been proposed for use in ICF capsules, to generate bright X-ray sources [1], for the study of equations of state [2], and more recently for efficient electron acceleration [3].
In the context of ICF, foams are considered as constituents of the outer layer of the capsule, because of their ability to smooth the inhomogeneities in the laser energy [4–7], to enhance the laser absorption efficiency [8, 9], and to increase the ablation loading on a substrate [10]. These features make a foam ablator a potentially important choice to increase the laser-target coupling, thus transferring the laser energy into compression of the inner layers of the capsule.
Most of the experimental data available regard plastic foams, whose internal structure is constituted by features in the micrometric scale, such as filaments and membranes. Typically, these foams are produced by chemical methods. However, the investigation of the behavior of porous materials constituted by different elements and materials is needed to explore other potential applications. Recent studies [11, 12] indicate the advantages of using mid-Z ablators, which can reduce the impact of laser plasma instabilities (LPIs), therefore increasing the performance of the target. Among them, high-density carbon (HDC), typically used in indirect-drive experiments, can be beneficial for enhancing the performances of capsules, not only for the indirect-drive scheme but also for the direct-drive approach. In this context, low-density carbon nanofoams could potentially combine the benefits coming from the use of an ablator made of mid-Z elements and from the features of plasma behavior deriving from its random internal structure.
In this work, we present the first numerical investigation of the behavior of nanostructured carbon foams produced with the pulsed laser deposition (PLD) technique [13–15] under the action of a high-power nanosecond-long laser beam, at conditions relevant to ICF. PLD nanofoams have already been studied in the framework of high-intensity laser-matter interaction as a near-critical layer in double-layer targets [16] for ion acceleration with ultrafast lasers (pulse duration from fs to 100s fs) [17, 18].
In this context, a solid theoretical understanding of the physical processes granting an enhanced acceleration has been developed thanks to particle-in-cell kinetic simulations [19–21], and the optimal nanofoam parameters (namely, thickness and density) have been determined as a function of laser intensity [22]. However, there are significant differences in the laser parameters between typical ultrafast laser-driven acceleration and ICF experiments, especially in terms of pulse duration (typically from 30 fs to 100s fs for the former, few ns for the latter) and irradiance (typically – W/cm2 and – W/cm2, respectively), and dedicated studies are thus required. Particle-in-cell codes are not suited to simulate dynamical processes evolving on the longer timescales relevant for ICF (up to hundreds of nanoseconds).
Here, we use the hydrodynamic MULTI-FM code [23], validated with experiments on microstructured plastic foams [24], to simulate the behavior of the PLD carbon nanofoams. The need for a dedicated code has been confirmed in the last years since modeling the foam as an equivalent medium of the same density in the simulations proved to be unsatisfactory [6, 23]. The MULTI-FM code features a model for reproducing the laser absorption in the foam and the behavior of the laser-generated plasma, by taking into account the randomly arranged internal structure of the material.
We consider a set of different nanofoam parameters, which depend on the manufacturing process used and in particular on the time duration of the laser employed. From the simulations, we find that the shockwave in all the foams is significantly slower than in a homogeneous medium with the same density as the nanofoam, as expected from the studies with the plastic foams [24, 25]. Moreover, we find the indication that the pressure at the front of the shockwave is significantly larger than the one obtained in the equivalent homogeneous medium.
2. Methods
2.1. Materials
The carbon nanofoam presented in this work is produced by pulsed laser deposition (PLD), a physical vapour deposition technique that employs a pulsed laser to remove material from a target and deposit it on a substrate in a controlled gas atmosphere. If the background gas pressure is sufficiently high with respect to the energy of the ablated species, they can be slowed to a diffusive regime and aggregate while in flight, giving rise to porous, nanostructured, low-density films. In the case of carbon, what can be obtained are carbon nanofoams, nanoparticle-assembled films with a fractal-like morphology, whose basic constituents are nanoparticles with dimensions in the order of 10 nm and densities as low as a few mg/cm3.
Two different deposition systems are employed, one exploiting a nanosecond duration laser (ns-PLD, conventional technique) and the other an ultrashort femtosecond laser (fs-PLD). The first is a Q-switched Nd : YAG laser (second harmonic = 532 nm), with 5–7 ns pulse duration, 1 J maximum pulse energy, 10 Hz repetition rate; the second is a Ti : sapphire-based CPA laser (Coherent Astrella, = 800 nm), 100 fs pulse duration, 5 mJ maximum pulse energy, 1 kHz repetition rate. As shown in previous work [15], the different laser pulse duration leads to very different ablation regimes and enables better control of the foam properties. Notably, carbon nanofoams with the same density and different spatial uniformity (i.e., different microscale homogeneity) can be obtained.
The PLD target is pyrolytic graphite, the laser incidence angle is 45°, and the distance between target and substrate is fixed at 7 cm. We set the fluence to 360 mJ/cm2 for both techniques, obtained with 8 mm spot size (top hat) and 260 mJ for ns-PLD and 0.8 mm spot size (FWHM) and 2.6 mJ for fs-PLD. By suitably tuning the vacuum chamber pressure in the range from 50 to 200 Pa (argon is used as an inert gas), carbon nanofoams of different densities are obtained. Combining the effect of the pressure with the effect of the different techniques (ns-PLD and fs-PLD), four kinds of representative carbon nanofoam are produced, as shown in Table 1.
The nanofoam morphology is characterized with a field emission scanning electron microscope (SEM, Zeiss Supra 40) with 3–10 kV accelerating voltage. The microscope electron is also exploited to perform EDXS (energy-dispersive X-rays spectroscopy) on the foam samples and to obtain information about the elemental composition and average density. This is done through the EDDIE method [26]: a theoretical model of the electron transport in the film and substrate allows to retrieve the film mass thickness starting from the ratio between the respective X-rays line intensities. The nanofoam thickness is measured from the SEM cross-sectional images, and the density is obtained as their ratio.
The determination of the nanofoam pore size , one of the parameters needed for the foam model implemented in the MULTI-FM code, is not straightforward: contrary to most of the chemically produced foams, pulsed laser-deposited nanofoams are cluster-assembled and fractal, with disordered voids without a clear shape. One possibility is to estimate the average pore size with the characteristic uniformity scale length of the nanofoam, that is, the average distance between two adjacent high-density regions of the foam (or equivalently the distance between two less-dense regions, i.e., the void dimension). The EDDIE method can be used point-by-point on the samples to build mass thickness maps of the foams, and by performing a Fourier transform analysis on the maps, the nanofoam uniformity scale length is obtained. In principle, it would be possible to perform the Fourier transform analysis directly on the SEM images, but they would be affected by electronic effects and arbitrary postprocessing (i.e., contrast and brightness), while the mass thickness maps cleanly convey the useful physical information. The structure factor is the radius-averaged squared Fourier transform of the mass thickness image, and since the spatial frequency —corresponding to the maximum value of —is the predominant spatial frequency in the image, its inverse is the uniformity scale length [15], and that value is used as an estimate for .
Four distinct categories of carbon nanofoams are produced and investigated in this work, as reported in Table 1: sample Figures 1(a) and 1(b) have a density of 6 mg/cm3, which is the characteristic lowest carbon nanofoam density achievable with the PLD technique; sample Figures 1(c) and 1(d) are three times as dense (18 mg/cm3), close to the highest density conditions in which the material retains a porous, foam-like structure. For each density, the samples produced with ns-PLD (Figures 1(a) and 1(c)) present a pore size that is half of the value obtained for corresponding fs-PLD samples (Figures 1(b) and 1(d)). This can be appreciated from Figure 1: Figures 1(a) and 1(d) are similar in the distance between two crests of the microscale foam structure, and Figure 1(b) has a greater spacing while Figure 1(c) is significantly more uniform. Figure 2 shows the respective EDXS mass thickness maps, highlighting the same nanofoam microstructure. From the Fourier transform image in the insets, the ring at can be appreciated (note the wider scale for the Figure 1(c) inset).
Since carbon nanofoams grow through the random stacking of aggregates which are fractal in nature, it is possible to relate the density of the nanofoams to the characteristics of the fractal aggregates [15]:where and are, respectively, the foam and the solid element (nanoparticle) density and k is a proportionality factor related to the packing of the aggregates in the tridimensional foam structure. , , and are the fractal aggregates properties: nanoparticle diameter, gyration radius, and fractal dimension. Since all the carbon nanofoams considered in this work have a fractal dimension [15],
The foam model implemented in the MULTI-FM code is still based on a fractal description but expresses the over ratio as a function of different parameters: the thickness of the solid elements , the pore size , and the fractal parameter .
By comparing the two equations, one could seek a way to express each term of (3) as a function of the parameters of the other, namely, , , and , since the left-hand term is the same in both equations. Given the number of parameters and their possible combinations, a complete term-by-term correspondence is not univocal. A very straightforward relationship between and can be derived, as they are both a measure of the fractality of the system, and it is reasonable to assume they should be independent of the other parameters. Under this hypothesis, equating the exponents of the two equations leads to . Then, since the pore dimension is obtained from the Fourier analysis of the mass thickness maps, the only remaining parameter can be calculated by inverting (3): . A value of 2.0 g/cm3 is taken as solid element density , coherent with the disordered sp2 structure of the carbon nanoparticles composing the nanofoam [13, 14, 27]. The main parameters for the different carbon nanofoams considered in this study (Series A, B, C, and D) are shown in Table 1: their density, pore dimension , and solid element dimension , along with the respective production technique (ns- or fs-PLD). It can be noted that the value of the solid element thickness , in the order of some tens of nanometers, is coherent with the foam microstructure as seen in high magnification SEM images. Since the nanoparticle diameter is around 10 nm, the carbon nanofoam solid element can be identified with the strands of connected nanoparticles that make up the short-range structure of the nanofoams, which delimit the voids—of dimension —from one another. Therefore, the foam model implemented in MULTI-FM is representative of the nanofoam microstructure, while it is not sensitive to the short-range nanostructure (namely, the nanoparticles composing the foam solid element).
2.2. MULTI-FM Simulations
We simulated the behavior of the laser-generated plasma with the 1D MULTI-FM code [23]. The MULTI-FM code is a modification of the MULTI code [28] specifically developed for reproducing the interaction of powerful laser pulses with porous materials. It is based on an effective model for laser absorption accounting for the effect of solid elements and empty spaces on the plasma behavior. In the model, the pore size , the thickness of the solid elements , the average density of the foam , and the density of the solid parts are related as in equation (3). In the code, the homogenization status of the plasma is determined through a parameter called ISFOAM, defined for each numerical cell at each time step as
In this equation,where is the characteristic homogenization time given as follows:where is the effective ionization, is the mass number, and is the temperature of the plasma. In the last step, equation (3) was substituted in place of to better illustrate the dependence of the homogenization time on pore dimension and density. Due to the low density of the nanofoams, , therefore, . When the homogenization is reached, is
As the plasma is heated and homogenized, the parameter ISFOAM changes from 1 for the cold, not yet irradiated foam, to 0, when the plasma in the cell is completely homogenized. In the model, the thermal conduction is limited by the same parameter ISFOAM, to account for the inhomogeneous distribution of free electrons in the pores during homogenization as follows:where is the heat flux, is the Spitzer conductivity, and is the electron temperature. In the same fashion, the force acting on the mesh nodes is limited by the homogenization aswhere is the node velocity, is the total pressure, and is the Lagrangian coordinates used in the calculations in the code.
The MULTI-FM code has been validated with an experimental campaign carried out at the ABC laser facility [24].
Here, the laser parameters are taken as the ones of the ABC laser facility, in view of an experimental campaign to be conducted there in the near future. The laser pulse has a temporal profile with a full-width half maximum (FWHM) ns, at the fundamental wavelength nm, with an intensity on target of W/cm2. The parameters for the carbon nanofoam samples used in the simulations are the ones indicated in Table 1. In all the simulations, a tabulated equation of state from the SESAME library has been used. The carbon ions are considered to be fully ionized, . For all the simulations, the thickness of the sample has been taken as 300 μm, the largest possible to obtain by PLD in a reasonable production time.
3. Results
We report a series of simulations performed with the MULTI-FM code. The foam parameters have been chosen as the ones reported in Table 1. We will refer to each simulated nanofoam sample with the corresponding series name reported in that table. As explained in the Methods section, the laser wavelength in the simulations is taken as nm. The corresponding critical mass density is, therefore, 3.3 mg/cm3, which means that all the foams we are considering are overcritical.
Figure 3 shows the results of the simulations for Series A and B, having both the density of 6 mg/cm3 and the pore size of 10 and 20 μm, respectively. The shockwave in the nanofoam is slower than in the homogeneous medium. This behavior is well known and theoretically and experimentally documented in the literature for plastic foams [5, 6, 23–25, 29]. It is due to the time required for the laser-generated plasma to homogenize and fill the empty spaces in the foam. The lack of free electrons and the inhomogeneity of the plasma lead to a reduced heat conduction and hydrodynamic motion, until the plasma is homogenized. At that point, the plasma behaves as an ordinary homogeneous plasma.
The characteristic time of homogenization is affected by the pore size as in (6). Because of this dependence, the homogenization time is increased in the case of the nanofoams of Series B, compared to the ones of Series A, and the shockwave is slower in the former than in the latter.
The plots on the right column show that the pressure of the plasma attained in the foam is noticeably higher than that one obtained in the homogeneous medium of the same density. As analytically demonstrated in Ref. [7], an overcritical foam is more efficient than a solid homogeneous material in converting the laser energy into ablation loading. The lower thermal conduction and the inhibited hydrodynamic motion of the cold part of the foam lower the fraction of the laser energy absorbed into the plasma converted into compression or heating of the material. Therefore, a large fraction of this energy is left in the plasma as thermal energy, increasing the pressure in the foam. This property of the foam plasma can be relevant for ICF because it will increase the efficiency of compression on the first inner layer in the fusion target. This feature of foams will be beneficial for ICF, especially considering also the other properties of foam materials, the most important being the high absorption efficiency of the laser energy, which can be as high as 90%.
Figure 4 shows the results of a simulation of a solid-density carbon sample obtained with the MULTI code and with the same equation of state used for the foams, and the density is equal to 2.0 g/cm3. The pressure at the shock front in the nanofoam Series A exceeds the pressure obtained in the diamond sample. This effect is even more evident in the case of the nanofoam Series B, where on average it exceeds 20 Mbar.
Figure 5 shows the results of the simulations performed with the nanofoams Series C and D. The oscillations which can be seen in the plots of the density and pressure, in the plasma behind the shock front, are due to the interplay between the homogenization of the plasma and the propagation of the shockwave during the simulation. They reflect the inhomogeneous nature of the plasma realized inside the foam, where some inhomogeneities on large scales are known to remain long after the homogenization time [30]. In nanofoams Series C and D, the shockwave in the foam is slower than in the case of Figure 3. The pore size is smaller, leading to a reduction in the homogenization time, but, at the same time, the density is increased, with the opposite effect (see (6)). The homogenization time is thus mostly unchanged, and the slower shockwave propagation can be attributed to the higher density of the material. The pressure at the shock front in the nanofoam Series C exceeds 20 Mbar on average, an effect even more noticeable in Series D, where it exceeds 30 Mbar.
These results show how a foam can improve the ablation loading, by giving a higher pressure at the front of the shock wave compared with a homogeneous medium. The foams show in all cases a first narrow peak of the pressure at the shock front, followed by a thicker zone of pressure which has generally the same value attained in the homogeneous material. This behavior is related to the limiters in the model, which act on the thermal conductivity and on the compression of the material, reducing both effects (see (8) and (9) in Ref. [23]). The fraction of the plasma energy transformed into the motion of the fluid, and the one transported by heat conduction is reduced in comparison with a homogeneous medium, leading to an increase in the pressure at the front of the shock, which travels at a lower speed.
The simulations can also give an estimate of the most effective choice for the thickness of the sample, to maximize the pressure at the shock front when it exits the ablator. While for Series A, the highest pressure is attained at the end of the laser pulse, in Series B, the peak pressure is reached at 1.7 ns, before the peak power of the laser. By looking at Series C and D, we see that the peak pressure is always reached at about half of the laser time, i.e., when the laser is at its maximum power. Therefore, by properly choosing the thickness of a carbon nanofoam ablator, depending on the specific laser temporal profile, one can increase the pressure at the shock front and maximize the ablation loading, analogously to what has been observed in Ref. [7] for plastic foams. When considering carbon nanofoams of 18 mg/cm3, for both Series C and D, a thickness of the sample of about 150 μm would be the optimal to improve the compression efficiency from the ablator. Samples with this thickness can easily be obtained with the PLD technique. We will perform dedicated experiments on the ABC facility in the near future to characterize the laser-generated plasma with these nanofoams and to test the results obtained with the MULTI-FM code described in this work.
4. Conclusion
In this work, we reported a first study on the use of carbon nanofoams obtained with the pulsed laser deposition technique as ablators for ICF experiments. We considered four kinds of samples obtained through the use of nanosecond and femtosecond pulsed lasers, reported in Table 1. Taking advantage of the flexibility provided by the two different PLD techniques, two pairs of carbon nanofoams samples with the same average density but different microscale structures (and thus pore size) have been produced: samples A and B with a density of 6 mg/cm3 and C and D of 18 mg/cm3. Samples A and C, produced with ns-PLD, have smaller pore size, while larger pore size is obtained for samples B and D, produced with fs-PLD. The nanofoams are characterized in terms of morphology, and a connection between the fractal nanofoam model [15] and the foam description implemented in the MULTI-FM code is proposed. We investigated the behavior of these samples under irradiation of a high-power laser by numerical simulations performed with the MULTI-FM code. In particular, we considered the laser parameters as one of the ABC laser facilities at ENEA Centro Ricerche Frascati. The simulations showed that the shockwave into the nanofoam propagates sensibly slower than in a homogeneous medium of the same average density because of the internal structure of the foam. This effect is enhanced in the case of the nanofoams of Series C and D with an average density of 18 mg/cm3 due to the larger average density compared to Series A and B, having an average density of 6 mg/cm3. In all the nanofoams, the pressure at the front of the shockwave is larger than the one obtained in the equivalent homogeneous medium and also higher than the one obtained in a compact carbon sample, whose behavior has been simulated with the same code and equation of state. These first simulations on compact carbon, while limited in scope, are preparatory for future development toward the study of HDC ablators, and the implementation of the proper equation of state into the MULTI-FM code is foreseen. This effect can be of use for realizing an ablator layer for a fusion capsule, increasing the ablation loading on a subsequent layer, thus increasing the compression efficiency. This can be important for direct-drive inertial fusion with carbon-based ablators, such as HDC ablators. Theoretical works showed that mid-Z ablators have a better overall performance than plastic ablators, also reducing the LPI detrimental effects [12]. The use of carbon nanofoams could improve the performances of fusion targets, combining the benefits that can be obtained from a mid-Z ablator with the ones that come from the use of a porous material, such as the reduction of laser imprint and the smoothing of the laser energy on the capsule. A potential issue related to this finding could reside in the ability of the model implemented in the MULTI-FM code, which presently does not include the possible effects deriving from the interaction of the laser with the peculiar nanoscale structure of the material. Dedicated experiments are foreseen in the near future to confirm the results of this work, through an experimental campaign to be held at the ABC laser facility.
Data Availability
The data that support the findings of this study are available from the corresponding authors upon reasonable request.
Disclosure
Current address of V. Ciardiello is Universitá degli Studi del Sannio - Department of Engineering (DING), Piazza Roma, 21 - 82100 Benevento, Italy. Current address of A. Formenti is Accelerator Technology and Applied Physics Division, Lawrence Berkeley National Laboratory, Livermore, CA, USA. Views and opinions expressed are, however, those of the author(s) only and do not necessarily reflect those of the European Union or the European Commission. Neither the European Union nor the European Commission can be held responsible for them. The involved teams have operated within the framework of the Enabling Research Project: ENR-IFE.01.CEA “Advancing shock ignition for direct-drive inertial fusion.”
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work has been carried out within the framework of the EUROfusion Consortium, funded by the European Union via the Euratom Research and Training Programme (Grant Agreement no. 101052200—EUROfusion).