3.1. Self-Organized Nanopatterning
Light manipulation at the nanoscale relies on the ability to control morphological parameters of nanostructures. In our case, the nanofabrication process exploits in the first steps a self-organized array of Au nanowires which acts as a stencil mask. Under grazing incidence ion irradiation, an initially flat Au film develops a periodic ripple undulation at the metal-vacuum interface (
Figure 1d). The ripple amplitude increases with ion dose while the residual thickness of the Au film decreases. When the bottom of the ripple valleys reaches the underlying glass, a disconnected array of nanowires is thus formed (
Figure 1e). It acts as a stencil mask for etching the underlying substrate at later stages. Representative AFM topographies of Au nanowire stencil masks prepared either with (
Figure 1a) or without (
Figure 1b) a Ti adhesion layer interposed at the Au/glass interface show strong differences in statistical height; RMS roughness (σ) decreases from σ = 17 ± 2 nm to σ = 11 ± 2 nm when a Ti layer is employed. One can also notice that the longitudinal coherence and the order of the ripples are improved in nanowires synthesized with Ti layer and their lateral size is narrower, as shown by the topographic profiles superimposed in
Figure 1c.
In order to quantitatively describe the role of titanium in the pattern formation process, we analyzed AFM morphological data in the framework of dynamic scaling theories, recurring to second order statistical quantities [
26,
27]. In the case of closely packed 1-dimensional nanostructures formed by IBS, one can easily derive the lateral scale of the structures from the lateral correlation length parameter, ξ, extracted from the Height-Height correlation function
Hcf(
r) = < [
h(
x+
r)-
h(
x)]
2 >.
Hcf(
r) has been systematically calculated for every AFM topographic line profile
h(
x), acquired along the fast scan direction which was kept orthogonal to ripple’s long axis. The data were cumulated and averaged over different AFM images acquired in equivalent areas in order to decrease statistical errors.
According to dynamic scaling theories for small values of the lateral displacement
r,
Hcf(
r) follows a power law, while for large displacement values,
Hcf(
r) tends to a saturation value 2σ
2, where σ is the RMS roughness [
26,
27]. The asymptotic break between the two regimes of
Hcf(
r) allows us to determine the statistical lateral correlation length, ξ, i.e., the distance r within which the surface height fluctuations are correlated. For distances
r > ξ, the height variations are instead spread randomly [
26,
27]. In
Figure 3b we report the lateral correlation length ξ extracted from the
Hcf(
r) as a function of the increasing ion fluence. In the case of Au nanowires mask it reads 95 ± 11 nm, while in presence of Ti this figure is reduced down to 65 ± 11 nm. The morphological differences introduced by Ti adhesion layer can be understood if we consider the role of Ti as an inhibitor of dewetting processes, either thermally activated or ion induced, thanks to the chemical and metallic bond of Ti with glass and Au respectively.
Once the Au nanowire array is formed, it acts as a stencil mask which guides selective etching of the underlying substrate in correspondence to the disconnected gaps. At this stage of the process the substrate is heated to 700 K and an auxiliary Au atom source is switched on (sketch in
Figure 2d). The vertical scale of projected nanostructures is enhanced since Au adatoms evaporated on the heated surface can diffuse for longer distances until they reach the top of ripples, where they enrich the stencil mask and prolong its lifetime during ion irradiation. At the end of the process, the samples were further exposed to the ion beam until complete removal of the gold stencil mask was achieved. Further details on the A-IPL process can be found in reference [
18].
In order to get high aspect ratio and tunable nanostructures capable of optical functionalization of the glass templates in terms of antireflection and light scattering effects [
28,
29], the A-IPL process has been applied to different stencil masks endowed with engineered morphologies presented in
Figure 1. To this purpose, different substrates were prepared by A-IPL at
T = 700 K and performed on the borosilicate glass with three different ion fluences:
F1 = 1.63 × 10
19 ions·cm
−2,
F2 = 2.08 × 10
19 ions·cm
−2, and
F3 = 2.97 × 10
19 ions cm
−2, corresponding to Sample 1, Sample 2, and Sample 3, respectively. For the sake of simplicity, the samples synthesized at the same fluences, but with a Ti layer interposed between the glass and Au film, are named Sample 1-T, Sample 2-T, and Sample 3-T, respectively.
In
Figure 2a,b the morphology of representative samples prepared respectively in the absence and in presence of titanium, for a given ion fluence (
F2), are compared. Significantly, we can notice that in both cases the pattern transferred from the Au NWs mask into the glass exhibits an almost ten-fold vertical amplification. RMS roughness goes from σ = 17 ± 2 nm (11 ± 2 nm) to an average σ = 136 ± 9 nm (σ = 94 ± 8 nm) respectively without (with) Ti. Moreover, nanostructures synthesized in presence of Ti adhesion layer appear more elongated in the longitudinal direction and present a narrower lateral size, as shown by the comparison of line profiles in
Figure 2c. In a more quantitative way, the lateral correlation length ξ reads 220 ± 11 nm for Sample 2, while we measure ξ = 160 ± 11 nm in presence of Ti (Sample 2-T).
In addition to the morphological tunability favored by Ti adhesion layer, we stress the broadband nature of the pattern presented in
Figure 2a,b ensured by the presence of multiscale nanostructures organized in a pseudo-periodic array. Since diffractive properties are dependent on the Fourier spectra of the spatial frequencies, we calculated one dimensional power spectral densities (1D-PSD plotted in
Figure 3a) from AFM data acquired on Sample 2 and 2-T, in similar manner to what was done for the
Hcf(
r) calculation.
Firstly, we highlight the broadband nature of the roughness distribution for both samples: a characteristic of paramount importance in view of optical broadband functionality required in PV applications; secondly, we also emphasize the peculiar behavior of Sample 2-T, which shows higher vertical amplification of surface features with large wavenumbers (small spatial wavelengths), while Sample 2 displays a strongly amplified roughness associated with small wavenumbers (large spatial wavelengths). This effect can be ascribed to the Ti layer presence which helps in preserving smaller structures on the sample and decreases coarsening effects.
Similar considerations also hold for the other couples of samples prepared at different ion fluences
F1 and
F3, which allow us to draw some conclusions on the general trend observed in the morphological data. In
Figure 3b, the averaged statistical parameters σ and ξ of all samples, extracted from their
Hcf(
r) functions, are summarized as a function of the increasing ion fluence. For what concerns the amplification of σ and ξ with increasing fluences, a clear trend can be inferred.
RMS roughness σ and lateral correlation length ξ both increase by prolonging the ion dose during the A-IPL process. This increase of roughness at increasing ion fluences can be understood considering that the A-IPL process involves the use of an auxiliary gold evaporation source during the ion bombardment in order to prolong the shielding action of the Au nanowire stencil mask; in this way, the glass substrate is etched at deeper levels. This is confirmed by the strong increase of σ evident in
Figure 3b, where points at a fluence lower than
F1 = 1.6 × 10
19 ions cm
−2 are extracted from Au NWs mask topographies. The increase in lateral correlation length at higher fluences can be ascribed to the coarsening of the smaller scale spatial modulations due to the ion and thermally activated glass relaxation. In this way, ion fluence is an experimental parameter which gives us the opportunity to easily tune the morphology of nanostructured glasses and hence their optical properties.
On the other side, we also highlight that, for each explored ion fluence, the Ti adhesion layer allows us to decrease the average lateral scale and vertical extension of the nanostructures, expressed respectively by the statistical parameters ξ and σ (red open squares and triangles in
Figure 3b). In particular, for the lowest ion fluence
F1, the lateral correlation length of Sample 1-T measures 105 ± 11 nm; this value is compatible with the one registered on the Au ripple pattern before its decomposition in nanowires [
18]. In the early stages of the A-IPL process the adhesion properties of the Ti layer are thus capable of stabilizing the Au nanowires stencil mask, preserving the small features present on the pristine Au ripple pattern.
The main physical mechanism is thus related to the inhibition of dewetting of the Au film once the continuous layer breaks up into an array of nanowires. The dewetting process is boosted by the hyperthermal energy deposition following ion impact which favours uphill mass transport of Au, thus exposing large portions of unshielded glass to harsh ion etching. This effect is manifested when the A-IPL process is prolonged at higher ion fluences F2 and F3, where the Au mask without Ti adhesion layer undergoes strong dewetting and allows the formation of large volume structures in the glass as a consequence of coarsening. In the presence of Ti instead, the pristine pattern of the Au continuous film is preserved and coarsening of the glass structures is prevented, resulting in more ordered and narrower glass ripples.
In this picture, the Ti adhesion layer thus represents a further experimental parameter (in addition to ion fluence), which allows to control the morphology of the nanopatterned substrates prepared by ion sputtering in the A-IPL configuration. The role of the engineered morphologies will be now explored from the optical point of view in order to elucidate their use in photon harvesting applications.
3.2. Optical Properties of Nanostructured Glasses
In order to characterize the optical performance of the nanostructured glass templates, we performed total integrated reflectance and haze measurements on the whole set of samples so far described (Samples 1, 1-T, 2, 2-T, 3, 3-T).
In
Figure 4a we study the differential variation of the total integrated reflectance in the Vis-NIR range (450–1000 nm) of the nanostructured glass substrates normalized to a flat reference borosilicate glass, according to the equation |Δ
R/
R| = |(
Rnano −
Rflat)/
Rflat|, where
Rflat ≈ 8%. All of the diffuse reflectance measurements were carried out by means of an integrating sphere setup. The Δ
R/
R quantity describes the amount of anti-reflective effects present on nanostructured samples. Sample 2 and 2-T (blue squares and open triangles, respectively) reflect about 12% and 20% less light than a flat glass averaged on the whole Vis-NIR range of wavelengths. To summarize the information relative to all samples synthesized at different ion fluences, in
Figure 4b (left scale) we plot the mean values of Δ
R/
R, integrated over the wavelength range Δλ = [450 nm; 1000 nm]. As the fluence is increased in the A-IPL process, the total reflectance raises towards the values of the flat glass and the antireflective effect vanishes. In particular, we observe an average reduction of reflectance from 23% (Sample 1-T) to 20% (Sample 2-T) to 6% (Sample 3-T-open red triangles). Correspondingly, for the samples without Ti, the average reduction of reflectance goes from 14% in Sample 1, to 5% for Sample 3 (blue squares). In other words, by modulating ion fluence one can control broadband anti-reflective effects, which are particularly enhanced in presence of the Ti layer.
Such behavior is in qualitative agreement with the “moth-eye” effect observed in bio-mimetic systems [
30]. In fact, high aspect ratio corrugations, whose amplitude is a significant fraction of the incident light wavelength and whose lateral size is below the diffraction limit, allow us to gradually match the refractive indexes of substrate and air, leading to a suppression of reflectance (index-grading effect [
29,
31]).
The morphological analysis carried out on the six samples (
Figure 3b) confirms that the Ti layer leads to the formation of more ordered patterns endowed with smaller scale non-diffractive features for photons belonging to the Vis-NIR spectrum. The observed index-grading anti-reflective behavior is predominant in Sample 1-T, which has the smallest correlation length in the range of 105 nm. We stress that these features can be useful for all kinds of applications in which a decrease of surface reflectivity is an important issue, as in de-trapping schemes employed in LED/OLED applications [
32].
At a higher ion fluence, the lateral correlation length of the spatial corrugations increases due to coarsening: the importance of the moth-eye antireflection effect decreases and large angle diffuse scattering becomes predominant [
21,
33]. In view of light trapping applications employing the self-organized glass gratings, we stress that the scattering efficiency can be of predominant importance with respect to the anti-reflection behavior; enhanced light absorption in a PV device can in fact be achieved when scattered photons penetrate the absorber layer at the grazing angle and eventually are coupled to wave-guided modes by total internal reflection [
34].
The scattering efficiency of the substrate, quantified as the ratio of diffuse reflectance to total reflectance (Haze), is determined experimentally by means of an integrating sphere. In
Figure 4a, the haze signals measured on sample 2 (blue line) and 2-T (red-dashed line) are reported as a function of light wavelength. Haze values decrease from 60% to 45% (43% to 28%) in the wavelength range 450–1000 nm; a comparison with conventional Asahi-U reference substrates, measured in the same experimental condition, which exhibit haze values in the range 30%–2% [
35] demonstrates the remarkable scattering efficiency of the novel substrates. The broadband scattering efficiency is high across the whole Vis-NIR spectrum and can be ascribed to the multiscale and pseudo-periodic nature of the self-organized structures highlighted in the reciprocal space spectral analysis (1D-PSD in
Figure 3a and AFM data shown in
Figure 2). The possibility to engineer the roughness spectrum of the nanostructured patterns, either by varying ion fluence or by adding a Ti adhesion layer, allows us to selectively enhance the scattering efficiency or the antireflective properties of the self-organized surfaces. The high haze of Sample 2 can be attributed to the presence of large scale spatial modulations which contribute to the 1D-PSD (
Figure 3a blue trace) in the small wavenumber region below 1.5 × 10
−3 nm
−1. On the contrary, the moth-eye antireflective effect measured on Sample 2-T can be attributed to the presence of small scale structures, as witnessed by the cross-over of the 1D-PSD (
Figure 3a red trace) observed in the large wavenumber region above 1.5 × 10
−3 nm
−1.
For what concerns the light scattering performance of the substrates prepared with different ion fluences, in
Figure 4b (right scale) we summarize the average values of haze integrated over the whole spectral range (full and open circles refer to samples with and without Ti, respectively). One can notice that (i) higher haze values are observed for samples exposed to larger ion fluences; (ii) at fixed ion dose, haze is higher on the sample prepared without Ti adhesion layer. These observations can be rationalized in terms of the morphological evolution of the substrates with ion fluence summarized in
Figure 3b. The RMS roughness and the lateral correlation length of the structures become larger (coarsening) as a consequence of dewetting of Au NWs mask, especially in absence of Ti. Due to coarsening, the dominant lateral size of the surface features grows at a point that scattering (haze) increases while index grading (moth-eye anti-reflective effect) decreases, as observed in the Δ
R/
R trace of
Figure 4b.
A further aspect, of paramount importance for light trapping amplification in opto-electronic devices, is represented by the possibility to couple incoming light to wave-guided modes which confine propagation in the high refractive index layers formed by the substrate and absorber film [
36,
37].
In
Figure 4d (see also sketch in
Figure 4e), we directly visualize the presence of waveguided light which propagates laterally via internal reflection through the slab after being scattered at the textured glass interface. A light beam is propagated through the nanostructured glass slab in transmission geometry (pattern analogous to Sample 2), while the directly transmitted component is blocked by a beam stopper. Red dashed lines indicate the glass borders, while white dashed lines highlight the limits of the patterned area. Pictures have been taken from the top side in absence of ambient light, excluding the direct beam in order to avoid saturation of the camera. For the flat glass,
Figure 4c, light is entirely transmitted in the direct beam (and thus blocked by the beam stopper). For the patterned glass,
Figure 4d, light scattered perpendicularly to the ripple axis is guided for tens of millimeters through the slab by total internal reflection before leaking out and propagating to the detector.
We stress that the scattering and antireflective efficiency highlighted above are expected to be even more important when a high refractive-index material (e.g., silicon) is deposited on the nanostructured glass. In fact, the glass surface, texturized with a pseudo-periodic pattern, works as diffraction grating that can couple incident light to optical modes supported by the high-refractive index semiconductor layer, strongly enhancing its effective absorption [
38].
3.3. Photon Harvesting in Thin Film Solar Cells
The optical performance of the patterned glass when applied as substrate for thin film solar cells has been investigated by fabricating p-i-n amorphous silicon solar cells and evaluating EQE and cell reflectance. The EQE measured for identical devices grown on flat and nanostructured substrates is shown in
Figure 5a for Sample 2 and 2-T, which exhibited both antireflective and scattering effects. The response measured for both the rough cells is significantly enhanced when compared to the flat reference cell (black line). The EQE enhancement in the long wavelength range is indicative of efficient light trapping by means of multiple passes of scattered light within the absorber layer. In addition, a reduction of the interference fringes is observed, due to scattering from the uncorrelated roughness of the internal interfaces conformed to the pattern of the substrate. Quantitatively,
JSC reaches 11.8 mA/cm
2 and 13.5 mA/cm
2 for Sample 2-T and 2, respectively, against 10.8 mA/cm
2 for the flat cell. The total reflectance of the same solar cells is reported in
Figure 5b. A significant reduction with respect to the flat case (black line) is measured for both the cells on textured glass and for Sample 2 in particular. This is in agreement with the results on the bare glass (
Figure 4a), where Sample 2 is characterized by higher haze with respect to 2-T, which translates into a more effective light trapping in the solar cell. Similarly to what has been done for the bare glass substrates, we evaluated the relative variation |Δ
R/
R|
cell = |
Rnano −
Rflat|/
Rflat of the integrated reflectance over the full spectral range for the cells on different nanostructured substrates versus the flat cell. The latter parameter is shown as blue bars in
Figure 5c,d, in order to quantify the efficiency of light handling within the cell. In general, a stronger reduction of reflectivity, i.e., improved light harvesting, is observed as the ion fluence increases, although for the samples without Ti adhesion layer (
Figure 5d), only two cells are available. A further observation is that, for fixed ion dose, the cell without Ti has a better light harvesting efficiency. These results compare well with the characterization of the substrates shown in the previous section; as expected the light trapping efficiency scales progressively with samples endowed with higher haze.
Figure 5c,d also reports the corresponding relative enhancement of
JSC measured for the same solar cells Δ
JSC/
JSC-flat. From
Figure 5c, for the glasses obtained with Ti layer, it appears that this parameter increases with ion fluence, reaching a maximum for the intermediate value (Sample 2-T). For the largest fluence, on the other hand, the relative enhancement of
JSC drops while the optical enhancement still increases. We attribute this divergent behavior to the monotonic increase of the RMS height of the nanopattern—as shown in
Figure 3b—which likely does not guarantee the full conformal growth of the semiconductor layers of the cell due to shadowing effects, with formation of defective regions. In fact, this issue is also found for the cell on Sample 3, without Ti layer, and is confirmed by the low shunt resistance and reduced open circuit voltage measured in both cases (data not shown). The severe deterioration of the electrical performance likely disturbs the charge collection with impact on the EQE and
JSC, thus masking the optical enhancement. It is known that high aspect ratio features, while valuable for an optical improvement, are usually detrimental to the electrical performance of the solar cells, as they can easily induce the formation of low-density faulty regions [
2,
38,
39]; a trade-off between optical and electrical properties must thus be found.
A post-growth refinement of the nanostructured glass morphology for improved electrical performance has been attempted via ion beam sputtering at normal incidence (normal-IBS), prior to solar cell deposition. The treatment exploits ballistic smoothening of the nanostructured surface [
40]. We selected Sample 2 as a reference, since this substrate produced the largest
JSC enhancement, while other relevant photovoltaic parameters where below the values obtained with the flat reference cell (e.g., open circuit voltage (V
OC) reads 776 mV versus 840 mV for the flat cell). We synthesized three nanostructured glasses with morphology equivalent to Sample 2 and followed the evolution of the morphological properties after 2 min and 4 min of normal-IBS process (see sketch in
Figure 6e). The pristine and the two treated glasses were finally used as substrates for co-deposited a-Si:H solar cells together, with a reference flat glass.
Figure 6 shows a morphological and optical comparison of the three samples. The high aspect ratio nanostructures typical of Sample 2 (see
Figure 2a) are heavily smoothed by 2 min and 4 min cycles of normal-IBS treatment (the AFM image in
Figure 6a refers to a 4 min cycle), and the V-shaped valleys typical of the pristine pattern are tailored to shallower U-shaped ones, as demonstrated by the superimposed AFM profiles in
Figure 6d. The RMS roughness σ decreases from the initial value of ~150 nm, to 120 nm after 2 min of normal-IBS, and down to only 56 nm when doubling the treatment time. In order to gain insight into the smoothing mechanisms active with normal incidence ion irradiation, an analysis of the slope distribution of the AFM topography has been carried out. The first derivatives of the AFM images were computed orthogonally to the ripple axis and converted into normalized histograms of the local slope distribution (
Figure 6b). For the pristine sample (blue line), the distribution is broad in the range of −60°/+60°. After 2 min of normal-IBS process (red line), a partial slope selection of the structures is observed, as evidenced by the two peaks in the distribution located approximately at −50° and +50°. A flattening effect in correspondence to the bottom of the U-shape valleys is evidenced by the growing peak around 0°. If the ion irradiation at normal incidence is protracted for 4 min (green line), the slope distribution is heavily impacted: the peak around 0° becomes dominant (strong smoothing effect) and two small peaks centered around ±15° are observed. The morphological reorganization is likely driven by the reduction of the surface free energy. It involves a smoothening of the corrugations caused mainly by ion induced relaxation mechanisms such as ion induced viscous flow, thermally activated surface self-diffusion, and ion induced surface diffusion [
40]. Normal incidence IBS thus appears to be an effective and fast process for tailoring the glass morphology, competitive with respect to other plasma treatments implemented on the transparent conductive oxide [
2,
39].
The effect of the treatment on the light scattering properties of the glasses is severe after a 4 min process, since the haze parameter drops to values in the range 40% ÷ 20% across the spectral interval (
Figure 6b). After a 2 min cycle, however, the haze maintains relevant values around 60% ÷ 50%, which qualifies the latter sample as a promising substrate for thin film solar cells.
Table 1 reports the photovoltaic parameters and relative variation of integrated reflectance of the co-deposited a-Si:H solar cells grown on the reference flat glass, on the pristine texture (0 min equivalent to Sample 2) and on treated nanostructured glass (2 min and 4 min). A relevant optical gain versus the flat reference is confirmed for the cells on pristine substrate and on the sample treated for 2 min (
JSC above 13 mA/cm
2 against 11.1 mA/cm
2 for the flat case), while as expected the gain is negligible for the cell on the glass treated for 4 min (Δ
JSC = 0.7 mA/cm
2). The lowered light harvesting capability, as determined by the total reflectance measurements on the cells, confirms the effect of the smoothing treatment. The relative variation of integrated reflectance |Δ
R|/
Rflat of the cells decreases dramatically from 50% for the cell on the pristine texture down to only 2% for the sample treated for 4 min. Remarkably, for the sample treated with a 2 min cycle, the |Δ
R|/
Rflat value remains almost unchanged at around 40%, while the electrical performance is positively impacted. The electrical parameters of the flat reference cell are modest since the p-layer/TCO contact has not been adjusted to the AZO layer here in use as front electrode; nevertheless, the improving trend after the smoothening treatment is clear. For the cell on the pristine texture, a severe drop is measured both for
VOC and FF (the latter caused by low shunt resistance), which can be ascribed to the V-shaped rough morphology of the substrate which prevents a conformal growth of the semiconducting layers. Conversely, the cells grown on treated substrates show a promising partial recovery of the electrical parameters. In particular Δ
VOC, evaluated versus the flat cell, is reduced to ~60 mV (for the 2 min substrate) and ~10 mV (for the 4 min substrate), against ~150 mV for the cell on the pristine texture. The relevant mechanism is the adjustment of the local slopes. It has been shown, at least for microcrystalline Si, that V-shaped valleys with unfavorable opening angles—corresponding to slopes higher than 35° in our reference system—are prone to developing voids and fissures during the silicon layer growth, with a detrimental impact on
VOC [
39]. The almost complete
VOC recovery achieved with the substrate treated for 4 min can thus be explained with the slope distribution fully constrained between ±35° (green line in
Figure 6b). Further work is in progress to optimize the proposed treatment.