1. Introduction
Tight (numerical aperture NA > 0.3) focusing of ultrashort (femto – or picosecond, fs/ps) laser pulses in bulk dielectrics was approved as a key-enabling process in laser inscription of functional embedded optical or microfluidic elements [
1,
2,
3]. This high-NA laser inscription regime provides in the focal region high enough peak laser intensities for the threshold-like optical breakdown (near-critical electron-hole plasma formation) and sufficient energy density deposition [
4]. The photogenerated near-critical electron-hole plasma (EHP) promptly reduces the real part of local effective dielectric function of the photoexcited dielectric material, Re{ε∗}, down to zero (by definition of the critical plasma frequency and density). This makes possible intra-pulse laser excitation of interfacial plasmon-polaritons at the interface between the photoexcited (dielectric function ε∗) and unexcited dielectric (dielectric function ε) for Re{ε∗} = ε
1* + iε
2*≤ −Re{ε} = ε
1 + iε
2 with ε
2 ≈ 0 [
5,
6]. The accompanying interference between the incident laser radiation and the induced interfacial electromagnetic waves strongly modulates the electric field and plasma density distributions in the interaction focal region; however, such near-or even sub-wavelength modulations of plasma density are hardly resolvable by optical imaging or diffraction methods, typically used for probing laser-generated plasma in dielectrics [
7,
8,
9]. Non-linear focusing of ultrashort laser pulses via Kerr and Raman-Kerr effects and their resulting channeling makes the plasma and related plasmonic paths considerably elongated in the pre-focal region and potentially more pronounced in the underlying laser propagation, spectral and structural effects, with the typical plasma densities in the near-critical range [
4,
6] due to NA-dependent scaling [
10].
Moreover, during the last decade fs/ps-laser damage tracks in dielectrics were intentionally characterized in their cross-sections at high magnification by electron [
11,
12,
13,
14] and atomic force microscopy [
6,
15]. Such cross-sectional analysis revealed complex, regular and hierarchical structure of the damage features, typically emerging on sub- and near-wavelength spatial scales and related to laser interference phenomena, accompanied by nanoplasmonic interactions [
6,
15,
16,
17,
18]. Though laser propagation regimes – geometrical (linear) focusing or filamentation (non-linear focusing) – were not identified in most of these previous studies, laser energy flow characteristics are rather monotonous, comparing to their modulation via interferential laser-plasmon interactions in the plasma channels [
4,
19,
20]. These promising potential indications of nanoplasmonic processes during laser-dielectric interactions call for new enlightening experimental and theoretical studies.
In this study, optical transmittance spectra of ultrashort (0.3 - 12 ps) laser pulses at 515- and 1030-nm wavelengths and variable pulse energy were acquired in the pre- and filamentation regimes in bulk natural and synthetic diamonds to search for the first time in situ for spectral signatures of dynamic electron-hole plasma density modulation related to non-local laser-plasmonic interactions.
2. Materials and Methods
Diamond samples used in this studies were cubes (dimensions – 2 × 2 × 2 mm) with six polished opposite facets. The first sample was a IaB-type natural diamond (total nitrogen content - 190 ppm), while the second one was a pure IIa-type synthetic diamond (total nitrogen content <1 ppm).
Their transmission spectra were examined, employing ultrashort laser pulses at the central 1030-nm (fundamental, FH) and 515-nm (second harmonic, SH) wavelengths of a Yb-laser Satsuma (Amplitude Systemes, St. Etienne, France), corresponding to the spectral full-widths at the half-maximum of 2.7 and 1.3 nm, respectively. The FH and SH ultrashort laser pulses were precisely focused within the crystal volume at the refractive-index corrected 400-μm depth , using a micro-objective with a numerical aperture NA = 0.25, to yield the focal spot radius at the 1/e-intensity level ≤ 4 and 2 μm, respectively. Their pulsewidths were varied at the maintained bandwidth shape via partial positive chirping (incomplete back-compression of pre-stretched pulses upon amplification) in the range τ=0.3 - 12 ps, measured by means of a single-pulse auto-correlator AA-10DD-12PS (Avesta Project Ltd., Moscow, Russia). In the transmission studies, laser pulse energy was adjusted, using a thin-film transmissive attenuator (Standa, Vilnius, Lithuania), in the ranges E=50 - 5000 nJ (FH) and E=50 - 800 nJ (SH). The transmitted FH or SH radiation was collected by a 0.2-NA quartz/fluorite microscope objective (LOMO, St. Petersburg, Russia) and directed to the entrance slit of a spectrometer ASP-IR-2.5 (spectral range - 0.9 - 2.5 μm) or ASP-190 (spectral range – 190 - 1100 nm), both manufactured by Avesta Project Ltd., Moscow, Russia (
Figure 1). Spectra were accumulated over a 10-second duration at the ultrashort laser pulse repetition rate of 10 kHz, while for the measurement at different pulse energies the samples were moved in steps by 50 μm, using a motorized translation stage for micro-positioning (Standa, Vilnius, Lithuania) after each spectrum acquisition.
Peak power and energy of ultrashort laser pulses were identified as critical parameters for filamentation in synthetic diamond samples. The onset of visible asymmetric elongation of glowing filamentation channels towards laser radiation was observed as a function of increasing laser pulse energy (peak power), as previously reported [
21,
22]. This observation suggested the formation of a nonlinear focus beyond the Rayleigh length (linear focus parameter). For linearly polarized ultrashort laser pulses with a wavelength of 515 nm and varying pulse durations, the threshold energy values were discovered to be in the range of ≈210 - 230 nJ [
21,
22].
4. Discussion
Previously, such strongly broadened (supercontinuum, SC [
19]) and modulated spectra were already observed for ultrashort – picosecond and femtosecond – laser pulses in transparent condensed media in the filamentation regime [
30,
31]. The periodical modulation, maximal at the both or one of shoulders of the self-phase modulation broadened spectra, was assigned to constructive or destructive interference of the new generated, phase-shifted spectral components [
30].
In our study, surprisingly, no dominating shoulder peaks were observed, with the sideband intensities monotonously decreasing for the larger spectral shifts from the central laser wavelength, while the overall spectral width was monotonously increased versus laser pulse energy. As a result, we were tempted to consider an alternative model of the ultramodulation in the SC spectra presented below.
First, we considered the observed filamentation of the ultrashort laser pulses within the “moving-foci” approach [
19] (
Figure 6), rather than their filamentary channeling. However, the plasma channels produced at the high-NA focusing and multi-TW/cm
2 laser intensities, contained near-critical or even supercritical EHP with the local density
Ncr(1030 nm) ~ 1×10
21 cm
-3 at the 1030-nm wavelength and
Ncr(515 nm) ~ 4×10
22 cm
-3 at the 515-nm wavelength for the refractive index of the plasma in the diamond n* ~ 1 (would be n
2-fold higher in unexcited diamond for its refractive index n(515,1030nm) ≈ 2.4 [
32]), appearing during the self-focusing beam collapse as highly reflective, rather than refractive, like in gases [
19].
Meanwhile, the supercritical EHP density (N
e > N
cr, Re{ε∗} < 0) in the filamentary plasma channels is required for laser excitation of plasmon-polaritons at the plasma channel interface with the surrounding less excited dielectric, according to the threshold condition Re{ε∗} ≤ − Re{ε} [
33] (
Figure 7). Owing to the uncertainty of the self-focusing collapse arrest by the EHP regarding the final laser incidence angle on the plasma channel and other details, the corresponding plasmon-polaritonic dispersion curves were simulated for the filamentary channel with the supercritical EHP at different tentative self-focusing angles (0, 30, 45 and 60
0), using the common equations for the plasmon-polariton wavenumber
q(ω) [
34] under the crude approximation of a flat conductive interface and p-polarized laser radiation [
33,
35]
where the corresponding expressions for the imaginary magnitudes ε∗ and ε with Re{ε∗,ε} ~ Im{ε∗,ε}, representing the lossy EHP in the dielectric, see in [
33]. In fact, the exact solutions for excitation of longitudinal plasmon-polaritons in subwavelength metallic nanowire-channels (radius aω/c«1) correspond to qualitatively similar dispersion curves h(ω) of surfaces plasmon-polaritons (specifically, with similar asymptotic dispersion trends at h → 0 and h → ∞ [
34]), some of which are overviewed in the book [
35].
This means that non-linear spectral interactions (SC and harmonics generation etc.) could occur not only in the self-focusing beam, but also in the plasmonic field along the plasma channel, thus becoming considerably enhanced in longer filaments, alike in “channeling” filamentation regime.
Second, prompt interference between the self-focusing laser radiation and interfacial plasmon-polaritons on the plasma channel interface provides their longitudinal interferential electric field patterns along the channels. Their potential periods are 1/q
P ~λ/10 (plasmon resonance in
Figure 7) or ~λ (plasmon-polaritonic curves near the light-cone line ω=cq/n), similarly to laser generation of periodical surface structures (LIPSS) on material surfaces [
18,
36]. Hence, an additional periodical modulation of plasma density (δN
e) and refractive index (δn*) emerges in the periodical electric-filed “hot spots” in the filamentary plasma channels, which could work as a longitudinal Bragg grating to filter spectral components at the wavelengths
Λfilt [
37]
for the effective refractive index in the filament
In the plasma channel with the near-critical or supercritical EHP density n* < n, n*+δn* ≤ n* with 2n
eff < n, thus, in many dielectrics 2n
eff ~ 1, though variation of n* near the critical EHP density is sharp and extensive (from n* ~ n till n* « 1 in the supercritical EHP). Comparing to the expected plasmon-polariton wavelength for the flat metal/dielectric interface [
18,
36] in Eq.(2), one could find q
P ≈ 1/λ and thus
Λfilt ≈ λ. This is consistent with our observations of the strongly modulated SC spectra around the central laser wavelength, which upon normalization to the SC intensity could be represented as the common Bragg-grating transmittance spectra (
Figure 8). The multiple modulation in this case results not from multiple orders, but from the limited longitudinal number of grating stripes along the plasma channels.
Counterintuitively, according to the common diffraction theory, the higher pulse energies in
Figure 5 result in the longer filaments and, potentially, the larger longitudinal numbers of regularly-spaced grating stripes along the plasma channels and more pronounced central minimum with less distinct side minima. However, in our spectral experiments at the higher pulse energies (pulsewidth – 0.6 ps,
Figure 8) we observed the increasing number of nearly equal, closely-spaced spectral modulation features, thus potentially indicating longitudinal non-uniformity of the Bragg grating periods. This could be a first indirect indication of non-local – e.g., plasmon-polariton propagation – effects along the plasma channels, accumulating plasmon-polariton energy deposition downstream the plasma channels and vice versa making the filaments/channels longitudinally non-uniform too.