Enhancement of X-ray emission from nanocolloidal gold suspensions under double-pulse excitation

  1. Wei-Hung Hsu1,
  2. Frances Camille P. Masim1,
  3. Armandas Balčytis2,
  4. Hsin-Hui Huang1,
  5. Tetsu Yonezawa3ORCID Logo,
  6. Aleksandr A. Kuchmizhak4,5,
  7. Saulius Juodkazis2,6 and
  8. Koji Hatanaka1,7,8ORCID Logo

1Research Center for Applied Sciences, Academia Sinica, Nankang, Taipei 11529, Taiwan
2Centre for Micro-Photonics, Swinburne University of Technology, Hawthorn, VIC 3122, Australia
3Division of Materials Science and Engineering, Faculty of Engineering, Hokkaido University, Sapporo, Hokkaido 0608628, Japan

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

  1. Corresponding author email

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 [1]. Femtosecond laser-based X-ray pulse sources [2] have been available on the basis of linear- and nonlinear-optical processes such as plasma formation [3-5] and ablation [6]. In practical applications for X-ray diffraction [7] 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 [15]. 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 [14] or terahertz wave emission from aqueous solutions [16].

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 [19]. In addition, there is another new theoretical approach on X-ray sources using plasmonic graphene [20]. Nanoparticles of silver [21] or aluminum [22] could be also considered. However, experimental approaches have been limited to flat samples [23] 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 [24]. This double-pulsed excitation is considered to be another way to efficiently augment the X-ray intensity [25]. 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 [28]. 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 [Graphic 1] = 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 [35].


Figure 1: (a) Schematics of experiment showing polarizations of the pre-pulse E1 and the main pulse E2 in the plane of incidence (xz-plane). (b) The absorption spectrum of the sample gold solution with an SEM image of gold nanoparticles. (c) Beam steering optimization in xy-projection for the highest yield of X-ray generation at 100 ps delay between the pre-pulse (energy E1 = 80 μJ) and the main pulse (E2 = 700 μJ). The distance between water film and Geiger counter was 15 cm, with an aperture of 2 mm. The orientation of the E1 polarization is shown for reference.

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 [25], 3–10 keV [27], and 3.8 × 1010 photons s−1 4πsr−1 [26], respectively.

Theory: Epsilon-Near-Zero (ENZ)-Material

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; [Graphic 2]). 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 [37]:


where ls = c/(κω) is the skin depth related to the imaginary part of the refractive index [Graphic 3], 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 [38]. 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 = ω2mem0/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 [Graphic 4] with use of Equation 2:

[Graphic 5]

Then AENZ = 4κ/[(κ + 1)2 + κ2] and the fluence threshold for the ENZ state is (Equation 3):

[Graphic 6]

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, [Graphic 7] = 3.1 × 1021 cm−3, [Graphic 8] = 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 [Graphic 9] = 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 [37]:

[Graphic 10]

For gold the following values are valid: εesc = 5.1 eV and [Graphic 11] = 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.

Results and Discussion

Two-pulse excitation for X-ray emission

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.


Figure 2: X-ray intensities from (a) a water film and (b) a film of a colloidal suspension of gold nanospheres irradiated by a pre-pulse at 80 μJ/pulse and a main pulse at 700 μJ/pulse with different time delays. The automatic positioning system was used to find the maximum X-ray intensity [35]. The laser irradiates the solution film from the right side in these figures. The width of the solution film hit by the laser pulse at an incidence angle of θ = 60° is 2h ≈ 40 μm (shown in panel (a)).


Figure 3: X-ray intensity at different positions of the water film (a) and the gold nanosphere colloidal solution film (b) at Δt = 6 ns. The laser irradiates the solution film from the right side in these figures. The fit constitutes a multi-Gaussian convolution. The intensity of the pre-pulse is 80 μJ. Under double-pulsed excitation, the spatial profiles of X-ray intensity are significantly different in the cases of water and the gold nanospheres colloidal solution.

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 [39] 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 [12].

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 [44]. 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 [45]. 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 [34], 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 [28] 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 [28], 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 [6].


Figure 4: Position of X-ray maximum generation (a) and its intensity (b) for distilled water films at short excitation times and at different pre-pulse energies E1, while the main pulse energy was fixed at E2 = 700 μJ. A positive z-shift corresponds to an upward movement of solution position towards the incoming laser pulse (see Figure 2a). The lines are guides to the eye.

Ponderomotive action and light field enhancement

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 = Δpt where the momentum in a medium of refractive index n of a laser pulse with energy E is p = nE/c[46]. 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 [46]:


where R1 = (n1n0)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 [47]. 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 [28] 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.


Figure 5: FDTD simulations of light-field enhancement under the conditions of the experiment: angle of incidence 60°, 800 nm wavelength, the concentration of gold nanoparticle yields a separation of ca. 1 μm between them. (a) Interaction of a water film with 20 nm-diameter gold nanoparticles and a laser pulse (modeled as a plane wave incidence). (b) The formation of a plasma around nanoparticles with permittivity near zero ε ≈ 0.86 + i0.86 (ENZ) in which light is efficiently absorbed: 20 nm wide shell around a 20 nm diameter nanoparticle. (c) Sphere of strongly ionized water around a gold nanoparticle. The absorption depth under ENZ plasma conditions is ls = 149 nm, which is taken as radius of the excited volume. Nanoroughening of water surface evolves at longer times. Disintegration of gold nanoparticles is expected inside the cavitation bubble. Intensity of incident beam: [Graphic 12].

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 [48]. 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 [17] 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 [50] 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 [51], while terahertz emission from various solutions under broad focusing was already known [52]. Simultaneous emission of X-ray and terahertz waves from gas clusters in vacuum [53,54], Al-coated glass substrates in vacuum [54], or a distilled water flow in air under double-pulse excitation [16] were independently reported. A simultaneous emission of X-ray and terahertz wave radiation is highly expected for novel pump–probe experiments [55]. 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.


  1. Sakdinawat, A.; Attwood, D. Nat. Photonics 2010, 4, 840–848. doi:10.1038/nphoton.2010.267
    Return to citation in text: [1]
  2. Mondal, S.; Chakraborty, I.; Ahmad, S.; Carvalho, D.; Singh, P.; Lad, A. D.; Narayanan, V.; Ayyub, P.; Kumar, G. R.; Zheng, J.; Sheng, Z. M. Phys. Rev. B 2011, 83, 035408. doi:10.1103/physrevb.83.035408
    Return to citation in text: [1]
  3. Kruer, W. L. The Physics of Laser-Plasma Interaction; Addison-Wesley: Boston, MA, U.S.A., 1988; pp 1–37. doi:10.1201/9780429502798
    Return to citation in text: [1] [2] [3] [4] [5]
  4. Kruer, W. L. The Physics of Laser-Plasma Interaction; Addison-Wesley: Boston, MA, U.S.A., 1988; pp 153–177. doi:10.1201/9780429502798
    Return to citation in text: [1] [2] [3] [4] [5]
  5. Liu, X.-L.; Lu, X.; Liu, X.; Xi, T.-T.; Liu, F.; Ma, J.-L.; Zhang, J. Opt. Express 2010, 18, 26007–26017. doi:10.1364/oe.18.026007
    Return to citation in text: [1]
  6. Gamaly, E. G. Phys. Rep. 2011, 508, 91–243. doi:10.1016/j.physrep.2011.07.002
    Return to citation in text: [1] [2]
  7. Bargheer, M.; Zhavoronkov, N.; Gritsai, Y.; Woo, J. C.; Kim, D. S.; Woerner, M.; Elsaesser, T. Science 2004, 306, 1771–1773. doi:10.1126/science.1104739
    Return to citation in text: [1]
  8. Lee, T.; Jiang, Y.; Rose-Petruck, C. G.; Benesch, F. J. Chem. Phys. 2005, 122, 084506. doi:10.1063/1.1852455
    Return to citation in text: [1]
  9. Hatanaka, K.; Fukumura, H. X-Ray Spectrom. 2012, 41, 195–200. doi:10.1002/xrs.2404
    Return to citation in text: [1]
  10. Bulanov, S. V.; Kovrizhnykh, L. M.; Sakharov, A. S. Phys. Rep. 1990, 186, 1. doi:10.1016/0370-1573(90)90002-j
    Return to citation in text: [1]
  11. Stearns, D. G.; Landen, O. L.; Campbell, E. M.; Scofield, J. H. Phys. Rev. A 1988, 37, 1684–1690. doi:10.1103/physreva.37.1684
    Return to citation in text: [1]
  12. Delfour, L.; Itina, T. E. J. Phys. Chem. C 2015, 119, 13893–13900. doi:10.1021/acs.jpcc.5b02084
    Return to citation in text: [1] [2] [3]
  13. Bauer, D.; Brabec, T.; Fehske, H.; Lochbrunner, S.; Meiwes-Broer, K.-H.; Redmer, R. New J. Phys. 2013, 15, 065015. doi:10.1088/1367-2630/15/6/065015
    Return to citation in text: [1] [2]
  14. Wopperer, P.; Dinh, P. M.; Reinhard, P.-G.; Suraud, E. Phys. Rep. 2015, 562, 1–68. doi:10.1016/j.physrep.2014.07.003
    Return to citation in text: [1] [2] [3]
  15. Link, S.; El-Sayed, M. A. J. Phys. Chem. B 1999, 103, 8410–8426. doi:10.1021/jp9917648
    Return to citation in text: [1]
  16. Huang, H.-h.; Nagashima, T.; Hsu, W.-h.; Juodkazis, S.; Hatanaka, K. Nanomaterials 2018, 8, 523–534. doi:10.3390/nano8070523
    Return to citation in text: [1] [2]
  17. Masim, F. C. P.; Porta, M.; Hsu, W.-H.; Nguyen, M. T.; Yonezawa, T.; Balčytis, A.; Juodkazis, S.; Hatanaka, K. ACS Photonics 2016, 3, 2184–2190. doi:10.1021/acsphotonics.6b00692
    Return to citation in text: [1] [2]
  18. Hatanaka, K.; Porta, M.; Masim, F. C. P.; Hsu, W.-H.; Nguyen, M. T.; Yonezawa, T.; Balčytis, A.; Juodkazis, S. arXiv 2016, No. 1604.07541.
    Return to citation in text: [1]
  19. Passig, J.; Irsig, R.; Truong, N. X.; Fennel, T.; Tiggesbäumker, J.; Meiwes-Broer, K. H. New J. Phys. 2012, 14, 085020. doi:10.1088/1367-2630/14/8/085020
    Return to citation in text: [1]
  20. Wong, L. J.; Kaminer, I.; Ilic, O.; Joannopoulos, J. D.; Soljačić, M. Nat. Photonics 2016, 10, 46–52. doi:10.1038/nphoton.2015.223
    Return to citation in text: [1]
  21. Grillet, N.; Manchon, D.; Cottancin, E.; Bertorelle, F.; Bonnet, C.; Broyer, M.; Lermé, J.; Pellarin, M. J. Phys. Chem. C 2013, 117, 2274–2282. doi:10.1021/jp311502h
    Return to citation in text: [1]
  22. Knight, M. W.; King, N. S.; Liu, L.; Everitt, H. O.; Nordlander, P.; Halas, N. J. ACS Nano 2014, 8, 834–840. doi:10.1021/nn405495q
    Return to citation in text: [1]
  23. Rousse, A.; Audebert, P.; Geindre, J. P.; Falliès, F.; Gauthier, J. C.; Mysyrowicz, A.; Grillon, G.; Antonetti, A. Phys. Rev. E 1994, 50, 2200–2207. doi:10.1103/physreve.50.2200
    Return to citation in text: [1]
  24. Winkler, T.; Haahr-Lillevang, L.; Sarpe, C.; Zielinski, B.; Götte, N.; Senftleben, A.; Balling, P.; Baumert, T. Nat. Phys. 2017, 14, 74–79. doi:10.1038/nphys4265
    Return to citation in text: [1]
  25. Hatanaka, K.; Miura, T.; Fukumura, H. Appl. Phys. Lett. 2002, 80, 3925–3927. doi:10.1063/1.1482135
    Return to citation in text: [1] [2]
  26. Hatanaka, K.; Ono, H.; Fukumura, H. Appl. Phys. Lett. 2008, 93, 064103. doi:10.1063/1.2967882
    Return to citation in text: [1] [2]
  27. Hatanaka, K.; Miura, T.; Fukumura, H. Chem. Phys. 2004, 299, 265–270. doi:10.1016/j.chemphys.2003.10.036
    Return to citation in text: [1] [2]
  28. Hsu, W.-H.; Masim, F. C. P.; Balčytis, A.; Juodkazis, S.; Hatanaka, K. Opt. Express 2017, 25, 24109–24118. doi:10.1364/oe.25.024109
    Return to citation in text: [1] [2] [3] [4] [5]
  29. Sandhu, A. S.; Kumar, G. R.; Sengupta, S.; Das, A.; Kaw, P. K. Phys. Rev. Lett. 2005, 95, 025005. doi:10.1103/physrevlett.95.025005
    Return to citation in text: [1] [2]
  30. Eder, D. C.; Pretzler, G.; Fill, E.; Eidmann, K.; Saemann, A. Appl. Phys. B: Lasers Opt. 2000, 70, 211–217. doi:10.1007/s003400050034
    Return to citation in text: [1] [2]
  31. Pisani, F.; Andiel, U.; Eidmann, K.; Witte, K.; Uschmann, I.; Morak, A.; Förster, E.; Sauerbrey, R. Appl. Phys. Lett. 2004, 84, 2772–2774. doi:10.1063/1.1695634
    Return to citation in text: [1]
  32. Bastús, N. G.; Comenge, J.; Puntes, V. Langmuir 2011, 27, 11098–11105. doi:10.1021/la201938u
    Return to citation in text: [1]
  33. Jana, N. R.; Gearheart, L.; Murphy, C. J. Adv. Mater. 2001, 13, 1389–1393. doi:10.1002/1521-4095(200109)13:18<1389::aid-adma1389>3.0.co;2-f
    Return to citation in text: [1]
  34. Masim, F. C. P.; Hsu, W.-H.; Liu, H.-L.; Yonezawa, T.; Balčytis, A.; Juodkazis, S.; Hatanaka, K. Opt. Express 2017, 25, 19497–19507. doi:10.1364/oe.25.019497
    Return to citation in text: [1] [2] [3]
  35. Hsu, W.-H.; Masim, F. C. P.; Porta, M.; Nguyen, M. T.; Yonezawa, T.; Balčytis, A.; Wang, X.; Rosa, L.; Juodkazis, S.; Hatanaka, K. Opt. Express 2016, 24, 19994. doi:10.1364/oe.24.019994
    Return to citation in text: [1] [2] [3]
  36. Gamaly, E. G.; Rode, A. V. Appl. Phys. A: Mater. Sci. Process. 2018, 124, 278. doi:10.1007/s00339-018-1693-3
    Return to citation in text: [1]
  37. Gamaly, E. G.; Rode, A. V.; Luther-Davies, B.; Tikhonchuk, V. T. Phys. Plasmas 2002, 9, 949–957. doi:10.1063/1.1447555
    Return to citation in text: [1] [2] [3]
  38. Linz, N.; Freidank, S.; Liang, X.-X.; Vogelmann, H.; Trickl, T.; Vogel, A. Phys. Rev. B 2015, 91, 134114. doi:10.1103/physrevb.91.134114
    Return to citation in text: [1]
  39. Wang, X. W.; Kuchmizhak, A. A.; Li, X.; Juodkazis, S.; Vitrik, O. B.; Kulchin, Y. N.; Zhakhovsky, V. V.; Danilov, P. A.; Ionin, A. A.; Kudryashov, S. I.; Rudenko, A. A.; Inogamov, N. A. Phys. Rev. Appl. 2017, 8, 044016. doi:10.1103/physrevapplied.8.044016
    Return to citation in text: [1]
  40. Baumgart, J.; Humbert, L.; Boulais, É.; Lachaine, R.; Lebrun, J.-J.; Meunier, M. Biomaterials 2012, 33, 2345–2350. doi:10.1016/j.biomaterials.2011.11.062
    Return to citation in text: [1]
  41. Boulais, É.; Lachaine, R.; Meunier, M. Nano Lett. 2012, 12, 4763–4769. doi:10.1021/nl302200w
    Return to citation in text: [1]
  42. Muto, H.; Miyajima, K.; Mafuné, F. J. Phys. Chem. C 2008, 112, 5810–5815. doi:10.1021/jp711353m
    Return to citation in text: [1]
  43. Pyatenko, A.; Yamaguchi, M.; Suzuki, M. J. Phys. Chem. C 2009, 113, 9078–9085. doi:10.1021/jp808300q
    Return to citation in text: [1]
  44. Boutopoulos, C.; Hatef, A.; Fortin-Deschênes, M.; Meunier, M. Nanoscale 2015, 7, 11758–11765. doi:10.1039/c5nr02721g
    Return to citation in text: [1]
  45. Nakajima, T.; Wang, X.; Chatterjee, S.; Sakka, T. Sci. Rep. 2016, 6, 28667. doi:10.1038/srep28667
    Return to citation in text: [1]
  46. Guck, J.; Ananthakrishnan, R.; Mahmood, H.; Moon, T. J.; Cunningham, C. C.; Käs, J. Biophys. J. 2001, 81, 767–784. doi:10.1016/s0006-3495(01)75740-2
    Return to citation in text: [1] [2]
  47. Shih, C.-Y.; Shugaev, M. V.; Wu, C.; Zhigilei, L. V. J. Phys. Chem. C 2017, 121, 16549–16567. doi:10.1021/acs.jpcc.7b02301
    Return to citation in text: [1]
  48. Wang, L.; Chen, Q.-D.; Cao, X.-W.; Buividas, R.; Wang, X.; Juodkazis, S.; Sun, H.-B. Light: Sci. Appl. 2017, 6, e17112. doi:10.1038/lsa.2017.112
    Return to citation in text: [1]
  49. Masim, F. C. P.; Liu, H.-L.; Porta, M.; Yonezawa, T.; Balčytis, A.; Juodkazis, S.; Hsu, W.-H.; Hatanaka, K. Opt. Express 2016, 24, 14781. doi:10.1364/oe.24.014781
    Return to citation in text: [1]
  50. Huang, H. J.; Yu, C.-p.; Chang, H. C.; Chiu, K. P.; Ming Chen, H.; Liu, R. S.; Tsai, D. P. Opt. Express 2007, 15, 7132–7139. doi:10.1364/oe.15.007132
    Return to citation in text: [1]
  51. Jin, Q.; E, Y.; Williams, K.; Dai, J.; Zhang, X.-C. Appl. Phys. Lett. 2017, 111, 071103. doi:10.1063/1.4990824
    Return to citation in text: [1]
  52. Dey, I.; Jana, K.; Fedorov, V. Y.; Koulouklidis, A. D.; Mondal, A.; Shaikh, M.; Sarkar, D.; Lad, A. D.; Tzortzakis, S.; Couairon, A.; Kumar, G. R. Nat. Commun. 2017, 8, 1184. doi:10.1038/s41467-017-01382-x
    Return to citation in text: [1]
  53. Balakin, A. V.; Dzhidzhoev, M. S.; Gordienko, V. M.; Esaulkov, M. N.; Zhvaniya, I. A.; Ivanov, K. A.; Kotelnikov, I. A.; Kuzechkin, N. A.; Ozheredov, I. A.; Panchenko, V. Y.; Savel’ev, A. B.; Smirnov, M. B.; Solyankin, P. M.; Shkurinov, A. P. IEEE Trans. Terahertz Sci. Technol. 2017, 7, 70–79. doi:10.1109/tthz.2016.2617867
    Return to citation in text: [1]
  54. Hamster, H.; Sullivan, A.; Gordon, S.; White, W.; Falcone, R. W. Phys. Rev. Lett. 1993, 71, 2725–2728. doi:10.1103/physrevlett.71.2725
    Return to citation in text: [1] [2]
  55. Zhang, X. C.; Shkurinov, A.; Zhang, Y. Nat. Photonics 2017, 11, 16–18. doi:10.1038/nphoton.2016.249
    Return to citation in text: [1]

Interesting articles

Imène Reghioua, Mattia Fanetti, Sylvain Girard, Diego Di Francesca, Simonpietro Agnello, Layla Martin-Samos, Marco Cannas, Matjaz Valant, Melanie Raine, Marc Gaillardin, Nicolas Richard, Philippe Paillet, Aziz Boukenter, Youcef Ouerdane and Antonino Alessi

Pavel Loiko, Tanujjal Bora, Josep Maria Serres, Haohai Yu, Magdalena Aguiló, Francesc Díaz, Uwe Griebner, Valentin Petrov, Xavier Mateos and Joydeep Dutta

Eduardo Aluicio-Sarduy, Simone Callegari, Diana Gisell Figueroa del Valle, Andrea Desii, Ilka Kriegel and Francesco Scotognella

© 2018 Hsu et al.; licensee Beilstein-Institut.
This is an Open Access article under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0). Please note that the reuse, redistribution and reproduction in particular requires that the authors and source are credited.
The license is subject to the Beilstein Journal of Nanotechnology terms and conditions: (https://www.beilstein-journals.org/bjnano)

Back to Article List

Other Beilstein-Institut Open Science Activities

Keep Informed

RSS Feed

Subscribe to our Latest Articles RSS Feed.


Follow the Beilstein-Institut


Twitter: @BeilsteinInst