Ultra-intense laser interaction with nanostructured near-critical plasmas

Near-critical plasmas irradiated at ultra-high laser intensities (I > 1018W/cm2) allow to improve the performances of laser-driven particle and radiation sources and to explore scenarios of great astrophysical interest. Near-critical plasmas with controlled properties can be obtained with nanostructured low-density materials. By means of 3D Particle-In-Cell simulations, we investigate how realistic nanostructures influence the interaction of an ultra-intense laser with a plasma having a near-critical average electron density. We find that the presence of a nanostructure strongly reduces the effect of pulse polarization and enhances the energy absorbed by the ion population, while generally leading to a significant decrease of the electron temperature with respect to a homogeneous near-critical plasma. We also observe an effect of the nanostructure morphology. These results are relevant both for a fundamental understanding and for the foreseen applications of laser-plasma interaction in the near-critical regime.

is the plasma frequency (m e is the electron mass, e is the elementary charge). Relativistic transparency 33 allows a laser to propagate in a plasma beyond the classical transparency threshold marked by n c , as long as < + n n C a 1 e c 0 2 (C = 0.5, 1, respectively for linear and circular polarization, = a eA m c / e 0 0 2 is the normalized laser amplitude and A 0 is the absolute value of the peak vector potential of the pulse). For the typical wavelengths of existing ultra-intense laser facilities 34 (λ ~ 0.8 -1 μm), NCPs still represent a considerable challenge from the targetry point of view 35 : even a moderately ionized solid has an electron density two orders of magnitude greater than n c , while typical gas-jets provide electron densities well below n c . Near-critical electron densities can be obtained irradiating a dense cryogenic gas-jet 36 or a pre-exploded target 4 . However, these techniques require complicated setups and lead to NCPs with non-finely controllable density gradients and compositions. This might be detrimental for schemes requiring a careful control of the electron density profile, such as laser-driven ion acceleration 37 with targets consisting in a near-critical layer attached to a solid foil 8,[10][11][12][13] or γ-ray sources 23 based on low-density channels embedded in a denser medium.
Low density nanostructured materials, such as nanotubes 38 , nanowires 39 or foams [40][41][42] allow to obtain NCPs without further complications for the experimental setup, potentially providing at the same time an enhanced control over the density profile of the plasma. These materials, while being near-critical on average, consist in alternating voids and solid-density structures. The typical scale-length of these density fluctuations might well be comparable with the laser wavelength (i.e. of the order of 1 μm) 41 . Modern-day laser facilities can provide laser pulses with temporal contrast well below 10 −10 at the few ps timescale and of 10 6 ∼ − at the 1 ps timescale [43][44][45] . Assuming a pre-heating of 100s eV, the characteristic expansion velocity of the plasma is 100 ∼ nm/ps, hence in ∼1 ps the plasma has not enough time to fill the micrometric voids typical of low density nanostructures. Thus, a nanostructure might survive long enough to influence the interaction with the laser, especially at lower laser intensities, since the pre-pulses are less intense. This is confirmed by several experimental works carried out irradiating nanostructured targets at relativistic intensities 3,44,46,47 .
Despite the interest in NCPs and the growing literature regarding experiments carried out using nanostructured low-density materials 5,10-13,46-48 , the role played by the nanostructure is still poorly understood. Particle-In-Cell simulations 49 provide a powerful tool to investigate laser-matter interaction at relativistic intensities, even with nanostructured targets. Nonetheless, laser interaction with nanostructured foams has been systematically investigated, very recently, only with 2D simulations 50,51 , which do not reliably allow to simulate realistic nanostructures and to study the effect of pulse polarization. In a 3D geometry realistic foams have been exclusively investigated in a specific configuration, in a work 12,13 focused on ion acceleration. Few more papers include numerical simulations of low-density wire targets [46][47][48][52][53][54][55] , though not in a near-critical regime and often in a 2D geometry.
In this work, by means of an extensive 2D and 3D Particle-In-Cell simulation campaign, we investigate how realistic nanostructures can influence laser-plasma interaction in the near-critical regime. We consider ultra-short laser pulses with intensities in the range a 0 = 5-45, which span the capabilities of current ultra-intense lasers, from table-top sub-100 TeraWatt systems up to state-of-the-art PetaWatt facilities. The specific case of a 0 = 5 is particularly interesting since it corresponds to a relatively moderate intensity I 5 10 W/cm 19 2 ∼ ⋅ , achievable in many facilities worldwide. Laser interaction with NCPs in such intensity regime could be appealing due to the increased absorption efficiency provided by the NCP, which could allow to observe or to enhance physical phenomena that would otherwise require higher laser intensities and, thus, larger and more expensive facilities. For this reason, for the lowest laser intensity case, we also investigate the effect of different nanostructure morphologies, considering two structures of wide experimental interest in the laser-plasma community and in materials science: foams 5,10,12,13 and nanowires [46][47][48]52,53,55,56 . The relevance of the investigation of how these kinds of nanostructured materials behave in extreme electromagnetic fields goes beyond the generation of NCPs, being of interest for both fundamental science 57 and specific applications 47,48 .
All the simulations have been carried out with an average density of 3 n c , which is both experimentally feasible with low density materials 41 and well within the near-critical regime for the relativistic intensities considered here. Similarly to what was done by Robinson et al. 30 , in this work we define a NCP as a plasma whose electron density n e satisfies . < < + n n n C a 0 1 1 , as in our work, the upper limit for relativistic transparency can be reformulated as S < 1, where S = n e /a 0 n c is the relativistic similarity parameter 58 .

Results
We simulate the interaction of an ultra-short (∼30 fs), relativistic (a 0 = 5-45), Ti:Sapphire (λ = 0.8 μm) laser pulse with the targets shown in Fig. 1. We consider foam targets made of nanospheres (radius of 40 nm) arranged in space according to two different models based on Diffusion Limited Aggregation 59 (see Methods section for details). Qualitatively these models provide fractal-like structures: an aggregate of random clusters of nanoparticles (DLCCA foam) and a tree-like morphology (DLA foam). Both these structures are representative of realistic foam-like plasmas (see Fig. 1(c)). As far as nanowires (radius of 96 nm) targets are of concern, both an ordered array and a random assembly of wires are simulated. In all these cases the average electron density is 3n c and the average filling factor is ~5%. Thus we explore scenarios ranging from near the transparency threshold (S = 0.6 for linear polarization) up to a moderately transparent regime (S = 0.067 for linear polarization).
We investigate the effect of the nanostructure morphologies on several quantities of experimental relevance, such as absorption efficiency, energy spectra and angular distribution of the emitted energetic particles.
We first discuss the effect of a foam-like nanostructure on the laser-matter coupling, comparing, at different laser intensities, some relevant quantities with those of a benchmark homogeneous plasma. In the second part we investigate different plasma nanostructures at fixed irradiation conditions. Lastly, we discuss the possible effects of the simulations dimensionality on the laser plasma-interaction with both homogeneous and nanostructured targets, which is of particular relevance for the numerical research in this field.
Effect of nanostructure on laser interaction with NCPs as a function of a 0 . In this section we investigate the main differences in laser-plasma coupling with either homogeneous and nanostructured plasmas. A first qualitative view of the interaction is given in Fig. 2(a and b), which show that the pulse propagation occurs differently in homogeneous and nanostructured plasmas. The nanostructured plasma allows for a stronger laser energy absorption and a much higher conversion efficiency of laser energy into ion kinetic energy. As shown in Fig. 2(c), the ion energy absorption remains well below 5% in all the homogeneous plasma cases, while it is remarkably enhanced with foams, up to 20% for a 0 = 5. This is a direct consequence of the explosion of the nanospheres constituting the nanostructured plasmas, whose signature can be found in the ion phase space, as shown in Fig. 2(f). Both these effects are more relevant at the lowest intensity. The enhancement of laser absorption at a 0 = 5 can be explained considering that the plasma is close to the transparency threshold. In the case of a uniform plasma, a denser layer is rapidly generated at the pulse front, hindering further propagation and increasing reflection. Conversely, a nanostructure lets the laser propagate through its voids deeper in the target, enhancing absorption. The higher absorption efficiency of laser energy into ion kinetic energy can be explained considering that laser energy is directly transferred to the electrons of the plasma. Indeed, due to their higher inertia, ions acquire energy only through the electric fields generated by charge separation. The local plasma density for the DLCCA foam target is ∼20 times higher than for the uniform plasma, so that significantly more intense electric fields are generated when electrons are blown away by the laser. Figure 2(d) shows that, at the lowest intensity, the nanostructured plasma leads to a doubled electron energy absorption if compared to the homogeneous one. This absorption enhancement decreases at a 0 = 15 and becomes negligible at a 0 = 45 (the curves of energy absorption over time are almost superimposed in this case). This is due to a faster homogenization of the nanostructures at higher laser intensities. Simulations show that at the lowest intensity, the nanostructures are still clearly distinguishable after 20 laser cycles, while at the highest intensity only a weak signature of the original structure is visible. In order to have a better insight of the electron heating we also analyzed the electron energy spectra (shown in Fig. 2(e)). As a general feature, the foam plasma allows for a more efficient acceleration of mildly energetic electrons (i.e. with a kinetic energy , where E cut-off is the cut-off energy). On the other hand, the plasma nanostructure has a detrimental effect on the acceleration of the high-energy electrons. At a 0 = 5 the effect is rather weak and the spectra are comparable, while at a 0 = 15 and a 0 = 45 they exhibit great differences, being the electron cut-off energies significantly reduced with foams. This also affects the mean kinetic energy of the most energetic component of the electron spectra T e . In Fig. 2(g) T e is shown as a function of the laser intensity a 0 . At a 0 = 5 the electron temperatures obtained with the two different plasmas are comparable, whereas at a 0 = 45 the homogeneous plasma T e can exceed three times the foam one (72 MeV vs. 23 MeV, C-pol).
Direct laser acceleration of the electrons in a plasma channel is a resonant process 1,60 . The particles in the tail of the electron energy distribution are those for which the resonance condition is met, so that they stay longer in phase with an accelerating electromagnetic field. It is reasonable to expect that perturbing the electromagnetic field distribution would result in a detuning of the resonance condition. Figure 2(h and i) show the electromagnetic field inside the plasma channel for a0 = 45, for a uniform plasma and for a DLCCA-plasma respectively. In the case of the nanostructured plasma, the diameter of the channel is larger and the field is less intense and significantly more disordered.
We also consider the effect of the laser polarization. Figure 2(c,d,e,g) show that with homogeneous plasmas there is a non negligible polarization dependence of the energy absorption and the electron spectra. For instance, the highest electron temperature difference is ∼14% for a 0 = 45 and the largest cut-off energy discrepancy is ∼30% for a 0 = 15. On the other hand, all these differences are almost completely removed with nanostructured plasmas. In a simplified picture, this could be attributed to the fact that a laser pulse interacting with a nanostructured target is focused onto a non-flat surface. A similar effect was also observed in previous works 12, 13 . We remark that the homogeneous plasma case for a 0 = 5, C-pol was not considered in this discussion, since it is opaque to the laser pulse so that the interaction is qualitatively different from the other cases: with C-pol and normal incidence, electron heating mechanisms are suppressed, leading to a more efficient formation of a denser layer at the vacuum-plasma boundary with respect to P-pol.
Effect of different nanostructure morphologies at low laser intensity. In order to assess the effect of different nanostructure morphologies on laser-plasma interaction, we performed three additional 3D simulations with laser intensity a 0 = 5 and linear polarization, considering a DLA-foam, a DLCCA-foam, an ordered nanowires array and a random nanowires assembly as targets (see Fig. 1(a)). Figure 3(a) shows the propagation of the laser pulse in all the simulated nanostructured targets. The first important observation is that the qualitative features of laser-plasma interaction do not significantly differ between each other. In particular, the penetration depth of the laser in the target is similar. For both nanowires targets a strong z-pinch 54 is observed, due to the large currents induced along the wires. This effect is not observed for the two foam targets.
As reported in the previous section, a nanostructured target leads to a strong increase of the absorption efficiency with respect to a homogeneous plasma (see Fig. 3(e and f)). However, the effect of the nanostructure morphology on the total absorption efficiency was observed to be surprisingly mild. At t = 35λ/c, the absorption efficiency of laser energy into particle kinetic energy is ∼75% for all the disordered nanostructures (random wires, DLA-foam and DLCCA-foam), while being slightly higher ( 85% ∼ ) for the ordered wires array. As shown in Fig. 3(e), this discrepancy reflects into a higher absorption efficiency of laser energy into ion kinetic energy for the ordered wires array (∼30% instead of 20% ∼ as for the other nanostructured targets). More significant discrepancies are observed in the electron energy spectra, as shown in Fig. 3(d), with differences in the absolute number of high-energy electrons greater than one order of magnitude. The highest number of particles in the tail of the spectrum was observed for the random wires array. The cut-off of the electron energy distribution was found to depend only mildly on the nanostructure morphology, being always in the range 15-20 MeV.
Target morphology was observed to significantly affect the angular distribution of both energetic electrons and ions. As shown in Fig. 3(b), foam and random wires targets lead to rather uniform angular distributions of the emitted electrons. On the contrary, for the ordered wires array a ring-like structure appears, with a minimum on the laser axis and a maximum around ∼°20 (a factor of two in the intensity is observed between the minimum and the maximum of the distribution). This effect is likely due to the characteristic configuration of the electromagnetic fields for the ordered nanowire array (see Fig. 3(g)). The return current flowing inside the nanowires drives a strong azimuthal magnetic field, while charge separation at the wire surface induces a strong, radially outward, electric field. The effect of the magnetic field is to expel forward-moving electrons from the wire, while the electrostatic field is attractive (as shown in Fig. 3(h) energetic particles are most likely found in an annular region surrounding each wire). These fields hinder the propagation of electrons along the nanowires axis and the cylindrical symmetry of the distribution is a result of the cylindrical symmetry of the system. Figure 3(c) shows the distribution of the emitted ions. The foam-like targets lead to an emission with a weak angular dependence, while a strong peak at 90 θ ∼° (the transverse direction with respect to laser axis) was observed for the ordered nanowires. A uniform angular distribution between ∼°±°90 30 was obtained with the random wires assembly (30° is the maximum tilt angle of the nanowires). This can be explained considering the cylindrical symmetry of plasma expansion for the nanowires target, which is radial for each one of them. A peak appears when they are all aligned along the laser axis, while the uniform angular distribution obtained with the random nanowires array reflects the uniform distribution of their tilt angle. In these conditions (relatively low laser intensities and relatively large nanoparticles) plasma expansion cannot be completely described as a Coulomb Explosion, showing also features better described with a fluid approach 61 . Effect of the dimensionality of the simulation. Numerical simulations are an important tool for the interpretation and the design of experimental activities. 3D simulations are usually necessary if a quantitative agreement between experimental and numerical results is desired. However, 3D simulations can be highly demanding in terms of computational resources. Understanding which features of laser-plasma interaction can be reliably reproduced in 2D could allow to perform numerical investigations with lower computational resources. This is especially beneficial if extensive parametric scans are required. The effect of the simulation dimensionality is particularly crucial for the case of a nanostructured target, since it might not be straightforward to reproduce in 2D an intrinsically 3D morphology. For this reason, for each 3D case, we performed a corresponding 2D simulation and we compared the main observables.
As far as a foam-like plasma is of concern, it is challenging to reproduce a 2D fractal-like structure with the same properties of its 3D analogous, keeping a similar filling factor and mean density. Thus, we used a random collection of nanospheres with the same radius, filling factor and density of those of the 3D case. Similar arguments apply for the nanowires targets, since in 2D it is not possible to reproduce at the same time the wire morphology (radius and spacing), the mean density and the local density. Thus we performed two simulations with different targets: one with the same spacing and wire radius of the 3D case but a higher density (∼12 n c ), the other one with the same mean density of the 3D case (3 n c ) but different geometry (thinner wires).
It should be highlighted that some physical processes such as charge separation, diffraction and self-focusing are inherently affected by the simulation dimensionality. For instance, in 2D geometry a charged particle is actually an infinite charged wire, thus the electrostatic potential is different from that of a charged particle in 3D. As a consequence, processes governed by charge separation, such as Coulomb explosion or Target Normal Sheath Acceleration, are expected to be heavily affected by the simulation dimensionality. We observed that the numerical results are indeed significantly affected by the simulation dimensionality. As a general feature, at fixed conditions, the pulse is able to propagate more deeply into the target in the 3D simulations. The laser absorption is rather similar (⪅10%) except for the a 0 = 5 homogeneous plasma case. In this last case, in 3D the ponderomotive force of the laser creates a clean density spike at the pulse erosion front, which reflects the incoming beam. On the other hand, in 2D the erosion front exhibits modulations which probably boost laser absorption with respect to the former case ( 40% ∼ in 2D, ∼25% in 3D). This effect is relevant only for plasmas close to the transparency threshold, since for highly transparent cases the laser penetrates deeply in the plasma and is efficiently absorbed regardless of the dimensionality.
We also report that the target morphology is crucial when studying nanowires in 2D. The plasma shaped to reproduce the 3D geometry resulted into an overdense material (∼12 n c ), reflecting the laser beam within few wavelengths. On the other hand the 2D wires with a 3 n c average density performed qualitatively similarly to the 3D case despite being much thinner.
Also the electron temperature is deeply influenced by the simulation dimensionality, being as far as twice the value in 2D with respect to the 3D case (the highest difference was found for the nanostructured plasma at a 0 = 45).
In this work we also studied the laser self-focusing within a near-critical plasma. At this regard we define the self-focusing parameter = | | + | | SF 0 where E and B are the electromagnetic fields and I 0 is the initial peak laser intensity. We always observed self-focusing except for a 0 = 5, due to the low transparency of this case. We report a low dependence of SF with respect to the simulation parameters. In 3D simulations SF was within the 4.5-6 range, regardless of the laser polarization or the target structure. As also reported in previous works 1 , the SF is lower in 2D simulation, being ∼3, since focusing can only occur along a single transverse direction instead of in a whole plane.

Discussion
In this work, we observed that the interaction of a relativistic laser pulse with a nanostructured NCP is significantly different than with a homogeneous one. Moreover, the morphology of the nanostructure was found to affect the interaction.
At moderate laser intensity (a 0 = 5), a strong enhancement of absorption efficiency with respect to the homogeneous plasma case is reported, regardless of the morphology. This suggests that the use of low-density nanostructured materials could be beneficial for applications strongly relying on absorption efficiency (e.g. laser-driven ion acceleration). In the literature, aligned nanowire arrays have been shown to lead to very high absorption efficiencies. Here we show that similar performances can be obtained also with foams and random nanowire assemblies. This enhancement was found to be milder for a 0 = 15 and negligible for a 0 = 45. However, in all the investigated cases, the presence of a nanostructure leads to a strong increase of the absorption efficiency of laser energy into ion kinetic energy, due to the explosion of the nanostructures.
At high laser-intensity (a 0 = 45), a marked detrimental effect of the nanostructure was observed in the tail of the electron energy distribution, with a ~50% reduction of the cut-off energy with respect to the homogeneous plasma. Negligible differences were obtained for the total absorption efficiency, however the target nanostructuring leads to an order of magnitude increase of the efficiency of energy transfer to the ion population. These results are relevant for experimental schemes aiming at the realization of high-energy synchrotron sources based on laser interaction with a NCP, since the energy of the emitted photons strongly depends on the energy of the electron population.
The morphology can strongly influence the angular distribution of energetic particles in laser interaction with nanostructured NCPs (this effect has been investigated only for a 0 = 5). In particular, using an ordered array of aligned nanowires instead of random nanostructures leads to a ring-like shape of the angular distribution of emitted electrons and to an ion emission strongly peaked near the transverse direction with respect to laser axis. The latter effect could find application as a possible diagnostic of the integrity of a nanowire array during an actual experiment. Indeed the observation of a peak in the angular distribution of emitted ions at 90° from laser axis would be a strong signature of the survival of the nanostructure for long enough to interact with the main pulse.
As an important remark, we have observed that, although 2D simulations are able to qualitatively reproduce the 3D results, the simulation dimensionality may deeply affect the physical processes at play. As an example, at low laser intensity the electron temperature is comparable between 2D and 3D cases despite a much different laser absorption. Moreover, reproducing the nanostructure morphology at lower dimensionality was found to be by no means straightforward.
In conclusion, our results show that including a realistic description of the nanostructure is essential to properly model physical processes of direct relevance both for a satisfactory understating and for the foreseen applications of laser interaction with nanostructured NCPs. Moreover, they suggest possible paths to guide the design of future experimental activities. The laser pulse had cos 2 temporal profile and a Gaussian transverse profile and it was focused at the vacuum-plasma boundary. The waist was 5λ, while the temporal duration was 15λ/c (full-width-half-maximum of the fields). The intensity was varied between a 0 = 5 and a 0 = 45 at fixed normal incidence. These parameters, if scaled to Ti:Sapphire lasers, correspond to a 40 fs beam focused into a 4 μm spot, providing a peak intensity in the range 5 ⋅ 10 19 W/cm 2 < I < 3.2 ⋅ 10 21 W/cm 2 . Corresponding simulations with different polarizations are performed keeping the pulse energy constant. The reported value of a 0 always refers to P-polarization and should be divided by 2 1/2 for C-polariazion.

Particle
For the 3D simulations with a uniform plasma, a 90λ × 50λ × 50λ box was used, with a resolution of 20 points per wavelength. A uniform plasma slab filled a 50λ × 50 λ × 50λ region with a density of 3n c . 12 macro-electrons per cell were used.
For the 3D simulations with a nanostructured plasma, a 80λ × 30λ × 30λ box was used, with a resolution of 50 points per wavelength. The plasma filled a 40λ × 30λ × 30λ region (30λ × 30λ × 30λ region-thick for the wire targets) with an average density of 3n c . The local density was n 60 c ∼ , with a filling factor of 5% ∼ (the local density was adjusted to have an average density of exactly 3n c ). 64 macro-electrons per cell were used.
We performed additional convergence tests to ensure that both the spatial resolution and the number of macro-electrons per cell are adequate for our physical scenario. The case a 0 = 45, P-polarization, homogeneous plasma was selected and simulations with resolution of 20, 30 and 40 points per wavelength were carried out. Negligible differences were observed for the electron energy spectra and for the absorption efficiency in these cases. Then, the case DLCCA-foam, P-polarization, a 0 = 45 was simulated in a reduced box (40λ × 6λ × 6λ) with a plane wave laser, using either 64 or 128 macro-electrons per cell. No significant difference was found in the observables relevant for this work.
In the 2D simulations with a uniform plasma, a 160λ × 70λ box was adopted with a resolution of 100 points per wavelength. The uniform plasma (100λ × 70λ with a 3n c density), was sampled with 4 macro-electrons per cells. The nanostructured plasma was obtained through a random ensemble of nanospheres with a radius of 0.05λ and density 60 n c . This corresponds to a filling factor of ∼5% and a volumetric mean density of 3n c . This plasma was sampled with 169 electrons per cell.
In all the simulated cases the plasma was fully pre-ionized and the the charge/mass ratio of the ions was 0.5 (e.g. C 6+ ). The electron population was initialized with a small temperature (few eVs) to avoid numerical artifacts. In order to verify that this choice does not affect our results, we performed an additional simulation with a high initial temperature (∼1 keV) in a reduced box. We observed only a small blurring of the nanostructures and negligible differences for the observables of interest in our work. Nanostructured foam models. The DLCCA-foam structure is generated through an extension of the Diffusion-Limited Cluster-Cluster Aggregation model (DLCCA), built to simulate the structure of fractal aggregates, such as colloids and soot 63 . In its original version, the DLCCA model starts with a collection of monodisperse spherical particles, randomly disposed into a cubic mesh. Then, these particles can move in random directions, to simulate the Brownian diffusing motion; once the spheres touch each other, they stick together to form a rigid dimer, which can undergo random motion as well. This mechanism is repeated until all the particles in the box aggregate together in a large fractal structure. In this work, the DLCCA model is used to produce aggregates with different particles numbers, which follow an exponential distribution. Thus, the aggregates are vertically deposited onto a substrate to form the foam structure. The spheres are generated on a regular lattice with a radius of 0.0375λ, barely touching each other. In order to obtain a more realistic structure, the spheres are inflated up to a radius of 0.05λ and randomly displaced by a small fraction of their radius.
The DLA foam is obtained with a simple Diffusion-Limited Aggregation (DLA) model on lattice 59 that works as follows. Particles are injected into a 3D uniform lattice from random sites located on the top of the lattice itself. Then, they are allowed to randomly walk in the lattice until they are halted because either they occupy a point adjacent to a pre-formed cluster or they ended up onto the bottom of the lattice, i.e. the substrate. Periodic conditions were used at the lateral boundaries. Using the same parameters (sphere radius and box size) as those used to generate the DLCCA foam we obtained DLA foams with much lower filling factors. For this reason, at first the DLA foam was generated with sphere particles with reduced radius (0.035λ) that then were inflated up to the desired radius (0.05λ). In this way we could obtain a filling factor value comparable to that of the DLCCA foam (~5%).
Both foams are generated in a 40λ × 6.4λ × 6.4λ box, which is then repeated periodically in the whole simulation box.
Nano-wire targets are modelled as a collection of cylinders, with a radius of 0.12λ (i.e. ∼100 nm for a Ti:Sapphire laser, since in this case λ = 800 nm). Nanowires with similar properties and with a length of tens of μm have been reported in the literature 39 . For the parallel wires target, the cylinders are aligned with the laser axis, while for the random wires target, their inclination angle is extracted randomly from a uniform distribution between −30° and +30° with respect to laser axis. The wires are truncated abruptly at the boundaries of the plasma region. Fig. 1(c) are produced through the Pulsed Laser Deposition (PLD) technique using the parameters described by Zani et al. 41 . To obtain the images, the Zeiss Supra 40 field emission Scanning Electron Microscope (SEM) is exploited with an accelerating voltage of 5 kV.

SEM images of low-density foam aggregates. Carbon Foams showed in
Data availability. The data supporting the findings of this work are available from the corresponding authors on request.