4School of Natural Sciences, Far Eastern Federal University (FEFU), Vladivostok 690041, Russia
5Institute of Automation and Control Processes, Far Eastern Branch, Russian Academy of Science, Vladivostok 690041, Russia
6Melbourne Centre for Nanofabrication, the Victorian Node of the Australian National Fabrication Facility, Clayton 3168 VIC, Australia
7College of Engineering, Chang Gung University , Guishan, Taoyuan 33302, Taiwan
8Department of Materials Science and Engineering, National Dong-Hwa University, Shoufeng, Hualien 97401, Taiwan
Associate Editor: P. Leiderer
Beilstein J. Nanotechnol. 2018, 9, 2609–2617. https://doi.org/10.3762/bjnano.9.242
Received 15 May 2018, Accepted 09 Sep 2018, Published 01 Oct 2018
Enhancement of X-ray emission was observed from a micro-jet of a nano-colloidal gold suspension in air under double-pulse excitation of ultrashort (40 fs) near-IR laser pulses. Temporal and spatial overlaps between the pre-pulse and the main pulse were optimized for the highest X-ray emission. The maximum X-ray intensity was obtained at a 1–7 ns delay of the main pulse irradiation after the pre-pulse irradiation with the micro-jet position shifted along the laser beam propagation. It was revealed that the volume around gold nanoparticles where the permittivity is near zero, ε ≈ 0, accounts for the strongest absorption, which leads to the effective enhancements of X-ray emission.
Keywords: double pulse; gold nanoparticles; intense femtosecond laser; plasma; water; X-ray
X-ray-related science and technology is one of the most powerful tools for material science, in particular regarding nanomaterials . Femtosecond laser-based X-ray pulse sources  have been available on the basis of linear- and nonlinear-optical processes such as plasma formation [3-5] and ablation . In practical applications for X-ray diffraction  and X-ray absorption fine structure (XAFS) measurements [8,9], or further nonlinear X-ray processes, a high flux of X-ray pulses is indispensable. X-ray intensity enhancements can be expected through an effective increase of the laser-absorption efficiency and an enhanced generation of high-energy electrons from a highly ionized plasma with high electron temperature [10,11] under intense laser irradiation.
For an efficient augmentation in X-ray intensity, femtosecond (fs)-laser excitation of metal nanoparticles and nanoclusters is of great interest, due to the efficient conversion of the absorbed energy into energetic ions, electrons and X-rays [12-14]. Plasmonic nanoparticles are expected to be highly useful for femtosecond laser-based X-ray emission due to high functionality, large absorption cross section and spectral selectivity based on surface plasmon resonance . An increase of the laser-absorption efficiency can be expected in plasmonic nanoparticles, which results in the efficient generation of highly ionized charge states [13,14]. It is expected that this characteristic interaction between intense femtosecond laser and plasmonic nanoparticles will contribute to a broad field of applications ranging from pulsed X-ray generation to energetic particle sources  or terahertz wave emission from aqueous solutions .
In recent experimental studies, intense fs-laser excitation of plasmonic gold nanoparticles with different sizes was applied for the efficient X-ray generation in an aqueous phase [17,18], and plasmonic silver nanoparticles were excited to accelerate electrons in vaccum . In addition, there is another new theoretical approach on X-ray sources using plasmonic graphene . Nanoparticles of silver  or aluminum  could be also considered. However, experimental approaches have been limited to flat samples  so far. Gold, with the higher atomic number and a high stability even in aqueous environments, should be considered a main target for fs-laser-induced X-ray experiments especially in aqueous solution.
Recently, it was reported that fs-double-pulsed excitation (including pump–probe type experiments) with dielectrics amplifies optical emission . This double-pulsed excitation is considered to be another way to efficiently augment the X-ray intensity . A significant enhancement in X-ray intensity from a water film was observed under double-pulse excitation with a pre-pulse ahead of the main pulse [26,27]. The pre-pulse irradiation with relatively low intensity induces a time-dependent laser-induced plasma (in the range of picoseconds) and ablation (in the range of nanoseconds). At a certain time delay, Δt, after the pre-pulse irradiation, the main-pulse irradiation of the solution surface results in enhanced X-ray generation. Investigation of X-ray generation intensity and spectrum in distilled water revealed a significant variation of the optimum position of the solution surface for different Δt up to 15 ns . Analyses of the relation between the time-dependent X-ray generation and the optimum laser irradiation position are still missing for other aqueous solutions, especially for experiments under ambient pressure. Previous studies, mostly with solid targets in vacuum, using metallic targets at high intensities of more than 1015 W/cm2 have shown strong X-ray generation enhancement when a pre-pulse was used. This was explained through the bremsstrahlung of hot electrons and plasma waves breaking at the peak X-ray intensities [29-31].
Here, we explore the optimization of enhancement in X-ray generation from an aqueous solution flow of plasmonic gold nanoparticles under fs-double-pulse excitation in air. The optimization of the solution film surface position with respect to the focal position of irradiation was systematically investigated for the most intense generation of hard X-rays.
Colloidal suspensions of gold nanospheres were prepared via kinetically controlled seeded-growth synthesis of citrate-stabilized nanoparticles [32-34]. In a 250 mL three-necked round-bottomed flask, 2.2 mM sodium citrate in deionized water (>18 MΩ, 150 mL) was heated at 115 °C for 15 min under vigorous stirring. After it reached the boiling point, 1 mL of HAuCl4(aq) (25 mM) was added. The color of the solution changed from yellow to bluish gray and to light pink over a period of 10 min. Immediately after the synthesis of the gold seed solution, the temperature was cooled down to 90 °C and the seeded growth of gold nanospheres was carried out. 1 mL of HAuCl4(aq) solution (25 mM) was injected in the reaction vessel. The reaction was finished after 30 min and the process was repeated twice. The sample was diluted by extracting 55 mL of the sample, adding 53 mL of deionized water and 2 mL of 60 mM sodium citrate. By changing the volume in each growth step, it is possible to control the particle concentration. The diameter of gold nanoparticles was 20 nm and the concentration was ca. 1.4 × 10−7 mol/cm3 or ρ ≈ 3.5 × 1011 nanoparticles/cm3. These values correspond to an intra-colloidal separation of = 1.4 μm or 1.07 × 102 nanoparticles in the focal volume approximated by a cylinder of volume πw2·2zR.
The micro-jet positioning system has been introduced in our earlier studies [28,35]. A beam-steering system for the pre-pulse irradiation was additionally installed to the setup used in this study. Figure 1 shows the principle and the optimization map in the lateral plane perpendicular to the laser beam. A precision of 1 μm for the axial displacement of the solution film with respect to the incoming laser beam was achieved within less than 1 s, which is satisfactory for practical applications of X-ray generation .
Transform-limited fs-laser pulses tp = 40 fs with central wavelength λ = 800 nm at 1 kHz repetition rate were horizontally polarized (Mantis, Legend Elite, HE USP, Coherent, Inc.). Two pulses were combined collinearly for the double-pulse irradiation using half-wave plates (65-906, Edmund Optics) and polarization beam splitters (47-048, Edmund Optics). Vertically and horizontally polarized pulses are defined as the pre-pulse and the main pulse, respectively. The vertically polarized pre-pulse intensity is below the threshold of X-ray generation, which is obtained only under double-pulse irradiation. The delay of the main pulse relative to the pre-pulse, Δt, was controlled within three time zones: 0–5 ns (I), 5–10 ns (II) and 10–15 ns (III), which were set into the optical path of the pre-pulse. The optical delay for each time zone is obtained with an 80 cm long variable delay line (SGSP 46-800, Sigma Koki). Change of the time zones from I to II (or III) was made with two independent fixed retroreflectors. The optical path of the main pulse had the same optical length as the longest optical path of the pre-pulse. After combination, the pulses were tightly focused onto a film from the solution jet with an off-axis parabolic mirror (effective focus length of f = 50.8 mm, reflection angle of 90°, and numerical aperture NA =0.25, 47-097, Edmund Optics) in air. Under these focusing conditions, the theoretical/geometrical diameter of the focal spot is 2w0 = 1.22λ/NA = 3.9 μm and the axial extent along the propagation estimated as two Rayleigh lengths 2zR = 2nλ/(NA2) = 25.5 μm where n = 1 is the refractive index of air. The focal region is longer than the spatial extent of the pulse ctp = 12 μm, c being the speed of light. A 17 ± 2 μm thick jet-flow of a gold nanocolloidal suspension was used as a target. The incident angle of the laser was set at θ = 60° to the normal of the surface of the solution flow. This corresponded to a twice larger effective thickness of the solution film h/cosθ ≡ 2h encountered by the laser beam. Usually, more than 120 mL of sample solution was being circulated by a pump during experiments. This amount is much larger than the volume of the laser irradiation at the focus, which can be estimated to be 1.5 × 10−4 mL. An apparent degradation of the solution absorption spectrum was not observed after experiments. It has been already reported, under similar experimental conditions, that the energy conversion efficiency of the laser pulse to hard X-rays, the electron temperatures, and X-ray photon numbers are estimated to be about 10−8 , 3–10 keV , and 3.8 × 1010 photons s−1 4πsr−1 , respectively.
A strong generation of X-rays is related to a large amount of absorbed energy and a high temperature [3,4]. In terms of relative permittivity (dielectric constant), εr, a transparent medium (water) with colloidal gold nanoparticles represents an effective medium with an overall positive real part of the permittivity. Upon excitation, the material is approaching dielectric breakdown and a state of εr ≈ 0 (epsilon-near-zero (ENZ) state; ). Under the ENZ condition, the strongest absorption takes place [36,37]:
Here, εi is the intrinsic permittivity of material, ν is the electron–ion collision rate, ω is the optical cyclic frequency at the wavelength of excitation, ne,a,cr are the electron, atom and critical densities, respectively. The imaginary part of the permittivity is given by:
The threshold fluence for water to reach the ENZ state can be calculated from the consideration that all absorbed energy density is converted to thermal energy of electrons. It can be calculated from the ablation threshold expression of a dielectric :
where ls = c/(κω) is the skin depth related to the imaginary part of the refractive index , c is the speed of light in vacuum, A is the absorbance, εb is the binding energy [eV/atom], and Ji is the ionization potential. The threshold fluence for the creation of the ENZ state is calculated for full ionization, i.e., with the bandgap energy Δg in Equation 3 instead of (εb + Ji). For water, Δg = 9.5 eV. However, the creation of a solvated electron out of the valence band requires only 6.2 eV, which was used in the estimations . The ENZ state is created when the electron density nENZ ≈ εincr (Equation 1); for water εi = n2 = 1.77, defined by the refractive index n = 1.33 at the excitation wavelength λ = 800 nm. The critical plasma density is ncr = ω2mem*ε0/e2, where me is the electron mass, m* = 1 is its effective mass in the plasma, e is the electron charge, ε0 is the permittivity of vacuum, and ω = 2πc/λ is the cyclic frequency of light. This expression can be used to calculate the corresponding plasma frequency ωp when the electron density, ne, is known. For instance, for gold ne = 5.9 × 1022 cm−3.
The refractive index under ENZ conditions can be calculated as with use of Equation 2:
Then AENZ = 4κ/[(κ + 1)2 + κ2] and the fluence threshold for the ENZ state is (Equation 3):
with Δg = 6.2 eV for a solvated electron and an electron collision rate for momentum exchange ν = 1015 s−1. The following numerical values are valid for the ENZ state in water: nENZ = κENZ = 0.86, AENZ = 0.82, = 3.1 × 1021 cm−3, = 42 mJ/cm2. This corresponds to the ionization threshold for water through the creation of solvated electrons. When the bandgap energy of water Δg = 9.5 eV is taken as the ionization threshold, a slightly larger fluence of = 60 mJ/cm2 will result.
For the metal (gold nanoparticles) we consider the ablation threshold expression equivalent to Equation 3 in which binding energy and electron work function (escape energy) need to be accounted for through (εb + εesc) together with the substitution A/ls = 4π/λ. For the ionization threshold, the electron work function is considered. Then, the ionization threshold reads :
For gold the following values are valid: εesc = 5.1 eV and = 230 mJ/cm2.
The above presented analysis takes into account the host material and the electron–ion plasma through the definition of permittivity, where the ENZ region defines the strongest light absorption. Due to the fast interaction, the permittivity conditions corresponding to the ENZ are changing in time and space.
Figure 2 shows the X-ray generation as a function of the time between the pre-pulse of 80 μJ and the main pulse of 700 μJ. Strong X-ray generation was observed at a delay of 1 ns in the case of the gold colloid. For the distilled water film, the required time was ca. 8 ns. Comparison of X-ray generation from water and colloidal gold solution at a time delay of Δt = 6 ns is shown in Figure 3. For the nanocolloidal gold suspension, a film shifted downwards by ca. 60 μm along the propagation of the laser beam (the negative side of the z-axis shown in Figure 3b) caused the strongest X-ray generation, which was approximately 5 to 6 times larger that during single-pulse exposure. This shift reflects the time integrated X-ray signal from the region of the strongest absorption due to ENZ conditions multiplied with the mass density of the solution. The shift of the solution film position was only ca. 10 μm (half of the film thickness) for distilled water. The strongest X-ray generation from pure water was about a third for the double-pulse irradiation. A long delay of 1–10 ns for intense X-rays from the colloidal solution film is consistent with the formation and growth of bubbles around gold nanoparticles and their ablation as discussed in the following.
The absorption of laser energy by the gold nanoparticles is expected to be relatively small considering their volume fraction and off-resonance excitation conditions. The plasmon resonance excitation for these NPs is expected to be near 520 nm. This provides a moderate enhancement of the optical near-fields as shown further by numerical simulations. The incident fluence of the main pulse used in the present study significantly exceeds both the ablation threshold for bulk gold  as well as the fragmentation thresholds for gold colloidal nanoparticles in water [40,41]. It should be stressed out that the increase of the effective size of the laser spot due to filamentation can significantly reduce the incident fluence [3,4]. However, it is reasonable to assume that under our experimental conditions the fragmentation threshold will be easily reached resulting in breaking up of the gold nanoparticles into smaller clusters. Fragmentation can proceed via explosive ablation at the characteristic time scales from few tens to few hundred of picoseconds [12,42,43]. The electrostatic scenario of rupture of the particles by the repulsive Coulomb forces resulting from ejection of the hot electrons absorbing the laser irradiation was shown to dominate for particles smaller than 60 nm .
In this study, the laser fluence reaching the gold nanoparticles in water is much higher than 0.2 J/cm2, which is also close to the threshold of melt-reshaping of nanoparticles and bubble initiation . Air breakdown at the close-to-film regions and filamentation of the laser pulse takes place, which obscures exact knowledge of the laser fluence at the nanoparticle location. The density of nanoparticles in the solution corresponds to the situation when a pressure build-up around the neighboring nanobubbles hinders their growth as observed from acoustic studies . The 2 ns transient was observed in photo-acoustic excitation experiments using the same density of nanoparticles in solution under two-pulse irradiation, while the ebergy off the main pulse was only 100 μJ , comparable with the energy used for the pre-pulse in this study.
It is reasonable to assume, that the maximum intensity of X-ray emission occurs when the temporal and spatial excitation of film solutions maximizes the volume of εr-near-zero (ENZ) regions and overlaps with the focal volume of the main pulse. Those ENZ regions are responsible for the most efficient energy deposition via absorption and define the transition of the propagating light to the ENZ surface (an optical near-field). This localization of light is affected by a complex pattern of surface nanoroughening of the liquid surface  and an internal pattern of nanobubbles. Figure 4a shows the temporal evolution of the position for maximum X-ray emission of the distilled water film consistent with a roughening of the liquid surface , where a positive z-shift corresponds to a better spatial overlap of the laser focus and the excited ENZ region on the film. A higher energy pre-pulse caused a faster evolution of maximum X-ray intensity at smaller z-shifts. The peak intensity of X-rays in water was on a picosecond time scale (Figure 4b) consistent with thermalization of electrons and ions .
The formation of surface waves via capillary driven instabilities occurring at longer times (nanoseconds) is affected by the linear momentum transferred to the film surface, which as be calculated from the force F = Δp/Δt where the momentum in a medium of refractive index n of a laser pulse with energy E is p = nE/c. When a light pulse of power P traverses from a low refractive index (n0 = 1, air) into a medium with larger refractive index (n1 = 1.33, water), momentum is gained and the force acting onto the front surface is :
where R1 = (n1 − n0)2/(n1 + n0)2 = 0.02 is the reflection coefficient. One finds for a pulse of Ep = 100 μJ of tp = 40 fs that Ffr = −2.32 N (along the z-axis), where the sign indicates that the force acting onto the surface is in the opposite direction to the laser pulse propagation (opposite to the z-axis in this study). The force acting on the rear surface of the solution jet can be estimated from:
which is Fb = 3.14 N, i.e., a pushing force along the light propagation direction. Considering a breakup of the incoming high intensity laser pulse into separate microfilaments, the transfer of momentum into the solution film surface should be one of the mechanisms [3,4] resulting in the formation of surface nanoroughening and surface waves leading to droplet formation at longer times. Different scenarios are encountered by the main pulse depending on the pre-pulse excitation of the solution film. When the pre-pulse energy increases, the main pulse encounters a strongly excited ENZ shell around the gold nanoparticles as well as, at longer delays, surface roughening of the solution film and cavitation bubbles below the surface . Such scenarios are simulated in the following using the finite-difference time-domain (FDTD) method (Lumerical).
Figure 5 shows the light-intensity distributions around gold nanoparticles modeling different stages of excitation from initial gold nanoparticles towards the formation of ENZ plasma around them. The absorption (skin) depth is defined by the imaginary part of permittivity and the ENZ conditions is ls = λ/(2πκ) ≈ 149 nm. A large volume in which the light intensity is enhanced by up to E2 ≈ 4-times is clearly distinguished. In contrast to the case without excitation (Figure 5a), the gold nanoparticles are coated with a 20 nm wide shell of ENZ material after excitation (Figure 5b). The formation of nanoscale roughening of the liquid surface occurrs at a later stage and affects the light intensity distribution  in subsurface regions of the colloidal suspension film (Figure 5c). The shell of ENZ material of radius equal to the skin depth ls shows that light is mostly reflected at the rim of this volume and is not reaching the interior of the volume (as expected). These FDTD simulations illustrate only qualitatively what the main pulse encounters in a solution film of colloidal gold nanoparticles and that a strong light-field enhancement of about four-times is expected at the gold–water interface, where the ENZ region is created.
New possibilities can be explored in applications by using the proposed hard X-ray generation under ambient conditions. Creation and localization in time and space of the high-temperature plasma as X-ray emitter is the key challenge. The most efficient energy deposition takes place in the ENZ regions. This is consistent with a high-energy plasma driven by a resonant absorption at near-critical plasma density regions [3,4,29,30]. The pump–probe concept leads to the evolution of a breakdown plasma, which is required for complex targets such as nanoparticle jets. The qualitative analysis outlined here needs detailed knowledge of the light filamentation and propagation in the regions of air breakdown close to the jet to better estimate the spatial and temporal evolution of the ENZ regions.
A strong enhancement of X-ray emission from a gold nanocolloidal suspension was achieved for double-pulse excitation and could be useful for the practical implementation of X-ray generation under ambient condition. High-precision beam steering allows for a stable tracking and X-ray intensity in real time. The enhancement of X-ray emission was discussed in terms of an optimized overlap of the focal region of the laser and the solution flow, which resulted in a nanorough surface and the generation of nanobubbles leading to a strong localization (ENZ regions) and absorption of light. The fast generation (ca. 1 ns) of X-rays from the gold nanocolloidal suspension compared to distilled water (7–8 ns) was observed. The presented analysis of energy deposition into the ENZ regions is important for a detailed description of X-ray generation as well as in other fs-laser processes . In this paper, spherical gold nanoparticles with a diameter of about 20 nm and a concentration of ca. 3.5 × 1011 nanoparticles/cm3 were used. Single-pulse excitation for X-ray emission  was carried out with spherical gold nanoparticles with different diameters from 10 to 100 nm and the results show that the appropriate size for the highest X-ray intensity is 40–50 nm. For ultrasound generation under fs-laser excitation [34,49], gold nanorod particles with more efficient surface plasmon resonance effects  were also used. Further enhancements of X-ray intensities are expected under double-pulse excitation with different sizes, shapes, or concentrations of nanoparticles.
It is clearly shown that the controlled spatial positioning of the solution flow at the laser focus with appropriate delay time induces effective X-ray emission. Recently, similar experiments with narrowly focussed fs-laser pulses for terahertz-wave emission from water flow were reported , while terahertz emission from various solutions under broad focusing was already known . Simultaneous emission of X-ray and terahertz waves from gas clusters in vacuum [53,54], Al-coated glass substrates in vacuum , or a distilled water flow in air under double-pulse excitation  were independently reported. A simultaneous emission of X-ray and terahertz wave radiation is highly expected for novel pump–probe experiments . The knowledge obtained in this paper will contribute to further advancements in this field.
T.Y. thanks the partial financial support from Kakenhi (B, 18H01820) from JSPS, Japan. S.J. is grateful for partial support via the ARC Discovery DP170100131 grant, A.A.K. for the RSF no.17-19-01325 grant. The authors acknowledge discussion of light-matter interaction aspects with Prof. Eugene G. Gamaly. FDTD simulations were performed on the swinSTAR supercomputer at Swinburne University of Technology.